Todo:

• use SlepC, which is built on Petsc, as it lets you choose different algorithms (e.g. Block Krylov-Schur is supposedly more efficient than Arnoldi -- not sure about Lanczos: cf http://slepc.upv.es/documentation/reports/str7.pdf )