Gaussian Process will return an (n,) array is the output is one dimensional, and an (n,m) array if the output is y dimensional. This makes it backwards compatible