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.

Gaussian Process Regression for Geometry Optimization

Global Optimization of Gaussian processes

Rectangularization of Gaussian process regression for optimization of hyperparameters

Using Gaussian process regression for efficient parameter reconstruction

Simultaneous and Meshfree Topology Optimization with Physics-informed Gaussian Processes

Scaling Gaussian Process Regression with Derivatives

Deep Gaussian process for multi-objective Bayesian optimization

Optimization of Gaussian Random Fields

A data-driven method for parametric PDE Eigenvalue Problems using Gaussian Process with different covariance functions

Optimistic Optimization of Gaussian Process Samples

Geometry-aware Bayesian Optimization in Robotics using Riemannian Matérn Kernels

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

Gaussian Processes and Bayesian Optimization for High Precision Experiments

Upgrading from Gaussian Processes to Student's-T Processes

Gaussian smoothing gradient descent for minimizing functions (GSmoothGD)

Applying Bayesian optimization with Gaussian process regression to computational fluid dynamics problems

Easy representation of multivariate functions with low-dimensional terms via Gaussian process regression kernel design: applications to machine learning of potential energy surfaces and kinetic energy densities from sparse data

GPfit: An R package for Gaussian Process Model Fitting using a New Optimization Algorithm

Efficient Two-Stage Gaussian Process Regression Via Automatic Kernel Search and Subsampling

Improving the Robustness of the Projected Gradient Descent Method for Nonlinear Constrained Optimization Problems in Topology Optimization

3DGS-LM: Faster Gaussian-Splatting Optimization with Levenberg-Marquardt