Algoritma konvolusi cepat & akurat (seperti FFT) untuk rentang dinamis tinggi?

Tampaknya konvolusi berbasis FFT menderita resolusi floating-point terbatas karena mengevaluasi segala sesuatu di sekitar akar persatuan, seperti yang dapat Anda lihat di 1014101410^{14}-faktor kesalahan dalam kode Python ini: from scipy.signal import convolve, fftconvolve a = [1.0, 1E-15] b =...