We show the use of the homogenous architecture in the parallel processing of long range interactions. We describe the implementation of a Monte Carlo algorithm for a two-dimensional Coulomb system on a parallel processor with hypercubic geometry (the 8-node concurrent processor at Caltech). Our res
Quasirandom Number Generators for Parallel Monte Carlo Algorithms
β Scribed by B.C. Bromley
- Publisher
- Elsevier Science
- Year
- 1996
- Tongue
- English
- Weight
- 191 KB
- Volume
- 38
- Category
- Article
- ISSN
- 0743-7315
No coin nor oath required. For personal study only.
β¦ Synopsis
A method for generating sequences of quasirandom numbers allows conventional serial Monte Carlo algorithms to be parallelized using a leapfrog scheme. Specifically, a Sobol' sequence can be broken up into interleaved subsets; with each processing node calculating a unique subset of the full sequence, all of the computational advantages of quasirandom Monte Carlo methods over pseudorandom algorithms or grid-based techniques are retained. Tests with several parallel supercomputers demonstrate that as many as 10 6 integration points (up to 6 dimensions) can be generated per second per node in the optimal case where the number of nodes is a power of 2. The speed of communication-free parallel Sobol' sequence generators and the rapid convergence properties of quasirandom Monte Carlo schemes indicate that the method described here may be gainfully applied to a wide range of problems.
π SIMILAR VOLUMES
## Abstract We have tested and compared several (pseudo) random number generators (RNGs) applied to a practical application, ground state energy calculations of molecules using variational and diffusion Monte Carlo metheds. A new multiple recursive generator with 8thβorder recursion (MRG8) and the
We describe the development of Metropolis Monte Carlo algorithms for a general network of multiple instruction multiple data (MIMD) parallel processors. The implementation of farm, event, and systolic parallel algorithms on transputer-based computers is detailed and their relative performance discus
I give a pedagogical introduction to the generalised Hybrid Monte Carlo and related algorithms. I shall explain why they work, how their performance depends upon the number of degrees of freedom and the correlation length, and how they can be tuned to reduce critical slowing down.