Skip to content

Changes to collision_terms_std #40

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

Merged
cgiovanetti merged 8 commits into from
Aug 31, 2026
Merged

Changes to collision_terms_std #40

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
aec4c19
Initial commit
graceannperk Jul 14, 2026
bf2fcb1
changes to collision_terms_std
graceannperk Aug 31, 2026
f3667ce
Merge branch 'main' into main
cgiovanetti Aug 31, 2026
e70acc0
Delete .vscode directory
cgiovanetti Aug 31, 2026
cc6ac60
Delete linx/chemical potential testing.ipynb
cgiovanetti Aug 31, 2026
9b61c1b
Delete linx/linx chemical potential testing 4-30.ipynb
cgiovanetti Aug 31, 2026
17ee3cc
Delete graphing 6.26.ipynb
cgiovanetti Aug 31, 2026
20bf5b8
revert to old imports in const
cgiovanetti Aug 31, 2026
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
14 changes: 13 additions & 1 deletion linx/background.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -242,7 +242,7 @@ def dY(self, t, Y, args):

H = thermo.Hubble(rho_EM + rho_nu + rho_extra)

C_rho_nue, C_rho_numu, _, _ = thermo.collision_terms_std(
C_rho_nue, C_rho_numu, C_n_nue, C_n_numu = thermo.collision_terms_std(
T_g, T_nu, T_nu, me=me, decoupled=self.decoupled, use_FD=self.use_FD, collision_me=self.collision_me
)

Expand All @@ -252,4 +252,16 @@ def dY(self, t, Y, args):
dT_g_dt = drho_EM_dt / drho_EM_dT_g
dT_nu_dt = drho_nu_dt / drho_nu_dT_nu

#P_e = thermo.rho_nue_std(T_g, me=me) / 3
#dndt_e = C_n_nue(T_g, T_nu, T_nu, me=me, decoupled=self.decoupled, use_FD=self.use_FD, collision_me=self.collision_me)

#dmue_dt = -(-3.* H *( (thermo.rho_nue_std+P_e) * thermo.dn_nue_dT_nue_std - thermo.n_nue_std * thermo.drho_nue_dT_nue_std) + thermo.dn_nue_dT_nue_std * thermo.drho_nue_dT_nue_std - thermo.drho_nue_dT_nue_std * dndt_e )/(thermo.dn_nue_dmu_nue_std * thermo.drho_nue_dT_nue_std - thermo.dn_nue_dT_nue_std * thermo.drho_nue_dmu_nue_std)


#P_mu = thermo.rho_numt_std(T_nu, me=me) / 3
#dndt_mu = C_n_numu(T_g, T_nu, T_nu, me=me, decoupled=self.decoupled, use_FD=self.use_FD, collision_me=self.collision_me)

#dmunu_dt = -(-3.*H*( (thermo.rho_numt_std+P_mu)* thermo.dn_numt_dT_numt_std - thermo.n_numt.std * thermo.drho_numt_dT_numt_std) + thermo.dn_numt_dT_numt_std * thermo.drho_numt_dT_numt_std - thermo.drho_numt_dT_numt_std * dndt_mu )/(thermo.dn_numt_dmu_numt_std * thermo.drho_numt_dT_numt_std - thermo.dn_numt_dT_numt_std * thermo.drho_numt_dmu_numt_std)

return H, dT_g_dt, dT_nu_dt
#return H, dT_g_dt, dT_nu_dt, dmue_dt, dmunu_dt
2 changes: 1 addition & 1 deletion linx/const.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -124,4 +124,4 @@
# T_Middle to T_switch: small network always sufficient.
T_switch = 1.25e9 * kB # MeV
# T_switch to T_end: turn on full network if desired.
T_end = 6.e7 * kB # MeV
T_end = 6.e7 * kB # MeV
Loading
Loading