English

A parallel non-uniform fast Fourier transform library based on an "exponential of semicircle" kernel

Numerical Analysis 2019-04-10 v2 Mathematical Software Numerical Analysis

Abstract

The nonuniform fast Fourier transform (NUFFT) generalizes the FFT to off-grid data. Its many applications include image reconstruction, data analysis, and the numerical solution of differential equations. We present FINUFFT, an efficient parallel library for type 1 (nonuiform to uniform), type 2 (uniform to nonuniform), or type 3 (nonuniform to nonuniform) transforms, in dimensions 1, 2, or 3. It uses minimal RAM, requires no precomputation or plan steps, and has a simple interface to several languages. We perform the expensive spreading/interpolation between nonuniform points and the fine grid via a simple new kernel---the `exponential of semicircle' eβ1x2e^{\beta \sqrt{1-x^2}} in x[1,1]x\in[-1,1]---in a cache-aware load-balanced multithreaded implementation. The deconvolution step requires the Fourier transform of the kernel, for which we propose efficient numerical quadrature. For types 1 and 2, rigorous error bounds asymptotic in the kernel width approach the fastest known exponential rate, namely that of the Kaiser--Bessel kernel. We benchmark against several popular CPU-based libraries, showing favorable speed and memory footprint, especially in three dimensions when high accuracy and/or clustered point distributions are desired.

Keywords

Cite

@article{arxiv.1808.06736,
  title  = {A parallel non-uniform fast Fourier transform library based on an "exponential of semicircle" kernel},
  author = {Alex H. Barnett and Jeremy F. Magland and Ludvig af Klinteberg},
  journal= {arXiv preprint arXiv:1808.06736},
  year   = {2019}
}

Comments

25 pages, 9 figures

R2 v1 2026-06-23T03:39:04.826Z