diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2018-11-14 10:46:51 -0600 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2018-11-14 10:46:51 -0600 |
| commit | 0ae23efdbba0168b6c32300c2f9881d2eddf8697 (patch) | |
| tree | 29a7a07cd86bd3179e3c212d58aee7be6658c23e | |
| parent | ef529416b2995889fa08c0de03f179b9ed3299a9 (diff) | |
| download | sddm-0ae23efdbba0168b6c32300c2f9881d2eddf8697.tar.gz sddm-0ae23efdbba0168b6c32300c2f9881d2eddf8697.tar.bz2 sddm-0ae23efdbba0168b6c32300c2f9881d2eddf8697.zip | |
speed things up by skipping zero values
| -rw-r--r-- | src/likelihood.c | 4 |
1 files changed, 4 insertions, 0 deletions
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]; } |