This method is implemented the same as goldstein method w.r.t. the overlap, blocksize etc. Algorithm per block (SIZE,SIZE) is to perform a 2D FFT of the block, and then to multiply pointwise with the kernel, which is padded with zeros. The kernel is centered around zero frequency.

Leijen 2009-04-14