FFT-BASED CONVOLUTION
The block convolutions in the overlap and add method can be more efficiently computed using a Fast Fourier Transform.

The gain in this method is approximately:

fast/slow = log2(N) / M
For N=1024, and M=100, we can achieve a factor of 10 improvement.