A Fast Analysis-Based Discrete Hankel Transform Using Asymptotic Expansions
The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
Citation Townsend, Alex. “A Fast Analysis-Based Discrete Hankel Transform Using Asymptotic Expansions.” SIAM J. Numer. Anal. 53, no. 4 (January 2015): 1897–1917. © 2015, Society for Industrial and Applied Mathematics
As Published http://dx.doi.org/10.1137/151003106
Publisher Society for Industrial and Applied Mathematics
Version Final published version
Citable link http://hdl.handle.net/1721.1/100550
Terms of Use Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use. SIAM J. NUMER. ANAL. c 2015 Society for Industrial and Applied Mathematics Vol. 53, No. 4, pp. 1897–1917
A FAST ANALYSIS-BASED DISCRETE HANKEL TRANSFORM USING ASYMPTOTIC EXPANSIONS ∗
ALEX TOWNSEND†
Abstract. A fast and numerically stable algorithm is described for computing the discrete Hankel transform of order 0 as well as evaluating Schl¨omilch and Fourier–Bessel expansions in O(N(log N)2/ loglog N) operations. The algorithm is based on an asymptotic expansion for Bessel functions of large arguments, the fast Fourier transform, and the Neumann addition formula. All the algorithmic parameters are selected from error bounds to achieve a near-optimal computational cost for any accuracy goal. Numerical results demonstrate theefficiency of the resulting algorithm.
Key words. Hankel transform, Bessel functions, asymptotic expansions, fast Fourier transform
AMS subject classifications. 65R10, 33C10
DOI. 10.1137/151003106
ν ≥ f , → 1. Introduction. The Hankel transform of order 0 of a function :[0 1] C is defined as [21, Chap. 9] 1 (1.1) Fν (ω)= f(r)Jν (rω)rdr, ω>0, 0 J ν ν where ν is the Bessel function of the first kind with parameter .When =0,(1.1) is essentially the two-dimensional Fourier transform of a radially symmetric function with support on the unit disc and appears in the Fourier–Hankel–Abel cycle [21, sect. 9.3]. Furthermore, the Hankel transform of integer order ν ≥ 1 is equivalent to the Fourier transform of functions of the form eiνθf(r), where θ ∈ [0, 2π). Throughout ν this paper we take to be an integer. Various semidiscrete and discrete variants of (1.1 ) are employed in optics [7], electromagnetics [17], medical imaging [13], and the numerical solution of partial differential equations [3]. At some stage these require the computation of the following sums: