message: Adapt FGMRESSolver print behavior
diff: @@ -802,7 +802,7 @@ void FGMRESSolver::Mult(const Vector &b, Vector &x) const
  return;
  }
  - if (print_level>=0)
  + if (print_level == 1)
  mfem::out << " Pass : " << setw(2) << 1
  << " Iteration : " << setw(3) << 0
  << " || r || = " << beta << endl;
  @@ -861,7 +861,7 @@ void FGMRESSolver::Mult(const Vector &b, Vector &x) const
  double resid = fabs(s(i+1));
  MFEM_ASSERT(IsFinite(resid), "resid = " << resid);
  - if (print_level >= 0)
  + if (print_level == 1)
  mfem::out << " Pass : " << setw(2) << (j-1)/m+1
  << " Iteration : " << setw(3) << j
  << " || r || = " << resid << endl;
  @@ -872,6 +872,12 @@ void FGMRESSolver::Mult(const Vector &b, Vector &x) const
  final_norm = resid;
  final_iter = j;
  converged = 1;
  +
  + if (print_level == 2)
  + {
  + mfem::out << "Number of FGMRES iterations: " << final_iter << endl;
  + }
  +
  for (i= 0; i<=m; i++)
  {
  if (v[i]) { delete v[i]; }
  @@ -881,7 +887,7 @@ void FGMRESSolver::Mult(const Vector &b, Vector &x) const
  }
  }
  - if (print_level>=0)
  + if (print_level == 1)
  {
  mfem::out << "Restarting..." << endl;
  }
  @@ -897,6 +903,12 @@ void FGMRESSolver::Mult(const Vector &b, Vector &x) const
  final_norm = beta;
  final_iter = j;
  converged = 1;
  +
  + if (print_level == 2)
  + {
  + mfem::out << "Number of FGMRES iterations: " << final_iter << endl;
  + }
  +
  for (i= 0; i<=m; i++)
  {
  if (v[i]) { delete v[i]; }
  @@ -912,8 +924,13 @@ void FGMRESSolver::Mult(const Vector &b, Vector &x) const
  if (z[i]) { delete z[i]; }
  }
  converged = 0;
  - return;
  + if (print_level >= 0)
  + {
  + mfem::out << "FGMRES: No convergence!" << endl;
  + }
  +
  + return;
  }
  