Toward large-scale Hybrid Monte Carlo simulations of the Hubbard model on graphics processing units