Set PETSc Solver Package for SLEPcEigenSolver object in python

Question

I would like to set the Matrix Solver Package to mumps for a SLEPcEigenSolver object, where the spectal transformation is set to "shift and invert". In the PETSc/SLEPc C-API this can be set with the following code

ST myst;
EPSGetST(eps_,&myst);
STSetType(myst,STSINVERT);

KSP myksp;
STGetKSP(myst, &myksp);
KSPSetType(myksp, KSPPREONLY);

PC mypc;
KSPGetPC(myksp,&mypc);
PCSetType(mypc,PCLU);
PCFactorSetMatSolverPackage(mypc,MATSOLVERMUMPS);

How can this be set from python?

Garth N. Wells · Answer 1 · Dec 23, 2014

Best answer

The DOLFIN development version now includes support for slepc4py, so I would suggest using slepc4py to get fine-grained control from Python.

answered Dec 23, 2014 by Garth N. Wells FEniCS Expert (35,930 points)
selected Dec 23, 2014 by thisch

So all i need to do is call the eps method, which returns a slepc4py object.

commented Dec 23, 2014 by thisch FEniCS Novice (580 points)

Yes. I haven't tested - it was implemented by Jan Blechta.

commented Dec 23, 2014 by Garth N. Wells FEniCS Expert (35,930 points)

This code works:

eps = eigensolver.eps()
st = eps.getST()
st.setType('sinvert')
ksp = st.getKSP()
ksp.setType('preonly')
pc = ksp.getPC()
pc.setType('lu')
pc.setFactorSolverPackage('mumps')

commented Dec 23, 2014 by thisch FEniCS Novice (580 points)

Set PETSc Solver Package for SLEPcEigenSolver object in python

1 Answer