Remove another mul or so.

diff --git a/.bzrignore b/.bzrignore
new file mode 100644 (file)
index 0000000..cbc5f51
--- /dev/null
+++ b/.bzrignore
@@ -0,0 +1 @@
+lost+found

diff --git a/foosrank.cpp b/foosrank.cpp
index 4b44706fe67fd77bcacadbfbe4093281a4c734b8..b0f30424ff7e6e7f7040f14356ff2651193e6211 100644 (file)
--- a/foosrank.cpp
+++ b/foosrank.cpp
@@ -122,11 +122,11 @@ static void compute_opponent_rating_pdf(int k, int a, double mu2, double sigma2,
                func1[i].real() = func1[i].imag() = func2[i].real() = func2[i].imag() = 0.0;
        }
 
-       double invsigma2 = 1.0 / sigma2;
+       double invsq2sigma2 = 1.0 / (sqrt(2.0) * sigma2);
        for (int i = 0; i < sz; ++i) {
                double x1 = 0.0 + h*i;
-               double z = (x1 - mu2) * invsigma2;
-               func1[i].real() = exp(-(z*z/2.0));
+               double z = (x1 - mu2) * invsq2sigma2;
+               func1[i].real() = exp(-z*z);
 
                double x2 = -3000.0 + h*i;
                func2[(i - sz/2 + sz*2)%(sz*2)].real() = prob_score_real(k, a, binomial_precompute, x2*winfac);
@@ -373,13 +373,13 @@ static void compute_new_double_rating(double mu1, double sigma1, double mu2, dou
                double r1 = i * h;
 
                // iterate over r2
-               double invsigma2 = 1.0 / sigma2;
+               double invsq2sigma2 = 1.0 / (sqrt(2.0) * sigma2);
                for (unsigned j = 0; j < curve.size(); ++j) {
                        double r1plusr2 = curve[j].first;
                        double r2 = r1plusr2 - r1;
 
-                       double z = (r2 - mu2) * invsigma2;
-                       double gaussian = exp(-(z*z/2.0));
+                       double z = (r2 - mu2) * invsq2sigma2;
+                       double gaussian = exp(-z*z);
                        sum += curve[j].second * gaussian;
                }