帮助辛积分器

Geo*_*rge 11 math physics scientific-computing numerical-integration differential-equations

我正在尝试开发物理模拟,我想实现一个四阶辛积分方法.问题是我必须弄错数学,因为在使用辛积分器时我的模拟根本不起作用(与模拟工作相当好的四阶Runge-Kutta积分器相比).我一直在谷歌搜索这个,我能找到的只是关于这个主题的科学文章.我试图改编文章中使用的方法,但我没有运气.我想知道是否有人有使用辛积分器的模拟的源代码,最好是模拟引力场,但任何辛积分器都可以.源代码的语言并不重要,但我会欣赏使用C风格语法的语言.谢谢!

Dr.*_*ius 7

至于你问的源代码:从这里你可以下载MATLAB和FORTRAN代码为哈密顿系统和可逆的问题对称方法辛方法.还有许多处理diff方程的其他方法.

而在纸,你可以找到的算法描述.

编辑

如果您使用Mathematica,也可能有所帮助.

  • @George如你所见,辛算法的源代码示例在网上很少见.完成后,请考虑将代码发布到某个地方,以帮助其他有需要的人. (2认同)