Composition constants for raising the orders of unconventional schemes for ordinary differential equations

Kahan, William; Li, Ren-Cang

doi:10.1090/S0025-5718-97-00873-9

Composition constants for raising the orders of unconventional schemes for ordinary differential equations
HTML articles powered by AMS MathViewer

by William Kahan and Ren-Cang Li PDF

Math. Comp. 66 (1997), 1089-1099 Request permission

Abstract:

Many models of physical and chemical processes give rise to ordinary differential equations with special structural properties that go unexploited by general-purpose software designed to solve numerically a wide range of differential equations. If those properties are to be exploited fully for the sake of better numerical stability, accuracy and/or speed, the differential equations may have to be solved by unconventional methods. This short paper is to publish composition constants obtained by the authors to increase efficiency of a family of mostly unconventional methods, called reflexive.

References

E. Hairer, S. P. Nørsett, and G. Wanner, Solving ordinary differential equations. I, 2nd ed., Springer Series in Computational Mathematics, vol. 8, Springer-Verlag, Berlin, 1993. Nonstiff problems. MR 1227985
E. Hairer and G. Wanner, Solving ordinary differential equations. II, Springer Series in Computational Mathematics, vol. 14, Springer-Verlag, Berlin, 1991. Stiff and differential-algebraic problems. MR 1111480, DOI 10.1007/978-3-662-09947-6
Nicholas J. Higham, The accuracy of floating point summation, SIAM J. Sci. Comput. 14 (1993), no. 4, 783–799. MR 1223274, DOI 10.1137/0914050
W. Kahan, Relaxation methods for solving systems of ordinary differential equations, manuscript, CS Division, Department of EECS, University of California at Berkeley, Oct. 1977.
W. Kahan, Analysis and application of simply compensated summation, work in progress, CS Division, Department of EECS, University of California at Berkeley, Oct. 1993.
W. Kahan, Unconventional numerical methods for trajectory calculations, lectures notes, CS Division, Department of EECS, University of California at Berkeley, Oct. 1993.
J. D. Lambert, Numerical methods for ordinary differential systems, John Wiley & Sons, Ltd., Chichester, 1991. The initial value problem. MR 1127425
R.-C. Li, Raising the Orders of Unconventional Schemes for Ordinary Differential Equations, PhD thesis, Department of Mathematics, University of California at Berkeley, CA, 1995.
Robert I. McLachlan, On the numerical integration of ordinary differential equations by symmetric composition methods, SIAM J. Sci. Comput. 16 (1995), no. 1, 151–168. MR 1311683, DOI 10.1137/0916010
J. M. Sanz-Serna and L. Abia, Order conditions for canonical Runge-Kutta schemes, SIAM J. Numer. Anal. 28 (1991), no. 4, 1081–1096. MR 1111455, DOI 10.1137/0728058
Masuo Suzuki, General theory of fractal path integrals with applications to many-body theories and statistical physics, J. Math. Phys. 32 (1991), no. 2, 400–407. MR 1088360, DOI 10.1063/1.529425
Masuo Suzuki, General theory of higher-order decomposition of exponential operators and symplectic integrators, Phys. Lett. A 165 (1992), no. 5-6, 387–395. MR 1169635, DOI 10.1016/0375-9601(92)90335-J
Haruo Yoshida, Construction of higher order symplectic integrators, Phys. Lett. A 150 (1990), no. 5-7, 262–268. MR 1078768, DOI 10.1016/0375-9601(90)90092-3