Skip to content

Grand neutrino #2

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.

Sign up for GitHub

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

Draft
EBHolm wants to merge 6 commits into
base: master
Choose a base branch
from
Draft

Grand neutrino #2

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
9f9a3af
Implemented collective_ncdm species
EBHolm Dec 19, 2023
83c3056
Bugfixes
EBHolm Dec 19, 2023
0cf184b
Fixed indexing issues
EBHolm Dec 19, 2023
344e82c
Final fix, now works
EBHolm Dec 19, 2023
e3e2c83
Added notebook
EBHolm Dec 19, 2023
8c32de4
Updated notebook
EBHolm Jan 3, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 6 additions & 0 deletions include/background.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -89,6 +89,10 @@ struct background
/* the following parameters help to define the analytical ncdm phase space distributions (p-s-d) */
double * ncdm_psd_parameters; /**< list of parameters for specifying/modifying ncdm p.s.d.'s, to be customized for given model
(could be e.g. mixing angles) */
int output_ncdm_binning;
int collective_ncdm; /**< regards all neutrinos as having the same distribution function,
thus solving only one hierarcy */
int collective_ncdm_N; /**< amount of collective ncdm species */
double * M_ncdm; /**< vector of masses of non-cold relic: dimensionless ratios m_ncdm/T_ncdm */
double * m_ncdm_in_eV; /**< list of ncdm masses in eV (inferred from M_ncdm and other parameters above) */
double * Omega0_ncdm, Omega0_ncdm_tot; /**< Omega0_ncdm for each species and for the total Omega0_ncdm */
Expand Down Expand Up @@ -189,6 +193,8 @@ struct background
int index_bg_rho_ncdm1; /**< density of first ncdm species (others contiguous) */
int index_bg_p_ncdm1; /**< pressure of first ncdm species (others contiguous) */
int index_bg_pseudo_p_ncdm1;/**< another statistical momentum useful in ncdma approximation */
int index_bg_q_ncdm1;
int index_bg_w_ncdm1;

int index_bg_rho_tot; /**< Total density */
int index_bg_p_tot; /**< Total pressure */
Expand Down
Loading