source: issm/oecreview/Archive/15392-16133/ISSM-15768-15769.diff

../trunk-jpl/src/c/classes/Elements/Penta.cpp

@@ -6875 +6875 @@
 ElementMatrix* Penta::CreateKMatrixCouplingSSAFSViscous(void){
 
         /*Constants*/
-        const int numdofm      = NDOF2 *NUMVERTICES2D;
-        const int numdofs      = NDOF4 *NUMVERTICES;
-        const int numdofstotal = NDOF4 *NUMVERTICES + NDOF3;
-        const int numdoftotal  = 2 *numdofm+numdofstotal;
+        const int numdofm     = NDOF2 *NUMVERTICES2D;
+        const int numdofs     = NDOF4 *NUMVERTICES + NDOF3;
+        const int numdoftotal = 2 *numdofm+numdofs;
 
         /*Intermediaries */
         int        i,j;
@@ -6886 +6885 @@
         IssmDouble viscosity,FSreconditioning; //viscosity
         IssmDouble epsilon[6]; /* epsilon=[exx,eyy,exy,exz,eyz];*/
         IssmDouble xyz_list[NUMVERTICES][3];
-        IssmDouble B[4][numdofstotal];
+        IssmDouble B[4][numdofs];
         IssmDouble Bprime[4][numdofm];
         IssmDouble B2[3][numdofm];
-        IssmDouble Bprime2[3][numdofstotal];
+        IssmDouble Bprime2[3][numdofs];
         IssmDouble D[4][4]={0.0};            // material matrix, simple scalar matrix.
         IssmDouble D2[3][3]={0.0};            // material matrix, simple scalar matrix.
         IssmDouble D_scalar;
         IssmDouble Ke_gg[numdofs][numdofm]={0.0}; //local element stiffness matrix 
         IssmDouble Ke_gg2[numdofm][numdofs]={0.0}; //local element stiffness matrix 
-        IssmDouble Ke_gg_gaussian[numdofstotal][numdofm]; //stiffness matrix evaluated at the gaussian point.
-        IssmDouble Ke_gg_gaussian2[numdofm][numdofstotal]; //stiffness matrix evaluated at the gaussian point.
+        IssmDouble Ke_gg_gaussian[numdofs][numdofm]; //stiffness matrix evaluated at the gaussian point.
+        IssmDouble Ke_gg_gaussian2[numdofm][numdofs]; //stiffness matrix evaluated at the gaussian point.
         GaussPenta *gauss=NULL;
         GaussTria  *gauss_tria=NULL;
         Node       *node_list[20];
@@ -6960 +6959 @@
                 D[3][3]=-gauss->weight*Jdet*FSreconditioning;
                 for (i=0;i<3;i++) D2[i][i]=D_scalar;
 
-                TripleMultiply( &B[0][0],4,numdofstotal,1,
+                TripleMultiply( &B[0][0],4,numdofs,1,
                                         &D[0][0],4,4,0,
                                         &Bprime[0][0],4,numdofm,0,
                                         &Ke_gg_gaussian[0][0],0);
 
                 TripleMultiply( &B2[0][0],3,numdofm,1,
                                         &D2[0][0],3,3,0,
-                                        &Bprime2[0][0],3,numdofstotal,0,
+                                        &Bprime2[0][0],3,numdofs,0,
                                         &Ke_gg_gaussian2[0][0],0);
 
-                for( i=0; i<numdofstotal; i++) for(j=0;j<numdofm; j++)      Ke_gg[i][j]+=Ke_gg_gaussian[i][j];
-                for( i=0; i<numdofm; i++)      for(j=0;j<numdofstotal; j++) Ke_gg2[i][j]+=Ke_gg_gaussian2[i][j];
+                for( i=0; i<numdofs; i++) for(j=0;j<numdofm; j++)      Ke_gg[i][j]+=Ke_gg_gaussian[i][j];
+                for( i=0; i<numdofm; i++)      for(j=0;j<numdofs; j++) Ke_gg2[i][j]+=Ke_gg_gaussian2[i][j];
         } 
-        for(i=0;i<numdofs;i++) for(j=0;j<numdofm;j++) Ke->values[(i+2*numdofm)*numdoftotal+j]+=Ke_gg[i][j];
-        for(i=0;i<numdofm;i++) for(j=0;j<numdofs;j++) Ke->values[i*numdoftotal+(j+2*numdofm)]+=Ke_gg2[i][j];
+        for(i=0;i<numdofs;i++) for(j=0;j<numdofm;j++)      Ke->values[(i+2*numdofm)*numdoftotal+j]+=Ke_gg[i][j];
+        for(i=0;i<numdofm;i++)      for(j=0;j<numdofs;j++) Ke->values[i*numdoftotal+(j+2*numdofm)]+=Ke_gg2[i][j];
 
         /*Transform Coordinate System*/
         TransformStiffnessMatrixCoord(Ke,node_list,numnodes,cs_list);