[Trilinos-Users] Printing Eigenvalues with LOCA
Casey Marks
casey at fivehut.com
Mon Feb 2 19:30:22 MST 2009
I'm updating some code that was written against an old (and no longer
available) version of Trilinos/LOCA. I've gotten it working with the
dev branch of Trilinos, but I have lost one important feature in the
conversion -- printing eigenvalues. I added my own
LOCA::SaveEigenData::AbstractStrategy, but its save() method is never
called. Any suggestions on how to get eigenvalues printed out on every
step?
Here is the code where I set the eigenvalue parameters:
locaStepperList.set("Compute Eigenvalues", true);
Teuchos::ParameterList& eList =
locaStepperList.sublist("Eigensolver");
eList.set("Method","Anasazi");
eList.set("Block Size", 1); // Size of blocks
eList.set("Num Blocks", 30); // Size of Arnoldi factorization
eList.set("Num Eigenvalues", 4); // Number of eigenvalues
eList.set("Convergence Tolerance", 2.0e-5); // Tolerance
eList.set("Step Size", 1); // How often to check convergence
eList.set("Maximum Restarts", 2);
eList.set("Sorting Order", "SM");
eList.set("Verbosity",
Anasazi::Errors +
Anasazi::Warnings +
Anasazi::FinalSummary); // Verbosity
eList.set("Save Eigen Data Method", "User-Defined");
eList.set("User-Defined Save Eigen Data Name", "CustomEigenOp");
Teuchos::RCP<LOCA::SaveEigenData::AbstractStrategy> sedm =
Teuchos::rcp(new EigenOp);
eList.set("CustomEigenOp", sedm);
Thank you for your time,
Casey Marks
More information about the Trilinos-Users
mailing list