wendy 0.1

概述

wendy解决一维引力n体问题 一种有效算法的机器精度[o（log n）]/ 粒子碰撞]。或者，它可以用 近似积分，但具有精确的力。

作者

jo bovy（多伦多大学）：bovy-at-astro-dot-utoronto- 点-ca

安装

克隆/fork/下载存储库并使用安装

sudo python setup.py install

或在本地使用

python setup.py install --user

用法

使用wendy.nbody初始化初始（x，v）的生成器对象 质量为m。然后，生成器返回系统的状态 等距时间间隔：

g= wendy.nbody(x,v,m,0.05) # delta t = 0.05
next_x, next_v= next(g) # at t=0.05
next_x, next_v= next(g) # at t=0.10
...

使用wendy.nbody初始化生成器时，可以选择（a） 精确地解决问题，或者不使用approx=，并且（b）包含 外部谐振子电势omega^2 x^2 / 2与 {TT6} $（对于精确解和近似解都）。

示例

您可以通过单击来运行这些，而无需安装wendy。 并导航到examples/目录。注意一些 电影可能无法在活页夹网页上呈现，因此您可能 在运行笔记本（或更改 subsamp输入。

• 一维的相位混合和强弛豫：example notebook（运行） 在本地观看电影，或在 nbviewer）
• 一个自引力的sech2磁盘：example notebook（本地运行到 看电影，或观看 nbviewer）
• 绝热收缩：example notebook（局部运行以查看 电影，或查看 nbviewer）
• 指数盘的绝热与非绝热能量注入： example notebook（运行 在本地观看电影，或在 nbviewer）
• wendy按粒子数缩放：example notebook（查看 nbviewer）