Scalable Bayesian Optimization Using Vecchia Approximations of Gaussian Processes
Not on the map. Why?
Abstract (arXiv)
Bayesian optimization is a technique for optimizing black-box target functions. At the core of Bayesian optimization is a surrogate model that predicts the output of the target function at previously unseen inputs to facilitate the selection of promising input values. Gaussian processes (GPs) are commonly used as surrogate models but are known to scale poorly with the number of observations. We adapt the Vecchia approximation, a popular GP approximation from spatial statistics, to enable scalable high-dimensional Bayesian optimization. We develop several improvements and extensions, including training warped GPs using mini-batch gradient descent, approximate neighbor search, and selecting multiple input values in parallel. We focus on the use of our warped Vecchia GP in trust-region Bayesian optimization via Thompson sampling. On several test functions and on two reinforcement-learning problems, our methods compared favorably to the state of the art.
Key facts. Cited by 11 papers in the map (data through 7 Oct 2026). Builds on 28 of its 54 references that are in the map. First posted on arXiv 2 Mar 2022.
Builds on
28 of 54 references.
References · 28
- 2017Billion-scale similarity search with GPUscited by 2,477
- 2016Automatic chemical design using a data-driven continuous representation of moleculescited by 762
- 2018A Tutorial on Bayesian Optimizationcited by 734
- 2018GPyTorch: Blackbox Matrix-Matrix Gaussian Process Inference with GPU Accelerationcited by 603
- 2015Deep Kernel Learningcited by 440
- 2013Gaussian Processes for Big Datacited by 377
- 2019Scalable Global Optimization via Local Bayesian Optimizationcited by 359
- 2018Neural Architecture Search with Bayesian Optimisation and Optimal Transportcited by 192
- 2021High-Dimensional Bayesian Optimization with Sparse Axis-Aligned Subspacescited by 139
- 2015High Dimensional Bayesian Optimisation and Bandits via Additive Modelscited by 137
- 2017Batched Large-scale Bayesian Optimization in High-dimensional Spacescited by 116
- 2018Maximizing acquisition functions for Bayesian optimizationcited by 113
16 more
- 2017Constrained Bayesian Optimization with Noisy Experimentscited by 103
- 2016The Parallel Knowledge Gradient Method for Batch Bayesian Optimizationcited by 73
- 2014Input Warping for Bayesian Optimization of Non-stationary Functionscited by 66
- 2014Generalized Product of Experts for Automatic and Principled Fusion of Gaussian Process Predictionscited by 66
- 2017A general framework for Vecchia approximations of Gaussian processescited by 66
- 2017The reparameterization trick for acquisition functionscited by 61
- 2016Parallel Bayesian Global Optimization of Expensive Functionscited by 43
- 2020Sparse Cholesky factorization by Kullback-Leibler minimizationcited by 40
- 2016Permutation and Grouping Methods for Sharpening Gaussian Process Approximationscited by 24
- 2021Conditioning Sparse Variational Gaussian Processes for Online Decision-makingcited by 23
- 2018Vecchia approximations of Gaussian-process predictionscited by 21
- 2020Sparse within Sparse Gaussian Processes using Neighbor Informationcited by 18
- 2020Scaled Vecchia approximation for fast computer-model emulationcited by 15
- 2021Correlation-based sparse inverse Cholesky factorization for fast Gaussian-process inferencecited by 12
- 2022Variational Nearest Neighbor Gaussian Processcited by 8
- 2021Correlated Product of Experts for Sparse Gaussian Process Regressioncited by 7
Built on by
Citing papers per year: 2022: 1; 2023: 2; 2024: 1; 2025: 4; 2026: 3.
11 papers.
- 2024Scalable Bayesian Optimization via Focalized Sparse Gaussian Processescited by 8
- 2025Collaborative Bayesian Optimization via Wasserstein Barycenterscited by 7
- 2023Vecchia Gaussian Process Ensembles on Internal Representations of Deep Neural Networkscited by 3
- 2023Scalable Model-Based Gaussian Process Clusteringcited by 1
- 2022Optimistic Optimization of Gaussian Process Samplescited by 0
- 2025Gradient-based Sample Selection for Faster Bayesian Optimizationcited by 0
- 2025Taking the GP Out of the Loopcited by 0
- 2025BOOST: A Data-Driven Framework for the Automated Joint Selection of Kernel and Acquisition Functions in Bayesian Optimizationcited by 0
- 2026Amortized Acquisition Optimization for Bayesian Optimization with Variational Mutual Informationcited by 0
- 2026Bayesian Optimization in Linear Timecited by 0
- 2026Scalable Derivative Gaussian Processes via Exact Gradient Reductioncited by 0
- 2026Scalable Derivative Gaussian Processes via Exact Gradient Reduction1 Jun 2026
- 2026Bayesian Optimization in Linear Time30 Apr 2026
- 2026Amortized Acquisition Optimization for Bayesian Optimization with Variational Mutual Information13 Jan 2026
- 2025BOOST: A Data-Driven Framework for the Automated Joint Selection of Kernel and Acquisition Functions in Bayesian Optimization4 Aug 2025
- 2025Taking the GP Out of the Loop15 Jun 2025
- 2025Collaborative Bayesian Optimization via Wasserstein Barycenters15 Apr 2025
- 2025Gradient-based Sample Selection for Faster Bayesian Optimization10 Apr 2025
- 2024Scalable Bayesian Optimization via Focalized Sparse Gaussian Processes29 Dec 2024
- 2023Scalable Model-Based Gaussian Process Clustering14 Sep 2023
- 2023Vecchia Gaussian Process Ensembles on Internal Representations of Deep Neural Networks26 May 2023
- 2022Optimistic Optimization of Gaussian Process Samples2 Sep 2022