277   PetscReal   dx, dy, a, epsilon, xc, xl, xr, yl, yr;  in RHSJacobian()  local
288   xc = -2.0 * epsilon * (1.0 / (dx * dx) + 1.0 / (dy * dy));  in RHSJacobian()
295   v[0]   = xc;  in RHSJacobian()
305   v[1]   = xc;  in RHSJacobian()
315     v[1]   = xc;  in RHSJacobian()
327     v[0]   = xc;  in RHSJacobian()
338     v[0]   = xc;  in RHSJacobian()
350       v[0]   = xc;  in RHSJacobian()
365   v[0]   = xc;  in RHSJacobian()
374   v[0]   = xc;  in RHSJacobian()
385     v[1]   = xc;  in RHSJacobian()
406   PetscReal          xc, xl, xr, yl, yr;  in RHSFunction()  local
425   xc = -2.0 * epsilon * (1.0 / (dx * dx) + 1.0 / (dy * dy));  in RHSFunction()
456   outptr[0]     = xc * inptr[0] + xr * inptr[1] + yr * inptr[m];  in RHSFunction()
457   outptr[m - 1] = 2.0 * xl * inptr[m - 2] + xc * inptr[m - 1] + yr * inptr[m - 1 + m];  in RHSFunction()
458 …for (i = 1; i < m - 1; i++) outptr[i] = xc * inptr[i] + xl * inptr[i - 1] + xr * inptr[i + 1] + yr…  in RHSFunction()
461 …outptr[j * m]         = xc * inptr[j * m] + xr * inptr[j * m + 1] + yl * inptr[j * m - m] + yr * i…  in RHSFunction()
462 …outptr[j * m + m - 1] = xc * inptr[j * m + m - 1] + 2.0 * xl * inptr[j * m + m - 1 - 1] + yl * inp…  in RHSFunction()
463 …for (i = 1; i < m - 1; i++) outptr[j * m + i] = xc * inptr[j * m + i] + xl * inptr[j * m + i - 1] …  in RHSFunction()
466   outptr[mn - m] = xc * inptr[mn - m] + xr * inptr[mn - m + 1] + 2.0 * yl * inptr[mn - m - m];  in RHSFunction()
467   outptr[mn - 1] = 2.0 * xl * inptr[mn - 2] + xc * inptr[mn - 1] + 2.0 * yl * inptr[mn - 1 - m];  in RHSFunction()
468 …for (i = 1; i < m - 1; i++) outptr[mn - m + i] = xc * inptr[mn - m + i] + xl * inptr[mn - m + i - …  in RHSFunction()