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

Gauss-Laguerre NCDM quadrature strategy results in segmentation fault #597

Open
robbieellis12 opened this issue Sep 27, 2024 · 0 comments
Open

Comments

@robbieellis12
Copy link

I am using:
N_ur=0
N_ncdm =3
m_ncdm = 0.04, 0.04, 0.04
T_ncdm = 0.71611, 0.71611, 0.71611
ncdm_quadrature_strategy =1, 1, 1
ncdm_N_momentum_bins =200,200, 200
(everything else matches explanatory.ini)

And I get:

"Reading input parameters
-> matched budget equations by adjusting Omega_Lambda = nan
Running CLASS version v3.2.3
Computing background
Omega0: nan: -> non-cold dark matter species with i=1 has m_i = 4.000000e-02 eV (so m_i / omega_i =nan eV)
-> ncdm species i=1 sampled with 200 (resp. 200) points for purpose of background (resp. perturbation) integration. In the relativistic limit it gives Delta N_eff = nan
..."

I tried the obvious stuff--I changed N_ur and added Omega_ncdm =0.000934,0.000934,0.000934 but I still get:

" -> matched budget equations by adjusting Omega_Lambda = 0.687261
Running CLASS version v3.2.3
Computing background
Omega0: 0.000934: -> non-cold dark matter species with i=1 has m_i = 4.000000e-02 eV (so m_i / omega_i =nan eV)
-> ncdm species i=1 sampled with 200 (resp. 200) points for purpose of background (resp. perturbation) integration. In the relativistic limit it gives Delta N_eff = nan"

I scoured input.c and background.c--I honestly can't tell where omega0_ncdm is even set when using the other quadrature methods (when it is being inferred from m_ncdm). I looked for every instance of omega0_ncdm in these modules and couldn't find anywhere it was being set. I don't get any error when using the trapezoidal methods (2/3). Am I missing something obvious or is there an issue with the code? Any way to resolve this?

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

1 participant