Optimal numerical integration for functions in fractional Gaussian Sobolev spaces
It addresses a theoretical challenge in numerical analysis for high-dimensional integration in Gaussian-weighted spaces, providing incremental advancements in quadrature methods.
This paper tackles the problem of numerical integration for functions in fractional Gaussian Sobolev spaces, constructing quadrature schemes on ℝ^d that achieve optimal asymptotic convergence rates, specifically for parameters 1 < p < ∞ and s > 1/p, and showing equivalence to Hermite spaces for p=2 with s > 1/2.
This paper investigates the numerical approximation of integrals for functions in fractional Gaussian Sobolev spaces $W^s_{p}(\mathbb{R}^d,γ)$ with dominating mixed smoothness defined via kernel related to the fractional Ornstein-Uhlenbeck operator. Building upon quadrature rules for fractional Sobolev spaces on the unit cube $[-\tfrac{1}{2}, \tfrac{1}{2}]^d$, we construct quadrature schemes on $\mathbb{R}^d$ that achieve the same rate of convergence. As a consequence, we establish the optimal asymptotic order of the integration error in the regime $1 < p < \infty$ and $s > \frac{1}{p}$. Furthermore, we show that the fractional Gaussian Sobolev spaces $W^s_{2}(\mathbb{R}^d,γ)$ coincide with Hermite spaces $\mathcal{H}^s(\mathbb{R}^d,γ)$ characterized by the weighted $\ell_2$-summability of their Fourier-Hermite coefficients. From this, we derive the optimal asymptotic order of the integration error for functions in these spaces for all $s > \frac{1}{2}$. We also establish the corresponding optimal asymptotic order for functions in fractional Sobolev spaces $W^s_{p,G}(\mathbb{R}^d,γ)$ defined via the Gagliardo seminorm.