45 #ifndef THYRA_BELOS_LINEAR_OP_WITH_SOLVE_DECL_HPP 46 #define THYRA_BELOS_LINEAR_OP_WITH_SOLVE_DECL_HPP 48 #include "Thyra_LinearOpWithSolveBase.hpp" 49 #include "Thyra_LinearOpSourceBase.hpp" 50 #include "BelosSolverManager.hpp" 52 #include "Teuchos_StandardMemberCompositionMacros.hpp" 65 template<
class Scalar>
74 typedef MultiVectorBase<Scalar>
MV_t;
76 typedef LinearOpBase<Scalar>
LO_t;
126 const RCP<Belos::LinearProblem<Scalar,MV_t,LO_t> > &lp,
127 const RCP<Teuchos::ParameterList> &solverPL,
128 const RCP<Belos::SolverManager<Scalar,MV_t,LO_t> > &iterativeSolver,
129 const RCP<
const LinearOpSourceBase<Scalar> > &fwdOpSrc,
130 const RCP<
const PreconditionerBase<Scalar> > &prec,
132 const RCP<
const LinearOpSourceBase<Scalar> > &approxFwdOpSrc,
134 const int convergenceTestFrequency
157 RCP<Belos::LinearProblem<Scalar,MV_t,LO_t> > *lp = NULL,
158 RCP<Teuchos::ParameterList> *solverPL = NULL,
159 RCP<Belos::SolverManager<Scalar,MV_t,LO_t> > *iterativeSolver = NULL,
160 RCP<
const LinearOpSourceBase<Scalar> > *fwdOpSrc = NULL,
161 RCP<
const PreconditionerBase<Scalar> > *prec = NULL,
163 RCP<
const LinearOpSourceBase<Scalar> > *approxFwdOpSrc = NULL,
172 RCP< const VectorSpaceBase<Scalar> >
range()
const;
174 RCP< const VectorSpaceBase<Scalar> >
domain()
const;
176 RCP<const LinearOpBase<Scalar> >
clone()
const;
185 Teuchos::FancyOStream &out,
186 const Teuchos::EVerbosityLevel verbLevel
214 const EOpTransp M_trans,
215 const MultiVectorBase<Scalar> &X,
216 const Ptr<MultiVectorBase<Scalar> > &Y,
228 const Ptr<
const SolveCriteria<Scalar> > solveCriteria)
const;
231 EOpTransp M_trans,
const SolveMeasureType& solveMeasureType
235 const EOpTransp transp,
236 const MultiVectorBase<Scalar> &B,
237 const Ptr<MultiVectorBase<Scalar> > &X,
238 const Ptr<
const SolveCriteria<Scalar> > solveCriteria
248 RCP<Belos::LinearProblem<Scalar,MV_t,LO_t> > lp_;
249 RCP<Teuchos::ParameterList> solverPL_;
250 RCP<Belos::SolverManager<Scalar,MV_t,LO_t> > iterativeSolver_;
251 int convergenceTestFrequency_;
253 RCP<const LinearOpSourceBase<Scalar> > fwdOpSrc_;
254 RCP<const PreconditionerBase<Scalar> > prec_;
255 bool isExternalPrec_;
256 RCP<const LinearOpSourceBase<Scalar> > approxFwdOpSrc_;
257 ESupportSolveUse supportSolveUse_;
259 typename Teuchos::ScalarTraits<Scalar>::magnitudeType defaultTol_;
261 void assertInitialized()
const;
269 #endif // THYRA_BELOS_LINEAR_OP_WITH_SOLVE_DECL_HPP RCP< const Teuchos::ParameterList > getParameterList() const
virtual bool solveSupportsNewImpl(EOpTransp transp, const Ptr< const SolveCriteria< Scalar > > solveCriteria) const
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
RCP< const LinearOpSourceBase< Scalar > > extract_approxFwdOpSrc()
void initialize(const RCP< Belos::LinearProblem< Scalar, MV_t, LO_t > > &lp, const RCP< Teuchos::ParameterList > &solverPL, const RCP< Belos::SolverManager< Scalar, MV_t, LO_t > > &iterativeSolver, const RCP< const LinearOpSourceBase< Scalar > > &fwdOpSrc, const RCP< const PreconditionerBase< Scalar > > &prec, const bool isExternalPrec, const RCP< const LinearOpSourceBase< Scalar > > &approxFwdOpSrc, const ESupportSolveUse &supportSolveUse, const int convergenceTestFrequency)
Initializes given precreated solver objects.
Concrete LinearOpWithSolveBase subclass in terms of Belos.
RCP< const PreconditionerBase< Scalar > > extract_prec()
RCP< const VectorSpaceBase< Scalar > > domain() const
RCP< const LinearOpBase< Scalar > > clone() const
std::string description() const
virtual bool solveSupportsSolveMeasureTypeImpl(EOpTransp M_trans, const SolveMeasureType &solveMeasureType) const
void setParameterList(RCP< Teuchos::ParameterList > const ¶mList)
Thyra specializations of MultiVecTraits and OperatorTraits.
LinearOpBase< Scalar > LO_t
bool isExternalPrec() const
virtual SolveStatus< Scalar > solveImpl(const EOpTransp transp, const MultiVectorBase< Scalar > &B, const Ptr< MultiVectorBase< Scalar > > &X, const Ptr< const SolveCriteria< Scalar > > solveCriteria) const
RCP< const Teuchos::ParameterList > getValidParameters() const
RCP< const VectorSpaceBase< Scalar > > range() const
virtual void applyImpl(const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha, const Scalar beta) const
RCP< const LinearOpSourceBase< Scalar > > extract_fwdOpSrc()
virtual bool solveSupportsImpl(EOpTransp M_trans) const
BelosLinearOpWithSolve()
Construct to unintialize.
MultiVectorBase< Scalar > MV_t
void uninitialize(RCP< Belos::LinearProblem< Scalar, MV_t, LO_t > > *lp=NULL, RCP< Teuchos::ParameterList > *solverPL=NULL, RCP< Belos::SolverManager< Scalar, MV_t, LO_t > > *iterativeSolver=NULL, RCP< const LinearOpSourceBase< Scalar > > *fwdOpSrc=NULL, RCP< const PreconditionerBase< Scalar > > *prec=NULL, bool *isExternalPrec=NULL, RCP< const LinearOpSourceBase< Scalar > > *approxFwdOpSrc=NULL, ESupportSolveUse *supportSolveUse=NULL)
Uninitializes and returns stored quantities.
ESupportSolveUse supportSolveUse() const
RCP< Teuchos::ParameterList > unsetParameterList()
virtual bool opSupportedImpl(EOpTransp M_trans) const
RCP< Teuchos::ParameterList > getNonconstParameterList()