Commit e28504e5 authored by Kipp Cannon's avatar Kipp Cannon

rate_estimation.c: address some review requests

- update comments
parent e001c2a0
...@@ -149,8 +149,9 @@ static double log_posterior(const double *ln_f_over_b, int n, double Rf, double ...@@ -149,8 +149,9 @@ static double log_posterior(const double *ln_f_over_b, int n, double Rf, double
if(i) { if(i) {
/* /*
* for these entries, compute the sum of log(Rf f / (Rb b) * for these entries, compute the sum of log(Rf f / (Rb b)
* + 1) by approximating it with a numerical integration * + 1) by approximating it with a numerical integration of
* of the exact addend * the addend (evaluated using the approximations described
* above)
*/ */
gsl_function _integrand = { gsl_function _integrand = {
.function = integrand, .function = integrand,
...@@ -163,9 +164,8 @@ static double log_posterior(const double *ln_f_over_b, int n, double Rf, double ...@@ -163,9 +164,8 @@ static double log_posterior(const double *ln_f_over_b, int n, double Rf, double
} }
/* /*
* now compute the sum of log(Rf f / (Rb b) + 1) for the remaining * now explicitly compute the sum of log(Rf f / (Rb b) + 1) for the
* entries with an explicit loop but using approximations * remaining entries
* (described above) for cases in which the addend is large
*/ */
for(; i < n; i++) { for(; i < n; i++) {
......
Markdown is supported
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