Changeset 16990

Timestamp:

12/03/13 09:21:07 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: SSAHO supposed to call SSA3d and not SSA

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (modified) (1 diff)

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

@@ -3199 +3199 @@
 
         /*compute all stiffness matrices for this element*/
-        ElementMatrix* Ke1=CreateKMatrixSSA(element);
+        ElementMatrix* Ke1=CreateKMatrixSSA3d(element);
         ElementMatrix* Ke2=CreateKMatrixHO(element);
         ElementMatrix* Ke3=CreateKMatrixCouplingSSAHO(element);
         ElementMatrix* Ke =new ElementMatrix(Ke1,Ke2,Ke3);
+
+        //Ke2=element->CreateKMatrix();
+        //for(int i=0;i<12;i++){
+        //      for(int j=0;j<12;j++){
+        //              if(Ke->values[i*12+j] != Ke2->values[i*12+j]){
+        //                      Ke->Echo();
+        //                      Ke2->Echo();
+        //                      printf("ID = %i\n",element->Id());
+        //                      _error_("S");
+        //              }
+        //      }
+        //}
 
         /*clean-up and return*/