45 #ifndef DUMMY_TEST_MODEL_EVALUATOR_DEF_HPP
46 #define DUMMY_TEST_MODEL_EVALUATOR_DEF_HPP
49 #include "Thyra_DummyTestModelEvaluator_decl.hpp"
50 #include "Thyra_SimpleDenseLinearOp.hpp"
51 #include "Thyra_DefaultSpmdVectorSpace.hpp"
52 #include "Thyra_DefaultSerialDenseLinearOpWithSolveFactory.hpp"
53 #include "Thyra_DefaultPreconditioner.hpp"
54 #include "Thyra_DetachedMultiVectorView.hpp"
55 #include "Thyra_DetachedVectorView.hpp"
56 #include "Thyra_MultiVectorStdOps.hpp"
57 #include "Thyra_VectorStdOps.hpp"
65 template<
class Scalar>
67 dummyTestModelEvaluator(
71 const bool supports_x_dot,
72 const bool supports_x_dot_dot,
73 const bool supports_extended_inargs,
74 const bool supports_extended_outargs,
75 const bool supports_derivatives
85 template<
class Scalar>
90 const bool supports_x_dot,
91 const bool supports_x_dot_dot,
92 const bool supports_extended_inargs,
93 const bool supports_extended_outargs,
94 const bool supports_derivatives
101 x_space_ = defaultSpmdVectorSpace<Scalar>(x_size);
103 p_space_.resize(p_sizes.
size());
105 p_space_[l] = defaultSpmdVectorSpace<Scalar>(p_sizes[l]);
110 g_space_.resize(g_sizes.
size());
112 g_space_[j] = defaultSpmdVectorSpace<Scalar>(g_sizes[j]);
115 W_factory_ = defaultSerialDenseLinearOpWithSolveFactory<Scalar>();
117 MEB::InArgsSetup<Scalar> inArgs;
118 inArgs.setModelEvalDescription(this->description());
119 inArgs.set_Np_Ng(p_space_.size(),g_sizes.
size());
120 inArgs.setSupports(MEB::IN_ARG_x);
122 inArgs.setSupports(MEB::IN_ARG_x_dot);
123 if (supports_x_dot_dot)
124 inArgs.setSupports(MEB::IN_ARG_x_dot_dot);
125 inArgs.setSupports(MEB::IN_ARG_step_size);
126 inArgs.setSupports(MEB::IN_ARG_stage_number);
127 inArgs.template setSupports<Thyra::MockExtendedInArgs<Scalar> >(
true);
129 if (!supports_extended_inargs)
130 inArgs.template setSupports<Thyra::MockExtendedInArgs<Scalar> >(
false);
131 prototypeInArgs_ = inArgs;
133 MEB::OutArgsSetup<Scalar> outArgs;
134 outArgs.setModelEvalDescription(this->description());
135 outArgs.set_Np_Ng(p_space_.size(), g_space_.size());
136 outArgs.setSupports(MEB::OUT_ARG_f);
137 outArgs.setSupports(MEB::OUT_ARG_W_op);
138 outArgs.setSupports(MEB::OUT_ARG_W_prec);
139 outArgs.template setSupports<Thyra::MockExtendedOutArgs<Scalar> >(
true);
141 if (!supports_extended_outargs)
142 outArgs.template setSupports<Thyra::MockExtendedOutArgs<Scalar> >(
false);
143 if (supports_derivatives)
145 #ifdef Thyra_BUILD_HESSIAN_SUPPORT
146 outArgs.setHessianSupports(
true);
149 MEB::DerivativeSupport derivativeSupport(MEB::DERIV_MV_BY_COL);
151 for (
int j=0; j<g_space_.size(); ++j)
152 outArgs.setSupports(MEB::OUT_ARG_DgDx,j,derivativeSupport);
154 prototypeOutArgs_ = outArgs;
156 nominalValues_ = inArgs;
158 V_S(x0.
ptr(), ST::zero());
159 nominalValues_.set_x(x0);
167 template<
class Scalar>
175 template<
class Scalar>
183 template<
class Scalar>
187 return Teuchos::null;
191 template<
class Scalar>
199 template<
class Scalar>
207 template<
class Scalar>
215 template<
class Scalar>
219 return nominalValues_;
223 template<
class Scalar>
231 template<
class Scalar>
239 template<
class Scalar>
243 return createNonconstSimpleDenseLinearOp<Scalar>(
244 createMembers<Scalar>(f_space_, x_space_->dim())
249 template<
class Scalar>
253 return nonconstUnspecifiedPrec<Scalar>(
254 createNonconstSimpleDenseLinearOp<Scalar>(
255 createMembers<Scalar>(f_space_, x_space_->dim())
261 template<
class Scalar>
269 template<
class Scalar>
273 return prototypeInArgs_;
277 template<
class Scalar>
286 template<
class Scalar>
293 MEB::InArgsSetup<Scalar> inArgs(prototypeInArgs_);
294 inArgs.set_Np(new_size);
295 prototypeInArgs_ = inArgs;
298 MEB::OutArgsSetup<Scalar> outArgs(prototypeOutArgs_);
299 outArgs.set_Np_Ng(new_size,g_space_.size());
300 prototypeOutArgs_ = outArgs;
305 template<
class Scalar>
307 DummyTestModelEvaluator<Scalar>::
308 change_p_size_correctly(
const Ordinal new_size)
310 this->change_p_size_incorrectly(new_size);
311 this->resetDefaultBase();
317 template<
class Scalar>
318 ModelEvaluatorBase::OutArgs<Scalar>
319 DummyTestModelEvaluator<Scalar>::createOutArgsImpl()
const
321 return prototypeOutArgs_;
325 template<
class Scalar>
326 void DummyTestModelEvaluator<Scalar>::evalModelImpl(
327 const ModelEvaluatorBase::InArgs<Scalar> &,
328 const ModelEvaluatorBase::OutArgs<Scalar> &
344 #define DUMMY_TEST_MODEL_EVALUATOR_INSTANT(SCALAR) \
346 template class DummyTestModelEvaluator<SCALAR >; \
348 template Teuchos::RCP<DummyTestModelEvaluator<SCALAR > > \
349 dummyTestModelEvaluator( \
350 const Ordinal x_size, \
351 const ArrayView<const Ordinal> &p_sizes, \
352 const ArrayView<const Ordinal> &g_sizes, \
353 const bool supports_x_dot, \
354 const bool supports_x_dot_dot, \
355 const bool supports_extended_inargs, \
356 const bool supports_extended_outargs, \
357 const bool supports_derivatives \
361 #endif // DUMMY_TEST_MODEL_EVALUATOR_DEF_HPP
ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
void reportFinalPoint(const ModelEvaluatorBase::InArgs< Scalar > &finalPoint, const bool wasSolved)
RCP< const VectorSpaceBase< Scalar > > get_p_space(int l) const
RCP< const Teuchos::Array< std::string > > get_p_names(int l) const
ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Test helper ModelEvaluator.
RCP< LinearOpBase< Scalar > > create_W_op() const
ModelEvaluatorBase::InArgs< Scalar > getUpperBounds() const
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Teuchos::Ordinal Ordinal
Type for the dimension of a vector space. `*.
RCP< const VectorSpaceBase< Scalar > > get_x_space() const
Base subclass for ModelEvaluator that defines some basic types.
Teuchos::ArrayView< const std::string > get_g_names(int j) const
RCP< const VectorSpaceBase< Scalar > > get_f_space() const
ModelEvaluatorBase::InArgs< Scalar > getLowerBounds() const
RCP< const VectorSpaceBase< Scalar > > get_g_space(int j) const
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)
DummyTestModelEvaluator(const Ordinal x_size, const ArrayView< const Ordinal > &p_sizes, const ArrayView< const Ordinal > &g_sizes, const bool supports_x_dot=false, const bool supports_x_dot_dot=false, const bool supports_extended_inargs=true, const bool supports_extended_outargs=true, const bool supports_derivatives=false)
RCP< const LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
RCP< PreconditionerBase< Scalar > > create_W_prec() const
Concrete aggregate class for all input arguments computable by a ModelEvaluator subclass object...