A very important optimization here is that
only needs to be recomputed within the finite-difference domain
of dependence of the
th grid point.