Parallel implementation of an adaptive and parameter-free N-body integrator

Published: 01-05-2011| Version 1 | DOI: 10.17632/ywgjysmp94.1
C. David Pruett,
William H. Ingham,
Ralph D. Herman


Abstract Previously, Pruett et al. (2003) [3] described an N-body integrator of arbitrarily high order M with an asymptotic operation count of O(M ^2 N ^2 ). The algorithm's structure lends itself readily to data parallelization, which we document and demonstrate here in the integration of point-mass systems subject to Newtonian gravitation. High order is shown to benefit parallel efficiency. The resulting N-body integrator is robust, parameter-free, hi... Title of program: PNB.f90 Catalogue Id: AEIK_v1_0 Nature of problem High accuracy numerical evalution of trajectories of N point masses each subject to Newtonian gravitation. Versions of this program held in the CPC repository in Mendeley Data AEIK_v1_0; PNB.f90; 10.1016/j.cpc.2011.01.014 This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2019)