Ignore:
Timestamp:
11/08/11 11:37:26 (13 years ago)
Author:
Mathieu Morlighem
Message:

Extended CS rotation to pattyn Stokes coupling

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/objects/Elements/Penta.cpp

    r10524 r10526  
    55475547
    55485548        /*Constants*/
    5549         const int    numdofp=NDOF2*NUMVERTICES;
    5550         const int    numdofs=NDOF4*NUMVERTICES;
    5551         const int    numdoftotal=(NDOF2+NDOF4)*NUMVERTICES;
    5552         int          i,j;
     5549        const int numdofp     = NDOF2 *NUMVERTICES;
     5550        const int numdofs     = NDOF4 *NUMVERTICES;
     5551        const int numdoftotal = (NDOF2+NDOF4) *NUMVERTICES;
     5552        int       i,j;
    55535553
    55545554        for(i=0;i<numdofs;i++) for(j=0;j<NUMVERTICES;j++){
     
    64696469        }
    64706470
     6471        /*Transform coordinate system*/
     6472        TransformLoadVectorCoord(pe,nodes,NUMVERTICES,XYZPEnum);
     6473
    64716474        /*Clean up and return*/
    64726475        delete gauss;
     
    65396542                }
    65406543        }
     6544
     6545        /*Transform coordinate system*/
     6546        TransformLoadVectorCoord(pe,nodes,NUMVERTICES,XYZPEnum);
    65416547
    65426548        /*Clean up and return*/
     
    76327638        for(i=0;i<numdofs;i++) stokes_values[i]=solution[doflists[i]];
    76337639
     7640        /*Transform solution in Cartesian Space*/
     7641        TransformSolutionCoord(&pattyn_values[0],this->nodes,NUMVERTICES,XYEnum);
     7642        TransformSolutionCoord(&stokes_values[0],this->nodes,NUMVERTICES,XYZPEnum);
     7643
    76347644        /*Ok, we have vx and vy in values, fill in vx and vy arrays: */
    76357645        for(i=0;i<NUMVERTICES;i++){
Note: See TracChangeset for help on using the changeset viewer.