The author presents a parallel implementation of the Vecchia approximation technique, utilizing batched matrix computations on contemporary GPUs to speed up evaluating the log-likelihood function.
Gaussian processes are efficiently approximated using GPU-accelerated Vecchia methods for geospatial data analysis.