Changeset 2007

Timestamp:

08/27/09 11:04:28 (16 years ago)

Author:

Mathieu Morlighem

Message:

minor leaks

Location:

issm/trunk/src/c/parallel

Files:

• steadystate_core.cpp (modified) (2 diffs)
• thermal_core_nonlinear.cpp (modified) (3 diffs)

issm/trunk/src/c/parallel/steadystate_core.cpp

@@ -111 +111 @@
                 //convergence? 
                 if(step>1){
-                        VecFree(&du_g);VecDuplicatePatch(&du_g,old_u_g);VecAYPX(du_g,-1.0,u_g);
-                        VecNorm(du_g,NORM_2,&ndu); VecNorm(old_u_g,NORM_2,&nu);
+                        VecDuplicatePatch(&du_g,old_u_g);VecAYPX(du_g,-1.0,u_g);
+                        VecNorm(du_g,NORM_2,&ndu); VecNorm(old_u_g,NORM_2,&nu); VecFree(&du_g);
 
-                        VecFree(&dt_g);VecDuplicatePatch(&dt_g,old_t_g); VecAYPX(dt_g,-1.0,t_g);
-                        VecNorm(dt_g,NORM_2,&normdt); VecNorm(old_t_g,NORM_2,&normt);
-
+                        VecDuplicatePatch(&dt_g,old_t_g); VecAYPX(dt_g,-1.0,t_g);
+                        VecNorm(dt_g,NORM_2,&normdt); VecNorm(old_t_g,NORM_2,&normt);VecFree(&dt_g);
                                         
                         if (debug) _printf_("%-60s%g\n                                     %s%g\n                                     %s%g%s\n",
@@ -135 +134 @@
         }
 
+        /*Free ressource*/
+        VecFree(&old_u_g);
+        VecFree(&old_t_g);
 
         /*Plug results into output dataset: */

issm/trunk/src/c/parallel/thermal_core_nonlinear.cpp

@@ -16 +16 @@
         Vec tg=NULL; 
         Vec tf=NULL; 
+        Vec tf_old=NULL; 
         double melting_offset;
 
@@ -100 +101 @@
                 /*Solve: */
                 if(debug)_printf_("%s\n","solving");
-                Solverx(&tf, Kff, pf, tf, solver_string);
+                Solverx(&tf, Kff, pf,tf_old, solver_string);
+                VecFree(&tf_old); VecDuplicatePatch(&tf_old,tf);
         
                 //no need for Kff and pf anymore
@@ -128 +130 @@
         VecFree(&pg_nopenalty);
         VecFree(&tf);
+        VecFree(&tf_old);
+        delete solver_string;
 
         /*Assign output pointers: */