transfer.NoWiggleEisensteinHu

Open adematti opened this issue 6 years ago • 1 comments

Hi,

I believe there is a mistake in gamma_eff l. 248 of cosmology/power/transfers.py: gamma_eff = self.Omh2 * (self.alpha_gamma + (1 - self.alpha_gamma) / (1 + (0.43*ks) ** 4))

One should use Omega_m * h = self.Omh2 / h here, as stated eq. 30 of https://arxiv.org/pdf/astro-ph/9709112.pdf.

Cheers

May 20 '19 15:05 adematti

Thanks for reporting this.

Could this be related to the units, e.g. of Mpc vs Mpc/h?

@nickhand

May 22 '19 03:05 rainwoodman