Gaussian process (GP) regression is a non-parametric, Bayesian framework to
approximate complex models. Standard GP regression can lead to an unbounded
model in which some points can take infeasible values. We introduce a new GP
method that enforces the physical constraints in a probab