Compute symmetric square root of A, assuming it is real, symmetric, positive definite
symm_square_root(A, prec = 15)
Useful links