#include <BelosOrthoManager.hpp>
Inheritance diagram for Belos::OrthoManager< ScalarType, MV >:
Public Member Functions | |
Constructor/Destructor | |
OrthoManager () | |
Default constructor. | |
virtual | ~OrthoManager () |
Destructor. | |
Orthogonalization methods | |
virtual void | innerProd (const MV &X, const MV &Y, Teuchos::SerialDenseMatrix< int, ScalarType > &Z) const =0 |
Provides the inner product defining the orthogonality concepts. | |
virtual void | norm (const MV &X, std::vector< typename Teuchos::ScalarTraits< ScalarType >::magnitudeType > normvec) const =0 |
Provides the norm induced by innerProd(). | |
virtual void | project (MV &X, Teuchos::Array< Teuchos::RCP< Teuchos::SerialDenseMatrix< int, ScalarType > > > C, Teuchos::Array< Teuchos::RCP< const MV > > Q) const =0 |
Given a list of (mutually and internally) orthonormal bases Q , this method takes a multivector X and projects it onto the space orthogonal to the individual Q[i] , optionally returning the coefficients of X for the individual Q[i] . All of this is done with respect to the inner product innerProd(). | |
virtual int | normalize (MV &X, Teuchos::RCP< Teuchos::SerialDenseMatrix< int, ScalarType > > B) const =0 |
This method takes a multivector X and attempts to compute an orthonormal basis for , with respect to innerProd(). | |
virtual int | projectAndNormalize (MV &X, Teuchos::Array< Teuchos::RCP< Teuchos::SerialDenseMatrix< int, ScalarType > > > C, Teuchos::RCP< Teuchos::SerialDenseMatrix< int, ScalarType > > B, Teuchos::Array< Teuchos::RCP< const MV > > Q) const =0 |
Given a set of bases Q[i] and a multivector X , this method computes an orthonormal basis for . | |
Error methods | |
virtual Teuchos::ScalarTraits< ScalarType >::magnitudeType | orthonormError (const MV &X) const =0 |
This method computes the error in orthonormality of a multivector. | |
virtual Teuchos::ScalarTraits< ScalarType >::magnitudeType | orthogError (const MV &X1, const MV &X2) const =0 |
This method computes the error in orthogonality of two multivectors. | |
Label methods | |
virtual void | setLabel (const std::string &label)=0 |
This method sets the label used by the timers in the orthogonalization manager. | |
virtual const std::string & | getLabel () const =0 |
This method returns the label being used by the timers in the orthogonalization manager. |
This class defines concepts of orthogonality through the definition of an inner product. It also provides computational routines for orthogonalization.
A concrete implementation of this class is necessary. The user can create their own implementation if those supplied are not suitable for their needs.
Definition at line 74 of file BelosOrthoManager.hpp.
Belos::OrthoManager< ScalarType, MV >::OrthoManager | ( | ) | [inline] |
virtual Belos::OrthoManager< ScalarType, MV >::~OrthoManager | ( | ) | [inline, virtual] |
virtual void Belos::OrthoManager< ScalarType, MV >::innerProd | ( | const MV & | X, | |
const MV & | Y, | |||
Teuchos::SerialDenseMatrix< int, ScalarType > & | Z | |||
) | const [pure virtual] |
Provides the inner product defining the orthogonality concepts.
All concepts of orthogonality discussed in this class are with respect to this inner product.
M
, then this might be the M
inner product (). Implemented in Belos::MatOrthoManager< ScalarType, MV, OP >.
virtual void Belos::OrthoManager< ScalarType, MV >::norm | ( | const MV & | X, | |
std::vector< typename Teuchos::ScalarTraits< ScalarType >::magnitudeType > | normvec | |||
) | const [pure virtual] |
Provides the norm induced by innerProd().
Implemented in Belos::MatOrthoManager< ScalarType, MV, OP >.
virtual void Belos::OrthoManager< ScalarType, MV >::project | ( | MV & | X, | |
Teuchos::Array< Teuchos::RCP< Teuchos::SerialDenseMatrix< int, ScalarType > > > | C, | |||
Teuchos::Array< Teuchos::RCP< const MV > > | Q | |||
) | const [pure virtual] |
Given a list of (mutually and internally) orthonormal bases Q
, this method takes a multivector X
and projects it onto the space orthogonal to the individual Q[i]
, optionally returning the coefficients of X
for the individual Q[i]
. All of this is done with respect to the inner product innerProd().
After calling this routine, X
will be orthogonal to each of the Q[i]
.
X | [in/out] The multivector to be modified. On output, X will be orthogonal to Q[i] with respect to innerProd(). | |
C | [out] The coefficients of X in the *Q [i], with respect to innerProd(). If C[i] is a non-null pointer and *C [i] matches the dimensions of X and *Q [i], then the coefficients computed during the orthogonalization routine will be stored in the matrix *C [i]. If C[i] is a non-null pointer whose size does not match the dimensions of X and *Q [i], then a std::invalid_argument std::exception will be thrown. Otherwise, if C.size() < i or C[i] is a null pointer, then the orthogonalization manager will declare storage for the coefficients and the user will not have access to them. | |
Q | [in] A list of multivector bases specifying the subspaces to be orthogonalized against. Each Q[i] is assumed to have orthonormal columns, and the Q[i] are assumed to be mutually orthogonal. |
Implemented in Belos::DGKSOrthoManager< ScalarType, MV, OP >, Belos::ICGSOrthoManager< ScalarType, MV, OP >, Belos::IMGSOrthoManager< ScalarType, MV, OP >, and Belos::MatOrthoManager< ScalarType, MV, OP >.
virtual int Belos::OrthoManager< ScalarType, MV >::normalize | ( | MV & | X, | |
Teuchos::RCP< Teuchos::SerialDenseMatrix< int, ScalarType > > | B | |||
) | const [pure virtual] |
This method takes a multivector X
and attempts to compute an orthonormal basis for , with respect to innerProd().
This routine returns an integer rank
stating the rank of the computed basis. If X
does not have full rank and the normalize() routine does not attempt to augment the subspace, then rank
may be smaller than the number of columns in X
. In this case, only the first rank
columns of output X
and first rank
rows of B
will be valid.
X | [in/out] The multivector to the modified. On output, X will have some number of orthonormal columns (with respect to innerProd()). | |
B | [out] The coefficients of the original X with respect to the computed basis. This matrix is not necessarily triangular; see the documentation for specific orthogonalization managers. |
Implemented in Belos::DGKSOrthoManager< ScalarType, MV, OP >, Belos::ICGSOrthoManager< ScalarType, MV, OP >, Belos::IMGSOrthoManager< ScalarType, MV, OP >, and Belos::MatOrthoManager< ScalarType, MV, OP >.
virtual int Belos::OrthoManager< ScalarType, MV >::projectAndNormalize | ( | MV & | X, | |
Teuchos::Array< Teuchos::RCP< Teuchos::SerialDenseMatrix< int, ScalarType > > > | C, | |||
Teuchos::RCP< Teuchos::SerialDenseMatrix< int, ScalarType > > | B, | |||
Teuchos::Array< Teuchos::RCP< const MV > > | Q | |||
) | const [pure virtual] |
Given a set of bases Q[i]
and a multivector X
, this method computes an orthonormal basis for .
This routine returns an integer rank
stating the rank of the computed basis. If the subspace does not have dimension as large as the number of columns of X
and the orthogonalization manager doe not attempt to augment the subspace, then rank
may be smaller than the number of columns of X
. In this case, only the first rank
columns of output X
and first rank
rows of B
will be valid.
Q[i]
as well as the orthonormality constraints. Therefore, this method is not necessarily equivalent to calling project() followed by a call to normalize(); see the documentation for specific orthogonalization managers.X | [in/out] The multivector to the modified. On output, the relevant rows of X will be orthogonal to the Q[i] and will have orthonormal columns (with respect to innerProd()). | |
C | [out] The coefficients of the original X in the *Q [i], with respect to innerProd(). If C[i] is a non-null pointer and *C [i] matches the dimensions of X and *Q [i], then the coefficients computed during the orthogonalization routine will be stored in the matrix *C [i]. If C[i] is a non-null pointer whose size does not match the dimensions of X and *Q [i], then a std::invalid_argument std::exception will be thrown. Otherwise, if C.size() < i or C[i] is a null pointer, then the orthogonalization manager will declare storage for the coefficients and the user will not have access to them. | |
B | [out] The coefficients of the original X with respect to the computed basis. This matrix is not necessarily triangular; see the documentation for specific orthogonalization managers. | |
Q | [in] A list of multivector bases specifying the subspaces to be orthogonalized against. Each Q[i] is assumed to have orthonormal columns, and the Q[i] are assumed to be mutually orthogonal. |
Implemented in Belos::DGKSOrthoManager< ScalarType, MV, OP >, Belos::ICGSOrthoManager< ScalarType, MV, OP >, Belos::IMGSOrthoManager< ScalarType, MV, OP >, and Belos::MatOrthoManager< ScalarType, MV, OP >.
virtual Teuchos::ScalarTraits< ScalarType >::magnitudeType Belos::OrthoManager< ScalarType, MV >::orthonormError | ( | const MV & | X | ) | const [pure virtual] |
This method computes the error in orthonormality of a multivector.
Implemented in Belos::DGKSOrthoManager< ScalarType, MV, OP >, Belos::ICGSOrthoManager< ScalarType, MV, OP >, Belos::IMGSOrthoManager< ScalarType, MV, OP >, and Belos::MatOrthoManager< ScalarType, MV, OP >.
virtual Teuchos::ScalarTraits<ScalarType>::magnitudeType Belos::OrthoManager< ScalarType, MV >::orthogError | ( | const MV & | X1, | |
const MV & | X2 | |||
) | const [pure virtual] |
This method computes the error in orthogonality of two multivectors.
Implemented in Belos::DGKSOrthoManager< ScalarType, MV, OP >, Belos::ICGSOrthoManager< ScalarType, MV, OP >, Belos::IMGSOrthoManager< ScalarType, MV, OP >, and Belos::MatOrthoManager< ScalarType, MV, OP >.
virtual void Belos::OrthoManager< ScalarType, MV >::setLabel | ( | const std::string & | label | ) | [pure virtual] |
This method sets the label used by the timers in the orthogonalization manager.
Implemented in Belos::DGKSOrthoManager< ScalarType, MV, OP >, Belos::ICGSOrthoManager< ScalarType, MV, OP >, and Belos::IMGSOrthoManager< ScalarType, MV, OP >.
virtual const std::string& Belos::OrthoManager< ScalarType, MV >::getLabel | ( | ) | const [pure virtual] |
This method returns the label being used by the timers in the orthogonalization manager.
Implemented in Belos::DGKSOrthoManager< ScalarType, MV, OP >, Belos::ICGSOrthoManager< ScalarType, MV, OP >, and Belos::IMGSOrthoManager< ScalarType, MV, OP >.