Changes
- The variance is now determined from using generalised least squares on a number of random samples from the Gaussian process, this has the benefit of being more accurate and efficient.
- As a result of the above there is no longer a "Sampling Likelihood" `tqdm` bar.
- There is now a `rtol` kwarg for the `bootstap_GLS` method, this controls the `rtol` in the call to the `scipy.linalg.pinvh` that is used to stabilise the large matrices.
- All documentation has been adapted to reflect this methodological change.