Co-simulation requires the ability to represent systems in a modular form, while guaranteeing that simulation algorithms rely exclusively on model external interface. These conditions enable models to be composed without exposing their internal state, a requisite needed for the co-simulation of cyber-physical systems. In this paper we provide a modular representation of geometric integrators, a type of integrator essential to the simulation of 2nd-order energy preserving systems. These integrators offer an alternative to the conventional decomposition of systems into 1st-order Ordinary Differential Equations (ODEs). This latter approach, although commonly used in nowadays M&S software, is not acceptable when long simulation runs are needed. Geometric integrators are represented in the Hybrid Systems Specification (HYFLOW), a modeling formalism to represent hybrid modular dynamic topology systems. We show that HYFLOW enables the composition of geometrical solvers, allowing the co-simulation of complex 2nd-order energy preserving systems.