| 42 Grav流水算法,结论 |
 
|
- if (( N mod 2) = 0 ) then !如果N是偶数 ,最后一个方法计算加速度
- Xc = cshift (Xc, dim=2, shift = -1)
- Mc = cshift (Mc, dim=2, shift = -1)
- Ac = cshift (Ac, dim=2, shift = -1)
- D = Xc-X
- R = sqrt (spread (sum (D*D, dim=1), dim=1, ncopies=3)))
- D = D/R**3
- A = A + Mc*D
- end if
- ! 对最后的结果组合加速度-在第i次计算粒子
- ! 相应于(i-(N-1)/2)粒子的位置
- Grav = A + cshift (Ac, dim=2, shift = (N-1)/2)
- 结束函数Grav
| Copyright: NPACT |
|