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

• Published in: Computer physics communications Vol. 182; no. 5; pp. 1187 - 1198
• Main Authors: Pruett, C. David, Ingham, William H., Herman, Ralph D.
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01.05.2011
• Subjects: Computation; Computer simulation; Gravitation; Initial-value problem; Integrators; Maclaurin series; Mathematical models; N-body problem; Ordinary differential equations; Parallel computation; Parallel processing; Power series; Processors; Summaries; UNIX (operating system); Ordinary differential equations; N-body problem; Initial-value problem; Parallel computation; Maclaurin series; Power series
• ISSN: 0010-4655; 1879-2944
• DOI: 10.1016/j.cpc.2011.01.014

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, highly accurate, and adaptive in both time-step and order. Moreover, it exhibits linear speedup on distributed parallel processors, provided that each processor is assigned at least a handful of bodies. Program title: PNB.f90 Catalogue identifier: AEIK_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEIK_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC license, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 3052 No. of bytes in distributed program, including test data, etc.: 68 600 Distribution format: tar.gz Programming language: Fortran 90 and OpenMPI Computer: All shared or distributed memory parallel processors Operating system: Unix/Linux Has the code been vectorized or parallelized?: The code has been parallelized but has not been explicitly vectorized. RAM: Dependent upon N Classification: 4.3, 4.12, 6.5 Nature of problem: High accuracy numerical evaluation of trajectories of N point masses each subject to Newtonian gravitation. Solution method: Parallel and adaptive extrapolation in time via power series of arbitrary degree. Running time: 5.1 s for the demo program supplied with the package.