为什么Math.round(0.4999999999999999994)返回1?

在下面的程序中,您可以看到每个略小于.5的值都是四舍五入,除了0.5

for (int i = 10; i >= 0; i--) {
long l = Double.doubleToLongBits(i + 0.5);
double x;
do {
x = Double.longBitsToDouble(l);
System.out.println(x + " rounded is " + Math.round(x));
l--;
} while (Math.round(x) > i);
}

打印

10.5 rounded is 11
10.499999999999998 rounded is 10
9.5 rounded is 10
9.499999999999998 rounded is 9
8.5 rounded is 9
8.499999999999998 rounded is 8
7.5 rounded is 8
7.499999999999999 rounded is 7
6.5 rounded is 7
6.499999999999999 rounded is 6
5.5 rounded is 6
5.499999999999999 rounded is 5
4.5 rounded is 5
4.499999999999999 rounded is 4
3.5 rounded is 4
3.4999999999999996 rounded is 3
2.5 rounded is 3
2.4999999999999996 rounded is 2
1.5 rounded is 2
1.4999999999999998 rounded is 1
0.5 rounded is 1
0.49999999999999994 rounded is 1
0.4999999999999999 rounded is 0

我使用Java 6更新31。

77717 次浏览

总结

在Java 6(可能更早)中,round(x)被实现为floor(x+0.5)。这是一个规范错误,因为这是一个病理病例。Java 7不再强制执行这个破碎的实现

这个问题

0.5+ 0.4999999999999999994正好是1的双精度:

static void print(double d) {
System.out.printf("%016x\n", Double.doubleToLongBits(d));
}


public static void main(String args[]) {
double a = 0.5;
double b = 0.49999999999999994;


print(a);      // 3fe0000000000000
print(b);      // 3fdfffffffffffff
print(a+b);    // 3ff0000000000000
print(1.0);    // 3ff0000000000000
}

这是因为0.4999999999999999994的指数小于0.5,所以当它们相加时,它的尾数会移位,ULP会变大。

解决方案

从Java 7开始,OpenJDK(例如)是这样实现的:4

public static long round(double a) {
if (a != 0x1.fffffffffffffp-2) // greatest double value less than 0.5
return (long)floor(a + 0.5d);
else
return 0;
}

< p > <子> 1. # EYZ0 < / sub > < / p > < p > <子> 2. http://bugs.java.com/bugdatabase/view_bug.do?bug_id=6430675(感谢@SimonNickerson找到这个) < / sub > < / p > < p > <子> 3.# EYZ0 < / sub > < / p > < p > <子> 4. # EYZ0 < / sub > < / p >

我在JDK 1.6 32位上得到了相同的结果,但在Java 7 64位上,我为0.4999999999999999994得到了0,四舍五入为0,最后一行不打印。这似乎是一个虚拟机问题,但是,使用浮点,您应该期望在不同的环境(CPU, 32位或64位模式)上的结果有所不同。

并且,当使用round或逆矩阵等时,这些可以产生巨大的差异。

x64输出:

10.5 rounded is 11
10.499999999999998 rounded is 10
9.5 rounded is 10
9.499999999999998 rounded is 9
8.5 rounded is 9
8.499999999999998 rounded is 8
7.5 rounded is 8
7.499999999999999 rounded is 7
6.5 rounded is 7
6.499999999999999 rounded is 6
5.5 rounded is 6
5.499999999999999 rounded is 5
4.5 rounded is 5
4.499999999999999 rounded is 4
3.5 rounded is 4
3.4999999999999996 rounded is 3
2.5 rounded is 3
2.4999999999999996 rounded is 2
1.5 rounded is 2
1.4999999999999998 rounded is 1
0.5 rounded is 1
0.49999999999999994 rounded is 0

这似乎是一个已知的错误(Java错误6430675:数学。Round对于0x1.fffffffffffffp-2具有惊人的行为),已在Java 7中修复。

JDK 6中的源代码:

public static long round(double a) {
return (long)Math.floor(a + 0.5d);
}

JDK 7中的源代码:

public static long round(double a) {
if (a != 0x1.fffffffffffffp-2) {
// a is not the greatest double value less than 0.5
return (long)Math.floor(a + 0.5d);
} else {
return 0;
}
}

当值为0.4999999999999999994d时,在JDK 6中,它将调用地板上并因此返回1,但在JDK 7中,if条件是检查该数字是否是小于0.5的最大双精度值。在本例中,该数字不是小于0.5的最大双精度值,因此else块返回0。

您可以尝试0.49999999999999999d,它将返回1,但不会返回0,因为这是小于0.5的最大双精度值。

下面的答案摘自甲骨文 错误报告6430675在。请访问报告以获得完整的解释。

方法{Math, StrictMath。轮在操作上定义为

(long)Math.floor(a + 0.5d)

对于双重论证。虽然这个定义通常如预期的那样工作,但对于0x1,它给出了令人惊讶的结果1,而不是0。fffffffffffffp-2(0.49999999999999994)。

值0.4999999999999999994是小于0.5的最大浮点值。作为十六进制浮点字面值,它的值是0x1。fffffffffffp-2,等于(2 -2 ^52)* 2^ 2。==(0.5 - 2^54)。因此,和的准确值

(0.5 - 2^54) + 0.5

是1 - 2^54。这是两个相邻浮点数(1 - 2^53)和1之间的中间值。在Java使用的IEEE 754算术舍入到最接近偶数舍入模式中,当一个浮点结果不精确时,必须返回两个可表示的浮点值中包含精确结果的较近的一个;如果两个值相等接近,则返回最后一位为0的值。在这种情况下,add的正确返回值是1,而不是小于1的最大值。

虽然方法是按定义操作的,但这个输入的行为是非常令人惊讶的;规范可以修改为“四舍五入到最接近的长,四舍五入绑定”,这将允许更改此输入的行为。