Optimal Monte Carlo integration on matrix manifolds
Abstract
This article develops optimal order (up to logarithmic factors) weighted cubature rules for Monte Carlo numerical integration in Sobolev reproducing kernel Hilbert spaces over common matrix manifolds. Specifically, we consider matrix models of the Stiefel, Grassmann, flag, and Lagrangian Grassmann manifolds each over either real numbers ℝ, complex numbers ℂ, or quaternions ℍ. The proposed cubature rules use evaluation points which are independent and identically distributed uniform on the manifold. The weights are chosen to be optimal with respect to the reproducing kernel Hilbert space worst-case error, under the constraint that the weights sum to one. Crucially, this constraint ensures optimal weights may be found without having to evaluate the kernel mean embedding (the integral of the kernel), which is often difficult to compute. We consider computable Matérn kernels of Frobenius distances between matrices which reproduce the Sobolev reproducing kernel Hilbert spaces with equivalent norms. Open-source numerical experiments validate our theoretical development, including an application to simultaneous approximation of intrinsic volumes in the Cauchy--Kubota formula. We also present preliminary empirical results beyond independent points including transformations of low-discrepancy quasi-Monte Carlo points from the unit cube and jointly optimized points and weights.
Related
- Slice Monte Carlo Integration
- Quasi-Monte Carlo based meshfree methods for parabolic equations on parameterized Riemannian manifolds
- Median randomized quasi-Monte Carlo for nested integration
- Optimal sampling strategies in event-chain Monte Carlo
- Sequential Entropic Optimal Transport for Bayesian Inverse Problems via Multilevel Monte Carlo
BibTeX
@misc{sorokin2026optimal,
title = {{Optimal Monte Carlo integration on matrix manifolds}},
author = {Aleksei G. Sorokin and Yihan Guo and Lek-Heng Lim},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.16965},
eprint = {2610.16965},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.16965},
}
RIS
TY - UNPB
TI - Optimal Monte Carlo integration on matrix manifolds
AU - Aleksei G. Sorokin
AU - Yihan Guo
AU - Lek-Heng Lim
PY - 2026
DA - 2026/10/12
PB - arXiv
DO - 10.48550/arxiv.2610.16965
UR - https://arxiv.org/abs/2610.16965
AB - This article develops optimal order (up to logarithmic factors) weighted cubature rules for Monte Carlo numerical integration in Sobolev reproducing kernel Hilbert spaces over common matrix manifolds. Specifically, we consider matrix models of the Stiefel, Grassmann, flag, and Lagrangian Grassmann manifolds each over either real numbers ℝ, complex numbers ℂ, or quaternions ℍ. The proposed cubature rules use evaluation points which are independent and identically distributed uniform on the manifold. The weights are chosen to be optimal with respect to the reproducing kernel Hilbert space worst-case error, under the constraint that the weights sum to one. Crucially, this constraint ensures optimal weights may be found without having to evaluate the kernel mean embedding (the integral of the kernel), which is often difficult to compute. We consider computable Matérn kernels of Frobenius distances between matrices which reproduce the Sobolev reproducing kernel Hilbert spaces with equivalent norms. Open-source numerical experiments validate our theoretical development, including an application to simultaneous approximation of intrinsic volumes in the Cauchy--Kubota formula. We also present preliminary empirical results beyond independent points including transformations of low-discrepancy quasi-Monte Carlo points from the unit cube and jointly optimized points and weights.
ER -