倾角函数的四精度计算试验

A Quad-precision Calculation Test for Inclination Function

倾角函数是天体力学分析理论中一种常用的函数.当把摄动方程展开成时间和根数的形式时需要用到.历史上提出了很多经典的倾角函数递推算法,并在双精度平台下开发了Fortran程序.进行了1次四精度计算倾角函数的试验,结果表明:L平面递推方法的四精度计算精度可达10^(-22),计算速度比双精度Jacobi方法快6倍.

In the analytical theories of celestial mechanics, the inclination function is a kind of frequently-used function when we expand the perturbation function as a form of time and orbit elements. A lot of classic recursion methods were brought out in history, and their Fortran-code were developed on the double-precision computer platform. While in this paper, a test of quad-precision calculation are reported. The results show that the quad-precision calculation precision of the L-plane recursion method can reach 10-22, and its computation speed is six times faster than the double-precision calculation of Jacobi method.