Multi-GPU Accelerated Multi-Spin Monte Carlo Simulations of the 2D Ising Model

Abstract

A modern graphics processing unit (GPU) is able to perform massively parallel scientific computations at low cost. We extend our implementation of the checkerboard algorithm for the two dimensional Ising model [T. Preis et al., J. Comp. Phys. 228, 4468 (2009)] in order to overcome the memory limitations of a single GPU which enables us to simulate significantly larger systems. Using multi-spin coding techniques, we are able to accelerate simulations on a single GPU by factors up to 35 compared to an optimized single Central Processor Unit (CPU) core implementation which employs multi-spin coding. By combining the Compute Unified Device Architecture (CUDA) with the Message Parsing Interface (MPI) on the CPU level, a single Ising lattice can be updated by a cluster of GPUs in parallel. For large systems, the computation time scales nearly linearly with the number of GPUs used. As proof of concept we reproduce the critical temperature of the 2D Ising model using finite size scaling techniques.

References in corpus (8)

• Implementation and evaluation of various demons deformable image registration algorithms on GPU
• Harvesting graphics power for MD simulations
• GPU-based ultra fast IMRT plan optimization
• GPU-based ultra fast dose calculation using a finite pencil beam model
• Air pollution modelling using a graphics processing unit with CUDA
• Phase diagram and critical behavior of the square-lattice Ising model with competing nearest- and next-nearest-neighbor interactions
• Parallel Algorithm for Solving Kepler's Equation on Graphics Processing Units: Application to Analysis of Doppler Exoplanet Searches
• Implementation and performance analysis of bridging Monte Carlo moves for off-lattice single chain polymers in globular states

Cited by in corpus (46)

• Massively Parallel Probabilistic Computing with Sparse Ising Machines
• Optimized simulated annealing for Ising spin glasses
• Speedup for quantum optimal control from automatic differentiation based on graphics processing units
• Highly accelerated simulations of glassy dynamics using GPUs: caveats on limited floating-point precision
• Roadmap for Unconventional Computing with Nanotechnology
• Large-scale Monte Carlo simulation of two-dimensional classical XY model using multiple GPUs
• Performance potential for simulating spin models on GPU
• Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units
• Pseudo-random number generators for Monte Carlo simulations on Graphics Processing Units
• GPU-based Swendsen-Wang multi-cluster algorithm for the simulation of two-dimensional classical spin systems
• Simulation of reaction-diffusion processes in three dimensions using CUDA
• Accelerated Quantum Monte Carlo with Probabilistic Computers
• q-State Potts model metastability study using optimized GPU-based Monte Carlo algorithms
• A Performance Study of the 2D Ising Model on GPUs
• Massively parallel multicanonical simulations
• Parallel Tempering Simulation of the three-dimensional Edwards-Anderson Model with Compact Asynchronous Multispin Coding on GPU
• Simulating Lattice Spin Models on Graphics Processing Units
• GPU-based single-cluster algorithm for the simulation of the Ising model
• Highly optimized simulations on single- and multi-GPU systems of 3D Ising spin glass
• CUDA programs for GPU computing of Swendsen-Wang multi-cluster spin flip algorithm: 2D and 3D Ising, Potts, and XY models
• Ageing at the Spin-Glass/Ferromagnet Transition: Monte Carlo Simulation using GPUs
• A GPU-based Large-scale Monte Carlo Simulation Method for Systems with Long-range Interactions
• Adaptive Multi-GPU Exchange Monte Carlo for the 3D Random Field Ising Model
• Monte Carlo simulation of the effects of higher order anisotropy on the spin reorientation transition in the two dimensional Heisenberg model with long range interactions
• Computation of magnetization, exchange stiffness, anisotropy, and susceptibilities in large-scale systems using GPU-accelerated atomistic parallel Monte Carlo algorithms
• Dynamic scaling in the 2D Ising spin glass with Gaussian couplings
• Multi-GPU-based Swendsen-Wang multi-cluster algorithm for the simulation of two-dimensional q-state Potts model
• Anisotropic interfacial tension, contact angles, and line tensions: A graphics-processing-unit-based Monte Carlo study of the Ising model
• GPU Based Parallel Ising Computing for Combinatorial Optimization Problems in VLSI Physical Design
• Pushing the Boundary of Quantum Advantage in Hard Combinatorial Optimization with Probabilistic Computers
• High Performance Monte Carlo Simulation of Ising Model on TPU Clusters
• Accelerating Binarized Neural Networks via Bit-Tensor-Cores in Turing GPUs
• Multipath Metropolis Simulation of Classical Heisenberg Model
• The Effect of Defects on Magnetic Droplet Nucleation
• Superdiffusion-like behavior in zero-temperature coarsening of the $d=3$ Ising model
• Fast GPU-based calculations in few-body quantum scattering
• Reaction-diffusion model Monte Carlo simulations on the GPU
• GPU acceleration of ab initio simulations of large-scale identical particles based on path integral molecular dynamics
• A differentiable programming framework for spin models
• Parallelization and implementation of multi-spin Monte Carlo simulation of 2D square Ising model using MPI and C++
• GPU Computing in Bayesian Inference of Realized Stochastic Volatility Model
• Fast algorithm for generating random bit strings and multispin coding for directed percolation
• Massive-Scale Simulations of 2D Ising and Blume-Capel Models on Rack-Scale Multi-GPU Systems
• Efficient sampling using Macrocanonical Monte Carlo and density of states mapping
• Metastability in the diluted parallel Ising model
• GPU-based Ising Computing for Solving Balanced Min-Cut Graph Partitioning Problem