44 includeXDot_ = includeXDot;
45 isInitialized_ =
false;
47 xDotIC_ = Scalar(lambda);
54 acceptModelParams_ =
false;
57 x_space_ = Thyra::defaultSpmdVectorSpace<Scalar>(dim_);
58 f_space_ = Thyra::defaultSpmdVectorSpace<Scalar>(dim_);
61 typedef Thyra::ModelEvaluatorBase MEB;
64 MEB::InArgsSetup<Scalar> inArgs;
65 inArgs.setModelEvalDescription(this->description());
66 inArgs.setSupports( MEB::IN_ARG_t );
67 inArgs.setSupports( MEB::IN_ARG_x );
68 inArgs.setSupports( MEB::IN_ARG_x_dot );
70 inArgs.setSupports( MEB::IN_ARG_beta );
71 inArgs.setSupports( MEB::IN_ARG_alpha );
72 if (acceptModelParams_) inArgs.set_Np( Np_ );
79 MEB::OutArgsSetup<Scalar> outArgs;
80 outArgs.setModelEvalDescription(this->description());
81 outArgs.setSupports( MEB::OUT_ARG_f );
83 outArgs.setSupports( MEB::OUT_ARG_W_op );
84 if (acceptModelParams_) {
85 outArgs.set_Np_Ng(Np_,Ng_);
86 outArgs.setSupports( MEB::OUT_ARG_DfDp,0,
87 MEB::DERIV_MV_JACOBIAN_FORM );
88 outArgs.setSupports( MEB::OUT_ARG_DgDp,0,0,
89 MEB::DERIV_MV_JACOBIAN_FORM );
90 outArgs.setSupports( MEB::OUT_ARG_DgDx,0,
91 MEB::DERIV_MV_GRADIENT_FORM );
97 nominalValues_ = inArgs_;
99 nominalValues_.set_t(Scalar(0.0));
100 const RCP<Thyra::VectorBase<Scalar> > x_ic = createMember(x_space_);
102 Thyra::DetachedVectorView<Scalar> x_ic_view( *x_ic );
105 nominalValues_.set_x(x_ic);
109 const RCP<Thyra::VectorBase<Scalar> > x_dot_ic = createMember(x_space_);
111 Thyra::DetachedVectorView<Scalar> x_dot_ic_view( *x_dot_ic );
112 x_dot_ic_view[0] = xDotIC_;
114 nominalValues_.set_x_dot(x_dot_ic);
117 isInitialized_ =
true;
126 Thyra::ModelEvaluatorBase::InArgs<Scalar> inArgs = inArgs_;
128 inArgs.set_t(exact_t);
131 Teuchos::RCP<Thyra::VectorBase<Scalar> > exact_x = createMember(x_space_);
133 Thyra::DetachedVectorView<Scalar> exact_x_view(*exact_x);
134 exact_x_view[0] = exp(lambda_*exact_t);
136 inArgs.set_x(exact_x);
140 Teuchos::RCP<Thyra::VectorBase<Scalar> > exact_x_dot = createMember(x_space_);
142 Thyra::DetachedVectorView<Scalar> exact_x_dot_view(*exact_x_dot);
143 exact_x_dot_view[0] = lambda_ * exp(lambda_*exact_t);
145 inArgs.set_x_dot(exact_x_dot);
201 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
202 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs
208 using Teuchos::rcp_dynamic_cast;
209 TEUCHOS_TEST_FOR_EXCEPTION( !isInitialized_, std::logic_error,
210 "Error, setupInOutArgs_ must be called first!\n");
212 const RCP<const VectorBase<Scalar> > x_in = inArgs.get_x().assert_not_null();
213 Thyra::ConstDetachedVectorView<Scalar> x_in_view( *x_in );
214 const RCP<VectorBase<Scalar> > f_out = outArgs.get_f();
215 const RCP<Thyra::LinearOpBase<Scalar> > W_out = outArgs.get_W_op();
217 if (inArgs.get_x_dot().is_null()) {
220 if (!is_null(f_out)) {
221 Thyra::DetachedVectorView<Scalar> f_out_view( *f_out );
222 f_out_view[0] = lambda_*x_in_view[0];
224 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
225 "Error -- Dahlquist Test Model requires f_out!\n");
228 if (!is_null(W_out)) {
229 RCP<Thyra::MultiVectorBase<Scalar> > matrix =
230 Teuchos::rcp_dynamic_cast<Thyra::MultiVectorBase<Scalar> >(W_out,
true);
231 Thyra::DetachedMultiVectorView<Scalar> matrix_view( *matrix );
232 matrix_view(0,0) = lambda_;
238 RCP<const VectorBase<Scalar> > x_dot_in;
239 x_dot_in = inArgs.get_x_dot().assert_not_null();
240 Scalar alpha = inArgs.get_alpha();
241 Scalar beta = inArgs.get_beta();
243 if (!is_null(f_out)) {
244 Thyra::DetachedVectorView<Scalar> f_out_view( *f_out );
245 Thyra::ConstDetachedVectorView<Scalar> x_dot_in_view( *x_dot_in );
246 f_out_view[0] = x_dot_in_view[0] - lambda_*x_in_view[0];
249 if (!is_null(W_out)) {
250 RCP<Thyra::MultiVectorBase<Scalar> > matrix =
251 Teuchos::rcp_dynamic_cast<Thyra::MultiVectorBase<Scalar> >(W_out,
true);
252 Thyra::DetachedMultiVectorView<Scalar> matrix_view( *matrix );
253 matrix_view(0,0) = alpha - beta*lambda_;
267 RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> > W_factory = this->get_W_factory();
268 RCP<Thyra::LinearOpBase<Scalar> > matrix = this->create_W_op();
270 RCP<Thyra::MultiVectorBase<Scalar> > multivec = Teuchos::rcp_dynamic_cast<Thyra::MultiVectorBase<Scalar> >(matrix,
true);
272 RCP<Thyra::VectorBase<Scalar> > vec = Thyra::createMember(x_space_);
274 Thyra::DetachedVectorView<Scalar> vec_view( *vec );
275 vec_view[0] = lambda_;
277 V_V(multivec->col(0).ptr(),*vec);
280 RCP<Thyra::LinearOpWithSolveBase<Scalar> > W =
281 Thyra::linearOpWithSolve<Scalar>(