From 17bc625e2d776589b2d4ba52d14f6bdfb56c9de7 Mon Sep 17 00:00:00 2001 From: Emilio Bellini Date: Tue, 12 May 2020 15:44:55 +0100 Subject: [PATCH] Corrected bug when integrating M2 --- source/background.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/background.c b/source/background.c index bd5dda03..b15d8c8e 100755 --- 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_) {