Abstract:Langevin Dynamics, Monte Carlo, and all-atom Molecular Dynamics simulations in implicit solvent, widely used to access the microscopic transitions in biomolecules, require a reliable source of random numbers. Here we present the two main approaches for implementation of random number generators (RNGs) on a GPU, which enable one to generate random numbers on the fly. In the one-RNG-per-thread approach, inherent in CPU-based calculations, one RNG produces a stream of random numbers in each thread of execution, whereas the one-RNG-for-all-threads approach builds on the ability of different threads to communicate, thus, sharing random seeds across the entire GPU device. We exemplify the use of these approaches through the development of Ran2, Hybrid Taus, and Lagged Fibonacci algorithms fully implemented on the GPU. As an application-based test of randomness, we carry out LD simulations of N independent harmonic oscillators coupled to a stochastic thermostat. This model allows us to assess statistical quality of random numbers by comparing the simulation output with the exact results that would be obtained with truly random numbers. We also profile the performance of these generators in terms of the computational time, memory usage, and the speedup factor (CPU/GPU time).

Batched Ranged Random Integer Generation

Batched ranged random integer generation

Bias-free True Random-Number Generator.

Algorithms for generating small random samples

A Class of Random Number Generators Based on Weyl Sequence

A fast normal random number generator

Fully scalable randomized benchmarking without motion reversal

Optimal Discrete Uniform Generation from Coin Flips, and Applications

An Alternative to Stride-Based RNG for Monte Carlo Transport

A Bias Free, Quantum Random Number Generator

Randomized Bit Vector

Software/hardware framework for generating parallel Gaussian random numbers based on the Monty Python method

A word-length optimized hardware Gaussian random number generator based on the Box-Muller method

A generalization of the Von Neumann extractor

Efficient pseudo-random number generators for biomolecular simulations on graphics processors

An asymptotically optimal algorithm for generating bin cardinalities

Implementing 64-bit Maximally Equidistributed $\mathbb{F}_2$-Linear Generators with Mersenne Prime Period

Restricted Randomized Benchmarking with Universal Gates of Fixed Sequence Length

Randen - fast backtracking-resistant random generator with AES+Feistel+Reverie

How Homogenize Chaos-Based Pseudo-random Sequences