[petsc-users] Restarting a SLEPC run to refine an eigenvector

[LUTSKO James]

Hello,
 I am a complete newbe so sorry if this is answered elsewhere or ill-posed but .. I am trying to get the smallest eigenvector of a  large matrix. Since the calculations are very long, I would like to be able to restart my code if I find that the "converged" eigenvector is not good enough without having to start from scratch. I am currently using the options

-eps_type jd  -eps_monitor -eps_smallest_real -eps_conv_abs -eps_tol 1e-8 -st_ksp_type gmres -st_pc_type jacobi -st_ksp_max_it 40

Is there any way to do this? I have tried using EPSSetInitialSubspace with the stored eigenvector but this does not seem to work.

thanks
jim

% James F. Lutsko, CNLPCS, Universite Libre de Bruxelles,
% Campus Plaine -- CP231  B-1050 Bruxelles, Belgium
% tel: +32-2-650-5997   email: jlutsko at ulb.ac.be
% fax: +32-2-650-5767
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230328/b1f106e8/attachment.html>