Hi,
I have noticed that the method growth_rate() in the class GrowthFactor, which claims to compute the linear growth rate $f$ following Eq. (4) of Hamilton (2001), is actually confusing $g(z)$ and $D(z)$.
In particular, in Eq. (4) of Hamilton (2001), the last term on the r.h.s. is divided by the $g(z)$. In growth_rate(), the return is
return (
-1
- self.cosmo.Om(z) / 2
+ self.cosmo.Ode(z)
+ 5 * self.cosmo.Om(z) / (2 * self.growth_factor(z))
)
with the growth_factor(z) actually returns $D(z)=g(z)/(1.0+z)$ by definition
return self._d_plus(z) / self._d_plus0
I do have further suggestions for these growth_rate() and growth_factor() methods, but I think this could be a quick fix.
Hi,
I have noticed that the method$f$ following Eq. (4) of Hamilton (2001), is actually confusing $g(z)$ and $D(z)$ .
growth_rate()in theclass GrowthFactor, which claims to compute the linear growth rateIn particular, in Eq. (4) of Hamilton (2001), the last term on the r.h.s. is divided by the$g(z)$ . In
growth_rate(), thereturniswith the$D(z)=g(z)/(1.0+z)$ by definition
growth_factor(z)actually returnsI do have further suggestions for these
growth_rate()andgrowth_factor()methods, but I think this could be a quick fix.