powered by
###An function to adjust
gmi(X, tol = sqrt(.Machine$double.eps))
a square matrix
a real value