from pyquante2.utils import dmat, geigh


def core_hamiltonian_rhf(H):
    name = 'Core hamiltonian guess'
    orbe, orbs = H.eigenv(H.h)
    return H.density(orbs)


def core_hamiltonian_uhf(H):
    name = 'Core hamiltonian guess'
    (orbea, orbsa), (orbeb, orbsb) = H.eigenv(H.h, H.h)
    return H.density(orbsa, orbsb)


def atomic_densities(H):
    name = 'Superposition of atomic densities guess'

@@ -314,24 +314,150 @@
 a, orbsb)%0A%0A%0A
+def core_hamiltonian_rohf(H):%0A    name = 'Core hamiltonian guess'%0A    orbe, orbs = H.eigenv(H.h)%0A    return H.density(orbs)%0A%0A%0A
 def atomic_d

implement ROHF