Fast calculation of Gaussian Process multiple-fold cross-validation residuals and their covariances