29 const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null)
32 using Teuchos::rcp_dynamic_cast;
34 RCP<SensitivityModelEvaluatorBase<Scalar> > wrapped_model;
37 RCP<const WrapperModelEvaluatorPairIMEX_Basic<Scalar> > modelPairIMEX =
38 rcp_dynamic_cast<const WrapperModelEvaluatorPairIMEX_Basic<Scalar> >(model);
39 RCP<const WrapperModelEvaluatorPairPartIMEX_Basic<Scalar> > modelPairPartIMEX =
40 rcp_dynamic_cast<const WrapperModelEvaluatorPairPartIMEX_Basic<Scalar> >(model);
44 if ((modelPairIMEX != Teuchos::null || modelPairPartIMEX != Teuchos::null) &&
45 (model.ptr() != sens_residual_model.ptr()) &&
46 (model.ptr() != sens_solve_model.ptr()))
47 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
"Unique model evaluators for state, sensitivity residual, and sensitivity solve is not supported for IMEX");
49 if (modelPairIMEX != Teuchos::null) {
51 rcp(
new WrapperModelEvaluatorPairIMEX_CombinedFSA<Scalar>(
52 modelPairIMEX, pList));
54 else if (modelPairPartIMEX != Teuchos::null) {
56 rcp(
new WrapperModelEvaluatorPairPartIMEX_CombinedFSA<Scalar>(
57 modelPairPartIMEX, pList));
61 rcp(
new CombinedForwardSensitivityModelEvaluator<Scalar>(
62 model, sens_residual_model, sens_solve_model, pList));
74 const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null)
76 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> > cmodel = model;
77 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> > csens_residual_model = sens_residual_model;
78 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> > csens_solve_model = sens_solve_model;
80 cmodel, csens_residual_model, csens_solve_model, pList);
Teuchos::RCP< SensitivityModelEvaluatorBase< Scalar > > wrapCombinedFSAModelEvaluator(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &model, const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &sens_residual_model, const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &sens_solve_model, const Teuchos::RCP< const Teuchos::ParameterList > &pList=Teuchos::null)