add {ma27,ma57,mumps}_print_level options

Author: svigerske

ChangeLog.md

@@ -7,6 +7,11 @@ More detailed information about incremental changes can be found in the
 
 ## 3.14
 
+### 3.14.8 (2022-07-13)
+
+- Added options ma27_print_level, ma57_print_level, and mumps_print_level
+  to enable output from these linear solvers.
+
 ### 3.14.7 (2022-06-24)
 
 - Fixed that ComputeSensitivityMatrix() of sIpopt assumed that there are

src/Algorithm/LinearSolvers/IpMa27TSolverInterface.cpp

@@ -84,6 +84,11 @@ void Ma27TSolverInterface::RegisterOptions(
    SmartPtr<RegisteredOptions> roptions
 )
 {
+   roptions->AddBoundedIntegerOption(
+      "ma27_print_level",
+      "Debug printing level for the linear solver MA27",
+      0, 4, 0,
+      "0: no printing; 1: Error messages only; 2: Error and warning messages; 3: Error and warning messages and terse monitoring; 4: All information.");
    roptions->AddBoundedNumberOption(
       "ma27_pivtol",
       "Pivot tolerance for the linear solver MA27.",
@@ -200,6 +205,8 @@ bool Ma27TSolverInterface::InitializeImpl(
       pivtolmax_ = Max(pivtolmax_, pivtol_);
    }
 
+   Index print_level;
+   options.GetIntegerValue("ma27_print_level", print_level, prefix);
    options.GetNumericValue("ma27_liw_init_factor", liw_init_factor_, prefix);
    options.GetNumericValue("ma27_la_init_factor", la_init_factor_, prefix);
    options.GetNumericValue("ma27_meminc_factor", meminc_factor_, prefix);
@@ -210,11 +217,13 @@ bool Ma27TSolverInterface::InitializeImpl(
 
    /* Set the default options for MA27 */
    ma27i(icntl_, cntl_);
-#if IPOPT_VERBOSITY == 0
 
-   icntl_[0] = 0;       // Suppress error messages
-   icntl_[1] = 0;       // Suppress diagnostic messages
-#endif
+   if( print_level == 0 )
+      icntl_[0] = 0;       // Suppress error messages
+   if( print_level <= 1 )
+      icntl_[1] = 0;       // Suppress warning messages
+   if( print_level >= 2 )
+      icntl_[2] = print_level - 2; // diagnostic messages level
 
    // Reset all private data
    initialized_ = false;

src/Algorithm/LinearSolvers/IpMa57TSolverInterface.cpp

@@ -192,6 +192,11 @@ void Ma57TSolverInterface::RegisterOptions(
    SmartPtr<RegisteredOptions> roptions
 )
 {
+   roptions->AddLowerBoundedIntegerOption(
+      "ma57_print_level",
+      "Debug printing level for the linear solver MA57",
+      0, 0,
+      "0: no printing; 1: Error messages only; 2: Error and warning messages; 3: Error and warning messages and terse monitoring; >=4: All information.");
    roptions->AddBoundedNumberOption(
       "ma57_pivtol",
       "Pivot tolerance for the linear solver MA57.",
@@ -321,6 +326,9 @@ bool Ma57TSolverInterface::InitializeImpl(
    DBG_ASSERT(ma57i != NULL);
 
    // Obtain the options settings
+   Index print_level;
+   options.GetIntegerValue("ma57_print_level", print_level, prefix);
+
    options.GetNumericValue("ma57_pivtol", pivtol_, prefix);
    if( options.GetNumericValue("ma57_pivtolmax", pivtolmax_, prefix) )
    {
@@ -361,7 +369,7 @@ bool Ma57TSolverInterface::InitializeImpl(
    wd_icntl_[1] = 0; /* Warning stream. */
 
    wd_icntl_[3] = 1; /* Print statistics.  NOT Used. */
-   wd_icntl_[4] = 0; /* Print error. */
+   wd_icntl_[4] = print_level; /* Print level. */
 
    wd_icntl_[5] = ma57_pivot_order; /* Pivoting order. */
 

src/Algorithm/LinearSolvers/IpMumpsSolverInterface.cpp

@@ -105,9 +105,8 @@ MumpsSolverInterface::MumpsSolverInterface()
 #endif
 
    mumps_c(mumps_);
-   mumps_->icntl[1] = 0;
-   mumps_->icntl[2] = 0; //QUIETLY!
-   mumps_->icntl[3] = 0;
+   mumps_->icntl[2] = 0;  // global info stream
+   mumps_->icntl[3] = 0;  // print level
    mumps_ptr_ = (void*) mumps_;
 }
 
@@ -131,6 +130,11 @@ void MumpsSolverInterface::RegisterOptions(
    SmartPtr<RegisteredOptions> roptions
 )
 {
+   roptions->AddLowerBoundedIntegerOption(
+      "mumps_print_level",
+      "Debug printing level for the linear solver MUMPS",
+      0, 0,
+      "0: no printing; 1: Error messages only; 2: Error, warning, and main statistic messages; 3: Error and warning messages and terse diagnostics; >=4: All information.");
    roptions->AddBoundedNumberOption(
       "mumps_pivtol",
       "Pivot tolerance for the linear solver MUMPS.",
@@ -190,6 +194,9 @@ bool MumpsSolverInterface::InitializeImpl(
    const std::string& prefix
 )
 {
+   Index print_level;
+   options.GetIntegerValue("mumps_print_level", print_level, prefix);
+
    options.GetNumericValue("mumps_pivtol", pivtol_, prefix);
    if( options.GetNumericValue("mumps_pivtolmax", pivtolmax_, prefix) )
    {
@@ -229,6 +236,12 @@ bool MumpsSolverInterface::InitializeImpl(
                        "MumpsSolverInterface called with warm_start_same_structure, but the problem is solved for the first time.");
    }
 
+   if( print_level > 0 )
+   {
+      mumps_->icntl[2] = 6;  // global info stream
+      mumps_->icntl[3] = print_level;  // print level
+   }
+
    return true;
 }
 
@@ -387,10 +400,6 @@ ESymSolverStatus MumpsSolverInterface::SymbolicFactorization()
 
    mumps_data->job = 1;      //symbolic ordering pass
 
-   //mumps_data->icntl[1] = 6;
-   //mumps_data->icntl[2] = 6;//QUIETLY!
-   //mumps_data->icntl[3] = 4;
-
    mumps_data->icntl[5] = mumps_permuting_scaling_;
    mumps_data->icntl[6] = mumps_pivot_order_;
    mumps_data->icntl[7] = mumps_scaling_;