The trick is I am setting gp.theta0 to the theta_ estimated from the full dataset (first call to `fit`), and setting thetaL and thetaU to None to force theta_ to remain the same (that's what's called the given case in the `fit`method, grep for it) during the cross-validation.