Conversation
… through class 2.9.3
theory/cosmo3D.c
Outdated
| @@ -529,11 +530,13 @@ double get_class_s8(struct file_content *fc, int *status){ | |||
| sprintf(fc->value[position_kmax],"%e",10.); | |||
| } | |||
| *status = run_class(fc,&ba,&th,&pt,&tr,&pm,&sp,&nl,&le); | |||
| double sigma8 = *nl.sigma8; | |||
There was a problem hiding this comment.
@elikrause @paulrogozenski,
Just let you know. I find that *nl.sigma8 gives you sigma8(cdm+baryon).
To get sigma8 (total matter), you should do nl.sigma8[nl.index_pk_total]
I have checked this by using python version of class (classy)
There was a problem hiding this comment.
In case you want to check yourself, here is an example
params = [
omega_m = 3.548971e-01
sigma_8 = 0.000000e+00
A_s = 2.028418e-09
n_s = 8.710132e-01
w0 = -1.000000e+00
wa = 0.000000e+00
omega_b = 5.310769e-02
omega_nuh2 = 4.166853e-03
h0 = 6.233805e-01
MGSigma = 0.000000e+00
MGmu = 0.000000e+00
]
From classy,
sigm8(total matter) = 0.6804764082044734
sigma8(cdm+b) = 0.692361342914594
There was a problem hiding this comment.
@paulrogozenski would you add this in your next commit?
…ions to use the total matter power spectra. Will look at implementing P_cluster for HOD contributions next
…rino modelling to change neff_standard to 3.044 with corresponding massive neutrino cases.
…rino modelling to change neff_standard to 3.044 with corresponding massive neutrino cases.
…rino modelling to change neff_standard to 3.044 with corresponding massive neutrino cases.
3 participants
No description provided.