Optimal quasi-Monte Carlo rules on order 2 digital nets for the numerical integration of multivariate periodic functions

Abstract

We investigate quasi-Monte Carlo rules for the numerical integration of multivariate periodic functions from Besov spaces $S_{p,q}^{r}B({\mathbb{T}}^d)$ with dominating mixed smoothness $1/p<r<2$. We show that order 2 digital nets achieve the optimal rate of convergence $N^{-r}(\log N{)}^{(d-1)(1-1/q)}$. The logarithmic term does not depend on $r$ and hence improves the known bound provided by J. Dick for the special case of Sobolev spaces $H_{\text{mix}}^r({\mathbb{T}}^d)$. Secondly, the rate of convergence is independent of the integrability $p$ of the Besov space, which allows for sacrificing integrability in order to gain Besov regularity. Our method combines characterizations of periodic Besov spaces with dominating mixed smoothness via Faber bases with sharp estimates of Haar coefficients for the discrepancy function of higher order digital nets. Moreover, we provide numerical computations which indicate that this bound also holds for the case $r=2$.

7 figures in Num. Mathem., online 2015

Cited by in corpus (7)

• The role of Frolov's cubature formula for functions with bounded mixed derivative
• Change of variable in spaces of mixed smoothness and numerical integration of multivariate functions on the unit cube
• Optimal order quasi-Monte Carlo integration in weighted Sobolev spaces of arbitrary smoothness
• An explicit construction of optimal order quasi-Monte Carlo rules for smooth integrands
• Optimal order quadrature error bounds for infinite-dimensional higher order digital sequences
• A higher order Faber spline basis for sampling discretization of functions
• Richardson extrapolation allows truncation of higher order digital nets and sequences