GPU-based Swendsen-Wang multi-cluster algorithm for the simulation of two-dimensional classical spin systems

Abstract

We present the GPU calculation with the common unified device architecture (CUDA) for the Swendsen-Wang multi-cluster algorithm of two-dimensional classical spin systems. We adjust the two connected component labeling algorithms recently proposed with CUDA for the assignment of the cluster in the Swendsen-Wang algorithm. Starting with the q-state Potts model, we extend our implementation to the system of vector spins, the q-state clock model, with the idea of embedded cluster. We test the performance, and the calculation time on GTX580 is obtained as 2.51 nano sec per a spin flip for the q=2 Potts model (Ising model) and 2.42 nano sec per a spin flip for the q=6 clock model with the linear size L=4096 at the critical temperature, respectively. The computational speed for the q=2 Potts model on GTX580 is 12.4 times as fast as the calculation speed on a current CPU core. That for the q=6 clock model on GTX580 is 35.6 times as fast as the calculation speed on a current CPU core.

accepted for publication in Comp. Phys. Commun

References in corpus (8)

• Multi-GPU Accelerated Multi-Spin Monte Carlo Simulations of the 2D Ising Model
• Probability-Changing Cluster Algorithm for Two-Dimensional XY and Clock Models
• Simulating spin models on GPU
• Finite-size Scaling of Correlation Ratio and Generalized Scheme for the Probability-Changing Cluster Algorithm
• Connected component identification and cluster update on GPU
• GPU-based single-cluster algorithm for the simulation of the Ising model
• Simulating Lattice Spin Models on Graphics Processing Units
• Comment on `Six-state clock model on the square lattice: Fisher zero approach with Wang-Landau sampling'

Cited by in corpus (19)

• Large-scale Monte Carlo simulation of two-dimensional classical XY model using multiple GPUs
• Berezinskii-Kosterlitz-Thouless transition on regular and Villain types of $q$-state clock models
• Entanglement entropy from non-equilibrium Monte Carlo simulations
• CUDA programs for GPU computing of Swendsen-Wang multi-cluster spin flip algorithm: 2D and 3D Ising, Potts, and XY models
• Parallelized event chain algorithm for dense hard sphere and polymer systems
• Adaptive Multi-GPU Exchange Monte Carlo for the 3D Random Field Ising Model
• TLBfind: a Thermal Lattice Boltzmann code for concentrated emulsions with FINite-size Droplets
• GPU Based Detection of Topological Changes in Voronoi Diagrams
• Duality transformations and the entanglement entropy of gauge theories
• Multi-GPU-based Swendsen-Wang multi-cluster algorithm for the simulation of two-dimensional q-state Potts model
• Large-scale calculation of ferromagnetic spin systems on the pyrochlore lattice
• Critical Nonequilibrium Cluster-flip Relaxations in Ising Models
• Measurement of entanglement entropy in the two-dimensional Potts model using wavelet analysis
• Parallelization and implementation of multi-spin Monte Carlo simulation of 2D square Ising model using MPI and C++
• The Nonclassical Regime of the Two-dimensional Long-range XY Model: a Comprehensive Monte Carlo Study
• Parallelized Stochastic Cutoff Method for Long-Range Interacting Systems
• Fast algorithm for generating random bit strings and multispin coding for directed percolation
• A generalized GPU-based connected component labeling algorithm
• Casimir effect in critical $\mathrm{O}(N)$ models from non-equilibrium Monte Carlo simulations