LOCA::Eigensolver::AnasaziStrategy Class Reference

Anasazi eigensolver strategy. More...

#include <LOCA_Eigensolver_AnasaziStrategy.H>

Inheritance diagram for LOCA::Eigensolver::AnasaziStrategy:

[legend]Collaboration diagram for LOCA::Eigensolver::AnasaziStrategy:[legend]List of all members.

Public Member Functions
	AnasaziStrategy (const Teuchos::RefCountPtr< LOCA::GlobalData > &global_data, const Teuchos::RefCountPtr< LOCA::Parameter::SublistParser > &topParams, const Teuchos::RefCountPtr< NOX::Parameter::List > &eigenParams)
	Constructor.
virtual	~AnasaziStrategy ()
	Destructor.
virtual NOX::Abstract::Group::ReturnType	computeEigenvalues (NOX::Abstract::Group &group, Teuchos::RefCountPtr< std::vector< double > > &evals_r, Teuchos::RefCountPtr< std::vector< double > > &evals_i, Teuchos::RefCountPtr< NOX::Abstract::MultiVector > &evecs_r, Teuchos::RefCountPtr< NOX::Abstract::MultiVector > &evecs_i)
	Compute eigenvalues/eigenvectors.
Protected Attributes
Teuchos::RefCountPtr< LOCA::GlobalData >	globalData
	Global data.

---

Detailed Description

Anasazi eigensolver strategy.

This class implements an eigensolver strategy using the generic Trilinos eigensolver package Anasazi. Since Anasazi is a generic solver, this strategy will work with any group implementation. This strategy references the following parameters passed through the eigenParams argument to the constructor:

• "Operator" -- [string] (default: "Jacobian inverse") Operator to compute eigenvalues of.
• "Block Size" -- [int] (default: 1) Block size
• "Arnoldi Size" -- [int] (default: 30) Maximum length of the Arnoldi factorization
• "NEV" -- [int] (default: 4) Number of requested eigenvalues
• "Tol" -- [double] (default: 1.0e-7) Tolerance for the converged eigenvalues
• "Convergence Check" -- [int] (default: 1) Checks convergence every so many steps
• "Restarts" -- [int] (default: 1) Number of restarts allowed
• "Debug Level" -- [int] (default: 1) Anasazi Debug level
• "Sorting Order" -- [string" (default: "LM") Sorting order of printed eigenvalues
• "Save Eigenvectors" -- [int] (default: 0) How many eigenvectors to save
• "Cayley Pole" -- [double] (default: 0.0) Cayley transformation pole
• "Cayley Zero" -- [double] (default: 0.0) Cayley tranformation zero

---

Constructor & Destructor Documentation

LOCA::Eigensolver::AnasaziStrategy::AnasaziStrategy (  const Teuchos::RefCountPtr< LOCA::GlobalData > &  global_data, const Teuchos::RefCountPtr< LOCA::Parameter::SublistParser > &  topParams, const Teuchos::RefCountPtr< NOX::Parameter::List > &  eigenParams )

Constructor. Parameters: global_data  [in] Global data object topParams  [in] Parsed top-level parameter list. Used to obtain linear-solver parameters and to instantiate sorting strategy. eigenParams  [in] Eigensolver parameters as described above. solver.

---

Member Function Documentation

NOX::Abstract::Group::ReturnType LOCA::Eigensolver::AnasaziStrategy::computeEigenvalues (  NOX::Abstract::Group &  group, Teuchos::RefCountPtr< std::vector< double > > &  evals_r, Teuchos::RefCountPtr< std::vector< double > > &  evals_i, Teuchos::RefCountPtr< NOX::Abstract::MultiVector > &  evecs_r, Teuchos::RefCountPtr< NOX::Abstract::MultiVector > &  evecs_i )  [virtual]

Compute eigenvalues/eigenvectors. The implementation here the sets up and calls the Anasazi BlockKrylovSchur solver for computing eigenvalues. Implements LOCA::Eigensolver::AbstractStrategy.

---

The documentation for this class was generated from the following files:

• LOCA_Eigensolver_AnasaziStrategy.H
• LOCA_Eigensolver_AnasaziStrategy.C

---