// integration step size
static const double step_size = 10.0;
+// rating constant (see below)
+static const double rating_constant = 455.0;
+
using namespace std;
double prob_score(double a, double rd);
//
double prob_score(double a, double rd)
{
- return prob_score_real(a, prodai(a), rd/455.0);
+ return prob_score_real(a, prodai(a), rd/rating_constant);
}
// Same, but takes in Product(a+i, i=1..9) as an argument in addition to a. Faster
{
double sum = 0.0;
double prodai_precompute = prodai(a);
- winfac /= 455.0;
+ winfac /= rating_constant;
for (double r2 = 0.0; r2 < 3000.0; r2 += step_size) {
double x = r2 + step_size*0.5;
double probscore = prob_score_real(a, prodai_precompute, (r1 - x)*winfac);