42 #ifndef THYRA_MULTI_VECTOR_LINEAR_OP_WITH_SOLVE_DECL_HPP 43 #define THYRA_MULTI_VECTOR_LINEAR_OP_WITH_SOLVE_DECL_HPP 46 #include "Thyra_LinearOpWithSolveBase.hpp" 47 #include "Thyra_DefaultDiagonalLinearOp.hpp" 48 #include "Thyra_LinearOpWithSolveBase.hpp" 49 #include "Thyra_DefaultMultiVectorProductVectorSpace.hpp" 50 #include "Teuchos_ConstNonconstObjectContainer.hpp" 61 template<
class Scalar>
88 RCP<LinearOpWithSolveBase<Scalar> >
92 RCP<const LinearOpWithSolveBase<Scalar> >
107 RCP< const VectorSpaceBase<Scalar> >
range()
const;
110 RCP< const VectorSpaceBase<Scalar> >
domain()
const;
113 RCP<const LinearOpBase<Scalar> >
clone()
const;
154 typedef Teuchos::ConstNonconstObjectContainer<LinearOpWithSolveBase<Scalar> > CNLOWS;
160 RCP<const DefaultMultiVectorProductVectorSpace<Scalar> > multiVecRange_;
161 RCP<const DefaultMultiVectorProductVectorSpace<Scalar> > multiVecDomain_;
166 static void validateInitialize(
180 template<
class Scalar>
181 RCP<DefaultMultiVectorLinearOpWithSolve<Scalar> >
192 template<
class Scalar>
193 RCP<DefaultMultiVectorLinearOpWithSolve<Scalar> >
200 RCP<DefaultMultiVectorLinearOpWithSolve<Scalar> >
202 mvlows->nonconstInitialize(lows,multiVecRange,multiVecDomain);
211 template<
class Scalar>
212 RCP<DefaultMultiVectorLinearOpWithSolve<Scalar> >
219 RCP<DefaultMultiVectorLinearOpWithSolve<Scalar> >
221 mvlows->initialize(lows,multiVecRange,multiVecDomain);
229 #endif // THYRA_MULTI_VECTOR_LINEAR_OP_WITH_SOLVE_DECL_HPP Implicit concrete LinearOpWithSolveBase subclass that takes a flattended out multi-vector and perform...
Base class for all linear operators that can support a high-level solve operation.
EOpTransp
Enumeration for determining how a linear operator is applied. `*.
bool solveSupportsImpl(EOpTransp M_trans) const
DefaultMultiVectorLinearOpWithSolve()
Construct to uninitialized.
bool opSupportedImpl(EOpTransp M_trans) const
bool solveSupportsSolveMeasureTypeImpl(EOpTransp M_trans, const SolveMeasureType &solveMeasureType) const
Interface for a collection of column vectors called a multi-vector.
void applyImpl(const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha, const Scalar beta) const
void nonconstInitialize(const RCP< LinearOpWithSolveBase< Scalar > > &lows, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
RCP< const LinearOpWithSolveBase< Scalar > > getLinearOpWithSolve() const
Simple struct for the return status from a solve.
Standard concrete implementation of a product vector space that creates product vectors fromed implic...
RCP< DefaultMultiVectorLinearOpWithSolve< Scalar > > nonconstMultiVectorLinearOpWithSolve(const RCP< LinearOpWithSolveBase< Scalar > > &lows, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
Nonmember constructor function.
RCP< const LinearOpBase< Scalar > > clone() const
RCP< DefaultMultiVectorLinearOpWithSolve< Scalar > > multiVectorLinearOpWithSolve()
Nonmember constructor function.
RCP< const VectorSpaceBase< Scalar > > domain() const
RCP< LinearOpWithSolveBase< Scalar > > getNonconstLinearOpWithSolve()
SolveStatus< Scalar > solveImpl(const EOpTransp transp, const MultiVectorBase< Scalar > &B, const Ptr< MultiVectorBase< Scalar > > &X, const Ptr< const SolveCriteria< Scalar > > solveCriteria) const
Simple struct that defines the requested solution criteria for a solve.
RCP< DefaultMultiVectorLinearOpWithSolve< Scalar > > multiVectorLinearOpWithSolve(const RCP< const LinearOpWithSolveBase< Scalar > > &lows, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
Nonmember constructor function.
void initialize(const RCP< const LinearOpWithSolveBase< Scalar > > &lows, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
RCP< const VectorSpaceBase< Scalar > > range() const