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.

A constructive approach to strong tractability using Quasi-Monte Carlo algorithms

Strong tractability of multivariate integration using quasi-Monte Carlo algorithms

Liberating The Weights

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

Strong tractability for multivariate integration in a subspace of the Wiener algebra

Linear Monte Carlo quadrature with optimal confidence intervals

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

Tractability of non-homogeneous tensor product problems in the worst case setting

Strong equivalences of approximation numbers and tractability of weighted anisotropic Sobolev embeddings

A Random Integration Algorithm for High-dimensional Function Spaces

Quasi-Monte Carlo Beyond Hardy-Krause

A new approach to strong convergence

Strong convergence of an explicit numerical approximation for n -dimensional superlinear SDEs with positive solutions

A universal median quasi-Monte Carlo integration

Quasi-Monte Carlo for unbounded integrands with importance sampling

On quantitative stability in infinite-dimensional optimization under uncertainty

Worst case tractability of $L_2$-approximation for weighted Korobov spaces

Good Lattice Rules in Weighted Korobov Spaces with General Weights

A rigorous integrator and global existence for higher-dimensional semilinear parabolic PDEs via semigroup theory

Strong convergence theorems by an extragradient-like approximation methods for quasi-variational inequalities