Adding MatDiagonalSet, MatAXPY and MatAYPX

Not sure why, but setting only […]

The st_ prefix is added when the KSP is created by SLEPc. If you create it yourself with set(...), the option should not be prefixed, unless you add the flag set(..., prefix = "st_").

[…] however, I prefer SLEPc to handle all that.

Would it be possible to ask you for a piece of code that calls PEPSolve() and where I could try the -ksp_reuse_preconditioner option? You can share that in private if it’s best for you.

[…] nonzero vectors, those are used as a guess […]

I’ve never seen tremendous improvements using this feature.

How does this work - I am struggling to figure this out on my own.

Please see Slepc ESPSetDeflationSpace function andd let me know if you need more details.