※ 引述《yellowfishie (喵喵喵喵~~~)》之铭言： : : 推 crazying:你都取 floor 还不一样就不会是 rounding 的问题了 = =a 03/17 23:37 我觉得不是rounding的问题, 因为都是IEEE 754的标准! 我不觉得Sun和AMD/intel会做错. (他们应该不想像Pentium一样赔到吐血) 所以我认为是software的问题居多. (也可能是gcc/glibc的问题, 因为float怪怪的) : solaris: : F = 2288000*0.576923 + 5250000*0.433333 = 3.595e+06 : ori: 3.595e+06 : (int) floor cast: 3595000 : linux: : F = 2288000*0.576923 + 5250000*0.433333 = 3.595e+06 : ori: 3.595e+06 : (int) floor cast: 3594999 orz是我机器的名字 XD 现在是ubuntu gwliao@orz:~/T$ cat a.c #include <stdio.h> int main(void) { double a,b,c,d,F; int I; a=2288000; b=0.576923; c=5250000; d=0.433333; F=a*b+c*d; I=F; printf("F=%2.30lf\n",F); printf("I=%-30d\n",I); } ---------------------------------------------- gwliao@orz:~/T$ gcc a.c;./a.out (这是intel+ubuntu) F=3594998.074000000022351741790771484375 I=3594998 [gwliao@edasun T]$ gcc a.c;./a.out (这是edasun) F=3594998.074000000022351741790771484375 I=3594998 vlsi-24 -r1- >./a.out (这也是Sun) F=3594998.074000000022351741790771484375 I=3594998 s917426@moon[/tmp/T]> ./a.out (这是AMD+FreeBSD) F=3594998.074000000022351741790771484375 I=3594998 ----------------------------------------------- 我不知道那里有问题. 因为我的答案是一样. Orz 不过用float的话, 结果的确不一样. [gwliao@orz:~/T$ gcc a.c;./a.out F=3594998.250000000000000000000000000000 I=3594998 [gwliao@edasun T]$ gcc a.c;./a.out F=3594998.000000000000000000000000000000 I=3594998 vlsi-24 -r1- >gcc a.c;./a.out F=3594998.000000000000000000000000000000 I=3594998 s854128@moon[/tmp/T]> ./a.out F=3594998.250000000000000000000000000000 I=3594998 ---------------------------------------------------- gcc有个OPTION是 -msoft-float, 我将它enable. float: s854128@moon[/tmp/T]> ./a.out F=3594998.000000000000000000000000000000 I=3594998 gwliao@orz:~/T$ ./a.out F=3594998.000000000000000000000000000000 I=3594998 Sun的机器少了有些东西, 所以没办法像ubuntu一样用软87. 但是ubuntu用软87後, 算出来的值跟Sun用硬体算出来的值一样. --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 140.112.230.125 ※ 编辑: gwliao 来自: 140.112.230.125 (03/18 12:50)