Changeset 27474

Timestamp:

12/21/22 13:03:32 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixing destructor crash

File:

• issm/trunk-jpl/src/c/classes/Loads/Friction.cpp (modified) (5 diffs)

issm/trunk-jpl/src/c/classes/Loads/Friction.cpp

@@ -19 +19 @@
 Friction::Friction(){/*{{{*/
         this->element=NULL;
-        this->law=0;
-        this->linearize=false;
-        this->apply_dim = 1.;
-        this->domaintype=-1;
+        this->law        = 0;
+        this->linearize  = 0;
+        this->apply_dim  = 1.;
+        this->domaintype = -1;
         this->vx_input=NULL;
         this->vy_input=NULL;
@@ -33 +33 @@
         /* Determine the dimension according to the domain type automatically. 
          * There are exceptions, e.g. HO, which needs the user to specify the dimension used in Friciton.*/
+
         this->element=element_in;
+        this->linearize  = 0;
 
         /* Load necessary parameters */
@@ -65 +67 @@
                 element_in->FindParam(&this->linearize,FrictionLinearizeEnum);
                 if(this->linearize){
-                        this->linearize = false; /*Change to make sure we do the calculation once*/
+                        this->linearize = 0; /*Change to make sure we do the calculation once*/
                         int numvertices = this->element->GetNumberOfVertices();
                         this->alpha2_list            = xNew<IssmDouble>(numvertices);
@@ -77 +79 @@
                                 this->GetAlphaComplement(&this->alpha2_complement_list[iv], gauss);
                         }
-                        this->linearize = true ; /*Change back, we are now all set!*/
+                        this->linearize = 1; /*Change back, we are now all set!*/
                 }
         }
@@ -92 +94 @@
 Friction::~Friction(){/*{{{*/
         if(this->linearize){
+                printf("linearize = %i\n",this->linearize);
                 xDelete<IssmDouble>(this->alpha2_list);
                 xDelete<IssmDouble>(this->alpha2_complement_list);