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

Jose E. Roman jroman at dsic.upv.es
Tue Mar 28 12:43:33 CDT 2023

What do you mean that EPSSetInitialSubspace() does not work? Doesn't it improve convergence with respect to not using it?

Is your smallest eigenvalue positive? And do you have many eigenvalues close to the smallest one?


> El 28 mar 2023, a las 17:32, LUTSKO James via petsc-users <petsc-users at mcs.anl.gov> escribió:
> 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

More information about the petsc-users mailing list