Good rotations
M. Henon, J-M. Petit
Abstract
Numerical integrations in celestial mechanics often involve the repeated computation of a rotation with a constant angle. A direct evaluation of these rotations yields a linear drift of the distance to the origin. This is due to roundoff in the representation of the sine s and cosine c of the angle theta. In a computer, one generally gets c2 + s2 <> 1, resulting in a mapping that is slightly contracting or expanding. In the present paper we present a method to find pairs of representable real numbers s and c such that c2 + s2 is as close to 1 as possible. We show that this results in a drastic decrease of the systematic error, making it negligible compared to the random error of other operations. We also verify that this approach gives good results in a realistic celestial mechanics integration.
Create a lesson
Related papers
Fast and Stable Arbitrary Precision Approximate Matrix Multiplication
Ioana Dumitriu, Olga Holtz, Abraham Hsu et al.
Analysis of an inexact Schwarz domain decomposition method in 3D
Nibedita Ghosh
Energy conserving compatible space--time FEM for deterministic and stochastic rotating shallow water
Maneesh Kumar Singh, Colin J. Cotter
Preconditioning proximal Galerkin: mesh and degree robust solvers for variational inequalities
Ioannis P. A. Papadopoulos, Pablo D. Brubeck
On a maximum principle in discretization with spherical harmonics
Arnold Reusken
Regularisation and Parameter Choice Rules for Composite Linear Inverse Problems
Christina Brandt, Patrick Horn