30   PetscReal ecc;    /* test problem parameter */  member
41 static PetscReal      p(PetscReal xi, PetscReal ecc);
68   user.ecc    = 0.1;  in main()
75   PetscCall(PetscOptionsGetReal(NULL, NULL, "-ecc", &user.ecc, &flg));  in main()
80 …x: %" PetscInt_FMT ",  my: %" PetscInt_FMT ",  ecc: %g \n\n", user.nx, user.ny, (double)user.ecc));  in main()
166 static PetscReal p(PetscReal xi, PetscReal ecc)  in p()  argument
168   PetscReal t = 1.0 + ecc * PetscCosScalar(xi);  in p()
179   PetscReal ecc = user->ecc;  in ComputeB()  local
186   ehxhy = ecc * hx * hy;  in ComputeB()
216   PetscReal  ecc = user->ecc, trule1, trule2, trule3, trule4, trule5, trule6;  in FormFunctionGradient()  local
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()
331   PetscReal ecc = user->ecc, trule1, trule2, trule3, trule4, trule5, trule6;  in FormHessian()  local
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()