作者gwliao (gwliao)
看板NTUGIEE_EDA
标题Re: [转录]Re: [问题]请问如何四舍五入...
时间Sat Mar 18 16:41:12 2006
※ 引述《yellowfishie (喵喵喵喵~~~)》之铭言：
: 测试程式 & 结果:
: cout << "ori: " << F << endl;
: cout << "cast_i: " << cast_i << endl;
: cout << "cast_floor_i: " << cast_floor_i << endl;
: cout << "cast_ceil_i: " << cast_ceil_i << endl;
printf("\n\n");
printf("ori: %2.30lf\n", F);
printf("cast_i: %-30d\n",cast_i);
printf("cast_floor_i: %-30d\n",cast_floor_i);
printf("cast_ceil_i: %-30d\n",cast_ceil_i);
: return 0;
: }
: linux:
: ori: 3.595e+06
: cast_i: 3594998
: cast_floor_i: 3594998
: cast_ceil_i: 3594999
: solaris:
: ori: 3.595e+06
: cast_i: 3594998
: cast_floor_i: 3594998
: cast_ceil_i: 3594998
正确答案是3594998.074 (这是我用整数算出来的)
float版本
Linux:(应该说x86)
ori: 3.595e+06
cast_i: 3594998
cast_floor_i: 3594998
cast_ceil_i: 3594999
ori: 3594998.250000000000000000000000000000
cast_i: 3594998
cast_floor_i: 3594998
cast_ceil_i: 3594999 (这个很正常, 因为x86就是比Sun多了0.25)
solaris:
ori: 3.595e+06
cast_i: 3594998
cast_floor_i: 3594998
cast_ceil_i: 3594998
ori: 3594998.000000000000000000000000000000
cast_i: 3594998
cast_floor_i: 3594998
cast_ceil_i: 3594998
double: x86跟Sun都一样
ori: 3.595e+06
cast_i: 3594998
cast_floor_i: 3594998
cast_ceil_i: 3594999
ori: 3594998.074000000022351741790771484375
cast_i: 3594998
cast_floor_i: 3594998
cast_ceil_i: 3594999
还有提醒你, cout预设的输出位数跟本不能用来讨论这问题.
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.230.125
※ 编辑: gwliao 来自: 140.112.230.125 (03/18 17:36)
1F：推 moonshade:能不能把ori 的值印出小数一下... 03/18 19:21
2F：→ moonshade:用float 的问题不是以前就发生过了 03/18 19:22
3F：→ moonshade:喔，我看到了XD 03/18 19:23
4F：→ moonshade:我就说了算出来答案就不一样了...做floor ceil 当然不一 03/18 19:23
5F：→ moonshade:样... 03/18 19:23
6F：推 moonshade:我觉得可以不用gcc再试一次，不然用sizeof 看一下 03/18 19:25
7F：→ moonshade:float 的大小，搞不好linux跟solaris 不一样 03/18 19:26
8F：推 Donnie:用 double 看看? 03/18 21:50
9F：推 gwliao:我有用double! Sun和x86一样 03/18 23:47
10F：→ gwliao:先做A=i1*f1和B=i2*f3, 再做F=A+B. 03/18 23:48
11F：→ gwliao:跟一次做完(i1*f1+i2*f2)的结果(x86的机器)不一样. XD 03/18 23:48
i1*f1+i2*f2 3594998.250000000000000000000000000000
A+B 3594998.000000000000000000000000000000
12F：推 moonshade:这个应该就是compiler处理的方式了XD 03/18 23:51
我觉得没问题 XD
秀asm code给你看.
f3=i1*f1;
f4=i2*f2;
f5=f3+f4;
f55=i1*f1+i2*f2;
fildl -72(%ebp)//i1*f1;
fmuls -80(%ebp)
fstps -40(%ebp)//f3 <=
fildl -68(%ebp)//i2*f2;
fmuls -76(%ebp)
fstps -36(%ebp)//f4 <=
flds -40(%ebp)//f3+f4;
fadds -36(%ebp)
fstps -32(%ebp)//f5 <=
fildl -72(%ebp)//i1*f1;
fmuls -80(%ebp)
fildl -68(%ebp)//i2*f2
fmuls -76(%ebp)
faddp %st, %st(1)//*
fstps -28(%ebp)//f55 <=
flds -32(%ebp)
subl $4, %esp
leal -8(%esp), %esp
fstpl (%esp)
pushl $.LC13//f5
call printf
addl $16, %esp
flds -28(%ebp)
subl $4, %esp
leal -8(%esp), %esp
fstpl (%esp)
pushl $.LC14
call printf/f55
addl $16, %esp
这一段asm code, 我只加空白行和注解.
13F：推 gwliao:我用gcc-4没有-O的flag, 不过一加-O就都一样. 03/18 23:56
14F：→ gwliao:前面的结果被forward到後面. 03/18 23:58
15F：→ gwliao:都跟Sun的结果一样 03/18 23:59
※ 编辑: gwliao 来自: 140.112.230.125 (03/19 00:15)
16F：推 gwliao:所以我觉得是X86的register-forwarding的实做有关. 03/19 00:15
17F：→ gwliao:眼睛瞎了faddp %st,%st(1)怎麽是乘, 应该是加的运算. 03/19 00:17
18F：推 moonshade:好伤眼睛的code啊Orz...不过我看不太懂CISC的asm 03/19 00:27
19F：→ moonshade:(在地上打滚) 03/19 00:28
20F：→ moonshade:新的x86 asm我已经看不懂了Orz 03/19 00:29
21F：推 gwliao:那是AT&T式的asm code, 跟MASM的不同. 那算x86而已. 03/19 00:47
22F：→ gwliao:没MMX/3D now/SSE....等, 看到就晕.... 03/19 00:48
23F：→ gwliao:不过.......刚刚那些就让我眼睛瞎过一次. Orz 03/19 00:48