Message242029
The next question is whether we care about that 1 ULP in the context of the random number calculations which already involve multi-step chains that aren't favorable to retaining precision:
s + (1.0 + s * s) ** 0.5
ainv = (2.0 * alpha - 1.0) ** 0.5
bbb = alpha - LOG4
ccc = alpha + ainv
g2rad = (-2.0 * _log(1.0 - random())) ** 0.5
z = _cos(x2pi) * g2rad
self.gauss_next = _sin(x2pi) * g2rad
stddev = (sqsum/n - avg*avg) ** 0.5 |
|
Date |
User |
Action |
Args |
2015-04-25 23:11:04 | rhettinger | set | recipients:
+ rhettinger, tim.peters, mark.dickinson, serhiy.storchaka |
2015-04-25 23:11:04 | rhettinger | set | messageid: <1430003464.08.0.153638572698.issue24059@psf.upfronthosting.co.za> |
2015-04-25 23:11:04 | rhettinger | link | issue24059 messages |
2015-04-25 23:11:03 | rhettinger | create | |
|