code: def geigh(H,S):
      "Solve the generalized eigensystem Hc = ESc"
      A = cholorth(S)
      E,U = np.linalg.eigh(simx(H,A))
      return E,np.dot(A,Uf geigh(H,S):
      "Solve the generalized eigensystem Hc = ESc"
      A = cholorth(S)
      E,U = np.linalg.eigh(simx(H,A))
      return E,np.dot(A,U)
docs: Solve the generalized eigensystem Hc = ESc