170   PetscReal   f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb;  in FormFunction()  local
218       d3 = (xc - xt);  in FormFunction()
223       d3 *= rhy;  in FormFunction()
227       f4 = PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunction()
236       d3 = (user->left[j - ys + 1] - user->left[j - ys + 2]) * rhy;  in FormFunction()
238       ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunction()
244       d3 = (user->bottom[i - xs + 1] - x[0][i]) * rhy;  in FormFunction()
245       ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunction()
304   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
362       d3 = (xc - xt);  in FormFunctionGradient()
371       df3dxc = d3 * hxdhy;  in FormFunctionGradient()
372       df4dxc = (d2 * hydhx + d3 * hxdhy);  in FormFunctionGradient()
378       d3 *= rhy;  in FormFunctionGradient()
387       f3 = PetscSqrtReal(1.0 + d3 * d3 + d8 * d8);  in FormFunctionGradient()
388       f4 = PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
408       d3 = (user->left[j - ys + 1] - user->left[j - ys + 2]) * rhy;  in FormFunctionGradient()
410       ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
416       d3 = (user->bottom[i - xs + 1] - x[0][i]) * rhy;  in FormFunctionGradient()
417       ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
508   PetscReal   f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb;  in QuadraticH()  local
576       d3 = (xc - xt) / hy;  in QuadraticH()
585       f3 = PetscSqrtReal(1.0 + d3 * d3 + d8 * d8);  in QuadraticH()
586       f4 = PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in QuadraticH()
591 …x * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) / (f4 * f4 …  in QuadraticH()
592 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) …  in QuadraticH()
596       htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3);  in QuadraticH()
598 …d4) / (f2 * f2 * f2) + (hxdhy * (1.0 + d2 * d2) + hydhx * (1.0 + d3 * d3) - 2 * d2 * d3) / (f4 * f…  in QuadraticH()