Changeset 25468

Timestamp:

08/25/20 13:57:17 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: working on adding spcthickness as control in MC

Location:

issm/trunk-jpl/src/c/classes/Constraints

Files:

• SpcStatic.cpp (modified) (1 diff)
• SpcStatic.h (modified) (2 diffs)

issm/trunk-jpl/src/c/classes/Constraints/SpcStatic.cpp

@@ -152 +152 @@
 /*}}}*/
 
+void SpcStatic::UpdateSpcThicknessAD(IssmDouble* vector,Nodes* nodes){/*{{{*/
+
+        /*Chase through nodes and find the node to which this SpcStatic applies: */
+        Node* node=(Node*)nodes->GetObjectById(NULL,nodeid);
+
+        /*Apply constraint: */
+        if(node){ //in case the spc is dealing with a node on another cpu
+                int sid = node->Sid();
+                this->value = vector[sid];
+                _assert_(!xIsNan<IssmDouble>(this->value)); 
+        }
+}
+/*}}}*/
+
 /*SpcStatic functions*/
 int        SpcStatic::GetDof(){/*{{{*/

issm/trunk-jpl/src/c/classes/Constraints/SpcStatic.h

@@ -14 +14 @@
 
         private: 
-                int        id;             /*! id, to track it */
-                int        nodeid;          /*!node id          */
-                int        dof;             /*!component        */
-                IssmDouble value;           /*value             */
+                int        id;
+                int        nodeid;
+                int        dof;
+                IssmDouble value;
                 int        analysis_type;
-                bool       penalty;         /*Is this a penalty constraint */
+                bool       penalty;
 
         public:
@@ -47 +47 @@
                 IssmDouble GetValue();
                 /*}}}*/
+                void UpdateSpcThicknessAD(IssmDouble* vector,Nodes* nodes);
 
 };