29#ifndef RYTHMOS_FORWARD_SENSITIVITY_MODEL_EVALUATOR_BASE_HPP
30#define RYTHMOS_FORWARD_SENSITIVITY_MODEL_EVALUATOR_BASE_HPP
33#include "Rythmos_IntegratorBase.hpp"
34#include "Thyra_ModelEvaluator.hpp"
35#include "Thyra_StateFuncModelEvaluatorBase.hpp"
36#include "Thyra_DefaultProductVectorSpace.hpp"
37#include "Thyra_PhysicallyBlockedLinearOpWithSolveBase.hpp"
38#include "Thyra_DefaultBlockedTriangularLinearOpWithSolve.hpp"
39#include "Thyra_ModelEvaluatorDelegatorBase.hpp"
40#include "Thyra_ModelEvaluatorHelpers.hpp"
41#include "Thyra_DefaultMultiVectorProductVectorSpace.hpp"
42#include "Thyra_DefaultMultiVectorProductVector.hpp"
43#include "Thyra_DefaultMultiVectorLinearOpWithSolve.hpp"
44#include "Teuchos_implicit_cast.hpp"
45#include "Teuchos_Assert.hpp"
62 :
virtual public Thyra::StateFuncModelEvaluatorBase<Scalar>
81 const RCP<
const Thyra::ModelEvaluator<Scalar> > &stateModel,
101 const RCP<
const Thyra::ModelEvaluator<Scalar> >& stateModel,
102 const RCP<
const Thyra::VectorSpaceBase<Scalar> >& p_space
106 virtual RCP<const Thyra::ModelEvaluator<Scalar> >
110 virtual RCP<Thyra::ModelEvaluator<Scalar> >
117 virtual RCP<const Thyra::DefaultMultiVectorProductVectorSpace<Scalar> >
125 Ptr<StepperBase<Scalar> > stateStepper,
135template<
class Scalar>
136RCP<const Thyra::VectorBase<Scalar> > create_s_bar_given_S(
137 const ForwardSensitivityModelEvaluatorBase<Scalar> &fwdSensModel,
138 const RCP<Thyra::MultiVectorBase<Scalar> > &S
141 return Thyra::multiVectorProductVector(fwdSensModel.get_s_bar_space(), S);
148template<
class Scalar>
149RCP<const Thyra::VectorBase<Scalar> > create_s_bar_given_S(
150 const ForwardSensitivityModelEvaluatorBase<Scalar> &fwdSensModel,
151 const RCP<
const Thyra::MultiVectorBase<Scalar> > &S
154 return Thyra::multiVectorProductVector(fwdSensModel.get_s_bar_space(), S);
Forward sensitivity transient ModelEvaluator node interface class.
virtual RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_sens_space() const =0
virtual int get_p_index() const =0
virtual RCP< const Thyra::DefaultMultiVectorProductVectorSpace< Scalar > > get_s_bar_space() const =0
virtual void initializePointState(Ptr< StepperBase< Scalar > > stateStepper, bool forceUpToDateW)=0
virtual RCP< const Thyra::ModelEvaluator< Scalar > > getStateModel() const =0
virtual RCP< Thyra::ModelEvaluator< Scalar > > getNonconstStateModel() const =0
virtual void initializeStructureInitCondOnly(const RCP< const Thyra::ModelEvaluator< Scalar > > &stateModel, const RCP< const Thyra::VectorSpaceBase< Scalar > > &p_space)=0
Initialize the structure of the model for an initial condition only sensitivity problem.
virtual void initializeStructure(const RCP< const Thyra::ModelEvaluator< Scalar > > &stateModel, const int p_index)=0
Initialize the structure of the model.