29 #ifndef ANASAZI_BASIC_EIGENPROBLEM_H 30 #define ANASAZI_BASIC_EIGENPROBLEM_H 47 template<
class ScalarType,
class MV,
class OP>
244 template <
class ScalarType,
class MV,
class OP>
254 template <
class ScalarType,
class MV,
class OP>
266 template <
class ScalarType,
class MV,
class OP>
280 template <
class ScalarType,
class MV,
class OP>
285 _Prec(Problem._Prec),
286 _InitVec(Problem._InitVec),
288 _isSym(Problem._isSym),
289 _isSet(Problem._isSet),
298 template <
class ScalarType,
class MV,
class OP>
305 if ( !_AOp.get() && !_Op.get() ) {
return false; }
308 if ( !_InitVec.get() ) {
return false; }
311 if (_nev == 0) {
return false; }
314 if (_AOp.get() && !_Op.get()) { _Op = _AOp; }
328 template <
class ScalarType,
class MV,
class OP>
334 "BasicEigenproblem::computeCurrResVec: This method is not currently " 335 "implemented for non-Hermitian problems. Sorry for any inconvenience.");
340 return Teuchos::null;
343 RCP<MV> X = sol.
Evecs;
344 std::vector<ScalarType> Lambda(sol.
numVecs);
345 for(
int i = 0; i < sol.
numVecs; i++)
346 Lambda[i] = sol.
Evals[i].realpart;
349 RCP<MV> AX = MVT::Clone(*X,sol.
numVecs);
350 if(_AOp != Teuchos::null)
351 OPT::Apply(*_AOp,*X,*AX);
353 OPT::Apply(*_Op,*X,*AX);
357 if(_MOp != Teuchos::null)
359 MX = MVT::Clone(*X,sol.
numVecs);
360 OPT::Apply(*_MOp,*X,*MX);
364 MX = MVT::CloneCopy(*X);
368 RCP<MV> R = MVT::Clone(*X,sol.
numVecs);
369 MVT::MvScale(*MX,Lambda);
370 MVT::MvAddMv(1.0,*AX,-1.0,*MX,*R);
bool setProblem()
Specify that this eigenproblem is fully defined.
std::vector< Value< ScalarType > > Evals
The computed eigenvalues.
MultiVecTraits< ScalarType, MV > MVT
Type-definition for the MultiVecTraits class corresponding to the MV type.
void setInitVec(const Teuchos::RCP< MV > &InitVec)
Set the initial guess.
This class defines the interface required by an eigensolver and status test class to compute solution...
bool isHermitian() const
Get the symmetry information for this eigenproblem.
Declaration of basic traits for the multivector type.
OperatorTraits< ScalarType, MV, OP > OPT
Type-definition for the OperatorTraits class corresponding to the OP type.
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
Virtual base class which defines basic traits for the operator type.
Teuchos::RCP< MV > Evecs
The computed eigenvectors.
int getNEV() const
Get the number of eigenvalues (NEV) that are required by this eigenproblem.
Teuchos::RCP< const MV > _AuxVecs
Reference-counted pointer for the auxiliary vector of the eigenproblem .
Teuchos::RCP< const OP > getA() const
Get a pointer to the operator A of the eigenproblem .
bool _isSym
Symmetry of the eigenvalue problem.
int _nev
Number of eigenvalues requested.
Teuchos::RCP< MV > _InitVec
Reference-counted pointer for the initial vector of the eigenproblem .
Teuchos::RCP< const OP > _Op
Reference-counted pointer for the operator of the eigenproblem .
Namespace Anasazi contains the classes, structs, enums and utilities used by the Anasazi package...
void setAuxVecs(const Teuchos::RCP< const MV > &AuxVecs)
Set auxiliary vectors.
int numVecs
The number of computed eigenpairs.
Teuchos::RCP< const OP > getOperator() const
Get a pointer to the operator for which eigenvalues will be computed.
Teuchos::RCP< const OP > getM() const
Get a pointer to the operator M of the eigenproblem .
Teuchos::RCP< const MV > computeCurrResVec() const
Returns the residual vector corresponding to the computed solution.
Abstract base class which defines the interface required by an eigensolver and status test class to c...
Traits class which defines basic operations on multivectors.
Virtual base class which defines basic traits for the operator type.
const Eigensolution< ScalarType, MV > & getSolution() const
Get the solution to the eigenproblem.
void setSolution(const Eigensolution< ScalarType, MV > &sol)
Set the solution to the eigenproblem.
virtual ~BasicEigenproblem()
Destructor.
Teuchos::RCP< const OP > _MOp
Reference-counted pointer for M of the eigenproblem .
Teuchos::RCP< const MV > getInitVec() const
Get a pointer to the initial vector.
Struct for storing an eigenproblem solution.
bool _isSet
Sanity Check Flag.
void setPrec(const Teuchos::RCP< const OP > &Prec)
Set the preconditioner for this eigenvalue problem .
bool isProblemSet() const
If the problem has been set, this method will return true.
void setHermitian(bool isSym)
Specify the symmetry of this eigenproblem.
This provides a basic implementation for defining standard or generalized eigenvalue problems...
Teuchos::RCP< const MV > getAuxVecs() const
Get a pointer to the auxiliary vector.
Teuchos::RCP< const OP > getPrec() const
Get a pointer to the preconditioner of the eigenproblem .
Teuchos::RCP< const OP > _Prec
Reference-counted pointer for the preconditioner of the eigenproblem .
void setM(const Teuchos::RCP< const OP > &M)
Set the operator M of the eigenvalue problem .
BasicEigenproblem()
Empty constructor - allows Anasazi::BasicEigenproblem to be described at a later time through "Set Me...
void setA(const Teuchos::RCP< const OP > &A)
Set the operator A of the eigenvalue problem .
Teuchos::RCP< const OP > _AOp
Reference-counted pointer for A of the eigenproblem .
void setNEV(int nev)
Specify the number of eigenvalues (NEV) that are requested.
Eigensolution< ScalarType, MV > _sol
Solution to problem.
void setOperator(const Teuchos::RCP< const OP > &Op)
Set the operator for which eigenvalues will be computed.