Compression of particle data from hierarchical approximate methods
ACM Transactions on Mathematical Software (TOMS)
Variable order revised binary treecode
Journal of Computational Physics
Journal of Computational Physics
Parallel multigrid summation for the N-body problem
Journal of Parallel and Distributed Computing
Improving boundary element methods for parasitic extraction
ASP-DAC '03 Proceedings of the 2003 Asia and South Pacific Design Automation Conference
| Hi-index | 0.01 |
Rapid evaluation of potentials in particle systems is an important, time-consuming step in many physical simulations. Over the past decade, the development of treecodes, such as the fast multipole method and Barnes--Hut method, has enabled large scale simulations in astrophysics, molecular dynamics, material science, etc. These methods use fixed-degree polynomial approximations of the potential at a point, due to a set of particles in a hierarchical representation of the particle system. In this paper, we present analysis and experiments to illustrate that fixed-degree multipole approximations can lead to large aggregate errors. We describe an alternate strategy based on careful selection of the multipole degree that leads to asymptotically lower errors while incurring minimal computation overhead. First, we estimate the error associated with each particle-cluster interaction and the aggregate error for each particle. Then we describe a technique for computing the multipole degree of each interaction with a view to reducing aggregate error, and establish the computational complexity of the new method. Following this, we discuss numerical experiments that demonstrate favorable error properties of the new method at the expense of marginal increase in computation. Finally, we report a parallel implementation on the SGI Origin 2000 that achieves excellent speedup for problems from domains such as astrophysics and boundary element solvers.