TY - JOUR
T1 - A cost optimal parallel algorithm for computing force field in N-body simulations on a CREW PRAM
AU - Xue, Guoliang
N1 - Funding Information:
Supported in part by the US Army Research OFce grant DAAH04-9610233 and by the National Science Foundation grants ASC-9409285 and OSR-9350540. E-mail address: xue@cs.uvm.edu (G. Xue).
PY - 2001
Y1 - 2001
N2 - We consider the following force field computation problem: given a cluster of n particles in three-dimensional space, compute the force exerted on each particle by the other particles. Depending on different applications, the pairwise interaction could be either gravitational or Lennard-Jones. In both cases, the force between two particles vanishes as the distance between them approaches to infinity. Since there are n(n-1)/2 pairs, direct method requires Θ(n2) time for force evaluation, which is very expensive for astronomical simulations. In 1985 and 1986, two famous Θ(nlogn) time hierarchical tree algorithms were published by Appel (SIAM J. Sci. Statist. Comput. 6 (1985) 85-103) and by Barnes and Hut (Nature, 324 (1980) 446 - 449), respectively. In a recent paper, we presented a linear time algorithm which builds the oct tree bottom-up and showed that Appel's algorithm can be implemented in Θ(n) sequential time. In this paper, we present an algorithm which computes the force field in Θ(logn) time on a Θ(n/log n) processor CREW PRAM. A key to this optimal parallel algorithm is replacing a recursive top-down force calculation procedure of Appel by an equivalent non-recursive bottom-up procedure. Our parallel algorithm also yields a new Θ(n) time-sequential algorithm for force field computation.
AB - We consider the following force field computation problem: given a cluster of n particles in three-dimensional space, compute the force exerted on each particle by the other particles. Depending on different applications, the pairwise interaction could be either gravitational or Lennard-Jones. In both cases, the force between two particles vanishes as the distance between them approaches to infinity. Since there are n(n-1)/2 pairs, direct method requires Θ(n2) time for force evaluation, which is very expensive for astronomical simulations. In 1985 and 1986, two famous Θ(nlogn) time hierarchical tree algorithms were published by Appel (SIAM J. Sci. Statist. Comput. 6 (1985) 85-103) and by Barnes and Hut (Nature, 324 (1980) 446 - 449), respectively. In a recent paper, we presented a linear time algorithm which builds the oct tree bottom-up and showed that Appel's algorithm can be implemented in Θ(n) sequential time. In this paper, we present an algorithm which computes the force field in Θ(logn) time on a Θ(n/log n) processor CREW PRAM. A key to this optimal parallel algorithm is replacing a recursive top-down force calculation procedure of Appel by an equivalent non-recursive bottom-up procedure. Our parallel algorithm also yields a new Θ(n) time-sequential algorithm for force field computation.
KW - Cost optimal algorithms
KW - Force field evaluation
KW - N-body simulations
KW - PRAM
KW - Spatial tree algorithms
UR - http://www.scopus.com/inward/record.url?scp=0034916340&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=0034916340&partnerID=8YFLogxK
U2 - 10.1016/S0304-3975(99)00326-6
DO - 10.1016/S0304-3975(99)00326-6
M3 - Article
AN - SCOPUS:0034916340
SN - 0304-3975
VL - 255
SP - 555
EP - 568
JO - Theoretical Computer Science
JF - Theoretical Computer Science
IS - 1-2
ER -