42 #ifndef THYRA_DEFAULT_INVERSE_LINEAR_OP_DECL_HPP 43 #define THYRA_DEFAULT_INVERSE_LINEAR_OP_DECL_HPP 45 #include "Thyra_InverseLinearOpBase.hpp" 46 #include "Teuchos_ConstNonconstObjectContainer.hpp" 85 template<
class Scalar>
223 RCP<LinearOpWithSolveBase<Scalar> >
226 RCP<const LinearOpWithSolveBase<Scalar> >
238 RCP< const VectorSpaceBase<Scalar> >
range()
const;
244 RCP< const VectorSpaceBase<Scalar> >
domain()
const;
247 RCP<const LinearOpBase<Scalar> >
clone()
const;
260 const Teuchos::EVerbosityLevel verbLevel
288 Teuchos::ConstNonconstObjectContainer<LinearOpWithSolveBase<Scalar> > lows_;
289 RCP<SolveCriteria<Scalar> > fwdSolveCriteria_;
291 RCP<SolveCriteria<Scalar> > adjSolveCriteria_;
294 void assertInitialized()
const;
298 const RCP<LOWS> &lows,
316 template<
class Scalar>
317 RCP<LinearOpBase<Scalar> >
331 template<
class Scalar>
332 RCP<LinearOpBase<Scalar> >
346 template<
class Scalar>
351 TEUCHOS_TEST_FOR_EXCEPT( !lows_.getConstObj().get() );
359 #endif // THYRA_DEFAULT_INVERSE_LINEAR_OP_DECL_HPP Base class for all linear operators that can support a high-level solve operation.
RCP< const VectorSpaceBase< Scalar > > domain() const
Returns this->getLows()->range() if <t>this->getLows().get()!=NULL and returns Teuchos::null otherwis...
EOpTransp
Enumeration for determining how a linear operator is applied. `*.
void uninitialize()
Set to uninitialized.
DefaultInverseLinearOp()
Constructs to uninitialized (see postconditions for uninitialize()).
bool opSupportedImpl(EOpTransp M_trans) const
Returns true only if all constituent operators support M_trans.
RCP< const LinearOpBase< Scalar > > clone() const
std::string description() const
RCP< const LinearOpWithSolveBase< Scalar > > getLows() const
RCP< const VectorSpaceBase< Scalar > > range() const
Returns this->getLows()->domain() if <t>this->getLows().get()!=NULL and returns Teuchos::null otherwi...
Throw an exception if a solve fails to converge.
Interface for a collection of column vectors called a multi-vector.
void describe(FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
void applyImpl(const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha, const Scalar beta) const
void initialize(const RCP< LinearOpWithSolveBase< Scalar > > &lows, const SolveCriteria< Scalar > *fwdSolveCriteria=NULL, const EThrowOnSolveFailure throwOnFwdSolveFailure=THROW_ON_SOLVE_FAILURE, const SolveCriteria< Scalar > *adjSolveCriteria=NULL, const EThrowOnSolveFailure throwOnAdjSolveFailure=THROW_ON_SOLVE_FAILURE)
Initialize given a non-const LinearOpWithSolveBase object and an optional .
Base interface for <ttLinearOpBase objects that are implemented in terms of the solve function on a L...
EThrowOnSolveFailure
Determines what to do if inverse solve fails.
Don't throw an exception if a solve fails to converge.
Simple struct that defines the requested solution criteria for a solve.
Concrete LinearOpBase subclass that creates an implicit LinearOpBase object using the inverse action ...
RCP< LinearOpWithSolveBase< Scalar > > getNonconstLows()