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. Jose
