Computing discrete convolutions with verified accuracy via Banach algebras and the FFT

[摘要] We introduce a method to compute rigorous component-wise enclosures of discrete convolutions using the fast Fourier transform, the properties of Banach algebras, and interval arithmetic. The purpose of this new approach is to improve the implementation and the applicability of computer-assisted proofs performed in weighed $\ell^1$ Banach algebras of Fourier/Chebyshev sequences, whose norms are known to be numerically unstable. We introduce some application examples, in particular a rigorous aposteriori error analysis for a steady state in the quintic Swift-Hohenberg PDE.

[学科分类] 应用数学

[关键词] discrete convolutions;Banach algebras;fast Fourier transform;interval arithmetic;rigorously verified numerics;quintic Swift-Hohenberg PDE