170   PetscReal   f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb;  in FormFunction()  local
217       d2 = (xc - xr);  in FormFunction()
222       d2 *= rhx;  in FormFunction()
227       f4 = PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunction()
237       d2 = (user->left[j - ys + 1] - x[j][0]) * rhx;  in FormFunction()
238       ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunction()
243       d2 = (user->bottom[i + 1 - xs] - user->bottom[i - xs + 2]) * rhx;  in FormFunction()
245       ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunction()
264     d2 = (user->bottom[0] - user->bottom[1]) * rhx;  in FormFunction()
265     ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2);  in FormFunction()
269     d2 = (user->top[xm + 1] - user->top[xm]) * rhx;  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
361       d2 = (xc - xr);  in FormFunctionGradient()
372       df4dxc = (d2 * hydhx + d3 * hxdhy);  in FormFunctionGradient()
373       df5dxc = d2 * hydhx;  in FormFunctionGradient()
377       d2 *= rhx;  in FormFunctionGradient()
388       f4 = PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
389       f5 = PetscSqrtReal(1.0 + d2 * d2 + d5 * d5);  in FormFunctionGradient()
409       d2 = (user->left[j - ys + 1] - x[j][0]) * rhx;  in FormFunctionGradient()
410       ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
415       d2 = (user->bottom[i + 1 - xs] - user->bottom[i - xs + 2]) * rhx;  in FormFunctionGradient()
417       ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
438     d2 = (user->bottom[0] - user->bottom[1]) * rhx;  in FormFunctionGradient()
439     ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2);  in FormFunctionGradient()
443     d2 = (user->top[xm + 1] - user->top[xm]) * rhx;  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
575       d2 = (xc - xr) / hx;  in QuadraticH()
586       f4 = PetscSqrtReal(1.0 + d3 * d3 + d2 * d2);  in QuadraticH()
587       f5 = PetscSqrtReal(1.0 + d2 * d2 + d5 * d5);  in QuadraticH()
591 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) …  in QuadraticH()
592 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) …  in QuadraticH()
595       hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6);  in QuadraticH()
598 … d4) - 2 * d1 * d4) / (f2 * f2 * f2) + (hxdhy * (1.0 + d2 * d2) + hydhx * (1.0 + d3 * d3) - 2 * d2…  in QuadraticH()