Mathematics – Numerical Analysis
Scientific paper
1998-05-20
Mathematics
Numerical Analysis
24 pages, 3 figures
Scientific paper
10.1006/jcph.1998.6066
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 c^2 + s^2 <> 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 c^2 + s^2 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.
Hénon Michel
Petit J-M
No associations
LandOfFree
Good rotations does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Good rotations, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Good rotations will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-498587