Changeset 18182

Timestamp:

06/25/14 09:52:47 (11 years ago)

Author:

seroussi

Message:

BUG: starting to debug LATaylorHood

Location:

issm/trunk-jpl/src/c

Files:

• analyses/StressbalanceAnalysis.cpp (modified) (3 diffs)
• analyses/UzawaPressureAnalysis.cpp (modified) (1 diff)
• classes/Elements/TriaRef.cpp (modified) (1 diff)

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

@@ -420 +420 @@
                                 case TaylorHoodEnum      : finiteelement = P2Enum;       break;
                                 case XTaylorHoodEnum     : finiteelement = P2Enum;       break;
+                                case LATaylorHoodEnum    : finiteelement = P2Enum;       break;
                                 case OneLayerP4zEnum     : finiteelement = P2xP4Enum;    break;
                                 case CrouzeixRaviartEnum : finiteelement = P2bubbleEnum; break;
@@ -2854 +2855 @@
         if(fe_FS==XTaylorHoodEnum)
          Ke1=CreateKMatrixFSViscousXTH(element);
-        else if(fe_FS==LATaylorHoodEnum)
+        else if(fe_FS==LATaylorHoodEnum){
          Ke1=CreateKMatrixFSViscousLATH(element);
+        }
         else
          Ke1=CreateKMatrixFSViscous(element);
@@ -3710 +3712 @@
                 
                 pressure_input->GetInputValue(&pressure, gauss);
-                element->NodalFunctionsDerivatives(dbasis,xyz_list,gauss);
+                element->NodalFunctionsDerivativesVelocity(dbasis,xyz_list,gauss);
 
                 for(i=0;i<numnodes;i++){

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

@@ -26 +26 @@
         }
 
-        iomodel->FetchDataToInput(elements,VxEnum);
-        iomodel->FetchDataToInput(elements,VyEnum);
-        if(iomodel->domaintype==Domain3DEnum) iomodel->FetchDataToInput(elements,VzEnum);
-        iomodel->FetchDataToInput(elements,PressureEnum);
+        iomodel->FetchDataToInput(elements,VxEnum,0.);
+        iomodel->FetchDataToInput(elements,VyEnum,0.);
+        if(iomodel->domaintype==Domain3DEnum) iomodel->FetchDataToInput(elements,VzEnum,0.);
+        iomodel->FetchDataToInput(elements,PressureEnum,0.);
 }/*}}}*/
 void UzawaPressureAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel){/*{{{*/

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

@@ -149 +149 @@
 
         switch(finiteelement){
+                case NoneEnum:
+                        return;
                 case P0Enum:
                         basis[0]=1.;