Changeset 3976

Timestamp:

05/27/10 15:26:16 (15 years ago)

Author:

seroussi

Message:

added UpdateInputsFromSolution for balancedthickness in Tria

Location:

issm/trunk/src/c/objects/Elements

Files:

• Penta.cpp (modified) (1 diff)
• Tria.cpp (modified) (1 diff)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -2116 +2116 @@
                 }
         }
-
 
         /*Add Ke_gg to global matrix Kgg: */

issm/trunk/src/c/objects/Elements/Tria.cpp

@@ -599 +599 @@
 /*FUNCTION Tria::UpdateInputsFromSolutionBalancedthickness {{{1*/
 void  Tria::UpdateInputsFromSolutionBalancedthickness(double* solution, int analysis_type, int sub_analysis_type){
-        ISSMERROR(" not supported yet!");
+
+        int i;
+
+        const int    numvertices=3;
+        const int    numdofpervertex=1;
+        const int    numdof=numdofpervertex*numvertices;
+
+        int          doflist[numdof];
+        double       values[numdof];
+        double       thickness[numvertices];
+
+        int          dummy;
+
+        /*Get dof list: */
+        GetDofList(&doflist[0],&dummy);
+
+        /*Use the dof list to index into the solution vector: */
+        for(i=0;i<numdof;i++){
+                values[i]=solution[doflist[i]];
+        }
+
+        /*Add thickness as inputs to the tria element: */
+        this->inputs->AddInput(new TriaVertexInput(ThicknessEnum,values));
 }
 /*}}}*/