409   const PetscScalar *inptr;  in RHSFunction()  local
451   PetscCall(VecGetArrayRead(tmp_in, &inptr));  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 …tr[j * m + m - 1] = xc * inptr[j * m + m - 1] + 2.0 * xl * inptr[j * m + m - 1 - 1] + yl * inptr[j…  in RHSFunction()
463 …r[j * m + i] = xc * inptr[j * m + i] + xl * inptr[j * m + i - 1] + xr * inptr[j * m + i + 1] + yl …  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 …++) outptr[mn - m + i] = xc * inptr[mn - m + i] + xl * inptr[mn - m + i - 1] + xr * inptr[mn - m +…  in RHSFunction()
470   PetscCall(VecRestoreArrayRead(tmp_in, &inptr));  in RHSFunction()