English

Random-prime--fixed-vector randomised lattice-based algorithm for high-dimensional integration

Numerical Analysis 2023-04-21 v1 Numerical Analysis

Abstract

We show that a very simple randomised algorithm for numerical integration can produce a near optimal rate of convergence for integrals of functions in the dd-dimensional weighted Korobov space. This algorithm uses a lattice rule with a fixed generating vector and the only random element is the choice of the number of function evaluations. For a given computational budget nn of a maximum allowed number of function evaluations, we uniformly pick a prime pp in the range n/2<pnn/2 < p \le n. We show error bounds for the randomised error, which is defined as the worst case expected error, of the form O(nα1/2+δ)O(n^{-\alpha - 1/2 + \delta}), with δ>0\delta > 0, for a Korobov space with smoothness α>1/2\alpha > 1/2 and general weights. The implied constant in the bound is dimension-independent given the usual conditions on the weights. We present an algorithm that can construct suitable generating vectors \emph{offline} ahead of time at cost O(dn4/lnn)O(d n^4 / \ln n) when the weight parameters defining the Korobov spaces are so-called product weights. For this case, numerical experiments confirm our theory that the new randomised algorithm achieves the near optimal rate of the randomised error.

Keywords

Cite

@article{arxiv.2304.10413,
  title  = {Random-prime--fixed-vector randomised lattice-based algorithm for high-dimensional integration},
  author = {Frances Y. Kuo and Dirk Nuyens and Laurence Wilkes},
  journal= {arXiv preprint arXiv:2304.10413},
  year   = {2023}
}

Comments

29 pages, 2 figures