170   PetscReal   f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb;  in FormFunction()  local
216       d1 = (xc - xl);  in FormFunction()
221       d1 *= rhx;  in FormFunction()
226       f2 = PetscSqrtReal(1.0 + d1 * d1 + d4 * d4);  in FormFunction()
250       d1 = (x[j][mx - 1] - user->right[j - ys + 1]) * rhx;  in FormFunction()
252       ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4);  in FormFunction()
257       d1 = (x[my - 1][i] - user->top[i - xs + 1]) * rhy;  in FormFunction()
259       ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4);  in FormFunction()
263     d1 = (user->left[0] - user->left[1]) / hy;  in FormFunction()
265     ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2);  in FormFunction()
268     d1 = (user->right[ym + 1] - user->right[ym]) * rhy;  in FormFunction()
270     ft += PetscSqrtReal(1.0 + d1 * d1 + 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
360       d1 = (xc - xl);  in FormFunctionGradient()
369       df1dxc = d1 * hydhx;  in FormFunctionGradient()
370       df2dxc = (d1 * hydhx + d4 * hxdhy);  in FormFunctionGradient()
376       d1 *= rhx;  in FormFunctionGradient()
385       f1 = PetscSqrtReal(1.0 + d1 * d1 + d7 * d7);  in FormFunctionGradient()
386       f2 = PetscSqrtReal(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
423       d1 = (x[j][mx - 1] - user->right[j - ys + 1]) * rhx;  in FormFunctionGradient()
425       ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
430       d1 = (x[my - 1][i] - user->top[i - xs + 1]) * rhy;  in FormFunctionGradient()
432       ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
437     d1 = (user->left[0] - user->left[1]) / hy;  in FormFunctionGradient()
439     ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2);  in FormFunctionGradient()
442     d1 = (user->right[ym + 1] - user->right[ym]) * rhy;  in FormFunctionGradient()
444     ft += PetscSqrtReal(1.0 + d1 * d1 + 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
574       d1 = (xc - xl) / hx;  in QuadraticH()
583       f1 = PetscSqrtReal(1.0 + d1 * d1 + d7 * d7);  in QuadraticH()
584       f2 = PetscSqrtReal(1.0 + d1 * d1 + d4 * d4);  in QuadraticH()
590 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) …  in QuadraticH()
593 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) …  in QuadraticH()
596       htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3);  in QuadraticH()
598 …y * (1.0 + d6 * d6) / (f6 * f6 * f6) + (hxdhy * (1.0 + d1 * d1) + hydhx * (1.0 + d4 * d4) - 2 * d1…  in QuadraticH()