207   PetscReal  f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb;  in FormFunctionGradient()  local
273       d3 = (xc - xt);  in FormFunctionGradient()
282       df3dxc = d3 * hxdhy;  in FormFunctionGradient()
283       df4dxc = (d2 * hydhx + d3 * hxdhy);  in FormFunctionGradient()
289       d3 *= rhy;  in FormFunctionGradient()
298       f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8);  in FormFunctionGradient()
299       f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
319       d3 = (left[j - ys + 1] - left[j - ys + 2]) * rhy;  in FormFunctionGradient()
321       ft = ft + PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
327       d3 = (bottom[i - xs + 1] - x[i - gxs]) * rhy;  in FormFunctionGradient()
328       ft = ft + PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
423   PetscReal  f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb;  in FormHessian()  local
496       d3 = (xc - xt) * rhy;  in FormHessian()
505       f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8);  in FormHessian()
506       f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2);  in FormHessian()
511 …x * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) / (f4 * f4 …  in FormHessian()
512 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) …  in FormHessian()
516       htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3);  in FormHessian()
518 …d4) / (f2 * f2 * f2) + (hxdhy * (1.0 + d2 * d2) + hydhx * (1.0 + d3 * d3) - 2 * d2 * d3) / (f4 * f…  in FormHessian()