# Gaussian process emulators for computer experiments with inequality constraints

arXiv: https://arxiv.org/abs/1606.01265 · HTML: https://mlsanity.com/paper/1606.01265/

## Key facts

- Cited by 17 papers in the map (data through 14 Oct 2026).
- First posted on arXiv 24 May 2016.

## Built on by (most cited first, 17 of 17)

- [Sensing Cox Processes via Posterior Sampling and Positive Bases](https://mlsanity.com/paper/2110.11181/) (2021)
- [High-dimensional additive Gaussian processes under monotonicity constraints](https://mlsanity.com/paper/2205.08528/) (2022)
- [Sampling from Gaussian Processes: A Tutorial and Applications in Global Sensitivity Analysis and Optimization](https://mlsanity.com/paper/2507.14746/) (2025)
- [Incorporating Sum Constraints into Multitask Gaussian Processes](https://mlsanity.com/paper/2202.01793/) (2022)
- [Relaxed Gaussian process interpolation: a goal-oriented approach to Bayesian optimization](https://mlsanity.com/paper/2206.03034/) (2022)
- [Adaptive Resolution for Finite-Rank Gaussian Processes](https://mlsanity.com/paper/2505.24066/) (2025)
- [Bayesian preference elicitation for decision support in multiobjective optimization](https://mlsanity.com/paper/2507.16999/) (2025)
- [LatentFlow: A General Framework for Conditioning Stochastic Processes](https://mlsanity.com/paper/2607.12922/) (2026)
- [Hierarchical Non-Stationary Temporal Gaussian Processes With L1-Regularization](https://mlsanity.com/paper/2105.09695/) (2021)
- [Gaussian Processes with State-Dependent Noise for Stochastic Control](https://mlsanity.com/paper/2305.16229/) (2023)
- [Spectrum Gaussian Processes Based On Tunable Basis Functions](https://mlsanity.com/paper/2107.06473/) (2021)
- [Beyond Surrogate Modeling: Learning the Local Volatility Via Shape Constraints](https://mlsanity.com/paper/2212.09957/) (2022)
- [Gaussian Process Regression with Soft Inequality and Monotonicity Constraints](https://mlsanity.com/paper/2404.02873/) (2024)
- [Determination of Particle-Size Distributions from Light-Scattering Measurement Using Constrained Gaussian Process Regression](https://mlsanity.com/paper/2507.03736/) (2025)
- [Fast Gaussian Processes under Monotonicity Constraints](https://mlsanity.com/paper/2507.06677/) (2025)
- [Multi-output Gaussian process prediction of physical fields under linear equality constraints](https://mlsanity.com/paper/2608.25709/) (2026)
- [Constrained Probabilistic Numerics: Gaussian Process Solvers for Linear PDEs under Linear Inequality Constraints](https://mlsanity.com/paper/2610.13804/) (2026)

## Abstract (arXiv)

Physical phenomena are observed in many fields (sciences and engineering) and are often studied by time-consuming computer codes. These codes are analyzed with statistical models, often called emulators. In many situations, the physical system (computer model output) may be known to satisfy inequality constraints with respect to some or all input variables. Our aim is to build a model capable of incorporating both data interpolation and inequality constraints into a Gaussian process emulator. By using a functional decomposition, we propose to approximate the original Gaussian process by a finite-dimensional Gaussian process such that all conditional simulations satisfy the inequality constraints in the whole domain. The mean, mode (maximum a posteriori) and prediction intervals (uncertainty quantification) of the conditional Gaussian process are calculated. To investigate the performance of the proposed model, some conditional simulations with inequality constraints such as boundary, monotonicity or convexity conditions are given. 1. Introduction. In the engineering activity, runs of a computer code can be expensive and time-consuming. One solution is to use a statistical surrogate for conditioning computer model outputs at some input locations (design points). Gaussian process (GP) emulator is one of the most popular choices [23]. The reason comes from the property of the GP that uncertainty quantification can be calculated. Furthermore, it has several nice properties. For example, the conditional GP at observation data (linear equality constraints) is still a GP [5]. Additionally, some inequality constraints (such as monotonicity and convexity) of output computer responses are related to partial derivatives. In such cases, the partial derivatives of the GP are also GPs. Incorporating an infinite number of linear inequality constraints into a GP emulator, the problem becomes more difficult. The reason is that the resulting conditional process is not a GP. In the literature of interpolation with inequality constraints, we find two types of meth-ods. The first one is deterministic and based on splines, which have the advantage that the inequality constraints are satisfied in the whole input domain (see e.g. [16], [24] and [25]). The second one is based on the simulation of the conditional GP by using the subdivision of the input set (see e.g. [1], [6] and [11]). In that case, the inequality constraints are satisfied in a finite number of input locations. Notice that the advantage of such method is that un-certainty quantification can be calculated. In previous work, some methodologies have been based on the knowledge of the derivatives of the GP at some input locations ([11], [21] and [26]). For monotonicity constraints with noisy data, a Bayesian approach was developed in [21]. In [11] the problem is to build a GP emulator by using the prior monotonicity