Local auxiliary domains are assigned to quadruples similarly to how they are assigned in the case of pairs and triples, to increase the efficiency of the computation of two-electron integrals. In density-fitting/resolution of the identity, two-electron integrals are approximated using an auxiliary basis set. Using this approximation, standard, four-center two-electron integrals are approximated as