Improved bounds on the gain coefficients for digital nets in prime base

07/11/2022
by   Takashi Goda, et al.
0

We study randomized quasi-Monte Carlo integration by scrambled nets. The scrambled net quadrature has long gained its popularity because it is an unbiased estimator of the true integral, allows for a practical error estimation, achieves a high order decay of the variance for smooth functions, and works even for L^p-functions with any p>1. The variance of the scrambled net quadrature for L^2-functions can be evaluated through the set of the so-called gain coefficients. In this paper, based on the system of Walsh functions and the concept of dual nets, we provide improved upper bounds on the gain coefficients for digital nets in general prime base. Our results explain the known bound by Owen (1997) for Faure sequences, the recently improved bound by Pan and Owen (2021) for digital nets in base 2 (including Sobol' sequences as a special case), and their finding that all the nonzero gain coefficients for digital nets in base 2 must be powers of two, all in a unified way.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
06/19/2021

The nonzero gain coefficients of Sobol's sequences are always powers of two

When a plain Monte Carlo estimate on n samples has variance σ^2/n, then ...
research
06/11/2020

Walsh functions, scrambled (0,m,s)-nets, and negative covariance: applying symbolic computation to quasi-Monte Carlo integration

We investigate base b Walsh functions for which the variance of the inte...
research
06/12/2023

Analysis and Synthesis of Digital Dyadic Sequences

We explore the space of matrix-generated (0, m, 2)-nets and (0, 2)-seque...
research
08/15/2023

Gain coefficients for scrambled Halton points

Randomized quasi-Monte Carlo, via certain scramblings of digital nets, p...
research
08/09/2022

Super-polynomial accuracy of multidimensional randomized nets using the median-of-means

We study approximate integration of a function f over [0,1]^s based on t...
research
09/08/2023

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

Let f:[0,1]^d→ℝ be a completely monotone integrand as defined by Aistlei...
research
10/05/2019

An algorithm to compute the t-value of a digital net and of its projections

Digital nets are among the most successful methods to construct low-disc...

Please sign up or login with your details

Forgot password? Click here to reset