42 #ifndef THYRA_EPETRA_MODEL_EVALUATOR_HPP 43 #define THYRA_EPETRA_MODEL_EVALUATOR_HPP 45 #include "Thyra_ModelEvaluatorDefaultBase.hpp" 47 #include "Thyra_LinearOpWithSolveFactoryBase.hpp" 48 #include "EpetraExt_ModelEvaluator.h" 49 #include "Epetra_Map.h" 50 #include "Teuchos_Array.hpp" 176 :
public ModelEvaluatorDefaultBase<double>,
177 virtual public Teuchos::ParameterListAcceptor
189 const RCP<const EpetraExt::ModelEvaluator> &epetraModel,
190 const RCP<LinearOpWithSolveFactoryBase<double> > &W_factory
195 const RCP<const EpetraExt::ModelEvaluator> &epetraModel,
196 const RCP<LinearOpWithSolveFactoryBase<double> > &W_factory
207 void setNominalValues(
const ModelEvaluatorBase::InArgs<double>& nominalValues );
217 const RCP<const Epetra_Vector> &stateVariableScalingVec
222 RCP<const Epetra_Vector>
227 RCP<const Epetra_Vector>
233 const RCP<const Epetra_Vector> &stateFunctionScalingVec
238 RCP<const Epetra_Vector>
243 RCP<const EpetraExt::ModelEvaluator> *epetraModel = NULL,
244 RCP<LinearOpWithSolveFactoryBase<double> > *W_factory = NULL
248 const ModelEvaluatorBase::InArgs<double>&
getFinalPoint()
const;
287 RCP<const VectorSpaceBase<double> >
get_x_space()
const;
289 RCP<const VectorSpaceBase<double> >
get_f_space()
const;
291 RCP<const VectorSpaceBase<double> >
get_p_space(
int l)
const;
293 RCP<const Teuchos::Array<std::string> >
get_p_names(
int l)
const;
295 RCP<const VectorSpaceBase<double> >
get_g_space(
int j)
const;
297 Teuchos::ArrayView<const std::string>
get_g_names(
int j)
const;
309 RCP<const LinearOpWithSolveFactoryBase<double> >
get_W_factory()
const;
311 ModelEvaluatorBase::InArgs<double>
createInArgs()
const;
314 const ModelEvaluatorBase::InArgs<double> &finalPoint
315 ,
const bool wasSolved
341 const ModelEvaluatorBase::InArgs<double> &inArgs,
342 const ModelEvaluatorBase::OutArgs<double> &outArgs
352 typedef Teuchos::Array<RCP<const Epetra_Map> >
p_map_t;
353 typedef Teuchos::Array<RCP<const Epetra_Map> >
g_map_t;
357 typedef Teuchos::Array<RCP<const VectorSpaceBase<double> > >
359 typedef Teuchos::Array<RCP<const VectorSpaceBase<double> > >
412 const EpetraExt::ModelEvaluator::InArgs &epetraInArgs,
413 ModelEvaluatorBase::InArgs<double> *inArgs
418 const ModelEvaluatorBase::InArgs<double> &inArgs,
419 EpetraExt::ModelEvaluator::InArgs *epetraInArgs
425 const ModelEvaluatorBase::OutArgs<double> &outArgs,
427 EpetraExt::ModelEvaluator::OutArgs *epetraUnscaledOutArgs,
429 RCP<LinearOpBase<double> > *W_op,
430 RCP<EpetraLinearOp> *efwdW,
432 RCP<Epetra_Operator> *eW
437 EpetraExt::ModelEvaluator::InArgs *epetraInArgs,
438 EpetraExt::ModelEvaluator::OutArgs *epetraUnscaledOutArgs,
439 const RCP<Teuchos::FancyOStream> &out,
440 const Teuchos::EVerbosityLevel verbLevel
445 const EpetraExt::ModelEvaluator::OutArgs &epetraUnscaledOutArgs,
446 const RCP<Teuchos::FancyOStream> &out,
447 const Teuchos::EVerbosityLevel verbLevel
452 const EpetraExt::ModelEvaluator::OutArgs &epetraOutArgs,
453 RCP<LinearOpBase<double> > &W_op,
454 RCP<EpetraLinearOp> &efwdW,
455 RCP<Epetra_Operator> &eW,
456 const ModelEvaluatorBase::OutArgs<double> &outArgs
482 RCP<EpetraModelEvaluator>
483 epetraModelEvaluator(
484 const RCP<const EpetraExt::ModelEvaluator> &epetraModel,
485 const RCP<LinearOpWithSolveFactoryBase<double> > &W_factory
492 ModelEvaluatorBase::EDerivativeMultiVectorOrientation
493 convert(
const EpetraExt::ModelEvaluator::EDerivativeMultiVectorOrientation &mvOrientation );
499 EpetraExt::ModelEvaluator::EDerivativeMultiVectorOrientation
500 convert(
const ModelEvaluatorBase::EDerivativeMultiVectorOrientation &mvOrientation );
506 ModelEvaluatorBase::DerivativeProperties
507 convert(
const EpetraExt::ModelEvaluator::DerivativeProperties &derivativeProperties );
513 ModelEvaluatorBase::DerivativeSupport
514 convert(
const EpetraExt::ModelEvaluator::DerivativeSupport &derivativeSupport );
520 EpetraExt::ModelEvaluator::Derivative
522 const ModelEvaluatorBase::Derivative<double> &derivative,
523 const RCP<const Epetra_Map> &fnc_map,
524 const RCP<const Epetra_Map> &var_map
527 EpetraExt::ModelEvaluator::MPDerivative
529 const ModelEvaluatorBase::MPDerivative &derivative,
530 const RCP<const Epetra_Map> &fnc_map,
531 const RCP<const Epetra_Map> &var_map
537 #endif // THYRA_EPETRA_MODEL_EVALUATOR_HPP
void preEvalScalingSetup(EpetraExt::ModelEvaluator::InArgs *epetraInArgs, EpetraExt::ModelEvaluator::OutArgs *epetraUnscaledOutArgs, const RCP< Teuchos::FancyOStream > &out, const Teuchos::EVerbosityLevel verbLevel) const
void setNominalValues(const ModelEvaluatorBase::InArgs< double > &nominalValues)
Set the nominal values.
RCP< const EpetraExt::ModelEvaluator > getEpetraModel() const
bool finalPointWasSolved() const
RCP< const VectorSpaceBase< double > > get_g_space(int j) const
RCP< Epetra_Vector > x_dot_unscaled_
RCP< const Epetra_Map > x_map_
RCP< EpetraLinearOp > create_epetra_W_op() const
void reportFinalPoint(const ModelEvaluatorBase::InArgs< double > &finalPoint, const bool wasSolved)
RCP< LinearOpBase< double > > create_W_op() const
void evalModelImpl(const ModelEvaluatorBase::InArgs< double > &inArgs, const ModelEvaluatorBase::OutArgs< double > &outArgs) const
void updateInArgsOutArgs() const
RCP< const VectorSpaceBase< double > > get_x_space() const
RCP< Teuchos::ParameterList > getNonconstParameterList()
p_map_is_local_t g_map_is_local_
RCP< Epetra_Vector > x_unscaled_
Teuchos::Array< RCP< const Epetra_Map > > p_map_t
RCP< LinearOpWithSolveFactoryBase< double > > W_factory_
EStateFunctionScaling stateFunctionScaling_
std::string description() const
Teuchos::Array< RCP< const Epetra_Map > > g_map_t
void setStateVariableScalingVec(const RCP< const Epetra_Vector > &stateVariableScalingVec)
Set the state variable scaling vector s_x (see above).
void convertInArgsFromThyraToEpetra(const ModelEvaluatorBase::InArgs< double > &inArgs, EpetraExt::ModelEvaluator::InArgs *epetraInArgs) const
RCP< const Teuchos::Array< std::string > > get_p_names(int l) const
RCP< const Teuchos::ParameterList > getParameterList() const
ModelEvaluatorBase::OutArgs< double > createOutArgsImpl() const
const ModelEvaluatorBase::InArgs< double > & getFinalPoint() const
Teuchos::Array< RCP< const VectorSpaceBase< double > > > g_space_t
void convertOutArgsFromThyraToEpetra(const ModelEvaluatorBase::OutArgs< double > &outArgs, EpetraExt::ModelEvaluator::OutArgs *epetraUnscaledOutArgs, RCP< LinearOpBase< double > > *W_op, RCP< EpetraLinearOp > *efwdW, RCP< Epetra_Operator > *eW) const
EpetraExt::ModelEvaluator::OutArgs epetraOutArgsScaling_
RCP< const VectorSpaceBase< double > > x_space_
ModelEvaluatorBase::OutArgs< double > prototypeOutArgs_
RCP< const Epetra_Vector > getStateVariableInvScalingVec() const
Get the state variable scaling vector s_x (see above).
ModelEvaluatorBase::InArgs< double > getNominalValues() const
RCP< LinearOpBase< double > > create_DgDx_dot_op_impl(int j) const
EpetraExt::ModelEvaluator::MPDerivative convert(const ModelEvaluatorBase::MPDerivative &derivative, const RCP< const Epetra_Map > &fnc_map, const RCP< const Epetra_Map > &var_map)
void finishConvertingOutArgsFromEpetraToThyra(const EpetraExt::ModelEvaluator::OutArgs &epetraOutArgs, RCP< LinearOpBase< double > > &W_op, RCP< EpetraLinearOp > &efwdW, RCP< Epetra_Operator > &eW, const ModelEvaluatorBase::OutArgs< double > &outArgs) const
RCP< const VectorSpaceBase< double > > f_space_
ModelEvaluatorBase::InArgs< double > finalPoint_
void setStateFunctionScalingVec(const RCP< const Epetra_Vector > &stateFunctionScalingVec)
Set the state function scaling vector s_f (see above).
void postEvalScalingSetup(const EpetraExt::ModelEvaluator::OutArgs &epetraUnscaledOutArgs, const RCP< Teuchos::FancyOStream > &out, const Teuchos::EVerbosityLevel verbLevel) const
ModelEvaluatorBase::InArgs< double > lowerBounds_
void initialize(const RCP< const EpetraExt::ModelEvaluator > &epetraModel, const RCP< LinearOpWithSolveFactoryBase< double > > &W_factory)
ModelEvaluatorBase::InArgs< double > getLowerBounds() const
Teuchos::Array< RCP< const VectorSpaceBase< double > > > p_space_t
EpetraExt::ModelEvaluator::InArgs epetraInArgsScaling_
bool currentInArgsOutArgs_
Concrete Adapter subclass that takes an EpetraExt::ModelEvaluator object and wraps it as a Thyra::Mod...
bool finalPointWasSolved_
void setParameterList(RCP< Teuchos::ParameterList > const ¶mList)
void updateNominalValuesAndBounds() const
RCP< const Epetra_Vector > stateFunctionScalingVec_
void uninitialize(RCP< const EpetraExt::ModelEvaluator > *epetraModel=NULL, RCP< LinearOpWithSolveFactoryBase< double > > *W_factory=NULL)
bool nominalValuesAndBoundsAreUpdated_
RCP< LinearOpBase< double > > create_DfDp_op_impl(int l) const
void convertInArgsFromEpetraToThyra(const EpetraExt::ModelEvaluator::InArgs &epetraInArgs, ModelEvaluatorBase::InArgs< double > *inArgs) const
ModelEvaluatorBase::InArgs< double > prototypeInArgs_
ModelEvaluatorBase::InArgs< double > getUpperBounds() const
RCP< const Epetra_Map > f_map_
RCP< LinearOpBase< double > > create_DgDp_op_impl(int j, int l) const
RCP< const Epetra_Vector > stateVariableScalingVec_
std::vector< bool > g_map_is_local_t
RCP< const EpetraExt::ModelEvaluator > epetraModel_
std::vector< bool > p_map_is_local_t
RCP< Teuchos::ParameterList > paramList_
ModelEvaluatorBase::InArgs< double > nominalValues_
RCP< PreconditionerBase< double > > create_W_prec() const
Returns null currently.
p_map_is_local_t p_map_is_local_
RCP< const Epetra_Vector > getStateFunctionScalingVec() const
Get the state function scaling vector s_f (see above).
RCP< const LinearOpWithSolveFactoryBase< double > > get_W_factory() const
Teuchos::ArrayView< const std::string > get_g_names(int j) const
RCP< LinearOpBase< double > > create_DgDx_op_impl(int j) const
ModelEvaluatorBase::InArgs< double > createInArgs() const
RCP< const VectorSpaceBase< double > > get_f_space() const
RCP< const VectorSpaceBase< double > > get_p_space(int l) const
RCP< const Epetra_Vector > getStateVariableScalingVec() const
Get the inverse state variable scaling vector inv_s_x (see above).
RCP< Teuchos::ParameterList > unsetParameterList()
RCP< const Teuchos::ParameterList > getValidParameters() const
RCP< const Epetra_Vector > invStateVariableScalingVec_
ModelEvaluatorBase::InArgs< double > upperBounds_