Testing for successfully compiled code, bug fixed with logBeta_tab

* Beta was accidentally being populated by Alpha values. Fixed in hyrec/hydrogen.c/interpolate_rates * Renamed files from RR to R to be compatible with previous naming convention. * Producing test_xe.dat to check that the recombination history is sensible.
hongwanliu · Sep 26, 2019 · e395201 · e395201
1 parent fe5301b
commit e395201
Show file tree

Hide file tree

Showing 6 changed files with 108,525 additions and 114 deletions.
diff --git a/hyrec/RR_2pto2s_BB_n_250.dat → hyrec/R_2pto2s_BB_n_250.dat b/hyrec/RR_2pto2s_BB_n_250.dat → hyrec/R_2pto2s_BB_n_250.dat
diff --git a/hyrec/RR_2sto2p_BB_n_250.dat → hyrec/R_2sto2p_BB_n_250.dat b/hyrec/RR_2sto2p_BB_n_250.dat → hyrec/R_2sto2p_BB_n_250.dat
diff --git a/hyrec/hydrogen.c b/hyrec/hydrogen.c
@@ -204,10 +204,10 @@ void interpolate_rates(double Alpha[2], double Beta[2], double *R2p2s, double *R
         + rate_table->logAlpha_tab[l][iTR+2]*coeff2[3]
       ); 
       Beta[l] = exp(
-        rate_table->logAlpha_tab[l][iTR-1]*coeff2[0]
-        + rate_table->logAlpha_tab[l][iTR]*coeff2[1]
-        + rate_table->logAlpha_tab[l][iTR+1]*coeff2[2]
-        + rate_table->logAlpha_tab[l][iTR+2]*coeff2[3]
+        rate_table->logBeta_tab[l][iTR-1]*coeff2[0]
+        + rate_table->logBeta_tab[l][iTR]*coeff2[1]
+        + rate_table->logBeta_tab[l][iTR+1]*coeff2[2]
+        + rate_table->logBeta_tab[l][iTR+2]*coeff2[3]
       ); 
     }
 

diff --git a/notebooks/recfast_fudge.ipynb b/notebooks/recfast_fudge.ipynb