Parallel P3M with exact calculation of short range forces

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 a $8k$ Connection Machine 2 with 8MHz clock. This decreases to an estimated 9.8 seconds on a full-sized $64k$ CM200.