QCD simulations with staggered fermions on GPUs