Abstract
A P3M (particle-particle, particle-mesh) algorithm to compute the gravitational force on a set of particles is described. The gravitational force is computed using fast Fourier transforms. This leads to an incorrect force when the distance between two particles is of the order of a grid cell. This incorrect force is subtracted exactly from all particles in parallel using convolution with the appropriate Green's function in real space in a time of order N T, irrespective of the degree of clustering of particles. Next, the correct 1/ r 2 force is added for all neighbouring particles in parallel, leading to an accurate algorithm which runs efficiently on a highly parallel computer. A full force calculation for 128K particles on a 128 3 grid in a mildly clustered situation requires approximately 196 seconds on an 8k Connection Machine 2 with 8MHz clock. This decreases to an estimated 9.8 seconds on a full-sized 64k CM200.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
