Add another way to compute S^-1/2 matrix