51 const Teuchos::RCP<EpetraExt::ModelEvaluator>& me_,
52 const Teuchos::Array<int>& mp_p_index_map_,
53 const Teuchos::Array<int>& mp_g_index_map_,
54 const Teuchos::Array< Teuchos::RCP<const Epetra_Map> >& base_g_maps_) :
56 mp_p_index_map(mp_p_index_map_),
57 mp_g_index_map(mp_g_index_map_),
58 base_g_maps(base_g_maps_),
61 num_p_mp(mp_p_index_map.size()),
62 num_g_mp(mp_g_index_map.size())
64 InArgs me_inargs =
me->createInArgs();
65 OutArgs me_outargs =
me->createOutArgs();
67 num_g = me_outargs.Ng();
69 TEUCHOS_TEST_FOR_EXCEPTION(
72 <<
"Error! Stokhos::MPInverseModelEvaluator::MPInverseModelEvaluator():"
73 <<
" Base response map array size does not match size of index array!");
153 OutArgsSetup outArgs;
154 OutArgs me_outargs = me->createOutArgs();
156 outArgs.setModelEvalDescription(this->description());
158 outArgs.set_Np_Ng(num_p, num_g);
159 for (
int i=0; i<num_g_mp; i++) {
160 outArgs.setSupports(OUT_ARG_g_mp, mp_g_index_map[i],
true);
161 for (
int j=0;
j<num_p;
j++)
162 outArgs.setSupports(OUT_ARG_DgDp_mp, mp_g_index_map[i],
j,
163 me_outargs.supports(OUT_ARG_DgDp,i,
j));
171 const OutArgs& outArgs)
const
174 InArgs me_inargs = me->createInArgs();
177 for (
int i=0; i<num_p; i++)
178 me_inargs.set_p(i, inArgs.get_p(i));
181 for (
int i=0; i<num_p_mp; i++) {
182 mp_const_vector_t p_mp = inArgs.get_p_mp(mp_p_index_map[i]);
183 if (p_mp != Teuchos::null) {
184 me_inargs.set_p(i+num_p, p_mp->getBlockVector());
189 OutArgs me_outargs = me->createOutArgs();
193 for (
int i=0; i<num_g_mp; i++) {
194 int ii = mp_g_index_map[i];
197 mp_vector_t g_mp = outArgs.get_g_mp(ii);
198 if (g_mp != Teuchos::null) {
199 me_outargs.set_g(i, Teuchos::rcp_dynamic_cast<Epetra_Vector>(g_mp->getBlockVector()));
203 for (
int j=0;
j<num_p;
j++) {
204 if (!outArgs.supports(OUT_ARG_DgDp_mp, ii,
j).none()) {
205 MPDerivative dgdp_mp = outArgs.get_DgDp_mp(ii,
j);
206 Teuchos::RCP<Stokhos::ProductEpetraMultiVector> dgdp_mp_mv =
207 dgdp_mp.getMultiVector();
208 Teuchos::RCP<Epetra_Operator> dgdp_mp_op =
209 dgdp_mp.getLinearOp();
210 if (dgdp_mp_mv != Teuchos::null) {
212 i,
j, Derivative(dgdp_mp_mv->getBlockMultiVector(),
213 dgdp_mp.getMultiVectorOrientation()));
215 else if (dgdp_mp_op != Teuchos::null) {
216 me_outargs.set_DgDp(i,
j, Derivative(dgdp_mp_op));
224 me->evalModel(me_inargs, me_outargs);