We present a reduced-space algorithm for solving the complex (damped) linear response equations
required to compute the complex linear response function for the hierarchy of methods: coupled cluster
singles, coupled cluster singles and iterative approximate doubles, and coupled cluster singles and
doubles. The solver is the keystone element for the development of damped coupled cluster response
methods for linear and nonlinear effects in resonant frequency regions