Hostname: page-component-76fb5796d-qxdb6 Total loading time: 0 Render date: 2024-04-27T13:24:57.995Z Has data issue: false hasContentIssue false
  • English
  • Français

An efficient integrator that uses Gauss-Radau spacings

Published online by Cambridge University Press:  12 April 2016

Edgar Everhart
Edgar Everhart*
Affiliation:
Physics Department and Chamberlin Observatory, University of Denver, Denver, Colorado 80208, USA

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the ‘Save PDF’ action button.

This describes our integrator RADAU, which has been used by several groups in the U.S.A., in Italy, and in the U.S.S.R. over the past 10 years in the numerical integration of orbits and other problems involving numerical solution of systems of ordinary differential equations. First- and second-order equations are solved directly, including the general second-order case. A self-starting integrator, RADAU proceeds by sequences within which the substeps are taken at Gauss-Radau spacings. This allows rather high orders of accuracy with relatively few function evaluations. After the first sequence the information from previous sequences is used to improve the accuracy. The integrator itself chooses the next sequence size. When a 64-bit double word is available in double precision, a 15th-order version is often appropriate, and the FORTRAN code for this case is included here. RADAU is at least comparable with the best of other integrators in speed and accuracy, and it is often superior, particularly at high accuracies.

Type
Section IV. Dynamics of Comets: Numerical Modelling
Information
International Astronomical Union Colloquium , Volume 83: Dynamics of Comets: Their Origin and Evolution , 1985 , pp. 185 - 202
Copyright
Copyright © Cambridge University Press 1985

References

Batrakov, Yu.N.: 1982, ‘Methods of Computation of the Perturbed Motion of Small Bodies in the Solar System’ in Sun and Planetary System, (Proceedings of the VI European Meeting in Astronomy), Fricke, W. and Teleki, G., editors, D. Reidel Publishing Co., Dordrecht. pp415–419.Google Scholar
Bettis, D.G. and Szebehely, V.: 1972, ‘Treatment of Close Approaches in the Numerical Integration of the Gravitational Problem of N-bodies’ Gravitational N-Body Problem, Lecar, M., editor, D. Reidel Publishing Co., Dordrecht, pp388–405. See p395.CrossRefGoogle Scholar
Bulirsch, R. and Stoer, J.: 1966 ‘Numerical Treatment of Ordinary Differential Equations by Extrapolatiom Methods’ Num. Math., 8, pp1–13.Google Scholar
Butcher, J.C.: 1964, ‘Integration Processes Based on Radau Quadrature Formulas’, Math. Comp. 18, pp233–344.CrossRefGoogle Scholar
Carusi, A., Kresak, L., Perozzi, E., and Valsecchi, G.B.: 1985 ‘The Long-Term Evolution Project’ in Dynamics of Comets; Their Origin and Evolution, Carusi, A. and Valsecchi, G.B., editors, D. Reidel Publishing Co., Dordrecht. (IAU Colloq. 83, Rome, June 1984) This volume.Google Scholar
Eckert, W.J., Brouwer, D., and Clemence, G.M.: 1951 ‘Coordinates of the Five Outer Planets, 1653-2060’ Astron. Papers Am. Ephemeris, 12.Google Scholar
Everhart, E.: 1974a, ‘Implicit Single-Sequence Methods for Integrating Orbits’, Celest. Mech. 10, pp35–55 Google Scholar
Everhart, E.: 197aa, ‘An Efficient Integrator of Very High Order and Accuracy’ Denver Res. Inst. Tech. Report, 1 July 1974 (unpublished).Google Scholar
Everhart, E. and Marsden, B.G.: 1983, ‘New Original and Future Cometary Orbits’ Astron. J., 88, pp135–137.Google Scholar
Fehlberg, E.: 1972, ‘Classical Eighth- and Lower Order Runge-Kutta-Nystrom Formulas with Stepsize Control for Special Second-Order Differential Equations’, NASA Tech. Rept. NASA TR R-381.Google Scholar
Gallaher, L.J. and Perlin, I.E.: 1966, ‘A comparison of Several Methods of Numerical Integration of Nonlinear Differential Equations’. Presented at the SIAM meeting, Univ. of Iowa, March. 1966 (unpublished). See Krogh, 1973.Google Scholar
House, F., Weiss, G., and Weigandt, R.: 1978, ‘Numerical Integration of Stellar Orbits’, Celest. Mech., 18, pp311–318.Google Scholar
Krogh, F.T.: 1973, ‘On Testing a Subroutine for Numerical Integration of Ordinary Differential Equations’, J. Assoc. Comput. Mach., 20, pp545–562.Google Scholar
Levi-Civita, T.: 1903, ‘Traiettorie singulari ed urti nel problema ristretto del tre corp’, Annali di Mat. 9 pp1–32.Google Scholar
Marsden, B.G., Sekanina, Z., and Everhart, E.: 1978, ‘New Osculating Orbits for 110 Comets and Analysis of Original Orbits for 200 Comets’, Astron. J., 83, pp64–71.Google Scholar
Newton, R.R.: 1959, ‘Periodic Orbits of a Planetoid’, Smithson. Contrib. Astrophys. 3, No. 7, pp69–78.Google Scholar
Papp, K.A., Innanen, K.A., and Patrick, A.T.: 1977, 1980, ‘A Comparison of Five Algolrithms for Numerical Orbit Computation in Galaxy Models’, Celest. Mech., 18, pp277–286, and 21, 337–349.Google Scholar
Shefer, V.A.: 1982, ‘Variational Equations of the Perturbed Two Body Problem in Regularized Form’, Institute of Theoretical Astronomy, Leningrad. Publication No. 37. An 11th-order version of RADAU is used.Google Scholar
Stroud, A.H. and Secrest, D.: 1966, Gaussian Quadrature Formulas Prentice Hall, Inc. Englewood Cliffs N.J., See Table 12 (Radau spacings).Google Scholar