Changeset 26107

Timestamp:

03/16/21 10:21:57 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added test on size to avoid segmentation faults

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (4 diffs)

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

@@ -5600 +5600 @@
         /*elastic green function:*/
         int index;
-        int         M;
+        int M,Mtest;
 
         /*Computational flags:*/
@@ -5625 +5625 @@
         DoubleVecParam* parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelchangeGRigidEnum)); _assert_(parameter);
         parameter->GetParameterValueByPointer((IssmDouble**)&G_rigid_precomputed,&M);
+        _assert_(M>0);
 
         /*allocate indices:*/
@@ -5632 +5633 @@
         if(computeelastic){
                 parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelchangeGElasticEnum)); _assert_(parameter);
-                parameter->GetParameterValueByPointer((IssmDouble**)&G_elastic_precomputed,&M);
+                parameter->GetParameterValueByPointer((IssmDouble**)&G_elastic_precomputed,&Mtest); _assert_(Mtest==M);
 
                 parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelchangeHElasticEnum)); _assert_(parameter);
-                parameter->GetParameterValueByPointer((IssmDouble**)&H_elastic_precomputed,&M);
+                parameter->GetParameterValueByPointer((IssmDouble**)&H_elastic_precomputed,&Mtest); _assert_(Mtest==M);
 
                 parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelchangeUElasticEnum)); _assert_(parameter);
-                parameter->GetParameterValueByPointer((IssmDouble**)&U_elastic_precomputed,&M);
+                parameter->GetParameterValueByPointer((IssmDouble**)&U_elastic_precomputed,&Mtest); _assert_(Mtest==M);
 
                 GU=xNewZeroInit<IssmDouble>(gsize);
@@ -5673 +5674 @@
                 indices[i]=alpha/M_PI*reCast<IssmDouble,int>(M-1);
                 index=reCast<int,IssmDouble>(indices[i]);
+                _assert_(index>=0); _assert_(index<M);
 
                 /*Rigid earth gravitational perturbation: */