9#ifndef Tempus_CombinedForwardSensitivityModelEvaluator_decl_hpp
10#define Tempus_CombinedForwardSensitivityModelEvaluator_decl_hpp
12#include "Tempus_config.hpp"
14#include "Thyra_StateFuncModelEvaluatorBase.hpp"
15#include "Thyra_DefaultMultiVectorProductVectorSpace.hpp"
16#include "Thyra_DefaultMultiVectorProductVector.hpp"
38template <
typename Scalar>
40 public Thyra::StateFuncModelEvaluatorBase<Scalar>,
77 const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null,
78 const Teuchos::RCP<MultiVector>& dxdp_init = Teuchos::null,
79 const Teuchos::RCP<MultiVector>& dx_dotdp_init = Teuchos::null,
80 const Teuchos::RCP<MultiVector>& dx_dotdot_dp_init = Teuchos::null);
97 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_p_space(
int p)
const;
99 Teuchos::RCP<const Teuchos::Array<std::string> >
get_p_names(
int p)
const;
101 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_x_space()
const;
103 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_f_space()
const;
105 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_g_space(
int j)
const;
107 Teuchos::ArrayView<const std::string>
get_g_names(
int j)
const;
109 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
create_W_op()
const;
113 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
create_DgDx_op(
int j)
const;
115 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
create_DgDp_op(
int j,
int l)
const;
117 Teuchos::RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> >
120 Thyra::ModelEvaluatorBase::InArgs<Scalar>
createInArgs()
const;
130 typedef Thyra::DefaultMultiVectorProductVectorSpace<Scalar>
DMVPVS;
135 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
136 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs)
const;
142 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
model_;
165 mutable Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
my_dfdx_;
168 mutable Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
my_dgdx_;
169 mutable Teuchos::RCP<Thyra::MultiVectorBase<Scalar> >
my_dgdx_mv_;
Transform a ModelEvaluator's sensitivity equations to its residual.
Thyra::MultiVectorBase< Scalar > MultiVector
bool use_dgdp_as_tangent_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getForwardModel() const
Get the underlying model 'f'.
Thyra::ModelEvaluatorBase::InArgs< Scalar > prototypeInArgs_
Teuchos::RCP< const DMVPVS > x_dxdp_space_
Teuchos::ArrayView< const std::string > get_g_names(int j) const
Teuchos::RCP< MultiVector > dx_dotdotdp_init_
Teuchos::RCP< MultiVector > dx_dotdp_init_
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int p) const
Teuchos::RCP< const DMVPVS > dgdp_space_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > sens_residual_model_
Thyra::ModelEvaluatorBase::OutArgs< Scalar > prototypeOutArgs_
Teuchos::RCP< MultiVector > dxdp_init_
Teuchos::RCP< const DMVPVS > f_dfdp_space_
Teuchos::RCP< const DMVPVS > dfdp_space_
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int p) const
Teuchos::RCP< const DMVPVS > dxdp_space_
Teuchos::RCP< Thyra::MultiVectorBase< Scalar > > my_dgdx_mv_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdxdotdot_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_g_space(int j) const
static Teuchos::RCP< const Teuchos::ParameterList > getValidParameters()
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdx_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_DgDx_op(int j) const
Thyra::DefaultMultiVectorProductVectorSpace< Scalar > DMVPVS
Thyra::VectorBase< Scalar > Vector
bool use_dfdp_as_tangent_
int xdotdot_tangent_index_
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_DgDp_op(int j, int l) const
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > sens_solve_model_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_DgDx_dot_op(int j) const
int getSensitivityParamIndex() const
Get sensitivity parameter index.
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdxdot_
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dgdx_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > model_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
A ModelEvaluator decorator for sensitivity analysis.