[Trilinos-Users] ansazi question

[Alicia Klinvex]

Hello David,

There is a way to retrieve the eigenpairs if convergence fails, but it's
not easy.  The issue is that you use the LOBPCGSolMgr, I assume, which does
not give you access to the solver.  The solver has a function to return the
Ritz values and vectors; the solver manager does not.  You would need to
create a custom debug status test (which allows you to grab the solver) and
provide it to the solver manager.  This is not something that we generally
recommend doing, but if you desperately want to do so, I can give you more
information about how to do this.

Best wishes,
Alicia

On Fri, Feb 13, 2015 at 6:14 PM, Hysom, David A. <hysom1 at llnl.gov> wrote:

>  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
>
> _______________________________________________
> Trilinos-Users mailing list
> Trilinos-Users at software.sandia.gov
> https://software.sandia.gov/mailman/listinfo/trilinos-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://software.sandia.gov/pipermail/trilinos-users/attachments/20150213/6c36dc5e/attachment.html>