round() は四捨五入ではないよという話（何回目?）

Python 公式ドキュメントに

> 注釈 浮動小数点数に対する round() の振る舞いは意外なものかもしれません: 例えば、 round(2.675, 2) は予想通りの 2.68 ではなく 2.67 を与えます。これはバグではありません: これはほとんどの小数が浮動小数点数で正確に表せないことの結果です。詳しくは 浮動小数点演算、その問題と制限 を参照してください。

と書いてあるそうだ。この説明は正確ではない。

2.5 も 3.5 も「浮動小数点数で正確に表せる」が，round の結果は異なる。

>>> round(2.5)
2
>>> round(3.5)
4

 

2.5 は 16 進表現で (2.8)₁₆ で，10 進数に変換するのは 2 * 16^0 + 8 * 16^(-1) = 2 + 0.5 である。

3.5 は 16 進表現で (3.8)₁₆ で，10 進数に変換するのは 3 * 16^0 + 8 * 16^(-1) = 2 + 0.5 である。

 

つまり，小数部が 1/2^1, 1/2^2, 1/2^3, 1/2^4, 1/2^5, 1/2^6... = 0.5, 0.25, 0.125, 0.0625, 0.03125, 0.015625... の任意の組み合わせで表せる場合は，「浮動小数点数で正確に表せる」。

 

0.1640625 = 1/2^3 + 1/2^5 + 1/2^7 であり，16進数で表せば (0.2A)₁₆，2進数で表せば (0.0010101)₂ つまり，2進数の小数点以下 3, 5, 7 が 1 の数であり，正確に表すことができるのである。

 

Python の公式文書は（翻訳者がミスしただけかもしれないが）ダメダメであるが，他の言語ではどうだろうか?

 

R: `? round` において

 

‘Details’ about “round to even” when rounding off a 5

 

Rounding to decimal digits in binary arithmetic is non-trivial (when digits != 0) and may be surprising. Be aware that most decimal fractions are not exactly representable in binary double precision. In R 4.0.0, the algorithm for round(x, d), for $d>0$, has been improved to measure and round “to nearest even”, contrary to earlier versions of R (or also to sprintf() or format() based rounding).

と書いているように，「四捨五入をしようとしたとき，結果（の末尾）が偶数になるように丸める」のである。

つまり，2.5 を「学校で習った四捨五入」をすると 3 になるのであるが，これは奇数なので 2 に丸める（2.5 は 2 と 3 のちょうど真ん中なので，2 に丸めても 3 に丸めてもよいのだけど結果が偶数になるように 2 に丸めるということ）。

3.5 を「学校で習った四捨五入」をすると 4 になるのであるが，これは偶数なので 4 のままにする。

Julia でも，being rounded to the nearest even integer と書いてある。

The RoundingMode r controls the direction of the rounding; the default is RoundNearest, which rounds to the nearest integer, with ties (fractional values of 0.5) being rounded to the nearest even integer. Note that round may give incorrect results if the global rounding mode is changed (see rounding).

こうなると，Python はどうなんだ?と思う。help(np.around) で表示される文書に，ちゃんと書いてある（np.round より np.around をつかえといわれる）。

For values exactly halfway between rounded decimal values, NumPy rounds to the nearest even value. Thus 1.5 and 2.5 round to 2.0, -0.5 and 0.5 round to 0.0, etc.

これに続いて更に詳しく書かれている。ぜひちゃんと読むべし。

冒頭に掲げた，寝言（の原文）はどこに書いてあるのかな?