Skip to content

Error initializing CLASS-PT for certain values of N_ur #13

@p-espenshade

Description

@p-espenshade

Hello, I am using the Python wrapper for CLASS-PT. I noticed that the following Python script works fine in initializing CLASS-PT:

from classy import Class

CPT = Class()
cosmology = {'h': 0.6736, 
      'n_s': 0.9649, 
      'omega_b': 0.02237, 
      'omega_cdm': 0.12, 
      'tau_reio': 0.0544, 
      'A_s': 2.0989031673191437e-09, 
      'z_pk': 1.5,
      'N_ur': 3.044, 
      'N_ncdm': 0
}
CPT.set(cosmology)
CPT.set({'output':'mPk', 
        'non linear':'PT',
        'IR resummation':'Yes', 
        'Bias tracers':'Yes'
})
CPT.compute()

However, when I change N_ur from 3.044 to 3.046 I get an error:

File "classy.pyx", line 393, in classy.Class.compute
classy.CosmoComputationError: 

Error in Class: spectra_init(L:5901) :error in spectra_pk(pba,ppt,ppm,pnlpt,pnl,psp);
=>spectra_pk(L:8564) :error in spectra_sigma(pba,ppm,psp,8./pba->h,0.,&(psp->sigma8));
=>spectra_sigma(L:9866) :error in spectra_pk_at_k_and_z(pba,ppm,psp,k,z,&pk,pk_ic);
=>spectra_pk_at_k_and_z(L:683) :error in spectra_pk_at_z(pba, psp, logarithmic, z, spectrum_at_z, spectrum_at_z_ic);
=>spectra_pk_at_z(L:408) :error in array_interpolate_spline(psp->ln_tau, psp->ln_tau_size, psp->ln_pk, psp->ddln_pk, psp->ln_k_size, ln_tau, &last_index, output_tot, psp->ln_k_size, psp->error_message);
=>array_interpolate_spline(L:1706) : x=9.559201e+00 > x_max=9.559201e+00

The reason I want to use N_ur = 3.046 is to compare the output with older datasets I have from when N_ur = 3.046 was the default value in CLASS.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions