GniCodes – Matlab programs for geometric numerical integration. Geometric numerical integration is synonymous with structure-preserving integration of ordinary differential equations. These notes, prepared for the Durham summer school 2002, are complementary to the monograph of E. Hairer, C. Lubich and G. Wanner [Geometric numerical integration. Structure-preserving algorithms for ordinary differential equations. (2002; Zbl 0994.65135)]. They give an introduction to the subject, and they discuss and explain the use of Matlab programs for experimenting with structure-preserving algorithms. We start by presenting some typical classes of problems having properties that are important and should be conserved by the discretization (Section 1). The flow of Hamiltonian differential equations is symplectic and possesses conserved quantities. Conservative systems have a time-reversible flow. Differential equations with first integrals and problems on manifolds are also considered. We then introduce in Section 2 simple symplectic and symmetric integrators, (partitioned) Runge-Kutta methods, composition and splitting methods, linear multistep methods, and algorithms for Hamiltonian problems on manifolds. We briefly discuss their symplecticity and symmetry. The improved performance of such geometric integrators is best understood with the help of a backward error analysis (Section 3). We explain some implications for the long-time integration of Hamiltonian systems and of completely integrable problems. Section 4 is devoted to a presentation and explanation of Matlab codes for implicit Runge-Kutta, composition, and multistep methods. The final Section 5 gives a comparison of the different methods and illustrates the use of these programs at some typical interesting situations: the computation of Poincaré sections, and the simulation of the motion of two bodies on a sphere. The Matlab codes as well as their Fortran 77 counterparts can be downloaded at http://www.unige.ch/math/folks/hairer under the item “software”.
Keywords for this software
References in zbMATH (referenced in 8 articles , 1 standard article )
Showing results 1 to 8 of 8.
- Zhang, Cheng; Huang, Jingfang; Wang, Cheng; Yue, Xingye: On the operator splitting and integral equation preconditioned deferred correction methods for the “good” Boussinesq equation (2018)
- Qu, Wenzhen; Brandon, Namdi; Chen, Dangxing; Huang, Jingfang; Kress, Tyler: A numerical framework for integrating deferred correction methods to solve high order collocation formulations of ODEs (2016)
- Console, Paola; Hairer, Ernst; Lubich, Christian: Symmetric multistep methods for constrained Hamiltonian systems (2013)
- Leung, Shingyu; Qian, Jianliang: Eulerian Gaussian beams for Schrödinger equations in the semi-classical regime (2009)
- Barrio, Roberto: Sensitivity tools vs. Poincaré sections (2005)
- Hairer, Ernst: Important aspects of geometric numerical integration (2005)
- Blowey, James F. (ed.); Craig, Alan W. (ed.); Shardlow, Tony (ed.): Frontiers in numerical analysis. 10th LMS-EPSRC numerical analysis summer school, Durham, UK, July 7--19, 2002 (2003)
- Hairer, Ernst; Hairer, Martin: GniCodes --Matlab programs for geometric numerical integration (2003)