In a resonant calculation, four sets of linear equations are solved and their respective Again, the computation of the property building blocks takes up a large part of the total wall time, amounting to 50% of the total wall time in the presented case. This is 20 percentage points less than observed for the quadratic response function and is a reflection of the substantial amount of additional linear equations that are required for the state-by-state approach, besides, of course, the lower computational cost and scaling of the densities.