Hostname: page-component-77c89778f8-cnmwb Total loading time: 0 Render date: 2024-07-21T20:16:25.438Z Has data issue: false hasContentIssue false

Parallel/Vector Integration Methods for Dynamical Astronomy

Published online by Cambridge University Press:  12 April 2016

Toshio Fukushima*
Affiliation:
National Astronomical Observatory 2-21-1, Ohsawa, Mitaka, Tokyo 181-8588, Japantoshio@nao.ac.jp

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 paper reviews three recent works on the numerical methods to integrate ordinary differential equations (ODE), which are specially designed for parallel, vector, and/or multi-processor-unit (PU) computers. The first is the Picard-Chebyshev method (Fukushima, 1997a). It obtains a global solution of ODE in the form of Chebyshev polynomial of large (> 1000) degree by applying the Picard iteration repeatedly. The iteration converges for smooth problems and/or perturbed dynamics. The method runs around 100-1000 times faster in the vector mode than in the scalar mode of a certain computer with vector processors (Fukushima, 1997b). The second is a parallelization of a symplectic integrator (Saha et al., 1997). It regards the implicit midpoint rules covering thousands of timesteps as large-scale nonlinear equations and solves them by the fixed-point iteration. The method is applicable to Hamiltonian systems and is expected to lead an acceleration factor of around 50 in parallel computers with more than 1000 PUs. The last is a parallelization of the extrapolation method (Ito and Fukushima, 1997). It performs trial integrations in parallel. Also the trial integrations are further accelerated by balancing computational load among PUs by the technique of folding. The method is all-purpose and achieves an acceleration factor of around 3.5 by using several PUs. Finally, we give a perspective on the parallelization of some implicit integrators which require multiple corrections in solving implicit formulas like the implicit Hermitian integrators (Makino and Aarseth, 1992), (Hut et al., 1995) or the implicit symmetric multistep methods (Fukushima, 1998), (Fukushima, 1999).

Type
Analytical and Numerical Tools
Copyright
Copyright © Kluwer 1999

References

Fukushima, T.: 1997 a, ‘Picard Iteration Method, Chebyshev Polynomial Approximation, and Global Numerical Integration of Dynamical Motions’, Astron. J., 113, 19091914.Google Scholar
Fukushima, T.: 1997b, ‘Vector Integration of Dynamical Motions by the Picard-Chebyshev Method’, Astron. J., 113, 23252328.CrossRefGoogle Scholar
Fukushima, T.: 1998, ‘Symmetric Multistep Methods Revisited’, in Proc. 30th Symp. on Cele. Mech. (Fukushima, et al. eds), 229247.Google Scholar
Fukushima, T.: 1999, ‘Symmetric Multistep Methods Revisited II: Numerical Experiments’, in Proc. IAU Coll. No. 173, to be printed.CrossRefGoogle Scholar
Hairer, E., Nøsett, S.P., and Wanner, G.: 1993, Solving Ordinary Differential Equations I (2nd ed.), Springer-Verlag, Berlin.Google Scholar
Hut, P., Makino, J., and McMillan, S.: 1995, Astrophys. J. Lett., 443, 93-.Google Scholar
Ito, T., and Fukushima, T.: 1997, ‘Parallelized Extrapolation Method and Its Application to the Orbital Dynamics’, Astron. J., 114, 12601267.Google Scholar
Kinoshita, H., Yoshida, H., and Nakai, H.: 1991, Cele. Mech. and Dyn. Astr., 50, 59.Google Scholar
Lambert, J.D., and Watson, I.A.: 1976, ‘Symmetric Multistep Methods for Periodic Initial Value Problems’, J. Inst. Maths Applics, 18, 189202.Google Scholar
Makino, J., and Aarseth, S.J.: 1992, Publ. Astron. Soc. Japan, 44, 141.Google Scholar
Miranker, W.L., and Liniger, W.: 1967, ‘Parallel Methods for the Numerical Integration of Ordinary Differential Equations’, Math. Comp., 21, 303320.Google Scholar
Quinlan, G.D., and Tremaine, S., 1990, ‘Symmetric Multistep Methods for the Numerical Integration of Planetary Orbits’, Astron. J., 100, 16941700.CrossRefGoogle Scholar
Saha, P., Stadel, J., and Tremaine, S.: 1997, ‘A Parallel Integration Method for Solar System Dynamics’, Astron. J., 114, 409415.Google Scholar
Wisdom, J., and Holman, M.: 1991, ‘Symplectic Maps for the N-body Problem’, Astron. J., 102, 15281538.CrossRefGoogle Scholar