Error estimates in the fast multipole method for scattering problems. Part 1: Truncation of the Jacobi-Anger series

Abstract

We perform a complete study of the truncation error of the Jacobi-Anger series. This series expands every plane wave eiš·υ→ in terms of spherical harmonics {Yℓ,m,(š)} |m|≤ℓ≤∞. We consider the truncated series where the summation is performed over the (ℓ,m)'s satisfying |m| ≤ ℓ ≤ L. We prove that if υ = |υ→| is large enough, the truncated series gives rise to an error lower than ε as soon as L satisfies L + 1/2 ≃ υ + CW2/3 (Kε-δυγ) υ1/3 where W is the Lambert function and C, K, δ, γ are pure positive constants. Numerical experiments show that this asymptotic is optimal. Those results are useful to provide sharp estimates for the error in the fast multipole method for scattering computation.