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.

Batch Active Learning in Gaussian Process Regression using Derivatives

Active Learning for Derivative-Based Global Sensitivity Analysis with Gaussian Processes

Scaling Gaussian Process Regression with Derivatives

Active Learning Methods with Deep Gaussian Processes.

DGP-LVM: Derivative Gaussian process latent variable models

Batched Gaussian Process Bandit Optimization via Determinantal Point Processes

Efficient Active Learning for Gaussian Process Classification by Error Reduction.

Black-Box Batch Active Learning for Regression

Batch-Mode Active Learning of Gaussian Process Regression with Maximum Model Change

Gaussian process learning of nonlinear dynamics

A Mini-Batch Method for Solving Nonlinear PDEs with Gaussian Processes

Active Learning for Abrupt Shifts Change-point Detection via Derivative-Aware Gaussian Processes

A global optimization of SVM batch active learning

Exploiting gradients and Hessians in Bayesian optimization and Bayesian quadrature

Active Learning for Gaussian Process Considering Uncertainties with Application to Shape Control of Composite Fuselage

Near-Optimal Active Learning of Multi-Output Gaussian Processes

Enhanced Gaussian Process Regression for Active Learning Model-based Predictive Control

Convex Batch Mode Active Sampling Via Α-Relative Pearson Divergence.

Epoch-Evolving Gaussian Process Guided Learning for Classification

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

Natural Gradients in Practice: Non-Conjugate Variational Inference in Gaussian Process Models