Large scale numerical simulation of the three-state Potts model