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" 60 template<
class Scalar>
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> >
121 virtual RCP<const Thyra::VectorSpaceBase<Scalar> >
get_p_sens_space()
const = 0;
135 template<
class Scalar>
136 RCP<const Thyra::VectorBase<Scalar> > create_s_bar_given_S(
138 const RCP<Thyra::MultiVectorBase<Scalar> > &S
141 return Thyra::multiVectorProductVector(fwdSensModel.
get_s_bar_space(), S);
148 template<
class Scalar>
149 RCP<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);
161 #endif // RYTHMOS_FORWARD_SENSITIVITY_MODEL_EVALUATOR_BASE_HPP Base class for defining stepper functionality.
virtual int get_p_index() const =0
virtual void initializePointState(Ptr< StepperBase< Scalar > > stateStepper, bool forceUpToDateW)=0
virtual void initializeStructure(const RCP< const Thyra::ModelEvaluator< Scalar > > &stateModel, const int p_index)=0
Initialize the structure of the model.
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.
Forward sensitivity transient ModelEvaluator node interface class.
virtual RCP< const Thyra::ModelEvaluator< Scalar > > getStateModel() const =0
virtual RCP< const Thyra::DefaultMultiVectorProductVectorSpace< Scalar > > get_s_bar_space() const =0
virtual RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_sens_space() const =0