source: issm/oecreview/Archive/24684-25833/ISSM-25727-25728.diff

../trunk-jpl/src/c/classes/Elements/Element.cpp

@@ -254 +254 @@
 
         /*Allocate arrays*/
         const int NUM_VERTICES = this->GetNumberOfVertices();
+        _assert_(NUM_VERTICES<=MAXVERTICES);
 
-        IssmDouble* eps_xx = xNew<IssmDouble>(NUM_VERTICES);
-        IssmDouble* eps_yy = xNew<IssmDouble>(NUM_VERTICES);
-        IssmDouble* eps_zz = xNew<IssmDouble>(NUM_VERTICES);
-        IssmDouble* eps_xy = xNew<IssmDouble>(NUM_VERTICES);
-        IssmDouble* eps_xz = xNew<IssmDouble>(NUM_VERTICES);
-        IssmDouble* eps_yz = xNew<IssmDouble>(NUM_VERTICES);
-        IssmDouble* eps_ef = xNew<IssmDouble>(NUM_VERTICES);
+        IssmDouble eps_xx[MAXVERTICES];
+        IssmDouble eps_yy[MAXVERTICES];
+        IssmDouble eps_zz[MAXVERTICES];
+        IssmDouble eps_xy[MAXVERTICES];
+        IssmDouble eps_xz[MAXVERTICES];
+        IssmDouble eps_yz[MAXVERTICES];
+        IssmDouble eps_ef[MAXVERTICES];
 
         /* Start looping on the number of vertices: */
         Gauss* gauss=this->NewGauss();
@@ -296 +297 @@
         }
 
         /*Add Stress tensor components into inputs*/
-        this->AddInput(StrainRatexxEnum,eps_xx,P1Enum);
-        this->AddInput(StrainRatexyEnum,eps_xy,P1Enum);
-        this->AddInput(StrainRatexzEnum,eps_xz,P1Enum);
-        this->AddInput(StrainRateyyEnum,eps_yy,P1Enum);
-        this->AddInput(StrainRateyzEnum,eps_yz,P1Enum);
-        this->AddInput(StrainRatezzEnum,eps_zz,P1Enum);
-        this->AddInput(StrainRateeffectiveEnum,eps_ef,P1Enum);
+        this->AddInput(StrainRatexxEnum,&eps_xx[0],P1Enum);
+        this->AddInput(StrainRatexyEnum,&eps_xy[0],P1Enum);
+        this->AddInput(StrainRatexzEnum,&eps_xz[0],P1Enum);
+        this->AddInput(StrainRateyyEnum,&eps_yy[0],P1Enum);
+        this->AddInput(StrainRateyzEnum,&eps_yz[0],P1Enum);
+        this->AddInput(StrainRatezzEnum,&eps_zz[0],P1Enum);
+        this->AddInput(StrainRateeffectiveEnum,&eps_ef[0],P1Enum);
 
         /*Clean up and return*/
         delete gauss;
         xDelete<IssmDouble>(xyz_list);
-        xDelete<IssmDouble>(eps_xx);
-        xDelete<IssmDouble>(eps_yy);
-        xDelete<IssmDouble>(eps_zz);
-        xDelete<IssmDouble>(eps_xy);
-        xDelete<IssmDouble>(eps_xz);
-        xDelete<IssmDouble>(eps_yz);
-        xDelete<IssmDouble>(eps_ef);
-
 }
 /*}}}*/
 void       Element::CoordinateSystemTransform(IssmDouble** ptransform,Node** nodes_list,int numnodes,int* cs_array){/*{{{*/