Corrected important memory leak

diff --git a/CDMATH/linearsolver/src/SparseMatrixPetsc.cxx b/CDMATH/linearsolver/src/SparseMatrixPetsc.cxx
old mode 100644 (file)
new mode 100755 (executable)
index 31a943c..651a744
--- a/CDMATH/linearsolver/src/SparseMatrixPetsc.cxx
+++ b/CDMATH/linearsolver/src/SparseMatrixPetsc.cxx
@@ -340,8 +340,16 @@ SparseMatrixPetsc::operator* (const Vector& vec) const
        MatAssemblyBegin(_mat, MAT_FINAL_ASSEMBLY);
        MatAssemblyEnd(_mat, MAT_FINAL_ASSEMBLY);
        MatMult(_mat,X,Y);
+ 
+       //Clean memory
+       VecDestroy(&X);
 
-       return vecToVector(Y);
+    Vector result=vecToVector(Y);
+
+       //Clean memory
+       VecDestroy(&Y);
+
+       return result;
 }
 
 SparseMatrixPetsc&