Changeset 15476

Timestamp:

07/10/13 10:01:45 (12 years ago)

Author:

seroussi

Message:

minor

Location:

issm/trunk-jpl/src/c/classes/Elements

Files:

• Penta.cpp (modified) (2 diffs)
• PentaRef.cpp (modified) (1 diff)

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

@@ -7029 +7029 @@
         /*Transform Coordinate System*/
         TransformStiffnessMatrixCoord(Ke,nodes,NUMVERTICES,XYZPEnum);
-        if(id==24){
-                //printarray(Ke->values,24,24);
-        }
 
         /*Clean up and return*/
@@ -7985 +7982 @@
                         for(p=0;p<6;p++){
                                 for(j=0;j<4;j++){
-                                        pe->values[p*4+j] += gauss->weight*Jdet*tau*rho_ice*forcex*SW[p][j][0];
-                                        pe->values[p*4+j] += gauss->weight*Jdet*tau*rho_ice*forcey*SW[p][j][1];
-                                        pe->values[p*4+j] += gauss->weight*Jdet*tau*rho_ice*(forcez-gravity)*SW[p][j][2];
+                                        pe->values[p*4+j] += gauss->weight*Jdet*tau*forcex*SW[p][j][0];
+                                        pe->values[p*4+j] += gauss->weight*Jdet*tau*forcey*SW[p][j][1];
+                                        pe->values[p*4+j] += gauss->weight*Jdet*tau*(forcez-rho_ice*gravity)*SW[p][j][2];
                                 }
                         }

issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp ¶

@@ -509 +509 @@
                 B_prime[(NDOF4*NUMNODESP1)*5+NDOF4*i+3] = 0.;
                 B_prime[(NDOF4*NUMNODESP1)*6+NDOF4*i+3] = 0.;
-                B_prime[(NDOF4*NUMNODESP1)*7+NDOF4*i+3] = basis[i];
+                B_prime[(NDOF4*NUMNODESP1)*7+NDOF4*i+3] = - basis[i];
         }