## Reaction-diffusion model Monte Carlo simulations on the GPU

Instituut-Lorentz, Leiden University, P.O. Box 9506, 2300 RA Leiden, The Netherlands

arXiv:1301.6082 [physics.comp-ph], (25 Jan 2013)

@article{2013arXiv1301.6082S,

author={Schram}, R.~D.},

title={"{Reaction-diffusion model Monte Carlo simulations on the GPU}"},

journal={ArXiv e-prints},

archivePrefix={"arXiv"},

eprint={1301.6082},

primaryClass={"physics.comp-ph"},

keywords={Physics – Computational Physics, Condensed Matter – Statistical Mechanics},

year={2013},

month={jan},

adsurl={http://adsabs.harvard.edu/abs/2013arXiv1301.6082S},

adsnote={Provided by the SAO/NASA Astrophysics Data System}

}

We created an efficient algorithm suitable for graphics processing units (GPUs) to perform Monte Carlo simulations of a subset of reaction-diffusion models. The algorithm uses techniques that are specific to GPU programming, and combines these with the multispin technique known from CPU programming 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. If we compare the performance of the GPU algorithm, between the GPU and CPU, we find a speed-up of about 130x.

January 28, 2013 by hgpu