pub trait SymmetricSqrt { type Output; fn ssqrt(&self, uplo: UPLO) -> Result<Self::Output>; }
Calculate symmetric square-root matrix using eigh
eigh