Hello,
I need to calculate the linear matter power spectrum for very small scales;
the output from camb (in output_matterpower.dat) is the same as the one in cosmomc
(given by the variable mpk_lin(i)=MatterPowerAt(Theory,k_scaled(i))/a_scl**3 in mpk.f90),
but they start being different on very small scales, starting from k~3, even though the cosmology used is the same. Any ideas? Maybe the variable "mpk_lin" is not the correct one to look at?
Thanks
p
camb vs cosmomc matter power spectrum
-
- Posts: 4
- Joined: February 26 2008
- Affiliation: University of California, Irvine
-
- Posts: 1943
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: camb vs cosmomc matter power spectrum
Possibly a k range issue - see the settings at the top of cmbtypes.f90 in cosmomc/source. cosmomc will do log-linear extrapolation if out of the computed range, which is obviously only OK if in small tails of window functions.