Fourier-Based Fast Multipole Method for the Helmholtz Equation

The fast multipole method (FMM) has had great success in reducing the computational complexity of solving the boundary integral form of the Helmholtz equation. We present a formulation of the Helmholtz FMM that uses Fourier basis functions rather than spherical harmonics. By modifying the transfer function in the precomputation stage of the FMM, time-critical stages of the algorithm are accelerated by causing the interpolation operators to become straightforward applications of fast Fourier transforms, retaining the diagonality of the transfer function, and providing a simplified error analysis. Using Fourier analysis, constructive algorithms are derived to a priori determine an integration quadrature for a given error tolerance. Sharp error bounds are derived and verified numerically. Various optimizations are considered to reduce the number of quadrature points and reduce the cost of computing the transfer function.

[1]  Francis Collino,et al.  Error estimates in the fast multipole method for scattering problems Part 1: Truncation of the Jacobi-Anger series , 2004 .

[2]  Weng Cho Chew,et al.  Error control of the translation operator in 3D MLFMA , 2003 .

[3]  Jussi Rahola,et al.  Diagonal forms of the translation operators in the fast multipole algorithm for scattering problems , 1995 .

[4]  Mark Tygert,et al.  Fast Algorithms for Spherical Harmonic Expansions , 2006, SIAM J. Sci. Comput..

[5]  Eric F Darve,et al.  A fast multipole method for Maxwell equations stable at all frequencies , 2004, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[6]  Vikram Jandhyala,et al.  Integration and interpolation based on fast spherical transforms for the multilevel fast multipole method , 2006 .

[7]  Jian-Ming Jin,et al.  Fast and Efficient Algorithms in Computational Electromagnetics , 2001 .

[8]  V. Rokhlin,et al.  The fast multipole method (FMM) for electromagnetic scattering problems , 1992 .

[9]  V. Rokhlin.  Diagonal Forms of Translation Operators for the Helmholtz Equation in Three Dimensions , 1993 .

[10]  N. Sneeuw,et al.  Global spherical harmonic computation by two-dimensional Fourier methods , 1996 .

[11]  Jukka Sarvas,et al.  Performing Interpolation and Anterpolation Entirely by Fast Fourier Transform in the 3-D Multilevel Fast Multipole Algorithm , 2003, SIAM J. Numer. Anal..

[12]  Jiming Song,et al.  Error Analysis for the Numerical Evaluation of the Diagonal Forms of the Scalar Spherical Addition Theorem , 1999 .

[13]  Eric Darve,et al.  The Fast Multipole Method , 2000 .

[14]  W. Marsden.  I and J , 2012 .

[15]  Rainer Kußmaul.  Clenshaw-Curtis quadrature with a weighting function , 2005, Computing.

[16]  Eric Darve,et al.  The Fast Multipole Method I: Error Analysis and Asymptotic Complexity , 2000, SIAM J. Numer. Anal..

[17]  V. Rokhlin.  Rapid Solution of Integral Equations of Scattering Theory , 1990 .

[18]  Francis Collino,et al.  ERROR ESTIMATES IN THE FAST MULTIPOLE METHOD FOR SCATTERING PROBLEMS PART 2 : TRUNCATION OF THE GEGENBAUER SERIES , 2005 .

[19]  Oscar P. Bruno,et al.  A fast, higher-order solver for scattering by penetrable bodies in three dimensions , 2005 .

[20]  Eric Darve,et al.  Efficient fast multipole method for low-frequency scattering , 2004 .

[21]  R. Coifman,et al.  The fast multipole method for the wave equation: a pedestrian prescription , 1993, IEEE Antennas and Propagation Magazine.

[22]  V. Rokhlin,et al.  Fast algorithms for polynomial interpolation, integration, and differentiation , 1996 .

[23]  Bradley K. Alpert,et al.  A Fast Spherical Filter with Uniform Resolution , 1997 .

[24]  A. D. McLaren,et al.  Optimal numerical integration on a sphere , 1963 .

[25]  Reiji Suda,et al.  A fast spherical harmonics transform algorithm , 2002, Math. Comput..

[26]  John J. Knab,et al.  Interpolation of band-limited functions using the approximate prolate series (Corresp.) , 1979, IEEE Trans. Inf. Theory.

[27]  D. Owen.  Handbook of Mathematical Functions with Formulas , 1965 .

[28]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[29]  D. Healy,et al.  Computing Fourier Transforms and Convolutions on the 2-Sphere , 1994 .

[30]  Eric F Darve.  The Fast Multipole Method , 2000 .

[31]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[32]  Sean S. B. Moore,et al.  FFTs for the 2-Sphere-Improvements and Variations , 1996 .

[33]  L. Gurel,et al.  Optimal Interpolation of Translation Operator in Multilevel Fast Multipole Algorithm , 2006, IEEE Transactions on Antennas and Propagation.