  In a recent publication [J. Chem. Phys. 142, 034115 (2015)] we have derived a
hierarchy of coupled differential equations in time domain to calculate the
linear optical properties of molecular aggregates. Here we provide details
about issues concerning the numerical implementation. In addition we present
the corresponding hierarchy in frequency domain.
