Alllow GP interface to return prediction variance

[iraikov]

Currently, the predict methods in the Gaussian Process interface module only return the mean at the requested points. However, the variance is also necessary to estimate the uncertainty at those points.