https://www.mdu.se/

mdu.sePublications
Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
A Parallel Nonuniform Fast Fourier Transform Library Based on an “Exponential of Semicircle" Kernel
Flatiron Institute, Simons Foundation, New York, USA.
Flatiron Institute, Simons Foundation, New York, USA.
Department of Mathematics, Simon Fraser University, Burnaby, BC, V5A 1S6, Canad.ORCID iD: 0000-0001-7425-8029
2019 (English)In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 41, no 5, p. C479-C504Article in journal (Refereed) Published
Abstract [en]

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 (nonuniform 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\beta \surd 1 - x 2 in 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. 

Place, publisher, year, edition, pages
2019. Vol. 41, no 5, p. C479-C504
National Category
Signal Processing
Identifiers
URN: urn:nbn:se:mdh:diva-65001DOI: 10.1137/18m120885xISI: 000493897100042Scopus ID: 2-s2.0-85074688530OAI: oai:DiVA.org:mdh-65001DiVA, id: diva2:1818748
Available from: 2023-12-12 Created: 2023-12-12 Last updated: 2023-12-12Bibliographically approved

Open Access in DiVA

No full text in DiVA

Other links

Publisher's full textScopus

Authority records

af Klinteberg, Ludvig

Search in DiVA

By author/editor
af Klinteberg, Ludvig
In the same journal
SIAM Journal on Scientific Computing
Signal Processing

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 16 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf