Gaussian Process Model
A Gaussian-process (GP) regression surrogate with per-point observation noise, suitable for the stochastic simulation setting: each observed point carries its own noise variance (the variance of the estimated mean). The GP uses a constant prior mean (the data mean by default) and a StationaryKernel.
The fit forms the covariance matrix K + diag(noiseVar) + jitter*I, factorizes it with a Cholesky decomposition (Hipparchus, already on the classpath), and caches the weight vector and the inverse for prediction. If the matrix is not numerically positive definite, the jitter is escalated and the fit retried.
Parameters
the problem definition (used for default kernel length scales)
the covariance kernel. Defaults to an RBFKernel with ARD length scales set to the problem's input ranges.
the constant prior mean. When null (the default), the mean of the observed responses is used (ordinary kriging).
a small positive value added to the diagonal for numerical stability. Defaults to defaultJitter.
Constructors
Functions
Fits the surrogate to the supplied observations.
The log marginal likelihood of the current fit, used for hyperparameter selection: -0.5 (y-m)^T K^{-1} (y-m) - 0.5 log|K| - (n/2) log(2π).
Predicts the posterior mean and variance at x. The model must have been fit first.