//**********************************************************************
PHX_EVALUATE_FIELDS(HessianTestEvaluator,workset)
{
// Grad = y * std::cos(x*y)
// = x * std::cos(x*y)-0.25*std::sin(y)
// Hess*v = dy * std::cos(x*y) - y * ( dx * sin(x*y) + dy * sin(x*y))
// = dx * std::cos(x*y) - x * ( dx * sin(x*y) + dy * sin(x*y)) - 0.25 * dy * std::cos(y)
//
for(int ip=0;ip<result.dimension_0();++ip)
result(ip) = std::sin(x(ip)*y(ip))+0.25*std::cos(y(ip));
}
