9 #ifndef Thyra_ImplicitAdjointModelEvaluator_hpp
10 #define Thyra_ImplicitAdjointModelEvaluator_hpp
12 #include "Thyra_ModelEvaluatorDelegatorBase.hpp"
13 #include "Thyra_DefaultScaledAdjointLinearOp.hpp"
14 #include "Thyra_DefaultAdjointLinearOpWithSolve.hpp"
24 template <
typename Scalar>
44 RCP<LinearOpWithSolveBase<Scalar> >
create_W()
const
46 return nonconstAdjointLows(this->getUnderlyingModel()->
create_W());
52 return nonconstAdjoint(this->getUnderlyingModel()->
create_W_op());
58 return nonconstAdjointPreconditioner(
65 return adjointLinearOpWithSolveFactory(
74 MEB::OutArgs<Scalar> model_outArgs =
75 this->getUnderlyingModel()->createOutArgs();
77 if (model_outArgs.supports(MEB::OUT_ARG_W) &&
78 outArgs.
get_W() != Teuchos::null) {
79 RCP<DefaultAdjointLinearOpWithSolve<Scalar> > adjoint_op =
81 outArgs.
get_W(),
true);
82 model_outArgs.set_W(adjoint_op->getNonconstOp());
85 if (model_outArgs.supports(MEB::OUT_ARG_W_op) &&
86 outArgs.
get_W_op() != Teuchos::null) {
87 RCP<DefaultScaledAdjointLinearOp<Scalar> > adjoint_op =
90 model_outArgs.set_W_op(adjoint_op->getNonconstOp());
93 if (model_outArgs.supports(MEB::OUT_ARG_W_prec) &&
95 RCP<AdjointPreconditioner<Scalar> > adjoint_op =
98 model_outArgs.set_W_prec(adjoint_op->getNonconstPreconditioner());
101 this->getUnderlyingModel()->evalModel(inArgs, model_outArgs);
105 template <
typename Scalar>
112 template <
typename Scalar>
Concrete PreconditionerBase subclass that wraps a preconditioner operator in MultiVectorLinearOp.
RCP< PreconditionerBase< Scalar > > create_W_prec() const
Create adjoint preconditioner.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
RCP< const LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
Get adjoint solver factory.
RCP< ImplicitAdjointModelEvaluator< Scalar > > implicitAdjointModelEvaluator(const RCP< const ModelEvaluator< Scalar > > &model)
RCP< LinearOpBase< Scalar > > create_W_op() const
Create adjoint op.
RCP< LinearOpWithSolveBase< Scalar > > get_W() const
virtual ~ImplicitAdjointModelEvaluator()=default
Destructor.
An implementation of AdjointModelEvaluatorBase that creates an implicit adjoint from the supplied mod...
RCP< PreconditionerBase< Scalar > > get_W_prec() const
RCP< LinearOpWithSolveBase< Scalar > > create_W() const
Create adjoint solver.
ImplicitAdjointModelEvaluator(const RCP< const ModelEvaluator< Scalar > > &model)
Constructor.
void evalModelImpl(const ModelEvaluatorBase::InArgs< Scalar > &inArgs, const ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
RCP< LinearOpBase< Scalar > > get_W_op() const
ImplicitAdjointModelEvaluator(const RCP< ModelEvaluator< Scalar > > &model)
Constructor.