On the Accuracy of Poisson's Formula Based N-Body Algorithms

We study the accuracy{cost tradeo s of a Poisson's formula based hierarchical N{body method. The parameters that control the degree of approximation of the computational elements and the separateness of interacting elements, govern both the arithmetic complexity and the accuracy of the method. Empirical models for predicting the execution time and the accuracy of the potential and force evaluations for three-dimensional problems are presented. We demonstrate how these models can be used to minimize the execution time for a prescribed error and verify the predictions through simulations on particle systems with up to one million particles. An interesting observation is that for a given error, de ning the near{ eld to consist of only nearest neighbor elements yields a lower computational complexity for a given error than the two{element separation recommended in the literature. We also show that the particle distribution may have a signi cant impact on the error.

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

[2]  L Greengard,et al.  On the Efficient Implementation of the Fast Multipole Algorithm. , 1988 .

[3]  Christopher R. Anderson,et al.  An Implementation of the Fast Multipole Method without Multipoles , 1992, SIAM J. Sci. Comput..

[4]  Michael S. Warren,et al.  A portable parallel particle program , 1995 .

[5]  V. Rokhlin,et al.  Rapid Evaluation of Potential Fields in Three Dimensions , 1988 .

[6]  Yu Hu,et al.  Implementing O(N) N-Body Algorithms Efficiently in Data-Parallel Languages , 1996, Sci. Program..

[7]  K. Esselink.  A comparison of algorithms for long-range interactions , 1995 .

[8]  J. CARRIERt,et al.  A FAST ADAPTIVE MULTIPOLE ALGORITHM FOR PARTICLE SIMULATIONS * , 2022 .

[9]  William Gropp,et al.  A Parallel Version of the Fast Multipole Method-Invited Talk , 1987, PPSC.

[10]  Piet Hut,et al.  A hierarchical O(N log N) force-calculation algorithm , 1986, Nature.

[11]  Jacob Katzenelson.  Computational structure of the N-body problem , 1989 .

[12]  Feng Zhao,et al.  The Parallel Multipole Method on the Connection Machine , 1991, SIAM J. Sci. Comput..

[13]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[14]  Yu Hu,et al.  A Data-Parallel Implementation of Hierarchical N-Body Methods , 1996, Int. J. High Perform. Comput. Appl..

[15]  K. Schmidt,et al.  Implementing the fast multipole method in three dimensions , 1991 .

[16]  Jiro Shimada,et al.  Performance of fast multipole methods for calculating electrostatic interactions in biomacromolecular simulations , 1994, J. Comput. Chem..