Abstract:Gaussian processes provide probabilistic surrogates for various applications including classification, uncertainty quantification, and optimization. Using a gradient-enhanced covariance matrix can be beneficial since it provides a more accurate surrogate relative to its gradient-free counterpart. An acute problem for Gaussian processes, particularly those that use gradients, is the ill-conditioning of their covariance matrices. Several methods have been developed to address this problem for gradient-enhanced Gaussian processes but they have various drawbacks such as limiting the data that can be used, imposing a minimum distance between evaluation points in the parameter space, or constraining the hyperparameters. In this paper a new method is presented that applies a diagonal preconditioner to the covariance matrix along with a modest nugget to ensure that the condition number of the covariance matrix is bounded, while avoiding the drawbacks listed above. Optimization results for a gradient-enhanced Bayesian optimizer with the Gaussian kernel are compared with the use of the new method, a baseline method that constrains the hyperparameters, and a rescaling method that increases the distance between evaluation points. The Bayesian optimizer with the new method converges the optimality, ie the $\ell_2$ norm of the gradient, an additional 5 to 9 orders of magnitude relative to when the baseline method is used and it does so in fewer iterations than with the rescaling method. The new method is available in the open source python library GpGradPy, which can be found at <a class="link-external link-https" href="https://github.com/marchildon/gpgradpy/tree/paper_precon" rel="external noopener nofollow">this https URL</a>. All of the figures in this paper can be reproduced with this library.

A Matrix-free Approach for Solving the Parametric Gaussian Process Maximum Likelihood Problem

A MATRIX-FREE APPROACH FOR SOLVING THE GAUSSIAN PROCESS MAXIMUM LIKELIHOOD PROBLEM

A Solution to the Ill-Conditioning of Gradient-Enhanced Covariance Matrices for Gaussian Processes

Representing Additive Gaussian Processes by Sparse Matrices

Vecchia Gaussian Processes: Probabilistic Properties, Minimax Rates and Methodological Developments

A computationally efficient symmetric diagonally dominant matrix projection-based Gaussian process approach

Bayesian and Maximum Likelihood Estimation for Gaussian Processes on an Incomplete Lattice

Contraction rates for conjugate gradient and Lanczos approximate posteriors in Gaussian process regression

MCMC methods for Gaussian process models using fast approximations for the likelihood

Iterative Methods for Vecchia-Laplace Approximations for Latent Gaussian Process Models

Composite Inference for Gaussian Processes

Approximate Implementation of Logarithm of Matrix Determinant in Gaussian Process Regression and Numerical Experiments

Iterative Methods for Full-Scale Gaussian Process Approximations for Large Spatial Data

Warm Start Marginal Likelihood Optimisation for Iterative Gaussian Processes

Scaling Gaussian Process Regression with Derivatives

Scalable Physics-based Maximum Likelihood Estimation using Hierarchical Matrices

Mixture of Gaussian Processes and its Applications

Efficient Gaussian process based on BFGS updating and logdet approximation

Inference for Gaussian Processes with Matérn Covariogram on Compact Riemannian Manifolds