Abstract:The classical approaches to numerically integrating a function $f$ are Monte Carlo (MC) and quasi-Monte Carlo (QMC) methods. MC methods use random samples to evaluate $f$ and have error $O(\sigma(f)/\sqrt{n})$, where $\sigma(f)$ is the standard deviation of $f$. QMC methods are based on evaluating $f$ at explicit point sets with low discrepancy, and as given by the classical Koksma-Hlawka inequality, they have error $\widetilde{O}(\sigma_{\mathsf{HK}}(f)/n)$, where $\sigma_{\mathsf{HK}}(f)$ is the variation of $f$ in the sense of Hardy and Krause. These two methods have distinctive advantages and shortcomings, and a fundamental question is to find a method that combines the advantages of both. In this work, we give a simple randomized algorithm that produces QMC point sets with the following desirable features: (1) It achieves substantially better error than given by the classical Koksma-Hlawka inequality. In particular, it has error $\widetilde{O}(\sigma_{\mathsf{SO}}(f)/n)$, where $\sigma_{\mathsf{SO}}(f)$ is a new measure of variation that we introduce, which is substantially smaller than the Hardy-Krause variation. (2) The algorithm only requires random samples from the underlying distribution, which makes it as flexible as MC. (3) It automatically achieves the best of both MC and QMC (and the above improvement over Hardy-Krause variation) in an optimal way. (4) The algorithm is extremely efficient, with an amortized $\widetilde{O}(1)$ runtime per sample. Our method is based on the classical transference principle in geometric discrepancy, combined with recent algorithmic innovations in combinatorial discrepancy that besides producing low-discrepancy colorings, also guarantee certain subgaussian properties. This allows us to bypass several limitations of previous works in bridging the gap between MC and QMC methods and go beyond the Hardy-Krause variation.

Randomized Halton Sequences

An improved Halton sequence for implementation in quasi-Monte Carlo methods

Scrambled Halton Subsequences and Inverse Star-Discrepancy

Randomization Using Quasigroups, Hadamard and Number Theoretic Transforms

Randomized Quasi-Monte Carlo Methods on Triangles: Extensible Lattices and Sequences

A Class of Orthogonal Sequences

Gain Coefficients for Scrambled Halton Points

Digital Nets and Sequences for Quasi-Monte Carlo Methods

Towards Derandomising Markov Chain Monte Carlo

Sampling with Hammersley and Halton Points

The importance of being scrambled: supercharged Quasi Monte Carlo

Low discrepancy sequences in high dimensions

On-Line Selection of Alternating Subsequences from a Random Sample

Quasi-Monte Carlo Beyond Hardy-Krause

Parallelized Midpoint Randomization for Langevin Monte Carlo

Sufficient Conditions for Central Limit Theorems and Confidence Intervals for Randomized Quasi-Monte Carlo Methods

A Universal Algorithm for Transforming Chaotic Sequences into Uniform Pseudo-Random Sequences

Outperforming the Best 1D Low-Discrepancy Constructions with a Greedy Algorithm

Randomness of D Sequences via Diehard Testing

Hyperchaotic bilateral random low-rank approximation random sequence generation method and its application on compressive ghost imaging