CLASS : #define _Z_PK_NUM_MAX_ modified not taken into account at execution

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Fabien Dournac
Posts: 20
Joined: March 27 2010
Affiliation: CERFACS
Contact:

CLASS : #define _Z_PK_NUM_MAX_ modified not taken into account at execution

Post by Fabien Dournac » October 06 2023

Hello,

I have modified in include/output.h the variable #define Z_PK_NUM_MAX 400 since I have 397 redshifts to compute for Pk. I am using current CLASS github version.

Compiles fine but at execution, I have only 57 output files.

Here my params.ini :

output=mPk, tCl, lCl
modes=s
ic=ad
tensor_method=
gauge=synchronous
h=0.67
YHe=BBN
recombination=HyRec
recfast_photoion_dependence=
reio_parametrization=reio_camb
z_reio=7.6711
reionization_exponent=1.5
reionization_width=0.5
helium_fullreio_redshift=3.5
helium_fullreio_width=0.5
binned_reio_num=3
binned_reio_z=8,12,16
binned_reio_xe=0.8,0.2,0.1
binned_reio_step_sharpness=0.3
many_tanh_num=2
many_tanh_z=3.5,11.3
many_tanh_xe=-2,-1
many_tanh_width=0.5
reio_inter_num=8
reio_inter_z=0, 3, 4, 8, 9, 10, 11, 12
reio_inter_xe=-2, -2, -1, -1, 0.9, 0.5, 0.1, 0
compute_damping_scale=
varying_fundamental_constants=none
varying_transition_redshift=
varying_alpha=1.
varying_me=1.
bbn_alpha_sensitivity=1.
T_cmb=2.7255
N_ur=3.044
N_ncdm=
ncdm_psd_filenames=psd_FD_single.dat
m_ncdm=0.06
T_ncdm=
ksi_ncdm=
deg_ncdm=
ncdm_quadrature_strategy=
ncdm_maximum_q=
ncdm_N_momentum_bins=
Omega_k=0.
Omega_dcdmdr=0.0
Gamma_dcdm=0.0
tau_dcdm=0.0
xi_idr=
stat_f_idr=0.875
Omega_idm_dr=0.
m_idm=1.0e11
a_idm_dr=0.
nindex_idm_dr=
idr_nature=
b_idr=
alpha_idm_dr=1.5
beta_idr=1.5
Omega_fld=0
Omega_scf=0
use_ppf=yes
c_gamma_over_c_fld=0.4
fluid_equation_of_state=EDE
scf_parameters=10.0, 0.0, 0.0, 0.0, 100.0, 0.0
attractor_ic_scf=yes
scf_tuning_index=0
scf_shooting_parameter=
DM_annihilation_efficiency=0.
DM_annihilation_variation=0.
DM_annihilation_z=1000
DM_annihilation_zmax=2500
DM_annihilation_zmin=30
DM_annihilation_f_halo=0
DM_annihilation_z_halo=8
DM_decay_fraction=0.
DM_decay_Gamma=0.
PBH_evaporation_fraction=0.
PBH_evaporation_mass=0.
PBH_accretion_fraction=0.
PBH_accretion_mass=0.
PBH_accretion_recipe=disk_accretion
PBH_accretion_relative_velocities=-1.
PBH_accretion_ADAF_delta=1.e-3
PBH_accretion_eigenvalue=0.1
f_eff_type=on_the_spot
chi_type=CK_2004
non_linear=none
nonlinear_min_k_max=200
pk_eq=
feedback model=
eta_0=
c_min=
Pk_ini_type=analytic_Pk
k_pivot=0.05
A_s=2.1260500000000005e-09
n_s=0.96
alpha_s=0.
f_bi=1.
n_bi=1.5
f_cdi=1.
f_nid=1.
n_nid=2.
alpha_nid=0.01
c_ad_bi=0.5
c_ad_cdi=-1.
c_bi_nid=1.
r=1.
n_t=scc
alpha_t=scc
potential=polynomial
V_0=1.e-13
V_1=-1.e-14
V_2=7.e-14
V_3=
V_4=
H_0=1.e-13
H_1=-1.e-14
H_2=7.e-14
H_3=
H_4=
phi_end=
full_potential=polynomial
Vparam0=
Vparam1=
Vparam2=
Vparam3=
Vparam4=
k1=0.002
k2=0.1
P_{RR}^1=2.3e-9
P_{RR}^2=2.3e-9
P_{II}^1=1.e-11
P_{II}^2=1.e-11
P_{RI}^1=-1.e-13
|P_{RI}^2|=1.e-13
special_iso=
command=cat external/external_Pk/Pk_example.dat
custom1=0.05 # In the example command: k_pivot
custom2=2.215e-9 # In the example command: A_s
custom3=0.9624 # In the example command: n_s
custom4=2e-10 # In the example (with tensors) command: A_t
custom5=-0.1 # In the example (with tensors) command: n_t
l_max_scalars=2500
l_max_tensors=500
selection=gaussian
selection_mean=0.98,0.99,1.0,1.1,1.2
selection_width=0.1
selection_bias=
selection_magnification_bias=
non_diagonal=4
dNdz_selection=
dNdz_evolution=
P_k_max_h/Mpc=1.
z_pk=3.731,3.721580808080808,3.712161616161616,3.702742424242424,3.693323232323232,3.68390404040404,3.6744848484848482,3.6650656565656563,3.6556464646464644,3.6462272727272724,3.6368080808080805,3.627388888888889,3.617969696969697,3.608550505050505,3.599131313131313,3.5897121212121212,3.5802929292929293,3.5708737373737374,3.5614545454545454,3.5520353535353535,3.5426161616161616,3.5331969696969696,3.5237777777777777,3.5143585858585857,3.504939393939394,3.495520202020202,3.48610101010101,3.476681818181818,3.467262626262626,3.457843434343434,3.448424242424242,3.4390050505050502,3.4295858585858583,3.420166666666667,3.410747474747475,3.4013282828282825,3.391909090909091,3.382489898989899,3.373070707070707,3.363651515151515,3.3542323232323232,3.3448131313131313,3.3353939393939394,3.3259747474747474,3.3165555555555555,3.3071363636363635,3.2977171717171716,3.2882979797979797,3.2788787878787877,3.269459595959596,3.260040404040404,3.250621212121212,3.24120202020202,3.231782828282828,3.222363636363636,3.212944444444444,3.203525252525252,3.1941060606060603,3.1846868686868683,3.1752676767676764,3.165848484848485,3.156429292929293,3.147010101010101,3.137590909090909,3.128171717171717,3.118752525252525,3.1093333333333333,3.0999141414141413,3.0904949494949494,3.0810757575757575,3.0716565656565655,3.0622373737373736,3.0528181818181817,3.0433989898989897,3.0339797979797978,3.024560606060606,3.015141414141414,3.005722222222222,2.99630303030303,2.986883838383838,2.977464646464646,2.968045454545454,2.9586262626262627,2.9492070707070708,2.939787878787879,2.930368686868687,2.920949494949495,2.911530303030303,2.902111111111111,2.892691919191919,2.883272727272727,2.8738535353535353,2.8644343434343433,2.8550151515151514,2.8455959595959595,2.8361767676767675,2.8267575757575756,2.8173383838383836,2.8079191919191917,2.7984999999999998,2.789080808080808,2.779661616161616,2.770242424242424,2.760823232323232,2.75140404040404,2.741984848484848,2.732565656565656,2.7231464646464647,2.7137272727272728,2.704308080808081,2.694888888888889,2.685469696969697,2.676050505050505,2.666631313131313,2.657212121212121,2.647792929292929,2.6383737373737373,2.6289545454545453,2.6195353535353534,2.6101161616161614,2.6006969696969695,2.5912777777777776,2.5818585858585856,2.5724393939393937,2.5630202020202018,2.55360101010101,2.544181818181818,2.534762626262626,2.525343434343434,2.5159242424242425,2.5065050505050506,2.4970858585858586,2.4876666666666667,2.4782474747474748,2.468828282828283,2.459409090909091,2.449989898989899,2.440570707070707,2.431151515151515,2.421732323232323,2.412313131313131,2.4028939393939392,2.3934747474747473,2.3840555555555554,2.3746363636363634,2.3652171717171715,2.3557979797979796,2.3463787878787876,2.3369595959595957,2.3275404040404037,2.318121212121212,2.30870202020202,2.2992828282828284,2.2898636363636364,2.2804444444444445,2.2710252525252526,2.2616060606060606,2.2521868686868687,2.2427676767676767,2.233348484848485,2.223929292929293,2.214510101010101,2.205090909090909,2.195671717171717,2.186252525252525,2.176833333333333,2.1674141414141412,2.1579949494949493,2.1485757575757574,2.1391565656565654,2.1297373737373735,2.1203181818181815,2.1108989898989896,2.1014797979797977,2.092060606060606,2.0826414141414142,2.0732222222222223,2.0638030303030304,2.0543838383838384,2.0449646464646465,2.0355454545454545,2.0261262626262626,2.0167070707070707,2.0072878787878787,1.9978686868686868,1.9884494949494949,1.979030303030303,1.969611111111111,1.960191919191919,1.950772727272727,1.9413535353535352,1.9319343434343432,1.9225151515151513,1.9130959595959593,1.9036767676767674,1.8942575757575757,1.8848383838383838,1.8754191919191918,1.8659999999999999,1.856580808080808,1.847161616161616,1.8377424242424243,1.8283232323232324,1.8189040404040404,1.8094848484848485,1.8000656565656565,1.7906464646464646,1.7812272727272727,1.7718080808080807,1.7623888888888888,1.7529696969696968,1.743550505050505,1.734131313131313,1.7247121212121213,1.715292929292929,1.7058737373737372,1.6964545454545454,1.6870353535353535,1.6776161616161616,1.6681969696969696,1.6587777777777777,1.6493585858585857,1.6399393939393938,1.6305202020202019,1.62110101010101,1.611681818181818,1.6022626262626263,1.5928434343434343,1.5834242424242424,1.5740050505050505,1.5645858585858585,1.5551666666666666,1.5457474747474746,1.5363282828282827,1.5269090909090908,1.5174898989898988,1.508070707070707,1.498651515151515,1.489232323232323,1.4798131313131313,1.4703939393939394,1.4609747474747474,1.4515555555555555,1.4421363636363635,1.4327171717171716,1.4232979797979797,1.4138787878787877,1.4044595959595958,1.3950404040404039,1.3856212121212121,1.3762020202020202,1.366782828282828,1.3573636363636363,1.3479444444444444,1.3385252525252525,1.3291060606060605,1.3196868686868686,1.3102676767676766,1.3008484848484847,1.2914292929292928,1.2820101010101008,1.2725909090909089,1.263171717171717,1.253752525252525,1.2443333333333333,1.2349141414141414,1.2254949494949494,1.2160757575757575,1.2066565656565655,1.1972373737373736,1.1878181818181817,1.1783989898989897,1.1689797979797978,1.1595606060606058,1.1501414141414141,1.1407222222222222,1.1313030303030303,1.1218838383838383,1.1124646464646464,1.1030454545454544,1.0936262626262625,1.0842070707070706,1.0747878787878786,1.0653686868686867,1.0559494949494947,1.0465303030303028,1.0371111111111109,1.0276919191919192,1.0182727272727272,1.0088535353535353,0.9994343434343433,0.9900151515151514,0.9805959595959596,0.9711767676767676,0.9617575757575758,0.9523383838383839,0.9429191919191919,0.9335,0.9240808080808081,0.9146616161616162,0.9052424242424243,0.8958232323232324,0.8864040404040404,0.8769848484848485,0.8675656565656565,0.8581464646464647,0.8487272727272728,0.8393080808080807,0.8298888888888889,0.820469696969697,0.811050505050505,0.8016313131313131,0.7922121212121213,0.7827929292929293,0.7733737373737374,0.7639545454545454,0.7545353535353535,0.7451161616161616,0.7356969696969697,0.7262777777777778,0.7168585858585859,0.7074393939393939,0.698020202020202,0.6886010101010102,0.6791818181818182,0.6697626262626263,0.6603434343434343,0.6509242424242424,0.6415050505050505,0.6320858585858586,0.6226666666666667,0.6132474747474748,0.6038282828282828,0.5944090909090909,0.584989898989899,0.575570707070707,0.5661515151515152,0.5567323232323232,0.5473131313131313,0.5378939393939394,0.5284747474747474,0.5190555555555555,0.5096363636363637,0.5002171717171717,0.4907979797979798,0.48137878787878785,0.47195959595959597,0.4625404040404041,0.45312121212121215,0.44370202020202026,0.43428282828282827,0.42486363636363633,0.41544444444444445,0.4060252525252525,0.39660606060606063,0.3871868686868687,0.37776767676767675,0.3683484848484848,0.35892929292929293,0.349510101010101,0.3400909090909091,0.33067171717171717,0.32125252525252523,0.3118333333333333,0.3024141414141414,0.2929949494949495,0.2835757575757576,0.27415656565656565,0.2647373737373737,0.2553181818181818,0.2458989898989899,0.23647979797979798,0.22706060606060607,0.21764141414141414,0.20822222222222223,0.1988030303030303,0.18938383838383838,0.17996464646464644,0.17054545454545456,0.16112626262626262,0.1517070707070707,0.1422878787878788,0.13286868686868686,0.12344949494949495,0.11403030303030304,0.10461111111111111,0.09519191919191919,0.08577272727272728,0.07635353535353535,0.06693434343434344,0.05751515151515152,0.048095959595959595,0.03867676767676768,0.02925757575757576,0.01983838383838384,0.01041919191919192,0.001,0.0
z_max_pk=10.
lensing=no
A_L=
sd_branching_approx=exact
sd_PCA_size=2
sd_detector_name=PIXIE
sd_only_exotic=no
sd_include_g_distortion=no
sd_add_y=0.
sd_add_mu=0.
include_SZ_effect=no
overwrite_root=no
headers=yes
format=class
write_background=no
write_thermodynamics=no
write_primordial=no
write_exotic_injection=no
write_distortions=no
write_parameters=yes
write_warnings=no
input_verbose=1
background_verbose=1
thermodynamics_verbose=1
perturbations_verbose=1
transfer_verbose=1
primordial_verbose=1
harmonic_verbose=1
fourier_verbose=1
lensing_verbose=1
distortions_verbose=1
output_verbose=1
Omega_b=0.05
Omega_m=0.32
output_root=test
Omega_EDE=0.68

