Bug in sigmasqr function in power.py

[ethlau]

The followling line in sigmasqr function in power.py,

y = romb(int_sigma, np.log10(kmin), np.log10(kmax), divmax=7)

Shouldn't the np.log10 function be np.log instead, since the integrand function int_sigma takes log of k, not log10 of k?

def int_sigma(logk):
    k = np.exp(logk)
    x = k * R
    w = 3.0 * (np.sin(x) - x * np.cos(x)) / (x * x * x)
    pk = transfer_fn(cosmo, k, **kwargs) ** 2 * primordial_matter_power(cosmo, k)
    return k * (k * w) ** 2 * pk

[EiffL]

Hi @ethlau, no I don't think that causes any bugs np.log10(kmin) and np.log10(kmax) are just rough integration bounds, they shouldn't have an impact on the result of the computation. Up to a redefinition of this kmin and kmax values, it shouldn't affect the result.