Hi,Preallocation of the matrix greatly improved the speed up of the eigensolver.
Is it possible to use the lapack routine 'DSPMV' in the external matrix-vector multiplication for the shellmatrix ? I couldn't find this routine in the petscblaslapack.h or slepcblaslapack.h. Am I missing something?
If there is a slepc/petsc matrix-vector multiplication routine that would also work.
Thanks, Tijo
