CCL: Integration error depending on IA amplitude
I’m seeing an integration error when calculating C_ell for the DES SV n(z) that depends on the value of the IA amplitude, which seems a bit odd.
I wanted to double check if I’m doing something wrong with ranges or input values or something else like that.
I’ve attached an example which uses the current CCL master branch. I had to call it issue.py.txt so that github would let me attach it. It loads the attached number densities.
If you use A=1.0 then it works, but if you use A=3.0 then it fails with
RuntimeError: Error CCL_ERROR_INTEG: ccl_cls.c: ccl_angular_cl(): error integrating over k
Can anyone advise?
About this issue
- Original URL
- State: closed
- Created 7 years ago
- Comments: 16 (2 by maintainers)
@elisachisari : yes, that sounds plausible. If the k->ell integral is all done in one step by first summing the relevant P(k) with their amplitudes and projection kernels, then this would make sense.
A=3 is indeed large for lensing sources, but not large for LRGs (which someone using CCL may want to consider).