807 research outputs found
A universal median quasi-Monte Carlo integration
We study quasi-Monte Carlo (QMC) integration over the multi-dimensional unit
cube in several weighted function spaces with different smoothness classes. We
consider approximating the integral of a function by the median of several
integral estimates under independent and random choices of the underlying QMC
point sets (either linearly scrambled digital nets or infinite-precision
polynomial lattice point sets). Even though our approach does not require any
information on the smoothness and weights of a target function space as an
input, we can prove a probabilistic upper bound on the worst-case error for the
respective weighted function space, where the failure probability converges to
0 exponentially fast as the number of estimates increases. Our obtained rates
of convergence are nearly optimal for function spaces with finite smoothness,
and we can attain a dimension-independent super-polynomial convergence for a
class of infinitely differentiable functions. This implies that our
median-based QMC rule is universal in the sense that it does not need to be
adjusted to the smoothness and the weights of the function spaces and yet
exhibits the nearly optimal rate of convergence. Numerical experiments support
our theoretical results.Comment: Major revision, 32 pages, 4 figure
Component-by-component construction of randomized rank-1 lattice rules achieving almost the optimal randomized error rate
We study a randomized quadrature algorithm to approximate the integral of
periodic functions defined over the high-dimensional unit cube. Recent work by
Kritzer, Kuo, Nuyens and Ullrich (2019) shows that rank-1 lattice rules with a
randomly chosen number of points and good generating vector achieve almost the
optimal order of the randomized error in weighted Korobov spaces, and moreover,
that the error is bounded independently of the dimension if the weight
parameters, , satisfy the summability condition
, where is a smoothness
parameter. The argument is based on the existence result that at least half of
the possible generating vectors yield almost the optimal order of the
worst-case error in the same function spaces.
In this paper we provide a component-by-component construction algorithm of
such randomized rank-1 lattice rules, without any need to check whether the
constructed generating vectors satisfy a desired worst-case error bound.
Similarly to the above-mentioned work, we prove that our algorithm achieves
almost the optimal order of the randomized error and that the error bound is
independent of the dimension if the same condition
holds. We also provide
analogous results for tent-transformed lattice rules for weighted half-period
cosine spaces and for polynomial lattice rules in weighted Walsh spaces,
respectively.Comment: major revision, 29 pages, 3 figure
- β¦