| 41 流水算法Grav--部分II |
 
|
- do k=1, (N-1)/2 ! 对计算拷贝的移动循环(Loop over shifts of circulating
copy)
- !右移计算数组Xc Mc Ac(Shift Circulating Arrays Xc Mc Ac to the right)
- Xc = cshift (Xc, dim=2, shift = -1)
- Mc = cshift (Mc, dim=2, shift = -1)
- Ac = cshift (Ac, dim=2, shift = -1)
- ! 对3维坐标计算距离R
- D= Xc-X
- R = sqrt (spread (sum (D*D, dim=1), dim=1, ncopies=3)))
- D = D/R**3
- A = A + Mc*D !确定的加速度( fixed acceleration)
- Ac = Ac - M*D ! 计算加速度
- end do
| Copyright: NPACT |
|