Classical simulations of Gaussian boson sampling