[numerics] add mumps verbosity if numerics verbose is set
@@ -3106,6 +3106,17 @@ int NM_gesv_expert(NumericsMatrix* A, double *b, unsigned keep)
/* the mumps instance is initialized (call with job=-1) */
NM_MUMPS_set_control_params(A);
NM_MUMPS(A, -1);
+ if ((NM_MUMPS_icntl(A, 1) == -1 ||
+ NM_MUMPS_icntl(A, 2) == -1 ||
+ NM_MUMPS_icntl(A, 3) == -1 ||
+ verbose) ||
+ (NM_MUMPS_icntl(A, 1) != -1 ||
+ NM_MUMPS_icntl(A, 2) != -1 ||
+ NM_MUMPS_icntl(A, 3) != -1 ||
+ !verbose))
+ {
+ NM_MUMPS_set_verbosity(A, verbose);
+ }
NM_MUMPS_set_icntl(A, 24, 1); // Null pivot row detection
NM_MUMPS_set_cntl(A, 5, 1.e20); // Fixation, recommended value
}
@@ -3476,6 +3487,18 @@ int NM_posv_expert(NumericsMatrix* A, double *b, unsigned keep)
NM_MUMPS_set_control_params(A);
NM_MUMPS_set_sym(A, 2); /* general symmetric */
NM_MUMPS(A, -1);
+ if ((NM_MUMPS_icntl(A, 1) == -1 ||
+ NM_MUMPS_icntl(A, 2) == -1 ||
+ NM_MUMPS_icntl(A, 3) == -1 ||
+ verbose) ||
+ (NM_MUMPS_icntl(A, 1) != -1 ||
+ NM_MUMPS_icntl(A, 2) != -1 ||
+ NM_MUMPS_icntl(A, 3) != -1 ||
+ !verbose))
+ {
+ NM_MUMPS_set_verbosity(A, verbose);
+ }
+
NM_MUMPS_set_icntl(A, 24, 1); // Null pivot row detection
NM_MUMPS_set_cntl(A, 5, 1.e20); // Fixation, recommended value
}
