Efficient Marginal Likelihood Computation for Gaussian Process Regression