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);
}
