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.

Spectral Preconditioning for Gradient Methods on Graded Non-convex Functions

Convergence of Projected Subgradient Method with Sparse or Low-Rank Constraints

Stochastic Gradient Methods with Preconditioned Updates

Spectral Method and Its Application to the Conjugate Gradient Method

A New Spectral Conjugate Subgradient Method with Application in Computed Tomography Image Reconstruction

Dual Space Preconditioning for Gradient Descent

A Spectral Conjugate Gradient Method for Nonlinear Inverse Problems

Spectral conjugate gradient methods for vector optimization problems

Some Nonlinear Conjugate Gradient Methods Based on Spectral Scaling Secant Equations

A family of spectral conjugate gradient method with strong convergence and its applications in image restoration and machine learning

Tight Convergence Rate Bounds for Optimization Under Power Law Spectral Conditions

An Improved Spectral Conjugate Gradient Algorithm for Nonconvex Unconstrained Optimization Problems

Convex quartic problems: homogenized gradient method and preconditioning

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

PROMISE: Preconditioned Stochastic Optimization Methods by Incorporating Scalable Curvature Estimates

Simba: A Scalable Bilevel Preconditioned Gradient Method for Fast Evasion of Flat Areas and Saddle Points

Preconditioning for Accelerated Gradient Descent Optimization and Regularization

Spectral analysis for preconditioning of multi-dimensional Riesz fractional diffusion equations

Gradient is All You Need?