get_symmetric_sqrt function

Calculate symmetric square root matrix of a positive definite covariance matrix