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