[numerics] release memory in gfc3d IPM method
@@ -1119,7 +1119,10 @@ void gfc3d_IPM(GlobalFrictionContactProblem* restrict problem, double* restrict
/* printIteresProbMatlabFile(iteration, globalVelocity, velocity, reaction, d, n, m, iterates); */
hasNotConverged = 0;
-
+ if (Qp)
+ Qp = NM_free(Qp);
+ if (Qpinv)
+ Qpinv = NM_free(Qpinv);
// numerics_printf_verbose(-1, "%9.2e %9.2e %9.2e\n", norm2VecDiff(tmpsol, globalVelocity, m), norm2VecDiff(tmpsol+m, velocity, nd), norm2VecDiff(tmpsol+m+nd, reaction, nd));
break;
}
