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.

Computable error bounds for quasi-Monte Carlo using points with non-negative local discrepancy

Digital nets with infinite digit expansions and construction of folded digital nets for quasi-Monte Carlo integration

Discrepancy bounds for low-dimensional point sets

Error Bounds for Quasi-Monte Carlo Integration for $${\fancyscript{l}_{\infty}}$$ with Uniform Point Sets

Error Bounds for Quasi-Monte Carlo Integration for \Mathscr{l}^{\infty}

Error bounds for quasi-Monte Carlo integration for ℒ^∞ with uniform point sets

Error bounds for quasi-Monte Carlo integration for L_∞ with uniform point sets

Error Bounds For Quasi-Monte Carlo Integration For L-Infinity With Uniform Point Sets

Low discrepancy sequences in high dimensions: How well are their projections distributed?

Quasi-Monte Carlo Beyond Hardy-Krause

Error bounds for quasi-Monte Carlo integration for <img src="/fulltext-image.asp?format=htmlnonpaginated&src=EW2351N2N3457633_html\605_2010_273_Article_IEq1.gif" border="0" alt="$${\fancyscript{L}_{\infty}}$$" /> with uniform point sets

Fast construction of higher order digital nets for numerical integration in weighted Sobolev spaces

The $b$-adic tent transformation for quasi-Monte Carlo integration using digital nets

Convergence of Uniformity Criteria and the Application in Numerical Integration

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

Achieving High Convergence Rates by Quasi-Monte Carlo and Importance Sampling for Unbounded Integrands

The $b$-adic symmetrization of digital nets for quasi-Monte Carlo integration

A universal median quasi-Monte Carlo integration

Construction of interlaced scrambled polynomial lattice rules of arbitrary high order

An efficient implementation algorithm for quasi-Monte Carlo approximations of high-dimensional integrals

On the upper bound of the $L_p$ discrepancy of Halton's sequence and the Central Limit Theorem for Hammersley's net