Changeset 18191

Timestamp:

06/26/14 16:00:59 (11 years ago)

Author:

seroussi

Message:

BUG: should not transform coordinates for augmented pressure

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

@@ -2942 +2942 @@
         }
 
+        /*Transform Coordinate System*/
+        element->TransformStiffnessMatrixCoord(Ke,cs_list);
+
+        /*The pressure augmentation should not be transformed*/
         MatrixMultiply(BtBUzawa,pnumdof,numdof,1,
                                 BtBUzawa,pnumdof,numdof,0,
                                 &Ke->values[0],1);
-
-        /*Transform Coordinate System*/
-        element->TransformStiffnessMatrixCoord(Ke,cs_list);
 
         /*Clean up and return*/
@@ -3723 +3724 @@
 
         /*Transform coordinate system*/
-        element->TransformLoadVectorCoord(pe,cs_list);
+        //element->TransformLoadVectorCoord(pe,cs_list); Do not transform pressure augmentation
 
         /*Clean up and return*/