Changeset 4479

Timestamp:

07/08/10 15:48:26 (15 years ago)

Author:

seroussi

Message:

InputUpdate for Balancedthickness in Penta

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (1 diff)

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

@@ -2556 +2556 @@
 /*FUNCTION Penta::InputUpdateFromSolutionBalancedthickness {{{1*/
 void  Penta::InputUpdateFromSolutionBalancedthickness(double* solution){
-        ISSMERROR(" not supported yet!");
+
+        int i;
+
+        const int    numvertices=6;
+        const int    numdofpervertex=1;
+        const int    numdof=numdofpervertex*numvertices;
+
+        int          doflist[numdof];
+        double       values[numdof];
+
+        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 PentaVertexInput(ThicknessEnum,values));
 }
 /*}}}*/