Arbitrarily High-Order Energy-Conserving Methods for Hamiltonian Problems with Quadratic Holonomic Constraints

In this paper, we define arbitrarily high-order energy-conserving methods for Hamiltonian systems with quadratic holonomic constraints. The derivation of the methods is made within the so-called line integral framework. Numerical tests to illustrate the theoretical findings are presented.