42#ifndef THYRA_MULTI_VECTOR_LINEAR_OP_WITH_SOLVE_HPP
43#define THYRA_MULTI_VECTOR_LINEAR_OP_WITH_SOLVE_HPP
45#include "Thyra_DefaultMultiVectorLinearOpWithSolve_decl.hpp"
46#include "Thyra_DefaultDiagonalLinearOp.hpp"
47#include "Thyra_LinearOpWithSolveBase.hpp"
48#include "Thyra_DefaultMultiVectorProductVectorSpace.hpp"
49#include "Thyra_DefaultMultiVectorProductVector.hpp"
50#include "Thyra_AssertOp.hpp"
51#include "Teuchos_dyn_cast.hpp"
72 validateInitialize(lows,multiVecRange,multiVecDomain);
74 multiVecRange_ = multiVecRange;
75 multiVecDomain_ = multiVecDomain;
86 validateInitialize(lows,multiVecRange,multiVecDomain);
88 multiVecRange_ = multiVecRange;
89 multiVecDomain_ = multiVecDomain;
97 return lows_.getNonconstObj();
101template<
class Scalar>
105 return lows_.getConstObj();
109template<
class Scalar>
112 lows_.uninitialize();
113 multiVecRange_ = Teuchos::null;
114 multiVecDomain_ = Teuchos::null;
121template<
class Scalar>
125 return multiVecRange_;
129template<
class Scalar>
133 return multiVecDomain_;
137template<
class Scalar>
141 return Teuchos::null;
151template<
class Scalar>
156 return Thyra::opSupported(*lows_.getConstObj(),M_trans);
160template<
class Scalar>
175 for (
Ordinal col_j = 0; col_j < numCols; ++col_j) {
185 Thyra::apply( *lows_.getConstObj(), M_trans, *X, Y.
ptr(), alpha, beta );
195template<
class Scalar>
201 return Thyra::solveSupports(*lows_.getConstObj(),M_trans);
205template<
class Scalar>
211 return Thyra::solveSupportsSolveMeasureType(
212 *lows_.getConstObj(),M_trans,solveMeasureType);
216template<
class Scalar>
227 using Teuchos::outArg;
228 using Teuchos::inOutArg;
234 accumulateSolveStatusInit(outArg(overallSolveStatus));
236 for (
Ordinal col_j = 0; col_j < numCols; ++col_j) {
247 Thyra::solve(*lows_.getConstObj(), transp, *B, X.
ptr(), solveCriteria);
249 accumulateSolveStatus(
251 solveStatus, inOutArg(overallSolveStatus) );
255 return overallSolveStatus;
263template<
class Scalar>
275 if (lows->range() != Teuchos::null)
277 "DefaultMultiVectorLinearOpWithSolve<Scalar>::initialize(lows,multiVecRange,multiVecDomain)",
278 *lows->range(), *multiVecRange->getBlock(0) );
279 if (lows->domain() != Teuchos::null)
281 "DefaultMultiVectorLinearOpWithSolve<Scalar>::initialize(lows,multiVecRange,multiVecDomain)",
282 *lows->domain(), *multiVecDomain->getBlock(0) );
286 (void)multiVecDomain;
const RCP< T > & assert_not_null() const
Implicit concrete LinearOpWithSolveBase subclass that takes a flattended out multi-vector and perform...
RCP< const LinearOpWithSolveBase< Scalar > > getLinearOpWithSolve() const
RCP< LinearOpWithSolveBase< Scalar > > getNonconstLinearOpWithSolve()
void applyImpl(const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha, const Scalar beta) const
DefaultMultiVectorLinearOpWithSolve()
Construct to uninitialized.
bool opSupportedImpl(EOpTransp M_trans) const
bool solveSupportsSolveMeasureTypeImpl(EOpTransp M_trans, const SolveMeasureType &solveMeasureType) const
RCP< const VectorSpaceBase< Scalar > > range() const
RCP< const LinearOpBase< Scalar > > clone() const
void nonconstInitialize(const RCP< LinearOpWithSolveBase< Scalar > > &lows, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
RCP< const VectorSpaceBase< Scalar > > domain() const
bool solveSupportsImpl(EOpTransp M_trans) const
SolveStatus< Scalar > solveImpl(const EOpTransp transp, const MultiVectorBase< Scalar > &B, const Ptr< MultiVectorBase< Scalar > > &X, const Ptr< const SolveCriteria< Scalar > > solveCriteria) const
void initialize(const RCP< const LinearOpWithSolveBase< Scalar > > &lows, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
Standard concrete implementation of a product vector space that creates product vectors fromed implic...
Concrete implementation of a product vector which is really composed out of the columns of a multi-ve...
virtual RCP< const VectorSpaceBase< Scalar > > domain() const =0
Return a smart pointer for the domain space for this operator.
Base class for all linear operators that can support a high-level solve operation.
Interface for a collection of column vectors called a multi-vector.
RCP< const VectorBase< Scalar > > col(Ordinal j) const
Calls colImpl().
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)
#define THYRA_ASSERT_VEC_SPACES(FUNC_NAME, VS1, VS2)
This is a very useful macro that should be used to validate that two vector spaces are compatible.
EOpTransp
Enumeration for determining how a linear operator is applied. `*.
Teuchos::Ordinal Ordinal
Type for the dimension of a vector space. `*.
T_To & dyn_cast(T_From &from)
Simple struct that defines the requested solution criteria for a solve.
Simple struct for the return status from a solve.