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

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)

摘要
此前，Pruett等人（2003）[3] 提出了一种任意高阶M的N体数值积分器（N-body integrator），其渐进运算复杂度为$O(M^2N^2)$。该算法的结构天然适配数据并行计算，本文针对受牛顿引力作用的质点系统积分场景，对该并行特性进行了阐述与验证。研究表明，高阶特性可有效提升并行计算效率。所得N体数值积分器具备鲁棒性、无参数依赖，且具有（原文截断）……

程序名称：PNB.f90
目录编号：AEIK_v1_0

问题属性
对每个受牛顿引力作用的N个质点的运动轨迹进行高精度数值求解。

本程序在Mendeley数据平台的CPC仓库中的版本信息：
AEIK_v1_0；PNB.f90；10.1016/j.cpc.2011.01.014

本程序源自贝尔法斯特女王大学托管的CPC程序库（1969-2019年）