Relaxed Gaussian Process Interpolation: a Goal-Oriented Approach to Bayesian Optimization (JMLR, 2026)