lesgourg / class_public

Public repository of the Cosmic Linear Anisotropy Solving System (master for the most recent version of the standard code; GW_CLASS to include Cosmic Gravitational Wave Background anisotropies; classnet branch for acceleration with neutral networks; ExoCLASS branch for exotic energy injection; class_matter branch for FFTlog)
225 stars 282 forks source link

Redefining the CLP #574

Open Kshitizpandey123 opened 6 months ago

Kshitizpandey123 commented 6 months ago

Sir, I want to implement a particular form of dark energy fluid, it has no extra parameter and is just a function of the scale factor (there is no w0_fld or wa_fld in the e.o.s expression). After that, I want to analyze the Planck2018TTTEEE using Montepython.

Here is what I understand so far: Since CLASS uses any one of these: LCDM (Omega_Lambda), Scalar Field (Omega_scf), or CPL (Omega_fld, by default), I need to set Omega_Lambda and Omega_scf to Zero and have to leave Omega_fld unspecified. 1.) My first doubt is: In which file do I set these values of Omega's? In the explanatory.ini? Or where exactly?

Now, knowing that by default the CLASS uses the CLP as the fluid model, I redefined CLP in background.c, set the values of w_fld, its derivative, and integral. 2.) I am not sure if I also have to make some changes in the input.c or not. I am particularly confused in lines 5817 to 5830: / ADDITIONAL SPECIES / Would you like me to change something here? Like setting Omega0_Lambda = Omega0_scf = 0
And, Omega0_fld = 1 - pba->Omega0_k-pba->Omega0_g . . . ? (closure)

Finally, do I need to make any change in the .param file in the montepython? I have already added this line in .param file: data.cosmo_arguments['Omega_Lambda'] = 0

483

Pls help!