Abstract
We propose a new rule of thumb for designing high-order composition methods for ODEs: instead of minimizing (some norm of) the principal error coefficients, simply set all the outer stages equal. This rule automatically produces families of minimum error 4th order and corrected 6th order methods, and very good standard 6th order methods, parameterized by the number of stages. Intriguingly, the most accurate methods (evaluated with the total work held fixed) have a very large number of stages.
Similar content being viewed by others
References
S. Blanes, High order numerical integrators for differential equations using composition and processing of low order methods, Appl. Numer. Math. 37 (2001) 289–306.
S. Blanes, F. Casas and J. Ros, Symplection integration with processing: A general study, SIAM J. Sci. Comput. 21 (1999) 711–727.
S. Blanes, F. Casas and J. Ros, Processing symplectic methods for near-integrable Hamiltonian systems, Celes. Mech. Dyn. Astron. 77 (2001) 17–35.
S. Blanes, F. Casas and J. Ros, High-order Runge-Kutta-Nyström geometric methods with processing, Appl. Numer. Math. 39 (2001) 245–259.
R.I. McLachlan, On the numerical integration of ordinary differential equations by symmetric composition methods, SIAM Sci. Comput. 16 (1995) 151–168.
R.I. McLachlan and P. Atela, The accuracy of symplectic integrators, Nonlinearity 5 (1991) 541–562.
M. Suzuki, Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simulations, Phys. Lett. A 146 (1990) 319–323.
H. Yoshida, Construction of higher order symplectic integrators, Phys. Lett. A 150 (1990) 262–269.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
McLachlan, R.I. Families of High-Order Composition Methods. Numerical Algorithms 31, 233–246 (2002). https://doi.org/10.1023/A:1021195019574
Issue Date:
DOI: https://doi.org/10.1023/A:1021195019574