我必须在电磁场的作用下解决带电粒子的运动方程。由于我必须处理速度超过精度,我无法使用自适应步长算法(如Runge-Kutta Cash-Karp),因为它们需要花费太多时间。我正在寻找一种算法,它既是辛(如Boris积分)又是指数拟合的(为了解决运动方程,即使方程是刚性的)。我找到了一种方法,但它适用于二阶微分方程:
https://www.math.purdue.edu/~xiaj/work/SEFRKN.pdf
后来我发现了一篇论文,它将描述一个四阶辛指数拟合Runge-Kutta:
http://users.ugent.be/~gvdbergh/files/publatex/annals1.pdf
由于我必须处理速度,我一直在寻找低阶算法。是否存在二阶辛指数拟合ODE算法?