Changeset 15544

Timestamp:

07/22/13 15:38:34 (12 years ago)

Author:

seroussi

Message:

FIX: Pvector front based on numnodes (not numvertices)

Location:

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

Files:

• Penta.cpp (modified) (2 diffs)
• Tria.cpp (modified) (1 diff)

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

@@ -8000 +8000 @@
 
         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-        ElementVector* pe=new ElementVector(nodes,NUMVERTICES,this->parameters,PattynApproximationEnum);
+        ElementVector* pe=new ElementVector(nodes,numnodes,this->parameters,PattynApproximationEnum);
         Input* surface_input=inputs->GetInput(SurfaceEnum); _assert_(surface_input);
         rho_water=matpar->GetRhoWater();
@@ -8205 +8205 @@
 
         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-        ElementVector* pe=new ElementVector(nodes,NUMVERTICES,this->parameters,StokesApproximationEnum);
+        ElementVector* pe=new ElementVector(nodes,numnodes,this->parameters,StokesApproximationEnum);
         rho_water=matpar->GetRhoWater();
         rho_ice  =matpar->GetRhoIce();

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

@@ -3181 +3181 @@
 
         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-        ElementVector* pe=new ElementVector(nodes,NUMVERTICES,this->parameters,MacAyealApproximationEnum);
+        ElementVector* pe=new ElementVector(nodes,numnodes,this->parameters,MacAyealApproximationEnum);
         Input* thickness_input=inputs->GetInput(ThicknessEnum); _assert_(thickness_input);
         Input* bed_input      =inputs->GetInput(BedEnum);       _assert_(bed_input);