1 Introduction Gaussian processes (GPs) are a preeminent framework for stochastic function approximation, statistical modeling of real-world measurements, non-parametric and nonlinear regression within machine learning (ML), and surrogate modeling. GPs are analytically tractable and can be fully specified in terms of a prior mean function and a kernel, also known as a covariance function.