Tables of maximally equidistributed combined LFSR generators
Mathematics of Computation
General purpose molecular dynamics simulations fully implemented on graphics processing units
Journal of Computational Physics
GPU accelerated Monte Carlo simulation of the 2D and 3D Ising model
Journal of Computational Physics
GPU implementations of the bond fluctuation model
Journal of Computational Physics
| Hi-index | 31.45 |
We created an efficient algorithm suitable for graphics processing units (GPUs) to perform Monte Carlo simulations on a subset of reaction-diffusion models. The set of reaction-diffusion models that the algorithm is applied to represents a seemingly simplistic set of problems on a one-dimensional lattice, where each site contains either a particle or is empty. However, these systems exhibit non-equilibrium phase transitions, with very large finite-time corrections, which mandates a fast algorithm to simulate them. The algorithm presented here uses techniques that are specific to GPU programming, and combines these with multispin coding to create one of the fastest algorithms for reaction-diffusion models. As an example, the algorithm is applied to the pair contact process with diffusion (PCPD). Compared to a simple algorithm on the CPU, our GPU algorithm is approximately 4000 times faster. The GPU algorithm is roughly 55 times faster than an optimized version for the CPU.