9#ifndef Tempus_ModelEvaluatorPairPartIMEX_StaggeredFSA_decl_hpp
10#define Tempus_ModelEvaluatorPairPartIMEX_StaggeredFSA_decl_hpp
12#include "Tempus_config.hpp"
14#include "Tempus_WrapperModelEvaluatorPairPartIMEX_Basic.hpp"
15#include "Tempus_StaggeredForwardSensitivityModelEvaluator.hpp"
17#include "Thyra_ProductMultiVectorBase.hpp"
18#include "Thyra_DefaultMultiVectorProductVectorSpace.hpp"
19#include "Thyra_DefaultMultiVectorProductVector.hpp"
31template <
typename Scalar>
41 const bool is_pseudotransient,
42 const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null);
53 virtual Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
61 virtual Teuchos::RCP<Thyra::VectorBase<Scalar> >
66 virtual Teuchos::RCP<const Thyra::VectorBase<Scalar> >
71 virtual Teuchos::RCP<Thyra::VectorBase<Scalar> >
76 virtual Teuchos::RCP<const Thyra::VectorBase<Scalar> >
86 virtual Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
99 const Teuchos::RCP<Thyra::NonlinearSolverBase<Scalar> >& solver,
100 const bool force_W_update);
107 virtual Thyra::ModelEvaluatorBase::InArgs<Scalar>
112 const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
113 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> & outArgs)
const;
125 typedef Thyra::DefaultMultiVectorProductVectorSpace<Scalar>
DMVPVS;
126 typedef Thyra::DefaultMultiVectorProductVector<Scalar>
DMVPV;
127 typedef Thyra::ProductMultiVectorBase<Scalar>
PMVB;
130 Teuchos::RCP<const WrapperModelEvaluatorPairPartIMEX_Basic<Scalar> >
forwardModel_;
142 Teuchos::RCP<const Tempus::SolutionHistory<Scalar> >
sh_;
149 mutable Teuchos::RCP<Thyra::MultiVectorBase<Scalar> >
my_dfdp_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.
ModelEvaluator pair for implicit and explicit (IMEX) evaulations.
Specialization of IMEX-Part ME for "combined" FSA method.
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters() const
Teuchos::RCP< const DMVPVS > imex_dxdp_prod_space_
Teuchos::RCP< const Tempus::SolutionState< Scalar > > implicit_x_state_
Teuchos::RCP< const Tempus::SolutionState< Scalar > > forward_state_
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::ModelEvaluator< Scalar > > appExplicitModel_
Teuchos::RCP< Thyra::MultiVectorBase< Scalar > > my_dfdp_mv_
Teuchos::RCP< Tempus::SolutionState< Scalar > > nc_forward_state_
Teuchos::RCP< FSAME > fsaExplicitModel_
virtual ~WrapperModelEvaluatorPairPartIMEX_StaggeredFSA()
Destructor.
virtual Teuchos::RCP< Thyra::VectorBase< Scalar > > getExplicitOnlyVector(const Teuchos::RCP< Thyra::VectorBase< Scalar > > &full) const
Extract explicit-only vector from a full solution vector.
Teuchos::RCP< const WrapperModelEvaluatorPairPartIMEX_Basic< Scalar > > forwardModel_
virtual void setForwardSolutionHistory(const Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > &sh)
Set solution history from forward state evaluation (for interpolation)
virtual void initialize()
Initialize after setting member data.
Teuchos::RCP< FSAME > fsaImplicitModel_
Teuchos::RCP< const Tempus::SolutionState< Scalar > > explicit_y_state_
virtual void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
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.
Thyra::DefaultMultiVectorProductVectorSpace< Scalar > DMVPVS
Thyra::DefaultMultiVectorProductVector< Scalar > DMVPV
virtual Teuchos::RCP< Thyra::VectorBase< Scalar > > getIMEXVector(const Teuchos::RCP< Thyra::VectorBase< Scalar > > &full) const
Extract IMEX vector from a full solution vector.
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > appImplicitModel_
WrapperModelEvaluatorPairPartIMEX_Basic< Scalar > Base
virtual Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getForwardModel() const
Get the underlying forward model.
virtual Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
virtual Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int i) const
Get the p space.
StaggeredForwardSensitivityModelEvaluator< Scalar > FSAME
Teuchos::RCP< const DMVPVS > explicit_dydp_prod_space_
void buildIMEXStates() const
Build implicit x and end explicit y states from forward_state_.
bool use_dfdp_as_tangent_
Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > sh_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > my_dfdp_op_
Thyra::ProductMultiVectorBase< Scalar > PMVB