From: Steinar H. Gunderson Date: Wed, 17 Oct 2007 16:36:41 +0000 (+0200) Subject: Fix another bad double bug... X-Git-Url: https://git.sesse.net/?p=foosball;a=commitdiff_plain;h=f74ca5032a740a1bb3c7a3d9edf8545c36cb0136 Fix another bad double bug... --- diff --git a/foosrank.cpp b/foosrank.cpp index b807084..b3eef1e 100644 --- a/foosrank.cpp +++ b/foosrank.cpp @@ -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))); } }