The numerical integration of the XC kernel benefits from a very efficient screening and, as made clear from the results, the ratio of computational efforts between the ERI and XC parts monotonically increases with system size. The same is not true for Fock matrix diagonalization in the SCF iterations, but even for the largest model with 11,934 basis functions, the associated wall time reaches a mere 11.3 seconds. The present implementation strategy, hardware, and math library make it affordable to treat systems up to and somewhat beyond 30,000 basis functions in the SCF optimization.