Index: /issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp	(revision 1658)
+++ /issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp	(revision 1659)
@@ -143,5 +143,5 @@
 			//compute KUoldF = KUold - F = K*Uold - F
 			VecDuplicate(uf,&KUold); MatMultPatch(Kff,old_uf,KUold);
-			VecDuplicate(KUold,&KUoldF);VecCopy(KUold,KUoldF); VecAYPX(KUoldF,-1.0,pf);
+			VecDuplicate(KUold,&KUoldF);VecCopy(KUold,KUoldF); VecAYPX(KUoldF,-1.0,pf); VecFree(&KUold);
 			//compute norm(KUF), norm(F) and residue
 			VecNorm(KUoldF,NORM_2,&nKUoldF); VecFree(&KUoldF);