code: def cholorth(S):
      "Cholesky orthogonalization"
      return np.linalg.inv(np.linalg.cholesky(S)).f cholorth(S):
      "Cholesky orthogonalization"
      return np.linalg.inv(np.linalg.cholesky(S)).T
docs: Cholesky orthogonalization