After execution, I get :

ls
output
params00_cl.dat
params00_parameters.ini
params00_unused_parameters
params00_z10_pk.dat
params00_z11_pk.dat
params00_z12_pk.dat
params00_z13_pk.dat
params00_z14_pk.dat
params00_z15_pk.dat
params00_z16_pk.dat
params00_z17_pk.dat
params00_z18_pk.dat
params00_z19_pk.dat
params00_z1_pk.dat
params00_z20_pk.dat
params00_z21_pk.dat
params00_z22_pk.dat
params00_z23_pk.dat
params00_z24_pk.dat
params00_z25_pk.dat
params00_z26_pk.dat
params00_z27_pk.dat
params00_z28_pk.dat
params00_z29_pk.dat
params00_z2_pk.dat
params00_z30_pk.dat
params00_z31_pk.dat
params00_z32_pk.dat
params00_z33_pk.dat
params00_z34_pk.dat
params00_z35_pk.dat
params00_z36_pk.dat
params00_z37_pk.dat
params00_z38_pk.dat
params00_z39_pk.dat
params00_z3_pk.dat
params00_z40_pk.dat
params00_z41_pk.dat
params00_z42_pk.dat
params00_z43_pk.dat
params00_z44_pk.dat
params00_z45_pk.dat
params00_z46_pk.dat
params00_z47_pk.dat
params00_z48_pk.dat
params00_z49_pk.dat
params00_z4_pk.dat
params00_z50_pk.dat
params00_z51_pk.dat
params00_z52_pk.dat
params00_z53_pk.dat
params00_z54_pk.dat
params00_z55_pk.dat
params00_z56_pk.dat
params00_z5_pk.dat
params00_z6_pk.dat
params00_z7_pk.dat
params00_z8_pk.dat
params00_z9_pk.dat

Why the modification of _Z_PK_NUM_MAX is not taken into account, I expect 397 output files. Any clue is welcome.

Regards

Post Reply