After localization, the resulting wavefunction is treated as a single determinant.
This causes the orbital occupation function to early return the mean field integer occupations instead of the actual one-RDM diagonal entries.