HOD_fast routine now returns the velocities of the galaxies