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/scattering.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/scattering.c')
| -rw-r--r-- | src/scattering.c | 8 |
1 files changed, 1 insertions, 7 deletions
diff --git a/src/scattering.c b/src/scattering.c index 25b17d4..2512cda 100644 --- a/src/scattering.c +++ b/src/scattering.c @@ -191,14 +191,8 @@ void init_interpolation(void) gsl_integration_cquad_workspace_free(w); } -double get_probability(double beta, double cos_theta, double theta0) +double get_probability(double beta, double cos_theta, double sin_theta, double theta0) { - double sin_theta; - - /* Technically this isn't defined up to a sign, but it doesn't matter since - * we are going to square it everywhere. */ - sin_theta = sqrt(1-pow(cos_theta,2)); - /* Make sure theta0 is less than MAX_THETA0, otherwise it's possible that * gsl_spline2d_eval() will quit. */ if (theta0 > MAX_THETA0) theta0 = MAX_THETA0; |