A Fast Algorithm for the Convolution of Functions with Compact Support Using Fourier Extensions

Journal Article · Thu Dec 21 00:00:00 UTC 2017 · SIAM Journal on Scientific Computing

DOI: https://doi.org/10.1137/17M1114764 · OSTI ID:1427516

Xu, Kuan ^[1]; Austin, Anthony P. ^[2]; Wei, Ke ^[3]

---

1. Univ. of Kent, Canterbury (United Kingdom). School of Mathematics, Statistics, and Actuarial Science
2. Argonne National Lab. (ANL), Argonne, IL (United States). Mathematics and Computer Science Division
3. Univ. of California, Davis, CA (United States). Dept. of Mathematics

In this paper, we present a new algorithm for computing the convolution of two compactly supported functions. The algorithm approximates the functions to be convolved using Fourier extensions and then uses the fast Fourier transform to efficiently compute Fourier extension approximations to the pieces of the result. Finally, the complexity of the algorithm is O(N(log N)²), where N is the number of degrees of freedom used in each of the Fourier extensions.

• Research Organization: Argonne National Lab. (ANL), Argonne, IL (United States); Univ. of California, Davis, CA (United States); Univ. of Kent, Canterbury (United Kingdom)
• Sponsoring Organization: National Science Foundation (NSF) (United States); Royal Society (United Kingdom); USDOE Office of Science (SC)
• Grant/Contract Number: AC02-06CH11357
• OSTI ID: 1427516
• Journal Information: SIAM Journal on Scientific Computing, Journal Name: SIAM Journal on Scientific Computing Journal Issue: 6 Vol. 39; ISSN 1064-8275
• Publisher: SIAMCopyright Statement
• Country of Publication: United States
• Language: English