Shifted lattice rules based on a general weighted discrepancy for integrals over Euclidean space

  • Authors:
  • Vasile Sinescu

  • Affiliations:
  • Département d'Informatique et de Recherche Opérationnelle, Université de Montréal, C.P. 6128 Succ. Centre-Ville, Montréal QC H3C 3J7, Canada

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2009

Quantified Score

Hi-index 7.30

Visualization

Abstract

We approximate weighted integrals over Euclidean space by using shifted rank-1 lattice rules with good bounds on the ''generalised weighted star discrepancy''. This version of the discrepancy corresponds to the classic L"~ weighted star discrepancy via a mapping to the unit cube. The weights here are general weights rather than the product weights considered in earlier works on integrals over R^d. Known methods based on an averaging argument are used to show the existence of these lattice rules, while the component-by-component technique is used to construct the generating vector of these shifted lattice rules. We prove that the bound on the weighted star discrepancy considered here is of order O(n^-^1^+^@d) for any @d0 and with the constant involved independent of the dimension. This convergence rate is better than the O(n^-^1^/^2) achieved so far for both Monte Carlo and quasi-Monte Carlo methods.