Setting SLEPcEigenSolver.parameters["verbose"] = True changes the solver type

Issue #265 resolved
Patrick Farrell created an issue

In SLEPcEigenSolver.cpp:144:

  if (parameters["verbose"])
  {
    KSP ksp;
    ST st;
    EPSMonitorSet(eps, EPSMonitorAll, PETSC_NULL, PETSC_NULL);
    EPSSetType(eps, EPSARPACK);
    EPSGetST(eps, &st);
    STGetKSP(st, &ksp);
    KSPMonitorSet(ksp, KSPMonitorDefault, PETSC_NULL, PETSC_NULL);
    EPSView(eps, PETSC_VIEWER_STDOUT_SELF);
  }

This is unexpected -- my version of SLEPc isn't compiled against ARPACK and was stalling with no error message or explanation because the EPSSetType(eps, EPSARPACK) call was failing.

Comments (5)

  1. Log in to comment