GPU implementation of a hybrid lattice Boltzmann method for non-isothermal flows