I would pass frequencies[i] to CrossSpectralDensity init and call np.atleast_1d in the CrossSpectralDensity init.