[Trilinos-Users] [EXTERNAL] Problem on Anassazi~

Thornquist, Heidi K hkthorn at sandia.gov
Wed Jan 2 21:08:19 MST 2013


Hi Sunghwan,

So how many eigenvalues did you request?  Also, there is usually a longer message that accompanies the throw, so could you send all the output that you see from the solver.

Thanks,
Heidi

________________________________
From: trilinos-users-bounces at software.sandia.gov [trilinos-users-bounces at software.sandia.gov] on behalf of SungHwan Choi [sunghwanchoi91 at gmail.com]
Sent: Wednesday, January 02, 2013 5:59 PM
To: trilinos-users at software.sandia.gov
Subject: [EXTERNAL] [Trilinos-Users] Problem on Anassazi~

Hi~,
I have problem with anassazi


Throw test that evaluated to true: ret != 0


When I tried to get eigen values, The code print above message. I think problem is on parameter setting but I don't know what parameters and how do I change them

    Teuchos::RCP< Anasazi::BasicEigenproblem<double,MV,OP> > MyProblem = Teuchos::rcp( new Anasazi::BasicEigenproblem<double,MV,OP>(matrix, ivec) );
    std::string which("SM");

    Teuchos::ParameterList MyPL;
    MyPL.set( "Verbosity", verbose );
    MyPL.set( "Which", which );
    MyPL.set( "Num Blocks", 20 );
    MyPL.set( "Maximum Restarts", 100 );
    MyPL.set( "Convergence Tolerance", 1.0e-8 );
    // Create the Block Krylov Schur solver
    // This takes as inputs the eigenvalue problem and the solver parameters
    Anasazi::BlockKrylovSchurSolMgr<double,MV,OP> MyBlockKrylovSchur(MyProblem, MyPL );
    // Solve the eigenvalue problem, and save the return code
    Anasazi::ReturnType solverreturn = MyBlockKrylovSchur.solve();

The size of matrix is 1000 by 1000

does anyone help me?

Happy new year~
Sunghwan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20130103/800dd77e/attachment.html 


More information about the Trilinos-Users mailing list