The SCF optimization involves the construction of a single Fock matrix at a time with a typical repetition of a few tens of times. In contrast, the response calculation involves the need to construct multiple Fock matrices with a number of matrices that can vary quite strongly over the course of the iterations in the iterative linear response equation solver. In the present work, we have not pursued any attempts to leverage the fact that the construction of multiple Fock matrices can be made in parallel in the sense that the set of ERIs can be computed only once but in conjunction with separate density matrix contractions and distributions, albeit with a penalty in the form of screening efficiency. The total number of Fock matrices constructed in the response part of the calculation typically amounts to several thousands in a spectrum simulation for a medium-or large-scale system, and the associated computational efforts therefore dwarfs those made in the SCF optimization.