diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2019-06-02 13:31:18 -0400 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2019-06-02 13:31:18 -0400 |
| commit | cc00800b0e55ab5a5d892b03e735302eb0c06d47 (patch) | |
| tree | 24f19717f241fc4393f227b06ab392dea5643749 /src/likelihood.c | |
| parent | 7fe82b13dc145c85dee8108299a1db54bb26b0c3 (diff) | |
| download | sddm-cc00800b0e55ab5a5d892b03e735302eb0c06d47.tar.gz sddm-cc00800b0e55ab5a5d892b03e735302eb0c06d47.tar.bz2 sddm-cc00800b0e55ab5a5d892b03e735302eb0c06d47.zip | |
update get_probability() to take sin(theta) as an argument
Since we already calculate sin(theta) in get_expected_charge() there's no
reason to calculate it again in get_probability(). This *may* already be
optimized out by the compiler.
Diffstat (limited to 'src/likelihood.c')
| -rw-r--r-- | src/likelihood.c | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/likelihood.c b/src/likelihood.c index 587349b..17cde5d 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -603,7 +603,7 @@ static void get_expected_charge(double beta, double theta0, double *pos, double * pretty low, this is expected to be a very small effect. */ prob_sct = 1.0 - get_fsct_d2o(l_d2o)*get_fsct_h2o(l_h2o)*get_fsct_acrylic(AV_THICKNESS); - charge = omega*FINE_STRUCTURE_CONSTANT*z*z*(1-(1/(beta*beta*n*n)))*get_probability(beta, cos_theta, theta0); + charge = omega*FINE_STRUCTURE_CONSTANT*z*z*(1-(1/(beta*beta*n*n)))*get_probability(beta, cos_theta, sin_theta, theta0); *reflected = (1.0-prob_abs)*(1.0-prob_sct)*f_reflec*charge + prob_sct*charge; |