9#ifndef TEMPUS_TEST_SINCOS_MODEL_IMPL_HPP
10#define TEMPUS_TEST_SINCOS_MODEL_IMPL_HPP
12#include "Teuchos_StandardParameterEntryValidators.hpp"
14#include "Thyra_DefaultSpmdVectorSpace.hpp"
15#include "Thyra_DetachedVectorView.hpp"
16#include "Thyra_DetachedMultiVectorView.hpp"
17#include "Thyra_DefaultSerialDenseLinearOpWithSolveFactory.hpp"
18#include "Thyra_DefaultMultiVectorLinearOpWithSolve.hpp"
19#include "Thyra_DefaultLinearOpSource.hpp"
20#include "Thyra_VectorStdOps.hpp"
21#include "Thyra_MultiVectorStdOps.hpp"
22#include "Thyra_DefaultMultiVectorProductVector.hpp"
31SinCosModel(Teuchos::RCP<Teuchos::ParameterList> pList_)
33 isInitialized_ =
false;
39 acceptModelParams_ =
false;
40 useDfDpAsTangent_ =
false;
50 x_space_ = Thyra::defaultSpmdVectorSpace<Scalar>(dim_);
51 f_space_ = Thyra::defaultSpmdVectorSpace<Scalar>(dim_);
53 p_space_ = Thyra::defaultSpmdVectorSpace<Scalar>(np_);
54 g_space_ = Thyra::defaultSpmdVectorSpace<Scalar>(ng_);
56 setParameterList(pList_);
59 DxDp_space_ = Thyra::multiVectorProductVectorSpace(x_space_, np_);
63Thyra::ModelEvaluatorBase::InArgs<Scalar>
67 TEUCHOS_TEST_FOR_EXCEPTION( !isInitialized_, std::logic_error,
68 "Error, setupInOutArgs_ must be called first!\n");
69 Thyra::ModelEvaluatorBase::InArgs<Scalar> inArgs = inArgs_;
71 inArgs.set_t(exact_t);
72 Teuchos::RCP<Thyra::VectorBase<Scalar> > exact_x = createMember(x_space_);
74 Thyra::DetachedVectorView<Scalar> exact_x_view(*exact_x);
75 exact_x_view[0] = a_+b_*sin((f_/L_)*t+phi_);
76 exact_x_view[1] = b_*(f_/L_)*cos((f_/L_)*t+phi_);
78 inArgs.set_x(exact_x);
79 Teuchos::RCP<Thyra::VectorBase<Scalar> > exact_x_dot = createMember(x_space_);
81 Thyra::DetachedVectorView<Scalar> exact_x_dot_view(*exact_x_dot);
82 exact_x_dot_view[0] = b_*(f_/L_)*cos((f_/L_)*t+phi_);
83 exact_x_dot_view[1] = -b_*(f_/L_)*(f_/L_)*sin((f_/L_)*t+phi_);
85 inArgs.set_x_dot(exact_x_dot);
96Thyra::ModelEvaluatorBase::InArgs<Scalar>
100 TEUCHOS_TEST_FOR_EXCEPTION( !isInitialized_, std::logic_error,
101 "Error, setupInOutArgs_ must be called first!\n");
102 Thyra::ModelEvaluatorBase::InArgs<Scalar> inArgs = inArgs_;
103 if (!acceptModelParams_) {
106 TEUCHOS_ASSERT_IN_RANGE_UPPER_EXCLUSIVE( j, 0, np_ );
112 inArgs.set_t(exact_t);
113 Teuchos::RCP<Thyra::VectorBase<Scalar> > exact_s = createMember(x_space_);
115 Thyra::DetachedVectorView<Scalar> exact_s_view(*exact_s);
117 exact_s_view[0] = 1.0;
118 exact_s_view[1] = 0.0;
120 exact_s_view[0] = (b/L)*t*cos((f/L)*t+phi);
121 exact_s_view[1] = (b/L)*cos((f/L)*t+phi)-(b*f*t/(L*L))*sin((f/L)*t+phi);
123 exact_s_view[0] = -(b*f*t/(L*L))*cos((f/L)*t+phi);
124 exact_s_view[1] = -(b*f/(L*L))*cos((f/L)*t+phi)+(b*f*f*t/(L*L*L))*sin((f/L)*t+phi);
127 inArgs.set_x(exact_s);
128 Teuchos::RCP<Thyra::VectorBase<Scalar> > exact_s_dot = createMember(x_space_);
130 Thyra::DetachedVectorView<Scalar> exact_s_dot_view(*exact_s_dot);
132 exact_s_dot_view[0] = 0.0;
133 exact_s_dot_view[1] = 0.0;
135 exact_s_dot_view[0] = (b/L)*cos((f/L)*t+phi)-(b*f*t/(L*L))*sin((f/L)*t+phi);
136 exact_s_dot_view[1] = -(2.0*b*f/(L*L))*sin((f/L)*t+phi)-(b*f*f*t/(L*L*L))*cos((f/L)*t+phi);
138 exact_s_dot_view[0] = -(b*f/(L*L))*cos((f/L)*t+phi)+(b*f*f*t/(L*L*L))*sin((f/L)*t+phi);
139 exact_s_dot_view[1] = (2.0*b*f*f/(L*L*L))*sin((f/L)*t+phi)+(b*f*f*f*t/(L*L*L*L))*cos((f/L)*t+phi);
142 inArgs.set_x_dot(exact_s_dot);
146template<
class Scalar>
147Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
155template<
class Scalar>
156Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
164template<
class Scalar>
165Thyra::ModelEvaluatorBase::InArgs<Scalar>
169 TEUCHOS_TEST_FOR_EXCEPTION( !isInitialized_, std::logic_error,
170 "Error, setupInOutArgs_ must be called first!\n");
171 return nominalValues_;
175template<
class Scalar>
176Teuchos::RCP<Thyra::LinearOpWithSolveBase<Scalar> >
181 RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> > W_factory = this->get_W_factory();
182 RCP<Thyra::LinearOpBase<Scalar> > matrix = this->create_W_op();
189 RCP<Thyra::MultiVectorBase<Scalar> > multivec = Teuchos::rcp_dynamic_cast<Thyra::MultiVectorBase<Scalar> >(matrix,
true);
191 RCP<Thyra::VectorBase<Scalar> > vec = Thyra::createMember(x_space_);
193 Thyra::DetachedVectorView<Scalar> vec_view( *vec );
197 V_V(multivec->col(0).ptr(),*vec);
199 Thyra::DetachedVectorView<Scalar> vec_view( *vec );
203 V_V(multivec->col(1).ptr(),*vec);
206 RCP<Thyra::LinearOpWithSolveBase<Scalar> > W =
207 Thyra::linearOpWithSolve<Scalar>(
221template<
class Scalar>
222Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
226 Teuchos::RCP<Thyra::MultiVectorBase<Scalar> > matrix = Thyra::createMembers(x_space_, dim_);
231template<
class Scalar>
232Teuchos::RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> >
236 Teuchos::RCP<Thyra::LinearOpWithSolveFactoryBase<Scalar> > W_factory =
237 Thyra::defaultSerialDenseLinearOpWithSolveFactory<Scalar>();
253template<
class Scalar>
254Thyra::ModelEvaluatorBase::InArgs<Scalar>
266template<
class Scalar>
267Thyra::ModelEvaluatorBase::OutArgs<Scalar>
276template<
class Scalar>
280 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
281 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs
284 typedef Thyra::DefaultMultiVectorProductVector<Scalar> DMVPV;
288 using Teuchos::rcp_dynamic_cast;
289 TEUCHOS_TEST_FOR_EXCEPTION( !isInitialized_, std::logic_error,
290 "Error, setupInOutArgs_ must be called first!\n");
292 const RCP<const VectorBase<Scalar> > x_in = inArgs.get_x().assert_not_null();
293 Thyra::ConstDetachedVectorView<Scalar> x_in_view( *x_in );
299 if (acceptModelParams_) {
300 const RCP<const VectorBase<Scalar> > p_in =
301 inArgs.get_p(0).assert_not_null();
302 Thyra::ConstDetachedVectorView<Scalar> p_in_view( *p_in );
308 RCP<const MultiVectorBase<Scalar> > DxDp_in, DxdotDp_in;
309 if (acceptModelParams_) {
310 if (inArgs.get_p(1) != Teuchos::null)
312 rcp_dynamic_cast<const DMVPV>(inArgs.get_p(1))->getMultiVector();
313 if (inArgs.get_p(2) != Teuchos::null)
315 rcp_dynamic_cast<const DMVPV>(inArgs.get_p(2))->getMultiVector();
318 Scalar beta = inArgs.get_beta();
320 const RCP<VectorBase<Scalar> > f_out = outArgs.get_f();
321 const RCP<Thyra::LinearOpBase<Scalar> > W_out = outArgs.get_W_op();
322 RCP<Thyra::MultiVectorBase<Scalar> > DfDp_out;
323 if (acceptModelParams_) {
324 Thyra::ModelEvaluatorBase::Derivative<Scalar> DfDp = outArgs.get_DfDp(0);
325 DfDp_out = DfDp.getMultiVector();
328 if (inArgs.get_x_dot().is_null()) {
331 if (!is_null(f_out)) {
332 Thyra::DetachedVectorView<Scalar> f_out_view( *f_out );
333 f_out_view[0] = x_in_view[1];
334 f_out_view[1] = (f/L)*(f/L)*(a-x_in_view[0]);
336 if (!is_null(W_out)) {
337 RCP<Thyra::MultiVectorBase<Scalar> > matrix =
338 Teuchos::rcp_dynamic_cast<Thyra::MultiVectorBase<Scalar> >(W_out,
true);
339 Thyra::DetachedMultiVectorView<Scalar> matrix_view( *matrix );
340 matrix_view(0,0) = 0.0;
341 matrix_view(0,1) = +beta;
342 matrix_view(1,0) = -beta*(f/L)*(f/L);
343 matrix_view(1,1) = 0.0;
346 if (!is_null(DfDp_out)) {
347 Thyra::DetachedMultiVectorView<Scalar> DfDp_out_view( *DfDp_out );
348 DfDp_out_view(0,0) = 0.0;
349 DfDp_out_view(0,1) = 0.0;
350 DfDp_out_view(0,2) = 0.0;
351 DfDp_out_view(1,0) = (f/L)*(f/L);
352 DfDp_out_view(1,1) = (2.0*f/(L*L))*(a-x_in_view[0]);
353 DfDp_out_view(1,2) = -(2.0*f*f/(L*L*L))*(a-x_in_view[0]);
356 if (useDfDpAsTangent_ && !is_null(DxDp_in)) {
357 Thyra::ConstDetachedMultiVectorView<Scalar> DxDp( *DxDp_in );
358 DfDp_out_view(0,0) += DxDp(1,0);
359 DfDp_out_view(0,1) += DxDp(1,1);
360 DfDp_out_view(0,2) += DxDp(1,2);
361 DfDp_out_view(1,0) += -(f/L)*(f/L) * DxDp(0,0);
362 DfDp_out_view(1,1) += -(f/L)*(f/L) * DxDp(0,1);
363 DfDp_out_view(1,2) += -(f/L)*(f/L) * DxDp(0,2);
369 RCP<const VectorBase<Scalar> > x_dot_in;
370 x_dot_in = inArgs.get_x_dot().assert_not_null();
371 Scalar alpha = inArgs.get_alpha();
372 if (!is_null(f_out)) {
373 Thyra::DetachedVectorView<Scalar> f_out_view( *f_out );
374 Thyra::ConstDetachedVectorView<Scalar> x_dot_in_view( *x_dot_in );
375 f_out_view[0] = x_dot_in_view[0] - x_in_view[1];
376 f_out_view[1] = x_dot_in_view[1] - (f/L)*(f/L)*(a-x_in_view[0]);
378 if (!is_null(W_out)) {
379 RCP<Thyra::MultiVectorBase<Scalar> > matrix =
380 Teuchos::rcp_dynamic_cast<Thyra::MultiVectorBase<Scalar> >(W_out,
true);
381 Thyra::DetachedMultiVectorView<Scalar> matrix_view( *matrix );
382 matrix_view(0,0) = alpha;
383 matrix_view(0,1) = -beta;
384 matrix_view(1,0) = +beta*(f/L)*(f/L);
385 matrix_view(1,1) = alpha;
388 if (!is_null(DfDp_out)) {
389 Thyra::DetachedMultiVectorView<Scalar> DfDp_out_view( *DfDp_out );
390 DfDp_out_view(0,0) = 0.0;
391 DfDp_out_view(0,1) = 0.0;
392 DfDp_out_view(0,2) = 0.0;
393 DfDp_out_view(1,0) = -(f/L)*(f/L);
394 DfDp_out_view(1,1) = -(2.0*f/(L*L))*(a-x_in_view[0]);
395 DfDp_out_view(1,2) = +(2.0*f*f/(L*L*L))*(a-x_in_view[0]);
398 if (useDfDpAsTangent_ && !is_null(DxdotDp_in)) {
399 Thyra::ConstDetachedMultiVectorView<Scalar> DxdotDp( *DxdotDp_in );
400 DfDp_out_view(0,0) += DxdotDp(0,0);
401 DfDp_out_view(0,1) += DxdotDp(0,1);
402 DfDp_out_view(0,2) += DxdotDp(0,2);
403 DfDp_out_view(1,0) += DxdotDp(1,0);
404 DfDp_out_view(1,1) += DxdotDp(1,1);
405 DfDp_out_view(1,2) += DxdotDp(1,2);
407 if (useDfDpAsTangent_ && !is_null(DxDp_in)) {
408 Thyra::ConstDetachedMultiVectorView<Scalar> DxDp( *DxDp_in );
409 DfDp_out_view(0,0) += -DxDp(1,0);
410 DfDp_out_view(0,1) += -DxDp(1,1);
411 DfDp_out_view(0,2) += -DxDp(1,2);
412 DfDp_out_view(1,0) += (f/L)*(f/L) * DxDp(0,0);
413 DfDp_out_view(1,1) += (f/L)*(f/L) * DxDp(0,1);
414 DfDp_out_view(1,2) += (f/L)*(f/L) * DxDp(0,2);
420 if (acceptModelParams_) {
421 RCP<VectorBase<Scalar> > g_out = outArgs.get_g(0);
422 if (g_out != Teuchos::null)
423 Thyra::assign(g_out.ptr(), *x_in);
425 RCP<Thyra::MultiVectorBase<Scalar> > DgDp_out =
426 outArgs.get_DgDp(0,0).getMultiVector();
427 if (DgDp_out != Teuchos::null)
428 Thyra::assign(DgDp_out.ptr(), Scalar(0.0));
430 RCP<Thyra::MultiVectorBase<Scalar> > DgDx_out =
431 outArgs.get_DgDx(0).getMultiVector();
432 if (DgDx_out != Teuchos::null) {
433 Thyra::DetachedMultiVectorView<Scalar> DgDx_out_view( *DgDx_out );
434 DgDx_out_view(0,0) = 1.0;
435 DgDx_out_view(0,1) = 0.0;
436 DgDx_out_view(1,0) = 0.0;
437 DgDx_out_view(1,1) = 1.0;
442template<
class Scalar>
443Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
447 if (!acceptModelParams_) {
448 return Teuchos::null;
450 TEUCHOS_ASSERT_IN_RANGE_UPPER_EXCLUSIVE( l, 0, Np_ );
453 else if (l == 1 || l == 2)
455 return Teuchos::null;
458template<
class Scalar>
459Teuchos::RCP<const Teuchos::Array<std::string> >
463 if (!acceptModelParams_) {
464 return Teuchos::null;
466 TEUCHOS_ASSERT_IN_RANGE_UPPER_EXCLUSIVE( l, 0, Np_ );
467 Teuchos::RCP<Teuchos::Array<std::string> > p_strings =
468 Teuchos::rcp(
new Teuchos::Array<std::string>());
470 p_strings->push_back(
"Model Coefficient: a");
471 p_strings->push_back(
"Model Coefficient: f");
472 p_strings->push_back(
"Model Coefficient: L");
475 p_strings->push_back(
"DxDp");
477 p_strings->push_back(
"Dx_dotDp");
481template<
class Scalar>
482Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
486 TEUCHOS_ASSERT_IN_RANGE_UPPER_EXCLUSIVE( j, 0, Ng_ );
492template<
class Scalar>
497 if (isInitialized_) {
502 typedef Thyra::ModelEvaluatorBase MEB;
505 MEB::InArgsSetup<Scalar> inArgs;
506 inArgs.setModelEvalDescription(this->description());
507 inArgs.setSupports( MEB::IN_ARG_t );
508 inArgs.setSupports( MEB::IN_ARG_x );
509 inArgs.setSupports( MEB::IN_ARG_beta );
510 inArgs.setSupports( MEB::IN_ARG_x_dot );
511 inArgs.setSupports( MEB::IN_ARG_alpha );
512 if (acceptModelParams_) {
520 MEB::OutArgsSetup<Scalar> outArgs;
521 outArgs.setModelEvalDescription(this->description());
522 outArgs.setSupports( MEB::OUT_ARG_f );
523 outArgs.setSupports( MEB::OUT_ARG_W_op );
524 if (acceptModelParams_) {
525 outArgs.set_Np_Ng(Np_,Ng_);
526 outArgs.setSupports( MEB::OUT_ARG_DfDp,0,
527 MEB::DERIV_MV_JACOBIAN_FORM );
528 outArgs.setSupports( MEB::OUT_ARG_DgDp,0,0,
529 MEB::DERIV_MV_JACOBIAN_FORM );
530 outArgs.setSupports( MEB::OUT_ARG_DgDx,0,
531 MEB::DERIV_MV_GRADIENT_FORM );
537 nominalValues_ = inArgs_;
540 nominalValues_.set_t(t0_ic_);
541 const RCP<Thyra::VectorBase<Scalar> > x_ic = createMember(x_space_);
543 Thyra::DetachedVectorView<Scalar> x_ic_view( *x_ic );
544 x_ic_view[0] = a_+b_*sin((f_/L_)*t0_ic_+phi_);
545 x_ic_view[1] = b_*(f_/L_)*cos((f_/L_)*t0_ic_+phi_);
547 nominalValues_.set_x(x_ic);
548 if (acceptModelParams_) {
549 const RCP<Thyra::VectorBase<Scalar> > p_ic = createMember(p_space_);
551 Thyra::DetachedVectorView<Scalar> p_ic_view( *p_ic );
556 nominalValues_.set_p(0,p_ic);
558 const RCP<Thyra::VectorBase<Scalar> > x_dot_ic = createMember(x_space_);
560 Thyra::DetachedVectorView<Scalar> x_dot_ic_view( *x_dot_ic );
561 x_dot_ic_view[0] = b_*(f_/L_)*cos((f_/L_)*t0_ic_+phi_);
562 x_dot_ic_view[1] = -b_*(f_/L_)*(f_/L_)*sin((f_/L_)*t0_ic_+phi_);
564 nominalValues_.set_x_dot(x_dot_ic);
567 isInitialized_ =
true;
571template<
class Scalar>
577 using Teuchos::ParameterList;
578 Teuchos::RCP<ParameterList> tmpPL = Teuchos::rcp(
new ParameterList(
"SinCosModel"));
579 if (paramList != Teuchos::null) tmpPL = paramList;
580 tmpPL->validateParametersAndSetDefaults(*this->getValidParameters());
581 this->setMyParamList(tmpPL);
582 Teuchos::RCP<ParameterList> pl = this->getMyNonconstParamList();
583 bool acceptModelParams = get<bool>(*pl,
"Accept model parameters");
584 bool haveIC = get<bool>(*pl,
"Provide nominal values");
585 bool useDfDpAsTangent = get<bool>(*pl,
"Use DfDp as Tangent");
586 if ( (acceptModelParams != acceptModelParams_) ||
589 isInitialized_ =
false;
591 acceptModelParams_ = acceptModelParams;
593 useDfDpAsTangent_ = useDfDpAsTangent;
594 a_ = get<Scalar>(*pl,
"Coeff a");
595 f_ = get<Scalar>(*pl,
"Coeff f");
596 L_ = get<Scalar>(*pl,
"Coeff L");
597 x0_ic_ = get<Scalar>(*pl,
"IC x0");
598 x1_ic_ = get<Scalar>(*pl,
"IC x1");
599 t0_ic_ = get<Scalar>(*pl,
"IC t0");
600 calculateCoeffFromIC_();
604template<
class Scalar>
605Teuchos::RCP<const Teuchos::ParameterList>
609 static Teuchos::RCP<const Teuchos::ParameterList> validPL;
610 if (is_null(validPL)) {
611 Teuchos::RCP<Teuchos::ParameterList> pl = Teuchos::parameterList();
612 pl->set(
"Accept model parameters",
false);
613 pl->set(
"Provide nominal values",
true);
614 pl->set(
"Use DfDp as Tangent",
false);
615 pl->set<std::string>(
"Output File Name",
"Tempus_BDF2_SinCos");
616 Teuchos::setDoubleParameter(
617 "Coeff a", 0.0,
"Coefficient a in model", &*pl);
618 Teuchos::setDoubleParameter(
619 "Coeff f", 1.0,
"Coefficient f in model", &*pl);
620 Teuchos::setDoubleParameter(
621 "Coeff L", 1.0,
"Coefficient L in model", &*pl);
622 Teuchos::setDoubleParameter(
623 "IC x0", 0.0,
"Initial Condition for x0", &*pl);
624 Teuchos::setDoubleParameter(
625 "IC x1", 1.0,
"Initial Condition for x1", &*pl);
626 Teuchos::setDoubleParameter(
627 "IC t0", 0.0,
"Initial time t0", &*pl);
628 Teuchos::setIntParameter(
629 "Number of Time Step Sizes", 1,
"Number time step sizes for convergence study", &*pl);
635template<
class Scalar>
640 phi_ = atan(((f_/L_)/x1_ic_)*(x0_ic_-a_))-(f_/L_)*t0_ic_;
641 b_ = x1_ic_/((f_/L_)*cos((f_/L_)*t0_ic_+phi_));
644template<
class Scalar>
645Teuchos::RCP<Thyra::LinearOpWithSolveBase<Scalar> >
650 RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> > W_factory = this->get_W_factory();
651 RCP<Thyra::LinearOpBase<Scalar> > matrix = this->create_W_op();
658 RCP<Thyra::MultiVectorBase<Scalar> > multivec = Teuchos::rcp_dynamic_cast<Thyra::MultiVectorBase<Scalar> >(matrix,
true);
660 RCP<Thyra::VectorBase<Scalar> > vec = Thyra::createMember(this->f_space_);
662 Thyra::DetachedVectorView<Scalar> vec_view( *vec );
666 V_V(multivec->col(0).ptr(),*vec);
668 Thyra::DetachedVectorView<Scalar> vec_view( *vec );
672 V_V(multivec->col(1).ptr(),*vec);
675 RCP<Thyra::LinearOpWithSolveBase<Scalar> > W =
676 Thyra::linearOpWithSolve<Scalar>(
683template<
class Scalar>
684Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
688 Teuchos::RCP<Thyra::MultiVectorBase<Scalar> > matrix = Thyra::createMembers(this->f_space_, this->dim_);
692template<
class Scalar>
693Thyra::ModelEvaluatorBase::InArgs<Scalar>
702 typedef Thyra::ModelEvaluatorBase MEB;
704 inArgs.setModelEvalDescription(this->description());
708template<
class Scalar>
709Thyra::ModelEvaluatorBase::OutArgs<Scalar>
713 typedef Thyra::ModelEvaluatorBase MEB;
714 MEB::OutArgsSetup<Scalar> outArgs;
715 outArgs.setModelEvalDescription(this->description());
716 outArgs.setSupports( MEB::OUT_ARG_f );
717 outArgs.setSupports( MEB::OUT_ARG_W_op );
718 outArgs.set_Np_Ng(this->Np_,0);
722template<
class Scalar>
726 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
727 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs
733 using Teuchos::rcp_dynamic_cast;
734 TEUCHOS_TEST_FOR_EXCEPTION( !this->isInitialized_, std::logic_error,
735 "Error, setupInOutArgs_ must be called first!\n");
737 const RCP<const VectorBase<Scalar> > x_in = inArgs.get_x().assert_not_null();
738 Thyra::ConstDetachedVectorView<Scalar> x_in_view( *x_in );
744 if (this->acceptModelParams_) {
745 const RCP<const VectorBase<Scalar> > p_in =
746 inArgs.get_p(0).assert_not_null();
747 Thyra::ConstDetachedVectorView<Scalar> p_in_view( *p_in );
753 Scalar beta = inArgs.get_beta();
755 const RCP<Thyra::LinearOpBase<Scalar> > W_out = outArgs.get_W_op();
756 if (inArgs.get_x_dot().is_null()) {
759 if (!is_null(W_out)) {
760 RCP<Thyra::MultiVectorBase<Scalar> > matrix =
761 Teuchos::rcp_dynamic_cast<Thyra::MultiVectorBase<Scalar> >(W_out,
true);
762 Thyra::DetachedMultiVectorView<Scalar> matrix_view( *matrix );
763 matrix_view(0,0) = 0.0;
764 matrix_view(1,0) = +beta;
765 matrix_view(0,1) = -beta*(f/L)*(f/L);
766 matrix_view(1,1) = 0.0;
772 RCP<const VectorBase<Scalar> > x_dot_in;
773 x_dot_in = inArgs.get_x_dot().assert_not_null();
774 Scalar alpha = inArgs.get_alpha();
775 if (!is_null(W_out)) {
776 RCP<Thyra::MultiVectorBase<Scalar> > matrix =
777 Teuchos::rcp_dynamic_cast<Thyra::MultiVectorBase<Scalar> >(W_out,
true);
778 Thyra::DetachedMultiVectorView<Scalar> matrix_view( *matrix );
779 matrix_view(0,0) = alpha;
780 matrix_view(1,0) = -beta;
781 matrix_view(0,1) = +beta*(f/L)*(f/L);
782 matrix_view(1,1) = alpha;
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs_bar, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs_bar) const
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
Teuchos::RCP< Thyra::LinearOpWithSolveBase< Scalar > > create_W() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int l) const
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > getExactSolution(double t) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters() const
Teuchos::RCP< Thyra::LinearOpWithSolveBase< Scalar > > create_W() const
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs_bar, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs_bar) const
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > getExactSensSolution(int j, double t) const
SinCosModel(Teuchos::RCP< Teuchos::ParameterList > pList=Teuchos::null)
void setParameterList(Teuchos::RCP< Teuchos::ParameterList > const ¶mList)
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_g_space(int j) const
void calculateCoeffFromIC_()
void setupInOutArgs_() const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int l) const