GPU-based single-cluster algorithm for the simulation of the Ising model

  • Authors:
  • Yukihiro Komura;Yutaka Okabe

  • Affiliations:
  • Department of Physics, Tokyo Metropolitan University, Hachioji, Tokyo 192-0397, Japan;Department of Physics, Tokyo Metropolitan University, Hachioji, Tokyo 192-0397, Japan

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2012

Quantified Score

Hi-index 31.45

Visualization

Abstract

We present the GPU calculation with the common unified device architecture (CUDA) for the Wolff single-cluster algorithm of the Ising model. Proposing an algorithm for a quasi-block synchronization, we realize the Wolff single-cluster Monte Carlo simulation with CUDA. We perform parallel computations for the newly added spins in the growing cluster. As a result, the GPU calculation speed for the two-dimensional Ising model at the critical temperature with the linear size L=4096 is 5.60 times as fast as the calculation speed on a current CPU core. For the three-dimensional Ising model with the linear size L=256, the GPU calculation speed is 7.90 times as fast as the CPU calculation speed. The idea of quasi-block synchronization can be used not only in the cluster algorithm but also in many fields where the synchronization of all threads is required.