Fix another bad double bug...

[foosball] / foosrank.cpp

diff --git a/foosrank.cpp b/foosrank.cpp
index 13e9a7914b93539b653a01d55ecdb5acb6019ba8..b3eef1eca712f6dc5887cdf3c5f76c3856387baa 100644 (file)
--- a/foosrank.cpp
+++ b/foosrank.cpp
@@ -92,7 +92,7 @@ double fac(int x)
 //   +inf
 //    /
 //    |
-//    | ProbScore[a] (r1-r2) Gaussian[mu2, sigma2] (dr2) dr2
+//    | ProbScore[a] (r1-r2) Gaussian[mu2, sigma2] (r2) dr2
 //    |
 //   /
 // -inf
@@ -421,13 +421,13 @@ void compute_new_double_rating(double mu1, double sigma1, double mu2, double sig
                for (double r1 = 0.0; r1 < 3000.0; r1 += pdf_step_size) {
                        double z = (r1 - mu1) / sigma1;
                        double gaussian = exp(-(z*z/2.0));
-                       curve.push_back(make_pair(r1, gaussian * simpson_integrate(OuterIntegralEvaluator(r1,mu2,sigma2,mu3,sigma3,mu4,sigma4,score1,score2,-0.5), 0.0, 3000.0, int_step_size)));
+                       curve.push_back(make_pair(r1, gaussian * simpson_integrate(OuterIntegralEvaluator(r1,mu2,sigma2,mu3,sigma3,mu4,sigma4,score1,score2,-1.0), 0.0, 3000.0, int_step_size)));
                }
        } else {
                for (double r1 = 0.0; r1 < 3000.0; r1 += pdf_step_size) {
                        double z = (r1 - mu1) / sigma1;
                        double gaussian = exp(-(z*z/2.0));
-                       curve.push_back(make_pair(r1, gaussian * simpson_integrate(OuterIntegralEvaluator(r1,mu2,sigma2,mu3,sigma3,mu4,sigma4,score2,score1,0.5), 0.0, 3000.0, int_step_size)));
+                       curve.push_back(make_pair(r1, gaussian * simpson_integrate(OuterIntegralEvaluator(r1,mu2,sigma2,mu3,sigma3,mu4,sigma4,score2,score1,1.0), 0.0, 3000.0, int_step_size)));
                }
        }