26   PetscReal epsilon; /* the diffusion coefficient     */  member
63   data.epsilon = 0.1;  in main()
277   PetscReal   dx, dy, a, epsilon, xc, xl, xr, yl, yr;  in RHSJacobian()  local
286   epsilon = data->epsilon;  in RHSJacobian()
288   xc = -2.0 * epsilon * (1.0 / (dx * dx) + 1.0 / (dy * dy));  in RHSJacobian()
289   xl = 0.5 * a / dx + epsilon / (dx * dx);  in RHSJacobian()
290   xr = -0.5 * a / dx + epsilon / (dx * dx);  in RHSJacobian()
291   yl = 0.5 * a / dy + epsilon / (dy * dy);  in RHSJacobian()
292   yr = -0.5 * a / dy + epsilon / (dy * dy);  in RHSJacobian()
407   PetscReal          a, epsilon;  in RHSFunction()  local
423   epsilon = data->epsilon;  in RHSFunction()
425   xc = -2.0 * epsilon * (1.0 / (dx * dx) + 1.0 / (dy * dy));  in RHSFunction()
426   xl = 0.5 * a / dx + epsilon / (dx * dx);  in RHSFunction()
427   xr = -0.5 * a / dx + epsilon / (dx * dx);  in RHSFunction()
428   yl = 0.5 * a / dy + epsilon / (dy * dy);  in RHSFunction()
429   yr = -0.5 * a / dy + epsilon / (dy * dy);  in RHSFunction()