-
Notifications
You must be signed in to change notification settings - Fork 0
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
Latest class #11
base: master
Are you sure you want to change the base?
Latest class #11
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.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@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.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
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.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@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.
No description provided.