214   PetscReal  hx, hy, hxhy, hxhx, hyhy;  in FormFunctionGradient()  local
227   hxhy = hx * hy;  in FormFunctionGradient()
248     trule1 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi, ecc)) / six;      /* L(i,j) */  in FormFunctionGradient()
249     trule2 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi, ecc)) / six;      /* U(i,j) */  in FormFunctionGradient()
250     trule3 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi + hx, ecc)) / six; /* U(i+1,j) */  in FormFunctionGradient()
251     trule4 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi - hx, ecc)) / six; /* L(i-1,j) */  in FormFunctionGradient()
329   PetscReal hx, hy, hxhy, hxhx, hyhy;  in FormHessian()  local
340   hxhy = hx * hy;  in FormHessian()
354     trule1 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi, ecc)) / six;      /* L(i,j) */  in FormHessian()
355     trule2 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi, ecc)) / six;      /* U(i,j) */  in FormHessian()
356     trule3 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi + hx, ecc)) / six; /* U(i+1,j) */  in FormHessian()
357     trule4 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi - hx, ecc)) / six; /* L(i-1,j) */  in FormHessian()