Gaussian process for arbitrary-dimensional output spaces

was #1611