“…In recent years, contact Hamiltonian systems have found many applications, first in the context of thermodynamics [9][10][11] and, more recently, in the context of the Hamiltonisation of several dissipative dynamical systems [12][13][14][15][16][17][18][19]. The large number of applications of contact systems that have appeared recently motivated research on geometric numerical integration [15,16,20,21]. Fortunately, contact flows possess geometric integrators (both variational and Hamiltonian) that precisely parallel their symplectic counterparts, and therefore they show remarkable numerical and analytical properties such as, e.g., increased stability, near-preservation of invariant quantities, and modified Hamiltonians.…”