Aliasing error of the exp$(eta sqrt{1-z^2})$ kernel in the nonuniform fast Fourier transform
Aliasing error of the exp$(eta sqrt{1-z^2})$ kernel in the nonuniform fast Fourier transform
复制标题
非均匀快速傅立叶变换中 exp$(eta sqrt{1-z^2})$ 核的混叠误差
DOI:
--
复制
发表时间:
2020
期刊:
影响因子:
--
通讯作者:
A. Barnett
中科院分区:
文献类型:
--
作者:
A. Barnett
The most popular algorithm for the nonuniform fast Fourier transform (NUFFT) uses the dilation of a kernel $phi$ to spread (or interpolate) between given nonuniform points and a uniform upsampled grid, combined with an FFT and diagonal scaling (deconvolution) in frequency space. The high performance of the recent FINUFFT library is in part due to its use of a new ``exponential of semicircle' kernel $phi(z)=e^{eta sqrt{1-z^2}}$, for $zin[-1,1]$, zero otherwise, whose Fourier transform $hatphi$ is unknown analytically. We place this kernel on a rigorous footing by proving an aliasing error estimate which bounds the error of the one-dimensional NUFFT of types 1 and 2 in exact arithmetic. Asymptotically in the kernel width measured in upsampled grid points, the error is shown to decrease with an exponential rate arbitrarily close to that of the popular Kaiser--Bessel kernel. This requires controlling a conditionally-convergent sum over the tails of $hatphi$, using steepest descent, other classical estimates on contour integrals, and a phased sinc sum. We also draw new connections between the above kernel, Kaiser--Bessel, and prolate spheroidal wavefunctions of order zero, which all appear to share an optimal exponential convergence rate.
影响因子:
1.7
作者:
F. Nestler
通讯作者:
F. Nestler
DOI:
10.1016/j.jcp.2014.12.052
发表时间:
2015-03
期刊:
J. Comput. Phys.
影响因子:
--
作者:
F. Nestler;Michael Pippig;D. Potts
通讯作者:
F. Nestler;Michael Pippig;D. Potts