Computes the generalized inverse of a matrix X. This function is used in
the m2LL.FPBK.nodet
functions in order
to estimate the spatial covariance parameters
mginv(X, tol = sqrt(.Machine$double.eps))
The matrix to be inverted
The tolerance of the estimation
The generalized inverse matrix