Return inverse square root of square positive-definite matrix.
invsqrt(A)
A
inverse square root of square positive-definite matrix.
require(Rsfar) X <- Bdiag(matrix(1:4,2,2), 3) invsqrt(t(X) %*% X)