[Trilinos-Users] ansazi question

[Hysom, David A.]

Hi,

We're using anasazi's LOBPCG, and have cases where convergence
fails, but we'd still like to retrieve the 2nd smallest eigenvector (we're
computing the smallesBasicEigenproblemt two). When we call:

Anasazi::Eigensolution<scalar_t, mvec_t> sol = problem->getSolution();
Teuchos::RCP<mvec_t> evecs = sol.Evecs;
cout << "evecs.size: " << evecs->getNumVectors() << endl;

we get a single vector. However, since anasazi can print out the
eigenvalues and residuals after every iteration, I'm guessing there
should be a way to retrieve the eigenpairs if convergence fails.

So, is there? If not, how easy would it be to modify the code to
do so, and where should I start?

thanks, David
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://software.sandia.gov/pipermail/trilinos-users/attachments/20150213/9b6ed93d/attachment.html>