Efficient implementation of Gauss collocation and Hamiltonian Boundary Value Methods

Abstract

In this paper we define an efficient implementation for the family of low-rank energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs), recently defined in the last years. The proposed implementation relies on the particular structure of the Butcher matrix defining such methods, for which we can derive an efficient splitting procedure. The very same procedure turns out to be automatically suited for the efficient implementation of Gauss-Legendre collocation methods, since these methods are a special instance of HBVMs. The linear convergence analysis of the splitting procedure exhibits excellent properties, which are confirmed by a few numerical tests.

19 pages, 7 tables, 2 figures

Cited by in corpus (14)

• Energy conservation issues in the numerical solution of the semilinear wave equation
• On the effectiveness of spectral methods for the numerical solution of multi-frequency highly-oscillatory Hamiltonian problems
• High-order energy-conserving Line Integral Methods for charged particle dynamics
• Analysis of Energy and QUadratic Invariant Preserving (EQUIP) methods
• Arbitrarily high-order energy-preserving methods for simulating the gyrocenter dynamics of charged particles
• Analysis of Spectral Hamiltonian Boundary Value Methods (SHBVMs) for the numerical solution of ODE problems
• Line Integral solution of Hamiltonian PDEs
• A new framework for polynomial approximation to differential equations
• A note on the continuous-stage Runge-Kutta-(Nyström) formulation of Hamiltonian Boundary Value Methods (HBVMs)
• Arbitrary high-order linearly implicit energy-preserving algorithms for Hamiltonian PDEs
• Arbitrarily high-order energy-conserving methods for Poisson problems
• Line Integral Solution of Hamiltonian Systems with Holonomic Constraints
• Spectrally accurate space-time solution of Manakov systems
• Arbitrary high-order methods for one-sided direct event location in discontinuous differential problems with nonlinear event function