Median QMC method for unbounded integrands over $\mathbb{R}^s$ in unanchored weighted Sobolev spaces

This paper investigates quasi-Monte Carlo (QMC) integration of Lebesgue integrable functions with respect to a density function over $\mathbb{R}^s$. We extend the construction-free median QMC rule proposed by Goda and L'ecuyer (SIAM J. Sci. Comput., 2022) to the weighted unanchored Sobolev space of functions defined over $\mathbb{R}^s$ introduced by Nichols and Kuo (J. Complexity, 2014). By taking the median of $k = \mathcal{O}(\log N)$ independent randomized QMC estimators, we prove that for any $ε\in (0,r-\frac{1}{2}]$, our method achieves a mean absolute error bound of $\mathcal{O}(N^{-r+ε})$, where $N$ is the number of points and $r>\frac{1}{2}$ is a parameter determined by the function space. This rate matches the rate of randomly shifted lattice rules obtained via a component-by-component (CBC) construction, while our approach requires no specific CBC constructions or prior knowledge of the space's weight structure. Numerical experiments demonstrate that our method attains an accuracy comparable to the CBC construction based method, and outperforms the Monte Carlo method.

Publication Details

Journal
ArXiv.org
Published
2025-03-07
Primary Topic
Economic Growth and Fiscal Policies
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Median QMC method for unbounded integrands over $\mathbb{R}^s$ in unanchored weighted Sobolev spaces

ArXiv.org
Economic Growth and Fiscal Policies
article

Median QMC method for unbounded integrands over $\mathbb{R}^s$ in unanchored weighted Sobolev spaces

article en

Abstract

This paper investigates quasi-Monte Carlo (QMC) integration of Lebesgue integrable functions with respect to a density function over $\mathbb{R}^s$. We extend the construction-free median QMC rule proposed by Goda and L'ecuyer (SIAM J. Sci. Comput., 2022) to the weighted unanchored Sobolev space of functions defined over $\mathbb{R}^s$ introduced by Nichols and Kuo (J. Complexity, 2014). By taking the median of $k = \mathcal{O}(\log N)$ independent randomized QMC estimators, we prove that for any $ε\in (0,r-\frac{1}{2}]$, our method achieves a mean absolute error bound of $\mathcal{O}(N^{-r+ε})$, where $N$ is the number of points and $r>\frac{1}{2}$ is a parameter determined by the function space. This rate matches the rate of randomly shifted lattice rules obtained via a component-by-component (CBC) construction, while our approach requires no specific CBC constructions or prior knowledge of the space's weight structure. Numerical experiments demonstrate that our method attains an accuracy comparable to the CBC construction based method, and outperforms the Monte Carlo method.

ArXiv.org
Openalex Percentile: Top 40%
Economic Growth and Fiscal Policies
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.