From 0ae23efdbba0168b6c32300c2f9881d2eddf8697 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Wed, 14 Nov 2018 10:46:51 -0600 Subject: speed things up by skipping zero values --- src/likelihood.c | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/src/likelihood.c b/src/likelihood.c index 97db51c..0855d08 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -266,10 +266,12 @@ double F(double t, double mu_noise, double mu_indirect, double *mu_direct, doubl mu_total = mu_noise + mu_indirect; for (i = 0; i < n; i++) { + if (mu_direct[i] == 0.0) continue; p += mu_direct[i]*norm_cdf(t,ts[i],sigma); mu_total += mu_direct[i]; } for (i = 0; i < n; i++) { + if (mu_shower[i] == 0.0) continue; p += mu_shower[i]*norm_cdf(t,ts_shower[i],ts_sigma[i]); mu_total += mu_shower[i]; } @@ -300,10 +302,12 @@ double f(double t, double mu_noise, double mu_indirect, double *mu_direct, doubl mu_total = mu_noise + mu_indirect; for (i = 0; i < n; i++) { + if (mu_direct[i] == 0.0) continue; p += mu_direct[i]*norm(t,ts[i],sigma); mu_total += mu_direct[i]; } for (i = 0; i < n; i++) { + if (mu_shower[i] == 0.0) continue; p += mu_shower[i]*norm(t,ts_shower[i],ts_sigma[i]); mu_total += mu_shower[i]; } -- cgit