From 7d9aeb6e2b711bd34a7e1d30ad0ef340257c1725 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Sun, 2 Jun 2019 16:22:38 -0400 Subject: fast_sqrt() -> sqrt() No reason to use fast_sqrt() in get_theta0_min() since the argument isn't guaranteed to be between 0 and 1. --- src/likelihood.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'src') diff --git a/src/likelihood.c b/src/likelihood.c index 950b646..696b52a 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -542,7 +542,7 @@ static void get_expected_charge_shower(particle *p, double *pos, double *dir, in * effectively averages across the face of a PMT. */ double get_theta0_min(double R, double r, double sin_theta_pmt) { - return fast_acos((R-r*sin_theta_pmt)/fast_sqrt(r*r + R*R - 2*r*R*sin_theta_pmt)); + return fast_acos((R-r*sin_theta_pmt)/sqrt(r*r + R*R - 2*r*R*sin_theta_pmt)); } static void get_expected_charge(double beta, double theta0, double *pos, double *dir, int pmt, double *q, double *reflected, double *t) -- cgit