[kernel] Multiple calls to NM_LU_solve with NM_preserve raises issues
@@ -55,10 +55,12 @@ int Equality::compute(double time)
double* z_ = _z->getArray();
for(size_t i = 0; i < _sizeOutput; ++i) z_[i] = -q_[i];
//info = NM_gesv(&*_M->numericsMatrix(), z_, true);
- info = NM_LU_solve(NM_preserve(&*_M->numericsMatrix()), z_, 1);
+ //info = NM_LU_solve(NM_preserve(&*_M->numericsMatrix()), z_, 1);
+ info = NM_LU_solve(&*_M->numericsMatrix(), z_, 1);
+
+
// --- Recovering of the desired variables from EQUALITY output ---
postCompute();
-
}
return info;
