Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Cobaya CCL connector: increase a_bg range to support sampling w0 and wa #270

Open
heatherprince opened this issue May 15, 2023 · 1 comment

Comments

@heatherprince
Copy link

When varying w0 and wa I got an error using the default minimum scale factor of 0.1, or even a smaller scale factor down to 5e-4 (pyccl.errors.CCLError: Comoving distance spline does not cover sufficiently high redshifts for HALOFIT. HALOFIT translates (w0, wa) to a w0_eff. This requires computing the comoving distance to the CMB, which requires the background splines being sampled to sufficiently high redshifts. If using the calculator mode, check the support of the background data.)

If I made the minimum scale factor too small, e.g. 1e-4, I got a different error (pyccl.errors.CCLError: Error CCL_ERROR_SPLINE: ccl_background.c: ccl_cosmology_compute_growth(): Error creating D(a) spline))

It ran without errors for a minimum a_bg of 2e-4, (zmax ~5000), so fixing it might be as simple as just changing the default a_bg range, but given the weird amount of tuning I had to do to get it to run without errors I'm not sure if the theory power spectra are reliable, and it would probably be good to check this with someone more familiar with ccl.

@vitenti
Copy link
Collaborator

vitenti commented Aug 28, 2023

Hello @damonge , when creating a pyccl cosmology using the calculator mode, we pass the arrays computed in the user defined intervals.

  • Is there a good way to probe pyccl to have the "good" intervals?
  • Or, maybe, is there a good way to report the issue back to the user with a good error message?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants