Implements a maximum likelihood estimation (MLE)
method for estimation and prediction in spatially varying coefficient
(SVC) models (Dambon et al. (2020) ). Covariance
tapering (Furrer et al. (2006) ) can be
applied such that the method scales to large data.