The linear density-density response function χ(r,r′) can be transformed from its Eigensystem representation into a computationally more efficient moment expansion representation using a suitable unitary transformation. Here, we propose an iterative approach for the direct calculation of this moment representation without resorting to either the direct-space or the conventional Eigensystem representation