This may be fastest, but have you done any performance tests comparing this to, e.g., `np.interp`?  Also, would it perhaps be faster to use Tom's `convolve`?   (`convolve_fft` will definitely be slower because the kernel size is tiny, so it should not be used)