Changeset 27059

Timestamp:

06/14/22 11:48:36 (3 years ago)

Author:

caronlam

Message:

BUG: memory management issue

Location:

issm/trunk-jpl/src/c

Files:

• classes/Elements/Tria.cpp (modified) (1 diff)
• cores/sealevelchange_core.cpp (modified) (4 diffs)

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

@@ -7482 +7482 @@
         }
 
+                cout << "DEBUG: Sealevelchange convo:";
         xDelete<IssmDouble>(G);
+                cout << " step1";
         xDelete<IssmDouble>(Grot);
+                cout << ", step2";
         for (int l=0;l<SLGEOM_NUMLOADS;l++) xDelete<IssmDouble>(Gsub[l]);
+                cout << ", step3";
         xDelete<IssmDouble*>(Gsub);
+                cout << ", step4\n";
 
         return;

issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp

@@ -340 +340 @@
                 }
 
+                cout << "DEBUG: Sealevelchange convolution cleared\n";
+
                 /*retrieve sea level average  and ocean area:*/
                 for(Object* & object : femmodel->elements->objects){
@@ -345 +347 @@
                         element->SealevelchangeOceanAverage(loads, oceanareas, subelementoceanareas, sealevelpercpu, slgeom);
                 }
+
+                cout << "DEBUG: Sealevelchange Oceanavaverage cleared\n";
 
                 loads->AssembleSealevelLoads();
@@ -360 +364 @@
                 ConserveOceanMass(femmodel,loads,barycontrib->Total()/totaloceanarea - oceanaverage,slgeom);
 
+                cout << "DEBUG: Sealevelchange Conserve ocean mass cleared\n";
+
                 //broadcast sea level loads 
                 loads->BroadcastSealevelLoads();
@@ -366 +372 @@
                 if(slcconvergence(loads->vsealevelloads,oldsealevelloads,eps_rel,eps_abs))break;
 
+                cout << "DEBUG: Sealevelchange convergence cleared\n";
                 //early return?
                 if(iterations>=max_nonlinear_iterations)break;