[petsc-users] request suggestions for most appropriate eigenvalue solver

Jose E. Roman jroman at dsic.upv.es
Wed Jul 18 10:28:22 CDT 2012

El 17/07/2012, a las 19:56, Giacomo Mulas escribió:

> Hello Jose.
> Some months ago, I wrote to the petsc-users mailing list asking if it would
> be possible to use an iterative solver in slepc which converges eigenvectors
> in order of maximum projection on a given basis vector.  Back then you told
> me you might look into it in some time and let me know.  Since some time
> passed, and you might quite understandably not remember about this, I am
> quoting that email exchange below.
> In the meanwhile, I found that there the Davidson algorithm and its derivatives (e.g. block Davidson) appears to behave as I would like, or close to it.
> You probably know them, but I send you a couple of references in any case:
> E. Davidson, J. Comput. Phys. 17, 87 ͑1975͒.
> E. Davidson, Comput. Phys. Commun. 53, 49 ͑1989͒.
> F. Ribeiro, C. Iung, and C. Leforestier, Chem. Phys. Lett. 362, 199
> ͑2002͒.
> F. Ribeiro, C. Iung, and C. Leforestier, J. Theor. Comput. Chem. 2, 609
> ͑2003͒.
> Any hope of implementing something like that (or some other algorithm to
> obtain the same behaviour) in Slepc/Petsc any time soon?
> Thanks in advance, bye
> Giacomo Mulas

Yes, we will try to have this included in the release. I will come back to you in a week or so.

More information about the petsc-users mailing list