#include <AnasaziBlockDavidsonSolMgr.hpp>
Inheritance diagram for Anasazi::BlockDavidsonSolMgr< ScalarType, MV, OP >:
Public Member Functions | |
Constructors/Destructor | |
BlockDavidsonSolMgr (const Teuchos::RCP< Eigenproblem< ScalarType, MV, OP > > &problem, Teuchos::ParameterList &pl) | |
Basic constructor for BlockDavidsonSolMgr. | |
virtual | ~BlockDavidsonSolMgr () |
Destructor. | |
Accessor methods | |
const Eigenproblem< ScalarType, MV, OP > & | getProblem () const |
Solver application methods | |
ReturnType | solve () |
This method performs possibly repeated calls to the underlying eigensolver's iterate() routine until the problem has been solved (as decided by the solver manager) or the solver manager decides to quit. |
This solver manager implements a hard-locking mechanism, whereby eigenpairs designated to be locked are moved from the eigensolver and placed in auxiliary storage. The eigensolver is then restarted and continues to iterate, always orthogonal to the locked eigenvectors.
Definition at line 79 of file AnasaziBlockDavidsonSolMgr.hpp.
Anasazi::BlockDavidsonSolMgr< ScalarType, MV, OP >::BlockDavidsonSolMgr | ( | const Teuchos::RCP< Eigenproblem< ScalarType, MV, OP > > & | problem, | |
Teuchos::ParameterList & | pl | |||
) |
Basic constructor for BlockDavidsonSolMgr.
This constructor accepts the Eigenproblem to be solved in addition to a parameter list of options for the solver manager. These options include the following:
string
specifying the desired eigenvalues: SM, LM, SR or LR. Default: "SR"int
specifying the block size to be used by the underlying block Davidson solver. Default: problem->getNEV()int
specifying the number of blocks allocated for the Krylov basis. Default: 2int
specifying the maximum number of restarts the underlying solver is allowed to perform. Default: 20MagnitudeType
specifying the level that residual norms must reach to decide convergence. Default: machine precision.bool
specifying whether residuals norms should be scaled by their eigenvalues for the purposing of deciding convergence. Default: truebool
specifying whether the algorithm should employ locking of converged eigenpairs. Default: falseint
specifying the maximum number of eigenpairs to be locked. Default: problem->getNEV()int
specifying the number of eigenpairs that must meet the locking criteria before locking actually occurs. Default: 1MagnitudeType
specifying the level that residual norms must reach to decide locking. Default: 0.1*convergence tolerancebool
specifying whether residuals norms should be scaled by their eigenvalues for the purposing of deciding locking. Default: true Definition at line 173 of file AnasaziBlockDavidsonSolMgr.hpp.
virtual Anasazi::BlockDavidsonSolMgr< ScalarType, MV, OP >::~BlockDavidsonSolMgr | ( | ) | [inline, virtual] |
ReturnType Anasazi::BlockDavidsonSolMgr< ScalarType, MV, OP >::solve | ( | ) | [virtual] |
This method performs possibly repeated calls to the underlying eigensolver's iterate() routine until the problem has been solved (as decided by the solver manager) or the solver manager decides to quit.
This method calls BlockDavidson::iterate(), which will return either because a specially constructed status test evaluates to Passed or an exception is thrown.
A return from BlockDavidson::iterate() signifies one of the following scenarios:
Implements Anasazi::SolverManager< ScalarType, MV, OP >.
Definition at line 277 of file AnasaziBlockDavidsonSolMgr.hpp.