ハイダイナミックレンジ用の高速かつ正確な畳み込みアルゴリズム(FFTなど)
FFTベースのたたみ込みでは、1の根の周りのすべてを評価するため、浮動小数点の解像度が制限されているようです。 1014101410^{14}このPythonコードの-factorエラー: from scipy.signal import convolve, fftconvolve a = [1.0, 1E-15] b = [1.0, 1E-15] convolve(a, b) # [ 1.00000000e+00, 2.00000000e-15, 1.00000000e-30] fftconvolve(a, b) # [ 1.00000000e+00, 2.11022302e-15, 1.10223025e-16] この問題の影響を受けない高速畳み込みアルゴリズムはありますか? または、直接(2次時間)畳み込みが正確な解を得る唯一の方法ですか? (そのような小さい数が切り落とさないように十分重要であるかどうかは私のポイントの外です。)