print PHASE and adding log for pardiso-mkl solve

coin-or · Dec 9, 2024 · 264b1c7 · 264b1c7
1 parent db0e648
commit 264b1c7
Showing 1 changed file with 5 additions and 2 deletions.
diff --git a/src/Algorithm/LinearSolvers/IpPardisoMKLSolverInterface.cpp b/src/Algorithm/LinearSolvers/IpPardisoMKLSolverInterface.cpp
@@ -471,7 +471,7 @@ ESymSolverStatus PardisoMKLSolverInterface::Factorization(
          PHASE = 11;
 
          Jnlst().Printf(J_DETAILED, J_LINEAR_ALGEBRA,
-                        "Calling Pardiso for symbolic factorization.\n");
+                        "Calling Pardiso for symbolic factorization (PHASE=%d).\n", PHASE);
          IPOPT_LAPACK_FUNC(pardiso, PARDISO)(PT_, &MAXFCT_, &MNUM_, &MTYPE_,
                                              &PHASE, &N, a_, ia, ja, &PERM,
                                              &NRHS, IPARM_, &MSGLVL_, &B, &X, &ERROR, DPARM_);
@@ -509,7 +509,7 @@ ESymSolverStatus PardisoMKLSolverInterface::Factorization(
          IpData().TimingStats().LinearSystemFactorization().Start();
       }
       Jnlst().Printf(J_MOREDETAILED, J_LINEAR_ALGEBRA,
-                     "Calling Pardiso for factorization.\n");
+                     "Calling Pardiso for factorization (PHASE=%d).\n", PHASE);
       // Dump matrix to file, and count number of solution steps.
       if( HaveIpData() )
       {
@@ -642,6 +642,9 @@ ESymSolverStatus PardisoMKLSolverInterface::Solve(
    Index NRHS = nrhs;
    Number* X = new Number[nrhs * dim_];
 
+   Jnlst().Printf(J_MOREDETAILED, J_LINEAR_ALGEBRA,
+                  "Calling Pardiso to solve (PHASE=%d).\n", PHASE);
+
    Number* ORIG_RHS = new Number[nrhs * dim_];
    Index ERROR;
    // Initialize solution with zero and save right hand side