Changeset 23206

Timestamp:

08/31/18 14:29:03 (7 years ago)

Author:

wester

Message:

CHG: PressureMassMatrix and SchurPreconditioner can now also be constructed in 2D

File:

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

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

@@ -3279 +3279 @@
                                 }
                         }
+                }else if(dim==2){
+                        /*Pressure mass matrix*/
+                        for(int k=0;k<pnumnodes;k++){
+                                for(int j=0;j<pnumnodes;j++){
+                                        Ke->values[(2*vnumnodes+k)*numdof+2*vnumnodes+j] += gauss->weight*Jdet*(pbasis[j]*pbasis[k]);
+                                }
+                        }
                 }
                 else{
@@ -3334 +3341 @@
                                 for(int j=0;j<pnumnodes;j++){
                                         Ke->values[(3*vnumnodes+k)*numdof+3*vnumnodes+j] += gauss->weight*1./viscosity*Jdet*(pbasis[j]*pbasis[k]);
+                                }
+                        }
+                }else if(dim==2){
+                        /*Pressure mass matrix*/
+                        for(int k=0;k<pnumnodes;k++){
+                                for(int j=0;j<pnumnodes;j++){
+                                        Ke->values[(2*vnumnodes+k)*numdof+2*vnumnodes+j] += gauss->weight*1./viscosity*Jdet*(pbasis[j]*pbasis[k]);
                                 }
                         }