Merge pull request #8 from emiliobellini/hi_class

Corrected bug when integrating M2 (thanks to Bill Wright)
miguelzuma · May 23, 2020 · 16ae0f6 · 16ae0f6
2 parents 78fcc86 + 17bc625
commit 16ae0f6
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/source/background.c b/source/background.c
@@ -3262,7 +3262,7 @@ int background_derivs(
     }
     /** - Planck mass equation (if parameterization in terms of alpha_m **/
     if (pba->M_pl_evolution_smg == _TRUE_)
-      dy[pba->index_bi_delta_M_pl_smg] = y[pba->index_bi_a]*pvecback[pba->index_bg_H]*pvecback[pba->index_bg_mpl_running_smg]*y[pba->index_bi_delta_M_pl_smg];   //in this case the running has to be integrated (eq 3.3 of 1404.3713 yields M2' = aH\alpha_M)
+      dy[pba->index_bi_delta_M_pl_smg] = y[pba->index_bi_a]*pvecback[pba->index_bg_H]*pvecback[pba->index_bg_mpl_running_smg]*(y[pba->index_bi_delta_M_pl_smg]+1.);   //in this case the running has to be integrated (eq 3.3 of 1404.3713 yields M2' = aH\alpha_M)
   }
 
   if (pba->has_fld == _TRUE_) {