Large scale random phase calculations for direct self-consistent field wavefunctions