Changeset 15925

Timestamp:

08/24/13 20:51:28 (12 years ago)

Author:

Mathieu Morlighem

Message:

BUG: IoToConstrain needs velocity interpolation only

File:

• issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateConstraintsStressbalance.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateConstraintsStressbalance.cpp

@@ -74 +74 @@
                 else if(isL1L2) finiteelement = P1Enum;
                 else if(isHO)   iomodel->Constant(&finiteelement,FlowequationFeHOEnum);
-                else if(isFS)   iomodel->Constant(&finiteelement,FlowequationFeFSEnum);
+                else if(isFS){  iomodel->Constant(&finiteelement,FlowequationFeFSEnum);
+                        /*Deduce velocity interpolation from finite element*/
+                        switch(finiteelement){
+                                case P1P1Enum          : finiteelement = P1Enum;       break;
+                                case P1P1GLSEnum       : finiteelement = P1Enum;       break;
+                                case MINIcondensedEnum : finiteelement = P1bubbleEnum; break;
+                                case MINIEnum          : finiteelement = P1bubbleEnum; break;
+                                case TaylorHoodEnum    : finiteelement = P2Enum;       break;
+                                default: _error_("finite element "<<finiteelement<<" not supported");
+                        }
+                }
 
                 IoModelToConstraintsx(constraints,iomodel,StressbalanceSpcvxEnum,StressbalanceAnalysisEnum,finiteelement,1);