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
271       d1 = (xc - xl);  in FormFunctionGradient()
280       df1dxc = d1 * hydhx;  in FormFunctionGradient()
281       df2dxc = (d1 * hydhx + d4 * hxdhy);  in FormFunctionGradient()
287       d1 *= rhx;  in FormFunctionGradient()
296       f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7);  in FormFunctionGradient()
297       f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
334       d1 = (x[(j + 1 - gys) * gxm - 1] - right[j - ys + 1]) * rhx;  in FormFunctionGradient()
336       ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
341       d1 = (x[(gym - 1) * gxm + i - gxs] - top[i - xs + 1]) * rhy;  in FormFunctionGradient()
343       ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
348     d1 = (left[0] - left[1]) * rhy;  in FormFunctionGradient()
350     ft += PetscSqrtScalar(1.0 + d1 * d1 + d2 * d2);  in FormFunctionGradient()
353     d1 = (right[ym + 1] - right[ym]) * rhy;  in FormFunctionGradient()
355     ft += PetscSqrtScalar(1.0 + d1 * d1 + 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
494       d1 = (xc - xl) * rhx;  in FormHessian()
503       f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7);  in FormHessian()
504       f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormHessian()
510 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) …  in FormHessian()
513 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) …  in FormHessian()
516       htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3);  in FormHessian()
518 …y * (1.0 + d6 * d6) / (f6 * f6 * f6) + (hxdhy * (1.0 + d1 * d1) + hydhx * (1.0 + d4 * d4) - 2 * d1…  in FormHessian()