double TestEvalMediator::evalDummyBasis(int m, const MultiIndex& mi) const { TEUCHOS_TEST_FOR_EXCEPTION(mi.order() > 1, std::runtime_error, "TestEvalMediator::evalDummyBasis found multiindex " "order > 1. The bad multiindex was " << mi.toString()); ADReal result = fields_[m].basis().evaluate(ADField::evalPoint()); SUNDANCE_MSG3(verb(), "basis.value() " << result.value()); SUNDANCE_MSG3(verb(), "basis.gradient() " << result.gradient()); if (mi.order()==0) { return result.value(); } else { return result.gradient()[mi.firstOrderDirection()]; } }