|
Computation of numerical method is inherently iterative: at each time step, the solution depends on the immediately preceding one.
|
At each time step, MPGrav is called (several times as using Runge Kutta):
-
For each particle i, one must sum over the forces due to all other particles j.
-
Computation is O(N2) - there is a potential for parallel computation in either i (normal case in standard Message Parallel) or j or both.
|
|
We will use 4th order Runge Kutta to integrate in time and the program is designed an overall routine looping over time with parallelism hidden in MPGrav routines
|
|
We analyse Message Parallel version and other foils discuss Data Parallel version
|
