We describe our implementation of the multivariate Matérn model for spatial datasets, using Vecchia’s approximation and a Fisher scoring optimization algorithm. consider various pararameterizations that have been proposed in literature ensuring validity, as well an unconstrained model. A strength study is code tested on many real-world datasets. use it to effect ordering conditioning restrictio...