Composition constants for raising the orders of unconventional schemes for ordinary differential equations
Authors:
William Kahan and RenCang Li
Journal:
Math. Comp. 66 (1997), 10891099
MSC (1991):
Primary 34A50, 65L05
MathSciNet review:
1423077
Fulltext PDF Free Access
Abstract 
References 
Similar Articles 
Additional Information
Abstract: Many models of physical and chemical processes give rise to ordinary differential equations with special structural properties that go unexploited by generalpurpose 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.
 1.
E.
Hairer, S.
P. Nørsett, and G.
Wanner, Solving ordinary differential equations. I, 2nd ed.,
Springer Series in Computational Mathematics, vol. 8, SpringerVerlag,
Berlin, 1993. Nonstiff problems. MR 1227985
(94c:65005)
 2.
E.
Hairer and G.
Wanner, Solving ordinary differential equations. II, Springer
Series in Computational Mathematics, vol. 14, SpringerVerlag, Berlin,
1991. Stiff and differentialalgebraic problems. MR 1111480
(92a:65016)
 3.
Nicholas
J. Higham, The accuracy of floating point summation, SIAM J.
Sci. Comput. 14 (1993), no. 4, 783–799. MR 1223274
(94a:65025), http://dx.doi.org/10.1137/0914050
 4.
W. Kahan, Relaxation methods for solving systems of ordinary differential equations, manuscript, CS Division, Department of EECS, University of California at Berkeley, Oct. 1977.
 5.
height 2pt depth 1.6pt width 23pt, Analysis and application of simply compensated summation, work in progress, CS Division, Department of EECS, University of California at Berkeley, Oct. 1993.
 6.
height 2pt depth 1.6pt width 23pt, Unconventional numerical methods for trajectory calculations, lectures notes, CS Division, Department of EECS, University of California at Berkeley, Oct. 1993.
 7.
J.
D. Lambert, Numerical methods for ordinary differential
systems, John Wiley & Sons Ltd., Chichester, 1991. The initial
value problem. MR 1127425
(92i:65114)
 8.
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.
 9.
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
(95j:65073), http://dx.doi.org/10.1137/0916010
 10.
J.
M. SanzSerna and L.
Abia, Order conditions for canonical RungeKutta schemes, SIAM
J. Numer. Anal. 28 (1991), no. 4, 1081–1096. MR 1111455
(92e:65103), http://dx.doi.org/10.1137/0728058
 11.
Masuo
Suzuki, General theory of fractal path integrals with applications
to manybody theories and statistical physics, J. Math. Phys.
32 (1991), no. 2, 400–407. MR 1088360
(92k:81096), http://dx.doi.org/10.1063/1.529425
 12.
Masuo
Suzuki, General theory of higherorder decomposition of exponential
operators and symplectic integrators, Phys. Lett. A
165 (1992), no. 56, 387–395. MR 1169635
(93b:22011), http://dx.doi.org/10.1016/03759601(92)90335J
 13.
Haruo
Yoshida, Construction of higher order symplectic integrators,
Phys. Lett. A 150 (1990), no. 57, 262–268. MR 1078768
(91h:70014), http://dx.doi.org/10.1016/03759601(90)900923
 1.
 E. Hairer, S. P. Nørsett, and G. Wanner, Solving Ordinary Differential Equations I, SpringerVerlag, New York, second ed., 1993. MR 94c:65005
 2.
 E. Hairer and G. Wanner, Solving Ordinary Differential Equations II, SpringerVerlag, New York, 1991. MR 92a:65016
 3.
 N. J. Higham, The accuracy of floating point summation, SIAM J. Sci. Comput., 14 (1993), pp. 783799. MR 94a:65025
 4.
 W. Kahan, Relaxation methods for solving systems of ordinary differential equations, manuscript, CS Division, Department of EECS, University of California at Berkeley, Oct. 1977.
 5.
 height 2pt depth 1.6pt width 23pt, Analysis and application of simply compensated summation, work in progress, CS Division, Department of EECS, University of California at Berkeley, Oct. 1993.
 6.
 height 2pt depth 1.6pt width 23pt, Unconventional numerical methods for trajectory calculations, lectures notes, CS Division, Department of EECS, University of California at Berkeley, Oct. 1993.
 7.
 J. D. Lambert, Numerical Methods for Ordinary Differential Systems, John Wiley & Sons, New York, 1991. MR 92i:65114
 8.
 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.
 9.
 R. I. McLachlan, On the numerical integration of ordinary differential equations by symmetric composition methods, SIAM Journal on Scientific Computing, 16 (1995), pp. 151168. MR 95j:65073
 10.
 J. M. SanzSerna and L. Abia, Order conditions for canonical RungeKutta schemes, SIAM Journal on Numerical Analysis, 28 (1991), pp. 10811096. MR 92e:65103
 11.
 M. Suzuki, General theory of fractal path integrals with applications to manybody theories and statistical physics, Journal of Mathematical Physics, 32 (1991), pp. 400407. MR 92k:81096
 12.
 height 2pt depth 1.6pt width 23pt, General theory of higherorder decomposition of exponential operators and symplectic integrators, Physics Letters A, 165 (1992), pp. 387395. MR 93b:22011
 13.
 H. Yoshida, Construction of higher order symplectic integrators, Physics Letters A, 150 (1990), pp. 262268. MR 91h:70014
Similar Articles
Retrieve articles in Mathematics of Computation of the American Mathematical Society
with MSC (1991):
34A50,
65L05
Retrieve articles in all journals
with MSC (1991):
34A50,
65L05
Additional Information
William Kahan
Affiliation:
Computer Science Division and Department of Mathematics, University of California at Berkeley, Berkeley, California 94720
Email:
wkahan@cs.berkeley.edu
RenCang Li
Affiliation:
Mathematical Science Section, Oak Ridge National Laboratory, P.O. Box 2008, Bldg. 6012, Oak Ridge, Tennessee 378316367
Email:
na.rcli@nanet.ornl.gov
DOI:
http://dx.doi.org/10.1090/S0025571897008739
PII:
S 00255718(97)008739
Keywords:
Ordinary differential equations,
reflexive methods,
composition schemes,
palindromic schemes
Received by editor(s):
June 10, 1996
Additional Notes:
The first author was supported in part by the Office of Naval Research contract N0001490J1372 and National Science Foundation contract ASC9005933.
The second author was supported in part by a Householder Fellowship in Scientific Computing at Oak Ridge National Laboratory, supported by the Applied Mathematical Sciences Research Program, Office of Energy Research, United States Department of Energy contract DEAC0596OR22464 with Lockheed Martin Energy Research Corporation.
Article copyright:
© Copyright 1997 American Mathematical Society
