Changeset 15420 for issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

Timestamp:

07/03/13 16:10:55 (12 years ago)

Author:

seroussi

Message:

NEW: starting stokes stabilization implementation

File:

• issm/trunk-jpl/src/c/classes/Elements/Penta.cpp (modified) (2 diffs)

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

@@ -6822 +6822 @@
 ElementMatrix* Penta::CreateKMatrixDiagnosticStokes(void){
 
-        /*compute all stiffness matrices for this element*/
-        ElementMatrix* Ke1=CreateKMatrixDiagnosticStokesViscous();
-        ElementMatrix* Ke2=CreateKMatrixDiagnosticStokesFriction();
-        ElementMatrix* Ke =new ElementMatrix(Ke1,Ke2);
+        int fe_stokes;
+        ElementMatrix* Ke1;
+        ElementMatrix* Ke2;
+        ElementMatrix* Ke;
+        parameters->FindParam(&fe_stokes,FlowequationFeStokesEnum);
+
+        switch(fe_stokes){
+                case 0:
+                        /*compute all stiffness matrices for this element*/
+                        Ke1=CreateKMatrixDiagnosticStokesViscous();
+                        Ke2=CreateKMatrixDiagnosticStokesFriction();
+                        Ke =new ElementMatrix(Ke1,Ke2);
+
+                case 1:
+                        /*compute all stiffness matrices for this element*/
+                        Ke1=CreateKMatrixDiagnosticStokesViscous();
+                        Ke2=CreateKMatrixDiagnosticStokesFriction();
+                        Ke =new ElementMatrix(Ke1,Ke2);
+
+                default:
+                        _error_("Finite element" << fe_stokes << " not supported yet");
+        }
 
         /*clean-up and return*/
@@ -6831 +6849 @@
         delete Ke2;
         return Ke;
+
 }
 /*}}}*/