43 #ifndef THYRA_DELAYED_LINEAR_OP_WITH_SOLVE_HPP 44 #define THYRA_DELAYED_LINEAR_OP_WITH_SOLVE_HPP 47 #include "Thyra_DelayedLinearOpWithSolve_decl.hpp" 48 #include "Thyra_LinearOpWithSolveFactoryBase.hpp" 49 #include "Thyra_LinearOpWithSolveBase.hpp" 50 #include "Thyra_SolveSupportTypes.hpp" 51 #include "Teuchos_VerboseObject.hpp" 60 template <
class Scalar>
66 template <
class Scalar>
76 TEUCHOS_TEST_FOR_EXCEPT(is_null(fwdOpSrc));
77 TEUCHOS_TEST_FOR_EXCEPT(!is_null(prec) && !is_null(approxFwdOpSrc));
78 TEUCHOS_TEST_FOR_EXCEPT(is_null(lowsf));
82 approxFwdOpSrc_ = approxFwdOpSrc;
84 supportSolveUse_ = supportSolveUse;
85 fwdOp_ = fwdOpSrc_->getOp().assert_not_null();
86 lows_is_valid_ =
false;
90 template <
class Scalar>
91 RCP<const LinearOpSourceBase<Scalar> >
98 template <
class Scalar>
99 RCP<const PreconditionerBase<Scalar> >
106 template <
class Scalar>
107 RCP<const LinearOpSourceBase<Scalar> >
110 return approxFwdOpSrc_;
114 template <
class Scalar>
118 return supportSolveUse_;
125 template<
class Scalar>
128 std::ostringstream oss;
129 oss << this->Teuchos::Describable::description()
132 if (!is_null(fwdOp_))
133 oss << fwdOp_->description();
138 oss << lows_->description();
149 template <
class Scalar>
150 RCP< const VectorSpaceBase<Scalar> >
153 if (!is_null(fwdOp_))
154 return fwdOp_->range();
155 return Teuchos::null;
159 template <
class Scalar>
160 RCP< const VectorSpaceBase<Scalar> >
163 if (!is_null(fwdOp_))
164 return fwdOp_->domain();
165 return Teuchos::null;
169 template <
class Scalar>
170 RCP<const LinearOpBase<Scalar> >
173 return Teuchos::null;
180 template<
class Scalar>
183 if (!is_null(lowsf_)) {
184 lowsf_->setVerbLevel(this->getVerbLevel());
185 lowsf_->setOStream(this->getOStream());
187 if (!is_null(lows_)) {
188 lows_->setVerbLevel(this->getVerbLevel());
189 lows_->setOStream(this->getOStream());
197 template <
class Scalar>
200 return Thyra::opSupported(*fwdOp_,M_trans);
204 template <
class Scalar>
213 Thyra::apply(*fwdOp_, M_trans, X, Y, alpha, beta);
220 template <
class Scalar>
224 return Thyra::solveSupports(*lows_, M_trans);
228 template <
class Scalar>
234 return Thyra::solveSupportsSolveMeasureType(*lows_, M_trans, solveMeasureType);
238 template <
class Scalar>
248 return Thyra::solve(*lows_, transp, B, X, solveCriteria);
255 template <
class Scalar>
259 lows_ = lowsf_->createOp();
260 if (!lows_is_valid_) {
262 lowsf_->initializePreconditionedOp(
263 fwdOpSrc_,prec_,&*lows_,supportSolveUse_);
264 else if (!is_null(approxFwdOpSrc_))
265 lowsf_->initializeApproxPreconditionedOp(
266 fwdOpSrc_,approxFwdOpSrc_,&*lows_,supportSolveUse_);
268 lowsf_->initializeOp(
269 fwdOpSrc_,&*lows_,supportSolveUse_);
270 lows_is_valid_ =
true;
278 #endif // THYRA_DELAYED_LINEAR_OP_WITH_SOLVE_HPP How the output LOWSB object will be useded for solves in unspecified.
EOpTransp
Enumeration for determining how a linear operator is applied. `*.
RCP< const VectorSpaceBase< Scalar > > range() const
RCP< const LinearOpSourceBase< Scalar > > getFwdOpSrc() const
virtual void applyImpl(const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha, const Scalar beta) const
SolveStatus< Scalar > solveImpl(const EOpTransp transp, const MultiVectorBase< Scalar > &B, const Ptr< MultiVectorBase< Scalar > > &X, const Ptr< const SolveCriteria< Scalar > > solveCriteria) const
std::string description() const
Simple interface class to access a precreated preconditioner as one or more linear operators objects ...
DelayedLinearOpWithSolve()
void informUpdatedVerbosityState() const
Overridden from Teuchos::VerboseObjectBase.
Interface for a collection of column vectors called a multi-vector.
Factory interface for creating LinearOpWithSolveBase objects from compatible LinearOpBase objects...
RCP< const LinearOpSourceBase< Scalar > > getApproxFwdOpSrc() const
void initialize(const RCP< const LinearOpSourceBase< Scalar > > &fwdOpSrc, const RCP< const PreconditionerBase< Scalar > > &prec, const RCP< const LinearOpSourceBase< Scalar > > &approxFwdOpSrc, const ESupportSolveUse supportSolveUse, const RCP< LinearOpWithSolveFactoryBase< Scalar > > &lowsf)
Simple struct for the return status from a solve.
Base interface for objects that can return a linear operator.
ESupportSolveUse
Enum that specifies how a LinearOpWithSolveBase object will be used for solves after it is constructe...
virtual bool opSupportedImpl(EOpTransp M_trans) const
virtual bool solveSupportsSolveMeasureTypeImpl(EOpTransp M_trans, const SolveMeasureType &solveMeasureType) const
virtual bool solveSupportsImpl(EOpTransp M_trans) const
RCP< const PreconditionerBase< Scalar > > getPrec() const
Simple struct that defines the requested solution criteria for a solve.
Delayed linear solver construction LinearOpWithSolveBase decorator class.
ESupportSolveUse getSupportSolveUse() const
RCP< const LinearOpBase< Scalar > > clone() const
RCP< const VectorSpaceBase< Scalar > > domain() const