openfermion.linalg.erpa_eom_hamiltonian

Evaluate sum{a,b,c,d}h{a, b, d, c}

openfermion.linalg.erpa_eom_hamiltonian(
    h_ijkl: numpy.ndarray, tpdm: numpy.ndarray, p: int, q: int, r: int, s: int
) -> Union[float, complex]

h_ijkl two-body integral tensors of full reduced Hamiltonian
tpdm two-RDM
p left creation op index
q left annihilation op index
r right creation op index
s right annihilation op index

float or complex number