177   PetscReal hx, hy, ehxhy;  in ComputeB()  local
184   hx    = two * pi / (nx + 1.0);  in ComputeB()
186   ehxhy = ecc * hx * hy;  in ComputeB()
197     temp = PetscSinScalar((i + 1) * hx);  in ComputeB()
214   PetscReal  hx, hy, hxhy, hxhx, hyhy;  in FormFunctionGradient()  local
225   hx   = two * pi / (nx + 1.0);  in FormFunctionGradient()
227   hxhy = hx * hy;  in FormFunctionGradient()
228   hxhx = one / (hx * hx);  in FormFunctionGradient()
247     xi     = (i + 1) * hx;  in FormFunctionGradient()
248     trule1 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi, ecc)) / six;      /* L(i,j) */  in FormFunctionGradient()
249     trule2 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi, ecc)) / six;      /* U(i,j) */  in FormFunctionGradient()
250     trule3 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi + hx, ecc)) / six; /* U(i+1,j) */  in FormFunctionGradient()
251     trule4 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi - hx, ecc)) / six; /* L(i-1,j) */  in FormFunctionGradient()
329   PetscReal hx, hy, hxhy, hxhx, hyhy;  in FormHessian()  local
338   hx   = two * pi / (nx + 1.0);  in FormHessian()
340   hxhy = hx * hy;  in FormHessian()
341   hxhx = one / (hx * hx);  in FormHessian()
353     xi     = (i + 1) * hx;  in FormHessian()
354     trule1 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi, ecc)) / six;      /* L(i,j) */  in FormHessian()
355     trule2 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi, ecc)) / six;      /* U(i,j) */  in FormHessian()
356     trule3 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi + hx, ecc)) / six; /* U(i+1,j) */  in FormHessian()
357     trule4 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi - hx, ecc)) / six; /* L(i-1,j) */  in FormHessian()