Accelerating Hybrid Monte Carlo simulations of the Hubbard model on the hexagonal lattice