Error estimates in the fast multipole method for scattering problems. Part 2 : truncation of the Gegenbauer series
ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Volume 39 (2005) no. 1, p. 183-221

We perform a complete study of the truncation error of the Gegenbauer series. This series yields an expansion of the Green kernel of the Helmholtz equation, $\frac{{\mathrm{e}}^{i|\stackrel{\to }{u}-\stackrel{\to }{v}|}}{4\pi i|\stackrel{\to }{u}-\stackrel{\to }{v}|}$, which is the core of the Fast Multipole Method for the integral equations. We consider the truncated series where the summation is performed over the indices $\ell \le L$. We prove that if $v=|\stackrel{\to }{v}|$ is large enough, the truncated series gives rise to an error lower than $ϵ$ as soon as $L$ satisfies $L+\frac{1}{2}\simeq v+C{W}^{\frac{2}{3}}\left(K\left(\alpha \right){ϵ}^{-\delta }{v}^{\gamma }\right)\phantom{\rule{0.166667em}{0ex}}{v}^{\frac{1}{3}}$ where $W$ is the Lambert function, $K\left(\alpha \right)$ depends only on $\alpha =\frac{|\stackrel{\to }{u}|}{|\stackrel{\to }{v}|}$ and $C\phantom{\rule{0.166667em}{0ex}},\delta ,\phantom{\rule{0.166667em}{0ex}}\gamma$ are pure positive constants. Numerical experiments show that this asymptotic is optimal. Those results are useful to provide sharp estimates of the error in the fast multipole method for scattering computation.

DOI : https://doi.org/10.1051/m2an:2005008
Classification:  33C10,  33C55,  41A80
Keywords: Gegenbauer, fast multipole method, truncation error
@article{M2AN_2005__39_1_183_0,
author = {Carayol, Quentin and Collino, Francis},
title = {Error estimates in the fast multipole method for scattering problems. Part 2 : truncation of the Gegenbauer series},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique},
publisher = {EDP-Sciences},
volume = {39},
number = {1},
year = {2005},
pages = {183-221},
doi = {10.1051/m2an:2005008},
zbl = {1087.33007},
mrnumber = {2136205},
language = {en},
url = {http://www.numdam.org/item/M2AN_2005__39_1_183_0}
}

Carayol, Quentin; Collino, Francis. Error estimates in the fast multipole method for scattering problems. Part 2 : truncation of the Gegenbauer series. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Volume 39 (2005) no. 1, pp. 183-221. doi : 10.1051/m2an:2005008. http://www.numdam.org/item/M2AN_2005__39_1_183_0/

[1] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions. Dover, New-York (1964).

[2] S. Amini and A. Profit, Analysis of the truncation errors in the fast multipole method for scattering problems. J. Comput. Appl. Math. 115 (2000) 23-33. | Zbl 0973.65092

[3] J.A. Barcelo, A. Ruiz and L. Vega, Weighted estimates for the Helmholtz equation and some applications. J. Funct. Anal. 150 (1997) 356-382. | Zbl 0890.35028

[4] H. Bateman, Higher transcendental Functions. McGraw-Hill (1953). | MR 58756

[5] Q. Carayol, Développement et analyse d'une méthode multipôle multiniveau pour l'électromagnétisme. Ph.D. thesis, Université Paris VI Pierre et Marie Curie, rue Jussieu 75005 Paris (2002).

[6] Q. Carayol and F. Collino, Error estimates in the fast multipole method for scattering problems. part 1: Truncation of the jacobi-anger series. ESAIM: M2AN 38 (2004) 371-394. | Numdam | Zbl 1077.41027

[7] T.M. Cherry, Uniform asymptotic formulae for functions with transition points. Trans. AMS 68 (1950) 224-257. | Zbl 0036.06102

[8] W.C. Chew, J.M. Jin, E. Michielssen and J.M. Song, Fast and Efficient Algorithms in Computational Electromagnetics. Artech House (2001).

[9] R. Coifman, V. Rokhlin and S. Greengard, The Fast Multipole Method for the wave equation: A pedestrian prescription. IEEE Antennas and Propagation Magazine 35 (1993) 7-12.

[10] D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory. Springer-Verlag 93 (1992). | MR 1183732 | Zbl 0760.35053

[11] E. Darve, The fast multipole method. I. Error analysis and asymptotic complexity. SIAM J. Numer. Anal. 38 (2000) 98-128 (electronic). | Zbl 0974.65033

[12] E. Darve, The fast multipole method: Numerical implementation. J. Comput. Physics 160 (2000) 196-240. | Zbl 0974.78012

[13] E. Darve and P. Havé, Efficient fast multipole method for low frequency scattering. J. Comput. Physics 197 (2004) 341-363. | Zbl 1073.65133

[14] B. Dembart and E. Yip, Accuracy of fast multipole methods for maxwell's equations. IEEE Comput. Sci. Engrg. 5 (1998) 48-56.

[15] M.A. Epton and B. Dembart, Multipole translation theory for the three-dimensional Laplace and Helmholtz equations. SIAM J. Sci. Comput. 16 (1995) 865-897. | Zbl 0852.31006

[16] I.S. Gradshteyn, I.M. Ryzhik, Table of integrals, series, and products, 5th edition. Academic Press (1994). | MR 1243179 | Zbl 0918.65002

[17] S. Koc, J. Song and W.C. Chew, Error analysis for the numerical evaluation of the diagonal forms of the scalar spherical addition theorem. SIAM J. Numer. Anal. 36 (1999) 906-921 (electronic). | Zbl 0924.65116

[18] L. Lorch, Alternative proof of a sharpened form of Bernstein's inequality for Legendre polynomials. Applicable Anal. 14 (1982/83) 237-240. | Zbl 0505.33007

[19] L. Lorch, Corrigendum: “Alternative proof of a sharpened form of Bernstein's inequality for Legendre polynomials” [Appl. Anal. 14 (1982/83) 237-240; MR 84k:26017]. Appl. Anal. 50 (1993) 47. | Zbl 0505.33007

[20] J.C. Nédélec, Acoustic and Electromagnetic Equation. Integral Representation for Harmonic Problems. Springer-Verlag 144 (2001). | MR 1822275 | Zbl 0981.35002

[21] J. Rahola, Diagonal forms of the translation operators in the fast multipole algorithm for scattering problems. BIT 36 (1996) 333-358. | Zbl 0854.65122

[22] G.N. Watson, A treatise on the theory of Bessel functions. Cambridge University Press (1966). | JFM 48.0412.02 | MR 1349110 | Zbl 0174.36202