Re: [请益] 如何直接判断浮点数运算时有误差(赠P币)

楼主: oopFoo (3d)   2018-12-02 21:35:39
你已经知道答案了,但你可能还不知道原理。
round "half" to nearest even. 又称gaussian/banker's rounding.
例子
round[{0.5, 1.5, 2.5, 3.5, .5}] =>{0, 2, 2, 4, 4}
为什么呢?因为四舍五入是biased. (.5)是刚好在(0,1)的一半(half)。
要公平,不累进error,其实要一半舍,一半入。最简单的方法就是odd或even入,有一半的机会。
banker's rounding 其实蛮常见的。
.net Math.round的 default 就是banker's rounding.
java也可以设定
https://docs.oracle.com/javase/8/docs/api/java/math/RoundingMode.html
floating point 的话就
010....(不到一半,舍)
101....(超过一半,入)
100....(刚好一半,用banker's rounding)看要round的bit是1或0.
硬件的floating point会有3个extra bits. guard, round, sticky. sticky bit 是or所有多出的bits.
2个bits不够,因为做完arithmatic还需要normalized 回来。
请参考
http://pages.cs.wisc.edu/~david/courses/cs552/S12/handouts/guardbits.pdf
金融业,应该都知道这个吧。
PS: "What Every Computer Scientist Should Know About Floating-Point Arithmetic"
https://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html
主要在在讨论rounding error。
作者: CorkiN (柯基)   2018-12-02 22:15:00
应该说,有念过计组的就会知道
作者: Muscovy (三分熟的闹钟)   2018-12-03 00:13:00
01234 舍, 56789 入, 一半一半看起来很公平啊.考虑连续区间也是, [0, 5) 舍, [5, 10) 入, 也是很公平.为什么会说“四舍五入是 biased”? 是有什么特别原因吗?路过纯好奇.
作者: guest0079 (SpongeBob SquarePants)   2018-12-03 03:47:00
楼上 哪有一半一半 零舍去得零算在01234那一半凑了五个 但0变0根本什么都没有舍去
作者: KeyFSN ( ~☼☽✩☁~ )   2018-12-03 16:29:00
专业
作者: v420746k (Tyrone_Huang)   2018-12-03 18:20:00
作者: Muscovy (三分熟的闹钟)   2018-12-04 19:49:00
呃, 0 包含 0~0.1 之间的所有可能啊, 不要只看 0.

Links booklink

Contact Us: admin [ a t ] ucptt.com