Commit 4e4795ed authored by Julia Stadler's avatar Julia Stadler
Browse files

functioning version with nuDM interactions

parent 8e4a910d
......@@ -51,6 +51,7 @@
-> the following commands are executed only if pth->has_coupling_nuDM is TRUE
line 6983: S_nuDM = 4./3. *pvecback[pba->index_bg_rho_ur]/pvecback[pba->index_bg_rho_cdm];
line 7238: dy[pv->index_pt_theta_cdm] += S_nuDM*pvecthermo[pth->index_th_dmu_nuDM]*(y[pv->index_pt_theta_ur]-Y[pv->index_pt_theta_cdm]);
-> additional condition ppw->approx[ppw->index_ap_rsa] == (int)rsa_off
line 7389: dy[pv->index_pt_theta_ur] += pvecthermo[pth->index_th_dmu_nuDM]*(y[pv->index_pt_theta_cdm]-y[pv->index_pt_theta_ur]);
line 7401: dy[pv->index_pt_shear_ur] += -9./10.*pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_shear_ur];
line 7408: dy[pv->index_pt_l3_ur] -= pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_l3_ur];
......
......@@ -47,7 +47,7 @@ write thermodynamics = no
input_verbose = 1
background_verbose = 1
thermodynamics_verbose = 1
perturbations_verbose = 10
perturbations_verbose = 1
transfer_verbose = 1
primordial_verbose = 1
spectra_verbose = 1
......
......@@ -7234,7 +7234,7 @@ int perturb_derivs(double tau,
dy[pv->index_pt_delta_cdm] = -(y[pv->index_pt_theta_cdm]+metric_continuity); /* cdm density */
dy[pv->index_pt_theta_cdm] = - a_prime_over_a*y[pv->index_pt_theta_cdm] + metric_euler; /* cdm velocity */
if(pth->has_coupling_nuDM==_TRUE_)
if(pth->has_coupling_nuDM==_TRUE_ && ppw->approx[ppw->index_ap_rsa] == (int)rsa_off)
dy[pv->index_pt_theta_cdm] += S_nuDM*pvecthermo[pth->index_th_dmu_nuDM]*(y[pv->index_pt_theta_ur]-y[pv->index_pt_theta_cdm]);
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment