[Trilinos-Users] ml with Anasazi

SungHwan Choi sunghwanchoi91 at gmail.com
Thu Jul 17 07:15:35 MDT 2014


Dear all,
Hi, I am an Anasazi user. I want to use a multilevel preconditioner. As far
as, I know AztecOO work well with ML package. However, I am not sure
multigrid like preconditioner can be applied on Anasazi.
The article titled "*A comparison of eigensolvers for large-scale 3D modal
analysis using AMG-preconditioned iterative methods*" said that ML package
can work with Anasazi solver but I couldn't find how ML can work with
Anasazi.

Thus, applying my pre-knowledge, I coded it but It does not work well.
Would you help me?


            Teuchos::RCP < Anasazi::BasicEigenproblem<double,MV,OP> >
MyProblem =Teuchos::rcp( new
Anasazi::BasicEigenproblem<double,MV,OP>(matrix, eigenvectors) );
            Teuchos::RCP<Epetra_Operator> PrecOp;
            Teuchos::ParameterList MLList;
            ML_Epetra::SetDefaults("DD",MLList);
            MLList.set("smoother: pre or post", "post");
            MLList.set("PDE equations", 1);
            MLList.set("smoother: type","IFPACK");
            MLList.set("smoother: ifpack type", "ILU");
            MLList.set("smoother: ifpack overlap", 1);
            MLList.sublist("smoother: ifpack list").set("fact:
level-of-fill", 5);
            ML_Epetra::MultiLevelPreconditioner* MLPrec =new
ML_Epetra::MultiLevelPreconditioner(*matrix.get(), MLList, true);
            PrecOp=Teuchos::rcp( MLPrec);
            MyProblem->setPrec(PrecOp);


Sunghwan Choi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://software.sandia.gov/pipermail/trilinos-users/attachments/20140717/b4dcd37b/attachment.html>


More information about the Trilinos-Users mailing list