9#ifndef Tempus_AdjointSensitivityModelEvaluator_decl_hpp
10#define Tempus_AdjointSensitivityModelEvaluator_decl_hpp
12#include "Thyra_StateFuncModelEvaluatorBase.hpp"
13#include "Thyra_DefaultMultiVectorProductVectorSpace.hpp"
14#include "Thyra_DefaultMultiVectorProductVector.hpp"
16#include "Tempus_config.hpp"
17#include "Tempus_SolutionHistory.hpp"
53template <
typename Scalar>
55 public Thyra::StateFuncModelEvaluatorBase<Scalar> {
84 const Scalar& t_final,
85 const bool is_pseudotransient,
86 const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null);
89 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
getModel()
const
110 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_p_space(
int p)
const;
112 Teuchos::RCP<const Teuchos::Array<std::string> >
get_p_names(
int p)
const;
114 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_x_space()
const;
116 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_f_space()
const;
118 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_g_space(
int j)
const;
120 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
create_W_op()
const;
122 Teuchos::RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> >
125 Thyra::ModelEvaluatorBase::InArgs<Scalar>
createInArgs()
const;
135 typedef Thyra::DefaultMultiVectorProductVectorSpace<Scalar>
DMVPVS;
136 typedef Thyra::DefaultMultiVectorProductVector<Scalar>
DMVPV;
141 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
142 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs)
const;
148 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
model_;
155 Teuchos::RCP<const Tempus::SolutionHistory<Scalar> >
sh_;
168 mutable Teuchos::RCP<Thyra::VectorBase<Scalar> >
my_x_dot_;
169 mutable Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
my_dfdx_;
172 mutable Teuchos::RCP<Thyra::MultiVectorBase<Scalar> >
my_dfdp_mv_;
173 mutable Teuchos::RCP<Thyra::MultiVectorBase<Scalar> >
my_dgdx_mv_;
ModelEvaluator for forming adjoint sensitivity equations.
bool response_gradient_is_computed_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getAdjointResidualModel() const
Get the underlying adjoint residual model.
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int p) const
Teuchos::RCP< Thyra::MultiVectorBase< Scalar > > my_dgdx_mv_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > adjoint_residual_model_
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
bool mass_matrix_is_identity_
Teuchos::RCP< Tempus::SolutionState< Scalar > > forward_state_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
Thyra::VectorBase< Scalar > Vector
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getAdjointSolveModel() const
Get the underlying adjoint solve model.
void setForwardSolutionHistory(const Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > &sh)
Set solution history from forward evaluation.
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
Thyra::MultiVectorBase< Scalar > MultiVector
bool mass_matrix_is_computed_
static Teuchos::RCP< const Teuchos::ParameterList > getValidParameters()
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > model_
Teuchos::RCP< const DMVPVS > residual_space_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getModel() const
Get the underlying model 'f'.
Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > sh_
Thyra::ModelEvaluatorBase::InArgs< Scalar > prototypeInArgs_
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
bool jacobian_matrix_is_computed_
Thyra::DefaultMultiVectorProductVector< Scalar > DMVPV
Thyra::DefaultMultiVectorProductVectorSpace< Scalar > DMVPVS
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int p) const
Teuchos::RCP< const DMVPVS > response_space_
Thyra::ModelEvaluatorBase::OutArgs< Scalar > prototypeOutArgs_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdp_op_
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdx_
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > adjoint_solve_model_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdxdot_
void setFinalTime(const Scalar t_final)
Set the final time from the forward evaluation.
bool mass_matrix_is_constant_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_g_space(int j) const
Teuchos::RCP< Thyra::MultiVectorBase< Scalar > > my_dfdp_mv_
Teuchos::RCP< Thyra::VectorBase< Scalar > > my_x_dot_
Teuchos::RCP< const DMVPVS > adjoint_space_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...