辛积分器 Mercury:错误报告
天体物理学
2008-08-05 v1
摘要
我们报告在天文学动力学问题中使用的混合辛积分器 MERCURY 中发现的一个问题。用于跟踪天体状态的变量未初始化,这可能导致天体以非物理的方式从模拟中消失。某些 FORTRAN 编译器会隐式初始化变量,从而避免模拟出现此问题。对于其他编译器,若最大天体数量参数值设置得足够大,模拟也不会受到影响。否则,该问题会在积分器启动后的第一个事件时显现,无论是从头开始还是继续之前停止的模拟。尽管该问题在某些条件下不会显现,但显式初始化该变量可以永久且无条件地解决此问题。
引用
@article{arxiv.0808.0483,
title = {Symplectic Integrator Mercury: Bug Report},
author = {K. de Souza Torres and D. R. Anderson},
journal= {arXiv preprint arXiv:0808.0483},
year = {2008}
}
备注
4 pages, 2 figures, 1 table