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