1 | 1)XdeltaA(i) = V(i) |
2 | Xuse(i) = X(i) running over all i in Processor |
3 | 2)VdeltaA(i) = MPGrav(i) involves communication |
4 | 3)XdeltaB(i) = V(i) + h*VdeltaA(i)/2 |
5 | Xuse(i) = X(i) +h*XdeltaA(i)/2 |
6 | 4)VdeltaB(i) = MPGrav(i) involves communication |
7 | 5)XdeltaC(i) = V(i) + h*VdeltaB(i)/2 |
8 | Xuse(i) = X(i) +h*XdeltaB(i)/2 |
9 | 6)VdeltaC(i) = MPGrav(i) involves communication |
10 | 7)XdeltaD(i) = V(i) + h*VdeltaC(i) |
11 | Xuse(i) = X(i) +h*XdeltaC(i) |
12 | 8)VdeltaD(i) = MPGrav(i) involves communication |
13 | 9)X(i) becomes (X(i) + h*(XdeltaA(i)+2*XdeltaB(i)+2*XdeltaC(i)+XdeltaD(i))/6 ) |
14 | V(i) becomes (V(i) + h*(VdeltaA(i)+2*VdeltaB(i)+2*VdeltaC(i)+VdeltaD(i))/6 ) |