131 …PetscReal          f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt…  in FormFunctionGradient()  local
183       d4 = (xc - xb);  in FormFunctionGradient()
190       df2dxc = (d1 * hydhx + d4 * hxdhy);  in FormFunctionGradient()
194       df6dxc = d4 * hxdhy;  in FormFunctionGradient()
199       d4 *= rhy;  in FormFunctionGradient()
206       f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
210       f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6);  in FormFunctionGradient()
239     d4 = (user->right[j] - user->right[j + 1]) * rhy;  in FormFunctionGradient()
240     ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
245     d4 = (user->top[i + 1] - user->top[i]) * rhx;  in FormFunctionGradient()
246     ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
310 …PetscReal          f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt…  in QuadraticH()  local
368       d4 = (xc - xb) * rhy;  in QuadraticH()
375       f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in QuadraticH()
379       f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6);  in QuadraticH()
381 …x * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) / (f2 * f2 …  in QuadraticH()
384 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) …  in QuadraticH()
386       hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6);  in QuadraticH()
389 …d6) / (f6 * f6 * f6) + (hxdhy * (1.0 + d1 * d1) + hydhx * (1.0 + d4 * d4) - 2 * d1 * d4) / (f2 * f…  in QuadraticH()