Reducing roundoff errors in numerical integration of planetary ephemeris

Abstract

Modern lunar-planetary ephemerides are numerically integrated on the observational timespan of more than 100 years (with the last 20 years having very precise astrometrical data). On such long timespans, not only finite difference approximation errors, but also the accumulating arithmetic roundoff errors become important because they exceed random errors of high-precision range observables of Moon, Mars, and Mercury. One way to tackle this problem is using extended-precision arithmetics available on x86 processors. Noting the drawbacks of this approach, we propose an alternative: using double-double arithmetics where appropriate. This will allow to use only double precision floating-point primitives which have ubiquitous support.

References in corpus (6)

• Lunar laser ranging in infrfared at hte Grasse laser station
• Masses of the Main Asteroid Belt and the Kuiper Belt from the Motions of Planets and Spacecraft
• Determining parameters of Moon's orbital and rotational motion from LLR observations using GRAIL and IERS-recommended models
• Mass of the Kuiper Belt
• Role of lunar laser ranging in realization of terrestrial, lunar, and ephemeris reference frames
• On the Extension of Adams--Bashforth--Moulton Methods for Numerical Integration of Delay Differential Equations and Application to the Moon's Orbit