Anasazi::Experimental::TraceMinDavidsonSolMgr< ScalarType, MV, OP > Class Template Reference

The Anasazi::TraceMinDavidsonSolMgr provides a flexible solver manager over the TraceMinDavidson eigensolver. More...

#include <AnasaziTraceMinDavidsonSolMgr.hpp>

Inheritance diagram for Anasazi::Experimental::TraceMinDavidsonSolMgr< ScalarType, MV, OP >:

List of all members.

Public Member Functions
Constructors
	TraceMinDavidsonSolMgr (const Teuchos::RCP< Eigenproblem< ScalarType, MV, OP > > &problem, Teuchos::ParameterList &pl)
	Basic constructor for TraceMinDavidsonSolMgr.

---

Detailed Description

template<class ScalarType, class MV, class OP>
class Anasazi::Experimental::TraceMinDavidsonSolMgr< ScalarType, MV, OP >

The Anasazi::TraceMinDavidsonSolMgr provides a flexible solver manager over the TraceMinDavidson eigensolver.

This solver manager implements a hard-locking mechanism, whereby eigenpairs designated to be locked are moved from the eigensolver and placed in auxilliary storage. The eigensolver is then restarted and continues to iterate, orthogonal to the locked eigenvectors.

The solver manager provides to the solver a StatusTestCombo object constructed as follows:
   combo = globaltest OR lockingtest OR debugtest
where

• globaltest terminates computation when global convergence has been detected.
  It is encapsulated in a StatusTestWithOrdering object, to ensure that computation is terminated only after the most significant eigenvalues/eigenvectors have met the convergence criteria.
  If not specified via setGlobalStatusTest(), globaltest is a StatusTestResNorm object which tests the 2-norms of the direct residuals relative to the Ritz values.
• lockingtest halts TraceMinDavidson::iterate() in order to deflate converged eigenpairs for locking.
  It will query the underlying TraceMinDavidson eigensolver to determine when eigenvectors should be locked.
  If not specified via setLockingStatusTest(), lockingtest is a StatusTestResNorm object.
• debugtest allows a user to specify additional monitoring of the iteration, encapsulated in a StatusTest object
  If not specified via setDebugStatusTest(), debugtest is ignored.
  In most cases, it should return Failed; if it returns Passed, solve() will throw an AnasaziError exception.

Additionally, the solver manager will terminate and restart solve() when the subspace is full, and it will terminate solve and return after a specified number of restarts.

Much of this behavior is controlled via parameters and options passed to the solver manager. For more information, see TraceMinDavidsonSolMgr().

Author:

Alicia Klinvex

Examples:

TraceMinDavidson/TraceMinDavidsonGeneralizedEx.cpp, TraceMinDavidson/TraceMinDavidsonLaplacianEx.cpp, TraceMinDavidson/TraceMinDavidsonSpecTransEx.cpp, and TraceMinDavidson/TraceMinDavidsonUserOpEx.cpp.

Definition at line 97 of file AnasaziTraceMinDavidsonSolMgr.hpp.

---

Constructor & Destructor Documentation

template<class ScalarType , class MV , class OP >

Anasazi::Experimental::TraceMinDavidsonSolMgr< ScalarType, MV, OP >::TraceMinDavidsonSolMgr	(	const Teuchos::RCP< Eigenproblem< ScalarType, MV, OP > > &	problem,
		Teuchos::ParameterList &	pl
	)

Basic constructor for TraceMinDavidsonSolMgr.

This constructor accepts the Eigenproblem to be solved in addition to a parameter list of options for the solver manager. Since this class inherits from TraceMinBaseSolMgr, it accepts the same options as TraceMinBaseSolMgr(), with a few additions:

• "Block Size" - an int specifying the block size to be used by the underlying solver. If the eigenvalues are clustered, you may want to use a larger block size to capture the correct multiplicity. Default: 1
• "Num Blocks" - an int specifying the maximum number of blocks in the subspace. After we compute this many blocks, we will restart.
• "Num Restart Blocks" - an int specifying how many blocks we keep when restarting. We will extract the most significant blocks from the previous restart and use them to reinitialize TraceMinDavidson. Default: NEV / block size
• "Maximum Restarts" - an int specifying the maximum number of restarts to be performed. Default: 50

A special note about Jacobi-Davidson:
Since Jacobi-Davidson is a special case of TraceMin-Davidson, it does not have its own solver class. If you want to run Jacobi-Davidson, all you need to do is set the following options:

• "When To Shift" = "Always" (which is the default)
• "How To Choose Shift" = "Ritz Values" (not recommended)

Choosing the Ritz shifts to be equal to the Ritz values destroys TraceMin-Davidson's guaranteed global convergence, so it is recommended that you consider TraceMin-Davidson (which is identical in every other way) with its default adjusted Ritz shifts instead.

Definition at line 166 of file AnasaziTraceMinDavidsonSolMgr.hpp.

---

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

• AnasaziTraceMinDavidsonSolMgr.hpp