It is clearly seen that the most time-consuming steps are the computation of the densities and matrix-vector contractions for building the quadratic response function: more than 70% the total wall time was spent by those operations in both the real and complex cases. From the timings it is observed that the RHSs and matrix-vector contractions take roughly twice as long for the complex response function, while the left-and right-hand transformations are faster when the frequency is complex. We note that the linear solver required nearly twice as many iterations for converging the response equations without a damping factor.