CosmoCoffee

Posted: **January 01 2016**

Pardon me if this has been addressed before. I'm using January 2015 version of CAMB, and wanted to use the 4th column (photon density perturbation) of transfer function outputs at various redshifts. I set
transfer_kmax = 20
transfer_k_per_logint = 100
in the parameter file.

With redshifts selected as in the figure file(1000, 203, 101), the transfer function has discontinuous points at z=203 and 101 which looks suspicious. I ran it both with intel-fortran compiled and gfortran-compiled ones with the same result.

As a test I went back to September 2013 version, and the same outputs for the same redshifts have null values for the transfer function at k>~0.02, so it used to be worse.

Is it a glitch from error accumulation of the ODE solver used?

[/b]

Posted: **January 04 2016**

By default it doesn't calculate the radiation multipoles accurately at late times, the glitches may be the sudden switch to the approximate solution. You can turn off DoLateRadTruncation and increase lAccuracyBoost if you want to get it more accurately, e.g. (in latest CAMB with python wrapper)

Code: Select all

options=&#91;&#123;&#125;,&#123;'DoLateRadTruncation'&#58;False&#125;,&#123;'DoLateRadTruncation'&#58;False, 'lAccuracyBoost'&#58;4&#125;&#93;
_, axs = plt.subplots&#40;1,3, figsize=&#40;12,4&#41;&#41;
for option, ax in zip&#40;options, axs&#41;&#58;
    pars = camb.CAMBparams&#40;&#41;
    pars.set_cosmology&#40;H0=67.5, ombh2=0.022, omch2=0.122&#41;
    pars.InitPower.set_params&#40;ns=0.965&#41;
    pars.set_matter_power&#40;redshifts=&#91;1000, 203, 10&#93;, kmax=0.2, k_per_logint=100&#41;
    pars.WantCls = False
    pars.NonLinear = model.NonLinear_none

    pars.set_accuracy&#40;**option&#41;

    results = camb.get_results&#40;pars&#41;

    M = results.get_matter_transfer_data&#40;&#41;
    for i in range&#40;3&#41;&#58;
        tz = M.transfer_z&#40;'delta_photon',i&#41;
        ax.plot&#40;M.q, tz*M.q**2&#41;
    ax.set_xlim&#40;&#91;0,0.1&#93;&#41;

gives

Posted: **January 05 2016**

Thank you so much, Antony. Not having the proper gfortran version for python wrapper, I just followed your instruction inside the usual input file. Now the figure is what I have with
do_late_rad_truncation = F
l_accuracy_boost = 3
which looks natural, and I'll test convergence of l_accuracy_boost later.

CosmoCoffee

CAMB: transfer function for gamma at z~100

CAMB: transfer function for gamma at z~100

Re: CAMB: transfer function for gamma at z~100

CAMB: transfer function for gamma at z~100