double PowerSpec_DM_2ndSpecies(double k) { /* at the moment, we simply call the Eistenstein & Hu spectrum * for the second DM species, but this could be replaced with * something more physical, say for neutrinos */ double power; power = Norm * k * pow(tk_eh(k), 2); return power; }
virtual double power(double k, int Type) { return k * pow(tk_eh(k), 2); }
double PowerSpec_EH(double k) /* Eisenstein & Hu */ { return Norm * k * pow(tk_eh(k), 2); }