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.

Using Quasi-Monte Carlo Method to Compute Volume of Point Set

A Quasi-Monte Carlo Method for Computing Areas of Point-Sampled Surfaces

Surface area estimation of digitized 3D objects using quasi-Monte Carlo methods

An Efficient Method to Mesh Point Cloud

Incremental Surface Approximation from Point Cloud

Volume Dispersion of Point Sets and Quasi-Monte Carlo Methods

Application of Quasi-Monte Carlo Method Based on Good Point Set in Tolerance Analysis

Blue Noise Sampling of Surfaces.

A Fast Progressive Surface Reconstruction Algorithm for Unorganized Points

3D Surface Segmentation from Point Clouds Via Quadric Fits Based on DBSCAN Clustering

The Use of Quasi-Symmetric Point Method in Probability Density Evolution Method

A High-Precision Reconstruction Algorithm of 3D Data Fields Based on Quasi-interpolation

Quasi-Monte Carlo Beyond Hardy-Krause

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

Mesh Reconstruction of Noisy Point Cloud Based on Intrinsic Property

A Point Cloud Segmentation Algorithm Using Local Convexity and Octree

Compute Quasi-Conformal Map on Point Cloud by Point Integral Method

Incremental Reconstruction of Water-Tight Surface from Unorganized Points

Robust Surface Area Measurement of Unorganized Point Clouds Based on Multiscale Supervoxel Segmentation

Novel Iso-surface Reconstruction Algorithm for Numerical Analysis from Medical Volume Data