Algorithme de convolution rapide et précis (comme FFT) pour une plage dynamique élevée?

Il semble que la convolution basée sur la FFT souffre d'une résolution limitée en virgule flottante en raison de tout évaluer autour des racines de l'unité, comme vous pouvez le voir dans l' erreur de facteur dans ce code Python:dix14101410^{14} from scipy.signal import convolve, fftconvolve a =...