Search code examples
pythonnumpyrandomnumbersscipy

Generating low discrepancy quasi-random sequences in python/numpy/scipy?

There is already a question on this but the answer contains a broken link, and being over two years old, I'm hoping there's a better solution now :)

Low discrepancy quasi-random sequences, e.g. Sobol sequences, fill a space more uniformly than uniformly random sequences. Is there a good/easy way to generate them in python?

Solution

  • I think the best alternative for Low Discrepancy sequences in Python is Sensitivity Analysis Library (SALib):

    https://github.com/SALib/SALib

    I think this is an active project and you can contact the author to check if the functionalities you need are already implemented. If that doesn't solve your problem, Corrado Chisari ported a SOBOL version made in Matlab (by John Burkardt) to Python, you can access it here:

    http://people.sc.fsu.edu/~jburkardt/py_src/sobol/sobol.html

    Someone cleaned up the comments in these sources and put them in the format of docstrings. It's much more readable and you can access it here:

    https://github.com/naught101/sobol_seq