//********************************************************************** 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)); }