Year: 2022
Author: Robert I. McLachlan
Communications in Computational Physics, Vol. 31 (2022), Iss. 3 : pp. 987–996
Abstract
Many applications in computational physics that use numerical integrators based on splitting and composition can benefit from the development of optimized algorithms and from choosing the best ordering of terms. The cost in programming and execution time is minimal, while the performance improvements can be large. In this note we report the influence of term ordering for random systems and for two systems from celestial mechanics that describe particle paths near black holes, quantifying its significance for both optimized and unoptimized methods. We also present a method for the computation of solutions of integrable monomial Hamiltonians that minimizes roundoff error and allows the effective use of compensation summation.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/10.4208/cicp.OA-2021-0154
Communications in Computational Physics, Vol. 31 (2022), Iss. 3 : pp. 987–996
Published online: 2022-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 10
Keywords: Symplectic integrators celestial mechanics splitting methods.