Large-scale Monte Carlo simulation of two-dimensional classical XY model using multiple GPUs