[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