[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