Scaled lattice rules for integration on ℝ^d achieving higher-order convergence with error analysis in terms of orthogonal projections onto periodic spaces

by   Dirk Nuyens, et al.

We show that to approximate the integral ∫_ℝ^d f(x) dx one can simply use scaled lattice rules from the unit cube [0,1]^d to properly sized boxes on ℝ^d, achieving higher-order convergence that matches the smoothness of the integrand function f in a certain Sobolev space of dominating mixed smoothness. Our method only assumes that we can evaluate the integrand function f and does not assume a particular density nor the ability to sample from it. In particular, for the analysis we show that the method of adding Bernoulli polynomials to a function to make it “periodic” on a box without changing its integral value over the box, is equivalent to an orthogonal projection from a well chosen Sobolev space of dominating mixed smoothness to an associated periodic Sobolev space of the same dominating mixed smoothness, which we call a Korobov space. We note that the Bernoulli polynomial method is often not used because of its computational complexity and also here we completely avoid applying it. Instead, we use it as a theoretical tool in the error analysis of applying a scaled lattice rule to increasing boxes in order to approximate integrals over the d-dimensional Euclidean space. Such a method would not work on the unit cube since then the committed error caused by non-periodicity of the integrand would be constant, but for integration on the Euclidean space we can use the certain decay towards zero when the boxes grow. Hence we can bound the truncation error as well as the projection error and show higher-order convergence in applying scaled lattice rules for integration on Euclidean space. We illustrate our theoretical analysis by numerical experiments which confirm our findings.


page 1

page 2

page 3

page 4


Construction-free median quasi-Monte Carlo rules for function spaces with unspecified smoothness and general weights

We study quasi-Monte Carlo (QMC) integration of smooth functions defined...

Digit-by-digit and component-by-component constructions of lattice rules for periodic functions with unknown smoothness

Lattice rules are among the most prominently studied quasi-Monte Carlo m...

A universal median quasi-Monte Carlo integration

We study quasi-Monte Carlo (QMC) integration over the multi-dimensional ...

Equivalence between Sobolev spaces of first-order dominating mixed smoothness and unanchored ANOVA spaces on ℝ^d

We prove that a variant of the classical Sobolev space of first-order do...

Weighted integration over a cube based on digital nets and sequences

Quasi-Monte Carlo (QMC) methods are equal weight quadrature rules to app...

Needlets Liberated

Spherical needlets were introduced by Narcowich, Petrushev, and Ward to ...