[Trilinos-Users] troubles on setting anasazi problem

Pavel Jiránek pavel.jiranek at gmail.com
Thu Dec 27 08:04:16 MST 2012


Hello SungHwan,

You need to set the number of requested eigenvalues using the method
setNEV() before calling setProblem().

Cheers,
Pavel Jiranek



2012/12/27 SungHwan Choi <sunghwanchoi91 at gmail.com>

> Hello,
> I have little trouble with anasazi. reffering example in didasko, I have
> written my own code but I have troubles. I think this is really easy to
> someone who are familiar to anasazi and Teuchos package. unfortunatly not
> me. Do you know how to solve the problem?
>
> here is my code, kinetic_matrix is defined as Epetra_CrsMatrix*
> kinetic_matrix
>
>
>
>
>
>     Teuchos::RCP<Epetra_CrsMatrix> matrix = Teuchos::rcp( kinetic_matrix ,
> false);
>     // Create an Anasazi MultiVector, based on Epetra MultiVector
>     const Epetra_Map * Map = &(kinetic_matrix->RowMap());
>     Teuchos::RCP<Epetra_MultiVector> ivec = Teuchos::rcp( new
> Epetra_MultiVector(*Map,1) );
>     //Teuchos::RCP<Epetra_MultiVector> ivec = Teuchos::rcp(  );
>     // Fill it with random numbers
>     ivec->Random();
>     // Setup the eigenproblem, with the matrix A and the initial vectors
> ivec
>     Teuchos::RCP< Anasazi::BasicEigenproblem<double,MV,OP> > MyProblem =
> Teuchos::rcp( new Anasazi::BasicEigenproblem<double,MV,OP>(matrix, ivec) );
>     MyProblem->setHermitian(true);
>     bool ierr = MyProblem->setProblem();
>     cout << ierr << endl;
>
> results is 0 which means that setting problem has problem. Do you have any
> hints to my trouble?
>
> Happy new year~
> _______________________________________________
> Trilinos-Users mailing list
> Trilinos-Users at software.sandia.gov
> http://software.sandia.gov/mailman/listinfo/trilinos-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20121227/b7c6dcf4/attachment.html 


More information about the Trilinos-Users mailing list