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 NT, irrespective of the degree of clustering of particles. Next, the correct 1/r2 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 1283 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.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…