From c9813b4e9d96a16129271afee59d47dd4220fd90 Mon Sep 17 00:00:00 2001 From: "Steinar H. Gunderson" Date: Sun, 21 Oct 2007 14:54:11 +0200 Subject: [PATCH] Remove another mul or so. --- .bzrignore | 1 + foosrank.cpp | 12 ++++++------ 2 files changed, 7 insertions(+), 6 deletions(-) create mode 100644 .bzrignore diff --git a/.bzrignore b/.bzrignore new file mode 100644 index 0000000..cbc5f51 --- /dev/null +++ b/.bzrignore @@ -0,0 +1 @@ +lost+found diff --git a/foosrank.cpp b/foosrank.cpp index 4b44706..b0f3042 100644 --- 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; } -- 2.39.2