To reduce memory costs, we compute our QNO overlap matrices "semi-directly", following the terminology of our prior work. In the semidirect algorithm, we first determine the extended PAO domain of ijkl by merging the PAO domain of ijkl with the PAO domain of all LMOs m such that mjkl, imkl, ijml, or ijkm are in the list of nonunique quadruplets.