English

Lattice Rules Meet Kernel Cubature

Numerical Analysis 2025-06-06 v2 Numerical Analysis Statistics Theory Statistics Theory

Abstract

Rank-1 lattice rules are a class of equally weighted quasi-Monte Carlo methods that achieve essentially linear convergence rates for functions in a reproducing kernel Hilbert space (RKHS) characterized by square-integrable first-order mixed partial derivatives. In this work, we explore the impact of replacing the equal weights in lattice rules with optimized cubature weights derived using the reproducing kernel. We establish a theoretical result demonstrating a doubled convergence rate in the one-dimensional case and provide numerical investigations of convergence rates in higher dimensions. We also present numerical results for an uncertainty quantification problem involving an elliptic partial differential equation with a random coefficient.

Keywords

Cite

@article{arxiv.2501.09500,
  title  = {Lattice Rules Meet Kernel Cubature},
  author = {Vesa Kaarnioja and Ilja Klebanov and Claudia Schillings and Yuya Suzuki},
  journal= {arXiv preprint arXiv:2501.09500},
  year   = {2025}
}

Comments

17 pages, 2 figures

R2 v1 2026-06-28T21:08:16.441Z