9#ifndef Tempus_AdjointAuxSensitivityModelEvaluator_decl_hpp
10#define Tempus_AdjointAuxSensitivityModelEvaluator_decl_hpp
12#include "Thyra_StateFuncModelEvaluatorBase.hpp"
13#include "Thyra_DefaultProductVectorSpace.hpp"
14#include "Thyra_DefaultProductVector.hpp"
15#include "Thyra_DefaultMultiVectorProductVectorSpace.hpp"
16#include "Thyra_DefaultMultiVectorProductVector.hpp"
18#include "Tempus_config.hpp"
19#include "Tempus_SolutionHistory.hpp"
56template <
typename Scalar>
58 public Thyra::StateFuncModelEvaluatorBase<Scalar> {
86 const Scalar& t_final,
87 const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null);
90 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
getModel()
const
107 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_p_space(
int p)
const;
109 Teuchos::RCP<const Teuchos::Array<std::string> >
get_p_names(
int p)
const;
111 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_x_space()
const;
113 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_f_space()
const;
115 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
create_W_op()
const;
117 Teuchos::RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> >
120 Thyra::ModelEvaluatorBase::InArgs<Scalar>
createInArgs()
const;
130 typedef Thyra::DefaultProductVectorSpace<Scalar>
DPVS;
131 typedef Thyra::DefaultProductVector<Scalar>
DPV;
132 typedef Thyra::DefaultMultiVectorProductVectorSpace<Scalar>
DMVPVS;
133 typedef Thyra::DefaultMultiVectorProductVector<Scalar>
DMVPV;
138 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
139 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs)
const;
145 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
model_;
153 Teuchos::RCP<const Tempus::SolutionHistory<Scalar> >
sh_;
163 mutable Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
my_dfdx_;
166 mutable Teuchos::RCP<Thyra::MultiVectorBase<Scalar> >
my_dfdp_mv_;
ModelEvaluator for forming adjoint sensitivity equations.
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdp_op_
Teuchos::RCP< const DMVPVS > residual_space_
Teuchos::RCP< const DPVS > x_prod_space_
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int p) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
Teuchos::RCP< const DMVPVS > adjoint_space_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > adjoint_model_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
Thyra::DefaultMultiVectorProductVectorSpace< Scalar > DMVPVS
Teuchos::RCP< const DPVS > f_prod_space_
Teuchos::RCP< const DMVPVS > response_space_
static Teuchos::RCP< const Teuchos::ParameterList > getValidParameters()
Teuchos::RCP< Thyra::MultiVectorBase< Scalar > > my_dfdp_mv_
bool mass_matrix_is_identity_
void setFinalTime(const Scalar t_final)
Set the final time from the forward evaluation.
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
Thyra::DefaultProductVector< Scalar > DPV
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getAdjointModel() const
Get the underlying adjoint model.
Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > sh_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getModel() const
Get the underlying model 'f'.
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdxdot_
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
void setForwardSolutionHistory(const Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > &sh)
Set solution history from forward evaluation.
Thyra::DefaultMultiVectorProductVector< Scalar > DMVPV
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > model_
Thyra::VectorBase< Scalar > Vector
Thyra::DefaultProductVectorSpace< Scalar > DPVS
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int p) const
bool mass_matrix_is_computed_
Thyra::ModelEvaluatorBase::OutArgs< Scalar > prototypeOutArgs_
Teuchos::RCP< Tempus::SolutionState< Scalar > > forward_state_
bool mass_matrix_is_constant_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdx_
Thyra::MultiVectorBase< Scalar > MultiVector
Thyra::ModelEvaluatorBase::InArgs< Scalar > prototypeInArgs_
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...