In my case it would be 2D convolution. So for example if the filter is 5x5 it would be applied to a 5x5 patch of the input and the two matrices would be multiplied using ordinary matrix multiplication to produce a 5x5 output matrix. Is there a nice way to implement this?