41 static PetscReal      p(PetscReal xi, PetscReal ecc);
166 static PetscReal p(PetscReal xi, PetscReal ecc)  in p()  argument
168   PetscReal t = 1.0 + ecc * PetscCosScalar(xi);  in p()
215   PetscReal  xi, v[5];  in FormFunctionGradient()  local
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()
330   PetscReal xi, v[5];  in FormHessian()  local
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()