9#ifndef Tempus_StaggeredForwardSensitivityModelEvaluator_decl_hpp
10#define Tempus_StaggeredForwardSensitivityModelEvaluator_decl_hpp
12#include "Thyra_StateFuncModelEvaluatorBase.hpp"
13#include "Thyra_DefaultMultiVectorProductVectorSpace.hpp"
16#include "Tempus_config.hpp"
17#include "Tempus_SolutionHistory.hpp"
48template <
typename Scalar>
50 public Thyra::StateFuncModelEvaluatorBase<Scalar>,
87 const bool is_pseudotransient,
88 const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null,
89 const Teuchos::RCP<MultiVector>& dxdp_init = Teuchos::null,
90 const Teuchos::RCP<MultiVector>& dx_dotdp_init = Teuchos::null,
91 const Teuchos::RCP<MultiVector>& dx_dotdot_dp_init = Teuchos::null);
110 const Teuchos::RCP<Thyra::NonlinearSolverBase<Scalar> >& solver,
111 const bool force_W_update) {
113 Teuchos::rcp_dynamic_cast<Thyra::NOXNonlinearSolver>(solver,
true);
114 lo_ = nox_solver->get_nonconst_W_op(force_W_update);
115 po_ = nox_solver->get_nonconst_prec_op();
123 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_p_space(
int p)
const;
125 Teuchos::RCP<const Teuchos::Array<std::string> >
get_p_names(
int p)
const;
127 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_x_space()
const;
129 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_f_space()
const;
131 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_g_space(
int j)
const;
133 Teuchos::ArrayView<const std::string>
get_g_names(
int j)
const;
135 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
create_W_op()
const;
139 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
create_DgDx_op(
int j)
const;
141 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
create_DgDp_op(
int j,
int l)
const;
143 Teuchos::RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> >
146 Thyra::ModelEvaluatorBase::InArgs<Scalar>
createInArgs()
const;
156 typedef Thyra::DefaultMultiVectorProductVectorSpace<Scalar>
DMVPVS;
161 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
162 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs)
const;
168 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
model_;
188 Teuchos::RCP<const Tempus::SolutionHistory<Scalar> >
sh_;
190 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
lo_;
191 Teuchos::RCP<Thyra::PreconditionerBase<Scalar> >
po_;
198 mutable Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
my_dfdx_;
201 mutable Teuchos::RCP<Thyra::MultiVectorBase<Scalar> >
my_dfdp_;
202 mutable Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
my_dgdx_;
203 mutable Teuchos::RCP<Thyra::MultiVectorBase<Scalar> >
my_dgdx_mv_;
A ModelEvaluator decorator for sensitivity analysis.
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
Solution state for integrators and steppers. SolutionState contains the metadata for solutions and th...
Transform a ModelEvaluator's sensitivity equations to its residual.
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
Teuchos::RCP< const DMVPVS > dxdp_space_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_DgDx_op(int j) const
Teuchos::RCP< Thyra::MultiVectorBase< Scalar > > my_dfdp_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > sens_solve_model_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
Teuchos::RCP< MultiVector > dx_dotdotdp_init_
static Teuchos::RCP< const Teuchos::ParameterList > getValidParameters()
Thyra::DefaultMultiVectorProductVectorSpace< Scalar > DMVPVS
Teuchos::RCP< MultiVector > dxdp_init_
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
Teuchos::ArrayView< const std::string > get_g_names(int j) const
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdxdotdot_
Teuchos::RCP< Thyra::MultiVectorBase< Scalar > > my_dgdx_mv_
Thyra::ModelEvaluatorBase::InArgs< Scalar > prototypeInArgs_
Teuchos::RCP< Tempus::SolutionState< Scalar > > nc_forward_state_
bool residual_sensitivity_is_computed_
Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > sh_
Teuchos::RCP< const DMVPVS > dgdp_space_
bool mass_matrix_is_computed_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_DgDx_dot_op(int j) const
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > sens_residual_model_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int p) const
Thyra::MultiVectorBase< Scalar > MultiVector
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getForwardModel() const
Get the underlying model 'f'.
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > model_
Teuchos::RCP< MultiVector > dx_dotdp_init_
Teuchos::RCP< const DMVPVS > dfdp_space_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
int xdotdot_tangent_index_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_DgDp_op(int j, int l) const
void setForwardSolutionHistory(const Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > &sh)
Set solution history from forward state evaluation (for interpolation)
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdxdot_
Thyra::ModelEvaluatorBase::OutArgs< Scalar > prototypeOutArgs_
bool use_dfdp_as_tangent_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dgdx_
Teuchos::RCP< Thyra::PreconditionerBase< Scalar > > po_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdx_
virtual void setSolver(const Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > &solver, const bool force_W_update)
Set the solver of the underlying model if you want to reuse it.
bool jacobian_matrix_is_computed_
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > lo_
Teuchos::RCP< const Tempus::SolutionState< Scalar > > forward_state_
Thyra::VectorBase< Scalar > Vector
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int p) const
bool use_dgdp_as_tangent_
bool acceleration_matrix_is_computed_
virtual void setForwardSolutionState(const Teuchos::RCP< const Tempus::SolutionState< Scalar > > &s)
Set solution state from forward state evaluation (for frozen state)
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_g_space(int j) const