[numerics] reduce verbose when solver fails
@@ -81,7 +81,7 @@ void compute_racines(double * Poly, int *nbRealRacines, double *Racines)
}
else
{
- printf("FC3D_unitary_enumerative: degre of polynom is 0.");
+ if (verbose) printf("FC3D_unitary_enumerative: degre of polynom is 0.");
degp1 = 1;
}
(*nbRealRacines) = 0;
@@ -244,9 +244,13 @@ int fc3d_unitary_enumerative_test_non_sliding(FrictionContactProblem* problem, d
}
/*sticking ? 0=MR+q*/
- int info = solv3x3(M, reaction, Q);
- if(info)
- numerics_warning("fc3d_unitary_enumerative_test_non_sliding", "NaN output in solv3x3");
+ //int info = solv3x3(M, reaction, Q);
+ reaction[0] = Q[0];
+ reaction[1] = Q[1];
+ reaction[2] = Q[2];
+ int info = solve_3x3_gepp(M, reaction);
+ if(info && (verbose > 0))
+ numerics_warning("fc3d_unitary_enumerative_test_non_sliding", "NaN output in solve_3x3_gepp");
M = M00;
reaction = reaction0;
Q = Q0;
@@ -749,11 +753,11 @@ int fc3d_unitary_enumerative_solve_poly_nu_sliding(FrictionContactProblem* probl
#ifdef FC3D_UE_DEBUG
#endif
- printf("fc3d_unitary_enumerative_poly_nu (FAILED) M:\n");
- print3x3(M);
+ if (verbose) printf("fc3d_unitary_enumerative_poly_nu (FAILED) M:\n");
+ if (verbose) print3x3(M);
M = M00;
- printf("fc3d_unitary_enumerative_poly_nu (FAILED) Q:\n");
- print3(Q);
+ if (verbose) printf("fc3d_unitary_enumerative_poly_nu (FAILED) Q:\n");
+ if (verbose) print3(Q);
Q = Q0;
return -1;
}
