Adaptive multi-resolution Gaussian processes: Scalable exact inference with naturally data-sparse covariance matrices
Not on the map. Why?
Abstract (arXiv)
Gaussian processes constitute a cornerstone of probabilistic machine learning, yet scaling them to large datasets typically forces a trade-off between computational efficiency and model fidelity. This work bridges this gap by presenting an adaptive multi-resolution Gaussian process framework that is both scalable and exact. Our key innovation is constructing a naturally data-sparse covariance matrix with adaptive multi-resolution basis functions. These basis functions are directly anchored to samples, eliminating the need for auxiliary points. By shrinking the support domains of multi-resolution basis, the matrix block sizes are limited, guaranteeing sparsity. The inverse of the data-sparse covariance matrix is computed exactly and efficiently via the sparse Cholesky inverse algorithm. To further improve predictive uncertainties, we construct an augmented basis function. Theoretical analysis and numerical experiments demonstrate that our model achieves exact inference with O(n log2 n) training cost and O(logd n) prediction cost, establishing a principled framework for scalable and high-fidelity Gaussian process regression.
Key facts. No citing papers in the map yet (data through 7 Oct 2026). Builds on 32 of its 69 references that are in the map. First posted on arXiv 24 Sep 2026.
Builds on
32 of 69 references.
References · 32
- 2013Gaussian Processes for Big Datacited by 377
- 2018When Gaussian Process Meets Big Data: A Review of Scalable GPscited by 205
- 2015Kernel Interpolation for Scalable Structured Gaussian Processes (KISS-GP)cited by 175
- 2017Machine Learning of Linear Differential Equations using Gaussian Processescited by 139
- 2019Exact Gaussian Processes on a Million Data Pointscited by 118
- 2017Numerical Gaussian Processes for Time-dependent and Non-linear Partial Differential Equationscited by 72
- 2014Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geostatistical Datasetscited by 71
- 2017A general framework for Vecchia approximations of Gaussian processescited by 66
- 2013Local Gaussian process approximation for large computer experimentscited by 59
- 2017A Case Study Competition Among Methods for Analyzing Large Spatial Datacited by 44
- 2020Sparse Cholesky factorization by Kullback-Leibler minimizationcited by 40
- 2017Scalable Log Determinants for Gaussian Process Kernel Learningcited by 34
20 more
- 2018Product Kernel Interpolation for Scalable Gaussian Processescited by 31
- 2023Sparse Cholesky Factorization for Solving Nonlinear PDEs via Gaussian Processescited by 31
- 2016Permutation and Grouping Methods for Sharpening Gaussian Process Approximationscited by 24
- 2017Compression, inversion, and approximate PCA of dense kernel matrices at near-linear computational complexitycited by 24
- 2022Posterior and Computational Uncertainty in Gaussian Processescited by 23
- 2015A multi-resolution approximation for massive spatial datasetscited by 20
- 2021Local approximate Gaussian process regression for data-driven constitutive laws: Development and comparison with neural networkscited by 17
- 2020Scaled Vecchia approximation for fast computer-model emulationcited by 15
- 2021Bayesian Numerical Methods for Nonlinear Partial Differential Equationscited by 15
- 2022Basis-Function Models in Spatial Statisticscited by 9
- 2022Exact Gaussian Processes for Massive Datasets via Non-Stationary Sparsity-Discovering Kernelscited by 9
- 2022Variational Nearest Neighbor Gaussian Processcited by 8
- 2016Fast spatial Gaussian process maximum likelihood estimation via skeletonization factorizationscited by 6
- 2022Gaussian Process Hydrodynamicscited by 4
- 2023Actually Sparse Variational Gaussian Processescited by 4
- 2017A class of multi-resolution approximations for large spatial datasetscited by 3
- 2025All Emulators are Wrong, Many are Useful, and Some are More Useful Than Others: A Reproducible Comparison of Computer Model Surrogatescited by 3
- 2018Conditionally Independent Multiresolution Gaussian Processescited by 1
- 2023Accurate and Scalable Stochastic Gaussian Process Regression via Learnable Coreset-based Variational Inferencecited by 1
- 2025Mind the Jumps: A Scalable Robust Local Gaussian Process for Multidimensional Response Surfaces with Discontinuitiescited by 1
Built on by
Not cited yet.