Changeset 10526 for issm/trunk/src/c/objects/Elements/Penta.cpp

Timestamp:

11/08/11 11:37:26 (13 years ago)

Author:

Mathieu Morlighem

Message:

Extended CS rotation to pattyn Stokes coupling

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (4 diffs)

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

@@ -5547 +5547 @@
 
         /*Constants*/
-        const int    numdofp=NDOF2*NUMVERTICES;
-        const int    numdofs=NDOF4*NUMVERTICES;
-        const int    numdoftotal=(NDOF2+NDOF4)*NUMVERTICES;
-        int          i,j;
+        const int numdofp     = NDOF2 *NUMVERTICES;
+        const int numdofs     = NDOF4 *NUMVERTICES;
+        const int numdoftotal = (NDOF2+NDOF4) *NUMVERTICES;
+        int       i,j;
 
         for(i=0;i<numdofs;i++) for(j=0;j<NUMVERTICES;j++){
@@ -6469 +6469 @@
         }
 
+        /*Transform coordinate system*/
+        TransformLoadVectorCoord(pe,nodes,NUMVERTICES,XYZPEnum);
+
         /*Clean up and return*/
         delete gauss;
@@ -6539 +6542 @@
                 }
         }
+
+        /*Transform coordinate system*/
+        TransformLoadVectorCoord(pe,nodes,NUMVERTICES,XYZPEnum);
 
         /*Clean up and return*/
@@ -7632 +7638 @@
         for(i=0;i<numdofs;i++) stokes_values[i]=solution[doflists[i]];
 
+        /*Transform solution in Cartesian Space*/
+        TransformSolutionCoord(&pattyn_values[0],this->nodes,NUMVERTICES,XYEnum);
+        TransformSolutionCoord(&stokes_values[0],this->nodes,NUMVERTICES,XYZPEnum);
+
         /*Ok, we have vx and vy in values, fill in vx and vy arrays: */
         for(i=0;i<NUMVERTICES;i++){