First step is to Circulate(shift) one position and calculate acelerations Fij and Fji in all index positions |
Shifting pipeline (N-1) times gives correct algorithm but does not save "Newton's factor of two". |
Just need (N-1)/2 steps when N is odd and N/2 steps when N is even which saves factor of two. |
Shift in blocks of J particles for message passing to block for good performance |