I want to get reproducible results from GP.sample. This function could take an optional parameter rng of type numpy.random.Generator. Rather than call np.random.multivariate_normal, we could then call ...