“…HBVMs are energy-conserving methods derived within the framework of (discrete) line integral methods, initially proposed in [54][55][56][57][58], and later refined in [22][23][24][29][30][31]. The approach has also been extended along several directions [10,14,18,24,25,27,28,39], including Hamiltonian BVPs [1], constrained Hamiltonian problems [15], highly-oscillatory problems [2,21,38], and Hamiltonian PDEs [3,13,16,17,21,40]. We also refer to the review paper [20] and to the monograph [19].…”