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.

Conditionally Gaussian PAC-Bayes

Improving Generalization of Complex Models under Unbounded Loss Using PAC-Bayes Bounds

Learning Gaussian Networks

Learning via Surrogate PAC-Bayes

Misclassification bounds for PAC-Bayesian sparse deep learning

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

Scalable PAC-Bayesian Meta-Learning via the PAC-Optimal Hyper-Posterior: From Theory to Practice

Scalable Bayesian Optimization Using Vecchia Approximations of Gaussian Processes

Controlling Multiple Errors Simultaneously with a PAC-Bayes Bound

Hybrid kernel approach to improving the numerical stability of machine learning for parametric equations with Gaussian processes in the noisy and noise-free data assumptions

Learning-to-Optimize with PAC-Bayesian Guarantees: Theoretical Considerations and Practical Implementation

Neural Conditional Probability for Inference

PAC-Bayesian Theory Meets Bayesian Inference

Streaming PAC-Bayes Gaussian process regression with a performance guarantee for online decision making

The haematopathological complications of gold therapy: effects of splenectomy and bal.

Better-than-KL PAC-Bayes Bounds

Predictive Coding beyond Gaussian Distributions

Iterative Construction of Gaussian Process Surrogate Models for Bayesian Inference

Probabilistic fine-tuning of pruning masks and PAC-Bayes self-bounded learning

A Bayesian Take on Gaussian Process Networks

Time-Varying Gaussian Process Bandits with Unknown Prior