Inside a tf.function, I am implementing a 2D IIR filter. Assuming my IIR filter is of shape (dim1, dim2), I need the (dim1, dim2) most recent outputs to compute the next iteration.