[scottdupoy]

That's exactly how it worked, hand rolled FFT and filtering following the method in "Numerical Recipes for C"

[gsliepen]

Oh, Numerical Recipes is nice but their algorithm implementations are not really state-of-the-art. I highly recommend using FFTW (https://fftw.org/) as it will likely give you a substantial performance improvement.