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