Random weights, robust lattice rules and the geometry of the cbc$r$c algorithm

TL;DR

This paper develops a robust method for constructing lattice rules with random weights in a reproducing kernel Hilbert space, introducing the cbc$r$c algorithm to ensure good performance across a range of weights, reducing the need for precise weight selection.

Contribution

It proposes the cbc$r$c algorithm that constructs lattice rules effective for all weights within a specified subspace, enhancing robustness and practical applicability.

Findings

Lattice rules generated by cbc$r$c perform well for all weights in the subspace.

The method reduces the need for precise weight determination in practical applications.

Error bounds increase by a factor of r, balancing robustness and computational cost.

Abstract

In this paper we study lattice rules which are cubature formulae to approximate integrands over the unit cube $[0,1]^s$ from a weighted reproducing kernel Hilbert space. We assume that the weights are independent random variables with a given mean and variance for two reasons stemming from practical applications: (i) It is usually not known in practice how to choose the weights. Thus by assuming that the weights are random variables, we obtain robust constructions (with respect to the weights) of lattice rules. This, to some extend, removes the necessity to carefully choose the weights. (ii) In practice it is convenient to use the same lattice rule for many different integrands. The best choice of weights for each integrand may vary to some degree, hence considering the weights random variables does justice to how lattice rules are used in applications. We also study a generalized version which uses $r$ constraints which we call the cbc$r$c (component-by-component with $r$ constraints) algorithm. We show that lattice rules generated by the cbc$r$c algorithm simultaneously work well for all weights in a subspace spanned by the chosen weights $\boldsymbol{\gamma}^{(1)},..., \boldsymbol{\gamma}^{(r)}$. Thus, in applications, instead of finding one set of weights, it is enough to find an $r$ dimensional convex polytope in which the optimal weights lie. The price for this method is a factor $r$ in the upper bound on the error and in the construction cost of the lattice rule. Thus the burden of determining one set of weights very precisely can be shifted to the construction of good lattice rules.