119   PetscScalar f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb;  in FormObjectiveLocal()  local
148       d1 = (xc - xl);  in FormObjectiveLocal()
153       d1 /= hx;  in FormObjectiveLocal()
158       f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormObjectiveLocal()
182       d1 = (x[j][mx - 1] - user->right[j + 1]) / hx;  in FormObjectiveLocal()
184       ft += PetscSqrtReal(1.0 + d1 * d1 + d4 * d4);  in FormObjectiveLocal()
189       d1 = (x[my - 1][i] - user->top[i + 1]) / hy;  in FormObjectiveLocal()
191       ft += PetscSqrtReal(1.0 + d1 * d1 + d4 * d4);  in FormObjectiveLocal()
195     d1 = (user->left[0] - user->left[1]) / hy;  in FormObjectiveLocal()
197     ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2);  in FormObjectiveLocal()
200     d1 = (user->right[ym + 1] - user->right[ym]) / hy;  in FormObjectiveLocal()
202     ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2);  in FormObjectiveLocal()
217   PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb;  in FormFunctionLocal()  local
255       d1 = (xc - xl);  in FormFunctionLocal()
264       df1dxc = d1 * hydhx;  in FormFunctionLocal()
265       df2dxc = (d1 * hydhx + d4 * hxdhy);  in FormFunctionLocal()
271       d1 /= hx;  in FormFunctionLocal()
280       f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7);  in FormFunctionLocal()
281       f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionLocal()
309   PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb;  in FormJacobianLocal()  local
355       d1 = (xc - xl) / hx;  in FormJacobianLocal()
364       f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7);  in FormJacobianLocal()
365       f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormJacobianLocal()
371 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) …  in FormJacobianLocal()
374 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) …  in FormJacobianLocal()
377       htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3);  in FormJacobianLocal()
379 … * (1.0 + d6 * d6) / (f6 * f6 * f6) + (hxdhy * (1.0 + d1 * d1) + hydhx * (1.0 + d4 * d4) - 2.0 * d…  in FormJacobianLocal()