Note that MPGrav(i) result depends on the array Xuse and fixed parameters such as mass
-
It happens to be independent of time and velocity although sucvh velocity and time dependent forces could be accomodated
|
|
The calculation involves an iteration over desired number of time steps increasing t by h each time.
|
|
Each time step involves 9 phases -- at each phase, one loops over all N particles (or rather over all N/Nproc particles stored in a given processor)
|
|
Some phases (4 of them) involve communication and computation; the others "just" SPMD computation.
|
|
The number of phases depends on ODE strategy used -- with the simple Euler method, there are 1 or 2 phases, depending on how one counts.
|
