2d convolution + cumsum