Reaction–diffusion model Monte Carlo simulations on the GPU
| Autor: | Raoul D. Schram |
|---|---|
| Rok vydání: | 2013 |
| Předmět: |
Numerical Analysis
Statistical Mechanics (cond-mat.stat-mech) Physics and Astronomy (miscellaneous) Computer science Applied Mathematics Monte Carlo method FOS: Physical sciences Computational Physics (physics.comp-ph) Computer Science Applications Computational science Hybrid Monte Carlo Computational Mathematics Computer Science::Graphics Modeling and Simulation Computer Science::Mathematical Software Dynamic Monte Carlo method Central processing unit Kinetic Monte Carlo General-purpose computing on graphics processing units Graphics Physics - Computational Physics Condensed Matter - Statistical Mechanics SIMPLE algorithm |
| Zdroj: | Journal of Computational Physics. 241:95-103 |
| ISSN: | 0021-9991 |
| DOI: | 10.1016/j.jcp.2013.01.041 |
| Popis: | 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. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|
Full Text from ScienceDirect