Sobol sequence generator
Project description
Sobol
Concise and fast implementation of the Sobol sequence for generating low-discrepancy quasi-random numbers in up to 1111 dimensions.
Install
pip install sobolsequence
Usage
import sobol
sobol.sample(dimension=3, n_points=5)
array([[0.5 , 0.5 , 0.5 ],
[0.75 , 0.25 , 0.75 ],
[0.25 , 0.75 , 0.25 ],
[0.375, 0.375, 0.625],
[0.875, 0.875, 0.125]])
Skip the first n
points:
sobol.sample(dimension=3, n_points=5, skip=10000)
Sample point by point using the underlying generator:
sob = sobol.generator(dimension=5)
for i in range(10):
print(next(sob))
References
This implementation is based on the Python version by Corrado Chisari available here.
- Antonov, Saleev, USSR Computational Mathematics and Mathematical Physics, Volume 19, 1980, pages 252 - 256.
- Paul Bratley, Bennett Fox, Algorithm 659: Implementing Sobol's Quasirandom Sequence Generator, ACM Transactions on Mathematical Software, Volume 14, Number 1, pages 88-100, 1988.
- Bennett Fox, Algorithm 647: Implementation and Relative Efficiency of Quasirandom Sequence Generators, ACM Transactions on Mathematical Software, Volume 12, Number 4, pages 362-376, 1986.
- Ilya Sobol, USSR Computational Mathematics and Mathematical Physics, Volume 16, pages 236-242, 1977.
- Ilya Sobol, Levitan, The Production of Points Uniformly Distributed in a Multidimensional Cube (in Russian), Preprint IPM Akad. Nauk SSSR, Number 40, Moscow 1976.
The direction numbers from Joe and Kuo are available here.
- Stephen Joe and Frances Kuo, Remark on Algorithm 659: Implementing Sobol's quasirandom sequence generator, ACM Trans. Math. Softw. 29, 49-57 (2003), http://doi.acm.org/10.1145/641876.641879
- Stephen Joe and Frances Kuo, Constructing Sobol sequences with better two-dimensional projections, SIAM J. Sci. Comput. 30, 2635-2654 (2008), https://doi.org/10.1137/070709359
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
SobolSequence-0.2.tar.gz
(23.4 kB
view hashes)
Built Distribution
Close
Hashes for SobolSequence-0.2-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | bfe90b4922335da900a119662c868a4cd0f2bd7944d30f38ae3ea88394e8b464 |
|
MD5 | 7b56618cb4d38845602925a8eea02f62 |
|
BLAKE2b-256 | 1c5498aad3dcf3289a145c09e564a7aa927b6b0d827d06c5a9163abcde099935 |