Statistics – Computation
Scientific paper
Sep 2002
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2002dda....33.0803c&link_type=abstract
American Astronomical Society, DDA Meeting #33, #08.03; Bulletin of the American Astronomical Society, Vol. 34, p.938
Statistics
Computation
Scientific paper
Symplectic integrators are widely used for addressing problems in celestial mechanics. To date, the most popular algorithms are the 2nd-order (leapfrog) and 4th-order methods, which require one and three force evaluations per timestep respectively. Of these, the 4th-order integrator generally gives better performance for a given amount of computational effort, but its performance is limited by the fact that some of its substeps have negative coefficients, i.e. they go backwards compared to the direction of the main integration. This implies that the other (positive) coefficients must be large to compensate. Here I describe a new approach to designing symplectic integrators that appears to have been overlooked to date. This approach makes it possible to design a 4th-order integrator with coefficients that all have magnitudes less than unity, substantially reducing the size of the leading error term compared to the conventional algorithm. I will also describe a new 3rd-order integrator with very small error terms which has no analogue among conventional symplectic algorithms. The performance of the new algorithms will be compared with the standard 2nd and 4th-order methods for several test problems.
No associations
LandOfFree
A New Direction for N-body Integrators 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 A New Direction for N-body Integrators, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and A New Direction for N-body Integrators will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1043266