Linearly constrained Gaussian processes