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.

Random Leja points

Adaptive Sparse Polynomial Chaos Expansions via Leja Interpolation

On the approximation properties of fast Leja points

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

Random points are good for universal discretization

Randomized least-squares with minimal oversampling and interpolation in general spaces

Quasi-Monte Carlo Beyond Hardy-Krause

Integer Feasibility of Random Polytopes

Proximal random reshuffling under local Lipschitz continuity

On Randomized Fictitious Play for Approximating Saddle Points Over Convex Sets

An explicit construction of optimized interpolation points on the 4-simplex

Random interpolation in the Nevanlinna and Smirnov classes and related spaces

Polynomial approximation and cubature at approximate Fekete and Leja points of the cylinder

Cubic Hermite interpolators on the space of probability measures

Bivariate Lagrange interpolation at the Padua points: the ideal theory approach

Random triangles in planar regions containing a fixed point

Randomized Iterative Methods for Linear Systems

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

Integral points on cubic surfaces: heuristics and numerics

Bravais Lattices for Euclidean Degree Efficient Polynomial Interpolation