English

On the Randomization of Frolov's Algorithm for Multivariate Integration

Numerical Analysis 2016-03-17 v2

Abstract

We are concerned with the numerical integration of functions from the Sobolev space Hr,mix([0,1]d)H^{r,\text{mix}}([0,1]^d) of dominating mixed smoothness rNr\in\mathbb{N} over the dd-dimensional unit cube. In 1976, K. K. Frolov introduced a deterministic quadrature rule whose worst case error has the order nr(logn)(d1)/2n^{-r} \, (\log n)^{(d-1)/2} with respect to the number nn of function evaluations. This is known to be optimal. 39 years later, Erich Novak and me introduced a randomized version of this algorithm using dd random dilations. We showed that its error is bounded above by a constant multiple of nr1/2(logn)(d1)/2n^{-r-1/2} \, (\log n)^{(d-1)/2} in expectation and by nr(logn)(d1)/2n^{-r} \, (\log n)^{(d-1)/2} almost surely. The main term nr1/2n^{-r-1/2} is again optimal and it turns out that the very same algorithm is also optimal for the isotropic Sobolev space Hs([0,1]d)H^s([0,1]^d) of smoothness s>d/2s>d/2. We also added a random shift to this algorithm to make it unbiased. Just recently, Mario Ullrich proved that the expected error of the resulting algorithm on Hr,mix([0,1]d)H^{r,\text{mix}}([0,1]^d) is even bounded above by nr1/2n^{-r-1/2}. This thesis is a review of the mentioned upper bounds and their proofs.

Cite

@article{arxiv.1603.04637,
  title  = {On the Randomization of Frolov's Algorithm for Multivariate Integration},
  author = {David Krieg},
  journal= {arXiv preprint arXiv:1603.04637},
  year   = {2016}
}

Comments

Master Thesis

R2 v1 2026-06-22T13:11:10.558Z