Page 1 of 1

Cls numerical accuracy in CosmoMC

Posted: March 17 2016
by Alex Zucca
Hi Antony,

I'm working on primordial magnetic fields and I have a CAMB code that works fine: compensated and passive modes are computed correctly. However as soon as I try to compute them in cosmomc, I face some issues. In particular the scalar compensated power spectrum appears to be full of spikes (I used the same parameters in CAMB and it produces a smooth power spectrum). I was wondering if cosmomc passes a different numerical accuracy to the camb subroutines (either CAMB_GetTransfers or CAMB_TransfersToPowers). If so, where are these parameters set?

Many thanks,
Alex Zucca

Re: Cls numerical accuracy in CosmoMC

Posted: March 17 2016
by Antony Lewis
Accuracy settings will vary depending on what data and settings you are using. I would write out the relevant things in CAMB and see how they differ when run with cosmomc.

Cls numerical accuracy in CosmoMC

Posted: March 21 2016
by Alex Zucca
Hi Antony,
thanks for your help. I managed to check all the parameters and the accuracy settings. In particular I checked the computation until the iCl_scalar array (before interpolation) and that is correct. In some way the interpolation process when called in CosmoMC is having some problems, but I cannot understand why.

Alex