| 1 |
MPGrav(i) uses Xuse(i) and mass M(i) and calculates a 3 vector force from 3 vector Xuse(i,j)
|
| 2 |
MPGrav(i) = S M(i)*M(j)*( Xuse(j)-Xuse(i))
|
| 3 |
where ri,j = | Xuse(i)-Xuse(j) | is distance between particle i and particle j
|
| 4 |
This calculation involves communication as N-Nproc of values of j (and parameters Xuse(j) M(j) ) are stored outside the processor holding the i'th processor
|
| 5 |
We describe algorithms with increasing complexity and efficiency!
|
| 6 |
Note in discussions, we are rather sloppy as to whether i,j are "local" (1...N/Nproc) or "global" (1...N) indices
|
