The surrogate. With observations y at points X, the Gaussian process prediction at a new point x* is a weighted average of the observed scores, weighted by proximity — plus an uncertainty that shrinks near observed points:
μ(x*) = k*ᵀ (K + σₙ²I)⁻¹ (y − ȳ) + ȳ
σ²(x*) = k(x*, x*) − k*ᵀ (K + σₙ²I)⁻¹ k*
K is the n × n grid of kernel values between observed points;
k* is the n-vector of kernel values from x* to each observation.
The inverse is the price of the answer — and it is tiny next to
training the model itself.
The acquisition. With best score f_best so far, EI is
z = (μ(x) − f_best) / σ(x)
EI(x) = (μ(x) − f_best)·Φ(z) + σ(x)·φ(z)
Φ and φ are the standard-normal CDF and PDF:
Φ(z) is "probability", φ(z) is "height of the bell".
In words: improvement × how likely, plus uncertainty × its upside.
Worked comparison. Best so far is 0.7. Candidate A has μ = 0.8, σ = 0.25; candidate B has μ = 0.6, σ = 0.6. Which does EI evaluate?
A: z = (0.8 − 0.7)/0.25 = 0.4
Φ(0.4) = 0.6554, φ(0.4) = 0.3683
EI = 0.1 × 0.6554 + 0.25 × 0.3683 = 0.1576
B: z = (0.6 − 0.7)/0.6 = −0.1667
Φ(−0.1667) = 0.4338, φ(−0.1667) = 0.3934
EI = −0.1 × 0.4338 + 0.6 × 0.3934 = 0.1927
EI picks B — lower mean, higher uncertainty.
UCB with κ = 2 agrees: A = 0.8 + 0.5 = 1.30,
B = 0.6 + 1.2 = 1.80.
The lab in numbers. The four starting points leave the best at 0.9218 (x = 0.55). EI at that high-mean point is only 0.0146, while the untouched edge x = 1.0 has EI 0.0967 — so the first step explores the edge (score 0.5033). The next step exploits: x = 0.653, score 0.9812, then x = 0.633, 0.9970. Five EI steps in, the best is 0.9989; ten steps reach 0.9992. Ten guided evaluations did what random search typically needs 25–50 draws to match — the 2–5× fewer evaluations the source reports.