The generalization of the implementation for the damped quadratic response function entailed generalizing the construction of the RHS for the response multiplier equations, the G-matrix contraction routines to handle two complex vectors, the B-matrix contraction routines to handle complex vectors from the left, the D F density to handle complex amplitudes, and the D A density to handle complex multipliers. The generalization of the latter also involved the on-the-fly evaluation of the doubles parts of the complex response vectors.