Typical computational time necessary to determine hyperparameters as a function of the number of data points used during training for a classical maximum likelihood method using Cholesky decomposition (red triangles) and when using the two-point correlation function computed by TreeCorr (blue circles). The typical number of training points for astrometric residuals modeled using GPs is between ∼10³ and ∼10⁴.