powered by
Same thing as 'expm::sqrtm', but faster.
sqrt_matrix(mat, symmetric)
Matrix to find square root matrix of
Is it symmetric? Passed to eigen.
Square root of mat
# NOT RUN { mat <- matrix(c(1,.1,.1,1), 2, 2) smat <- sqrt_matrix(mat=mat, symmetric=TRUE) smat %*% smat # }
Run the code above in your browser using DataLab