[Trilinos-Users] Anasazi problem

[SungHwan Choi]

I have question;;;
I think I have no problem and actually I used this codes before,,,,, after
I used Teuchos::Array It doesn't work;;;  please let me know the solution
here is pseudo code;



Teuchos::Array<Teuchos::RCP <Epetra_CrsMatrix > > crs_array ;
for (int i =0;i<2;i++){

       crs_array.push_back(Teuchos::rcp(new Epetra_CrsMatrix(copy,map,0)) );
}

functionA(crs_array[0]);



int functionA(Teuchos::RCP<Epetra_CrsMatrix> matrix )
    cout << *matrix << endl;   //Epetra_CrsMatrix is printed well

    Teuchos::RCP < Anasazi::BasicEigenproblem<double,MV,OP> > MyProblem
=Teuchos::rcp( new Anasazi::BasicEigenproblem<double,MV,OP>() );
    MyProblem->setNEV( 5);
    MyProblem->setHermitian(true);
    MyProblem->setOperator(matrix);
    MyProblem->setInitVec(initial_eigenvector);
    MyProblem->setProblem();

    Teuchos::ParameterList MyPL;
    MyPL.set( "Verbosity", verb );
    MyPL.set( "Which", "SR" );
    MyPL.set( "Block Size", 20);
    MyPL.set( "Num Blocks", 24);


    Anasazi::BlockKrylovSchurSolMgr<double,MV,OP>
MyBlockKrylovSchur(MyProblem, MyPL );
    cout <<"111111111111111" <<endl;  //it is not printed
    Anasazi::ReturnType solverreturn = MyBlockKrylovSchur.solve();

terminate called after throwing an instance of 'std::invalid_argument'
  what():
/appl/trilinos/libs/include/AnasaziBlockKrylovSchurSolMgr.hpp:266:
Throw number = 1
Throw test that evaluated to true: !_problem->isProblemSet()


I already spent more than a day';;;  please what is the problem and the
solution
Sunghwan Choi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20130321/f21cb96b/attachment.html