Changeset 26134

Timestamp:

03/23/21 12:53:17 (4 years ago)

Author:

tsantos

Message:

CHG: working on mono layer HO

File:

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

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

@@ -2727 +2727 @@
         if(!element->IsOnBase() || element->IsFloating()) return NULL;
         Element* basalelement = element->SpawnBasalElement();
-        ElementMatrix* Ke = CreateKMatrixSSAFriction(basalelement);
+        ElementMatrix* Ke    = basalelement->NewElementMatrix(MLHOApproximationEnum);
+        ElementMatrix* KeSSA = CreateKMatrixSSAFriction(basalelement); //only to get K11 and K33
+        
+        /*Fetch number of nodes and dof for this finite element*/
+        int numnodes = basalelement->GetNumberOfNodes();
+
+        for(int i=0;i<numnodes;i++){
+      for(int j=0;j<numnodes;j++){
+         Ke->values[(i+0)*2*2*numnodes+j+0] = KeSSA->values[2*i*2*numnodes+2*j]; //K11
+         Ke->values[(i+6)*2*2*numnodes+j+6] = KeSSA->values[(2*i+1)*2*numnodes+2*j+1]; //K33
+      }  
+   }
 
         /*clean-up and return*/
         if(basalelement->IsSpawnedElement()){basalelement->DeleteMaterials(); delete basalelement;};
+        delete KeSSA;
         return Ke;
 }/*}}}*/