LOCA::Eigensolver::AnasaziStrategy Class Reference

Anasazi eigensolver strategy. More...

#include <LOCA_Eigensolver_AnasaziStrategy.H>

Inheritance diagram for LOCA::Eigensolver::AnasaziStrategy:

Collaboration diagram for LOCA::Eigensolver::AnasaziStrategy:

List of all members.

Public Member Functions
	AnasaziStrategy (const Teuchos::RCP< LOCA::GlobalData > &global_data, const Teuchos::RCP< LOCA::Parameter::SublistParser > &topParams, const Teuchos::RCP< Teuchos::ParameterList > &eigenParams)
	Constructor.
virtual	~AnasaziStrategy ()
	Destructor.
virtual NOX::Abstract::Group::ReturnType	computeEigenvalues (NOX::Abstract::Group &group, Teuchos::RCP< std::vector< double > > &evals_r, Teuchos::RCP< std::vector< double > > &evals_i, Teuchos::RCP< NOX::Abstract::MultiVector > &evecs_r, Teuchos::RCP< NOX::Abstract::MultiVector > &evecs_i)
	Compute eigenvalues/eigenvectors.
Protected Types
typedef NOX::Abstract::MultiVector	MV
	Typedef to simplify templating.
typedef LOCA::AnasaziOperator::AbstractStrategy	OP
	Typedef to simplify templating.
Protected Attributes
Teuchos::RCP< LOCA::GlobalData >	globalData
	Global data.
Teuchos::RCP < LOCA::Parameter::SublistParser >	topParams
	Parsed parameter list.
Teuchos::RCP < Teuchos::ParameterList >	eigenParams
	Eigensolver parameter list.
Teuchos::RCP < Teuchos::ParameterList >	solverParams
	Linear solver parameter list.
int	blksz
	The block size.
int	nev
	The number of requested eigenvalues.
bool	isSymmetric
	Is problem symmetric.

---

Detailed Description

Anasazi eigensolver strategy.

This class implements an eigensolver strategy using the generic Trilinos eigensolver package Anasazi. In particular, this strategy uses the Anasazi::BlockKrylovSchur solver. 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 (this list is passed directly to the Anasazi::BlockKrylovSchulSolMgr solver manager):

• "Operator" -- [string] (default: "Jacobian inverse") Operator to compute eigenvalues of.
• "Block Size" -- [int] (default: 1) Block size
• "Num Blocks" -- [int] (default: 30) Maximum number of blocks (equals the maximum length of the Arnoldi factorization
• "Num Eigenvalues" -- [int] (default: 4) Number of requested eigenvalues
• "Convergence Tolerance" -- [double] (default: 1.0e-7) Tolerance for the converged eigenvalues
• "Step Size" -- [int] (default: 1) Checks convergence every so many steps
• "Maximum Restarts" -- [int] (default: 1) Number of restarts allowed
• "Symmetric" -- [bool] (default: false) Is the operator symmetric
• "Verbosity" -- [Anasazi::MsgType] (default: Anasazi::Errors+Anasazi::Warnings+Anasazi::FinalSummary) Verbosity level
• "Sorting Order" -- [string" (default: "LM") Sorting order of printed eigenvalues

---

Constructor & Destructor Documentation

LOCA::Eigensolver::AnasaziStrategy::AnasaziStrategy	(	const Teuchos::RCP< LOCA::GlobalData > &	global_data,
		const Teuchos::RCP< LOCA::Parameter::SublistParser > &	topParams,
		const Teuchos::RCP< Teuchos::ParameterList > &	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.

References blksz, eigenParams, Teuchos::ParameterList::get(), LOCA::Parameter::SublistParser::getSublist(), Teuchos::ParameterList::isParameter(), isSymmetric, nev, Teuchos::rcp(), Teuchos::ParameterList::set(), solverParams, and topParams.

---

Member Function Documentation

NOX::Abstract::Group::ReturnType LOCA::Eigensolver::AnasaziStrategy::computeEigenvalues	(	NOX::Abstract::Group &	group,
		Teuchos::RCP< std::vector< double > > &	evals_r,
		Teuchos::RCP< std::vector< double > > &	evals_i,
		Teuchos::RCP< NOX::Abstract::MultiVector > &	evecs_r,
		Teuchos::RCP< 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.

References LOCA::AnasaziOperator::AbstractStrategy::beginPostProcessing(), NOX::Abstract::Vector::createMultiVector(), NOX::Abstract::Group::getX(), LOCA::AnasaziOperator::AbstractStrategy::label(), NOX::Abstract::Group::NotConverged, NOX::Abstract::Group::Ok, LOCA::AnasaziOperator::AbstractStrategy::preProcessSeedVector(), LOCA::AnasaziOperator::AbstractStrategy::rayleighQuotient(), Teuchos::rcp(), NOX::Utils::StepperIteration, and LOCA::AnasaziOperator::AbstractStrategy::transformEigenvalue().

---

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

• LOCA_Eigensolver_AnasaziStrategy.H
• LOCA_Eigensolver_AnasaziStrategy.C