Treffer: A multipole based treecode using spherical harmonics for potentials of the form r-λ
CC BY 4.0
Sauf mention contraire ci-dessus, le contenu de cette notice bibliographique peut être utilisé dans le cadre d’une licence CC BY 4.0 Inist-CNRS / Unless otherwise stated above, the content of this bibliographic record may be used under a CC BY 4.0 licence by Inist-CNRS / A menos que se haya señalado antes, el contenido de este registro bibliográfico puede ser utilizado al amparo de una licencia CC BY 4.0 Inist-CNRS
Weitere Informationen
In this paper we describe an efficient algorithm for computing the potentials of the form r-λ where λ > 1. This treecode algorithm uses spherical harmonics to compute multipole coefficients that are used to evaluate these potentials. The key idea in this algorithm is the use of Gegenbauer polynomials to represent r-λ in a manner analogous to the use of Legendre polynomials for the expansion of the Coulomb potential r-1. We exploit the relationship between Gegenbauer and Legendre polynomials to come up with a natural generalization of the multipole expansion theorem used in the classical fast multipole algorithm [2]. This theorem is used with a hierarchical scheme to compute the potentials. The resulting algorithm has known error bounds and can be easily implemented with modification to the existing fast multipole algorithm. The complexity of the algorithm is O(p3N log N) and has several advantages over the existing Cartesian coordinates based expansion schemes.