[numerics] gfc3d_admm: set mumps solver if WITH_MUMPS defined
@@ -608,11 +608,15 @@ void gfc3d_ADMM(GlobalFrictionContactProblem* restrict problem, double* restrict
{
NM_gesv_expert(W,v,NM_KEEP_FACTORS);
}
+
else
{
- /* NM_gesv_expert(W,v,NM_KEEP_FACTORS); */
NSM_linear_solver_params* p = NSM_linearSolverParams(W);
+#ifdef WITH_MUMPS
+ p->solver = NSM_MUMPS;
+#else
p->solver = NSM_CS_CHOLSOL;
+#endif
NM_posv_expert(W,v,NM_KEEP_FACTORS);
}
