Changeset 18143

Timestamp:

06/12/14 11:59:42 (11 years ago)

Author:

Mathieu Morlighem

Message:

NEW: some functionalities required for Nathan's algorithm in 3d

Location:

issm/trunk-jpl/src/c/classes/Elements

Files:

• Tetra.cpp (modified) (1 diff)
• TetraRef.cpp (modified) (4 diffs)

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

@@ -844 +844 @@
                         break;
                 case TaylorHoodEnum:
+                case XTaylorHoodEnum:
                         numnodes        = 14;
                         tetra_node_ids  = xNew<int>(numnodes);

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

@@ -354 +354 @@
                 case MINIEnum:              return NUMNODESP1b+NUMNODESP1;
                 case TaylorHoodEnum:        return NUMNODESP2+NUMNODESP1;
+                case XTaylorHoodEnum:       return NUMNODESP2+NUMNODESP1;
                 default: _error_("Element type "<<EnumToStringx(finiteelement)<<" not supported yet");
         }
@@ -368 +369 @@
                 case MINIEnum:          return P1bubbleEnum;
                 case TaylorHoodEnum:    return P2Enum;
+                case XTaylorHoodEnum:   return P2Enum;
                 default:       _error_("Element type "<<EnumToStringx(fe_stokes)<<" not supported yet");
         }
@@ -374 +376 @@
 }
 /*}}}*/
-int TetraRef::PressureInterpolation(int fe_stokes){/*{{{*/
+int  TetraRef::PressureInterpolation(int fe_stokes){/*{{{*/
 
         switch(fe_stokes){
@@ -382 +384 @@
                 case MINIEnum:          return P1Enum;
                 case TaylorHoodEnum:    return P1Enum;
+                case XTaylorHoodEnum:   return P1Enum;
                 default:       _error_("Element type "<<EnumToStringx(fe_stokes)<<" not supported yet");
         }