We would like to extract the two-electron integrals from Psi4. This piece of code:
mints = MintsHelper(wfn.basisset())
TWO = mints.ao_eri()
prints the entire array, but can we get a folded list of ERIs (i.e. only those that cannot be obtained from symmetry)?