@ChristopherWu
2015-02-16T11:23:44.000000Z
字数 1118
阅读 1280
未分类
亲爱的程序猿们,你们肯定都用过printf吧?你们知道,看起来理所当然的简单的printf,实际上是个难倒众多计算机科学家的难题吗?直到1971年,才有我们的毒师Mr.White科学家(Jon White)解决这个问题,直到1990年,我们的Mr.White才正式发表这算法的最终版本,Dragon4,
在随后到最近的几十年来,语言上的各种浮点数打印字符串都是用Dragon4算法,其他人的研究都只是对这个算法修修补补,直到Grisu算法的出现。Grisu算法由Florian Loitsch发表,64位的浮点数可以用它表示,但是有0.6%的依然要用Dragon4算法来表示。
因此,要是想做轮子,无论如何都要懂Dragon4算法!!! ————引言
为了修复我在wine中发现的bug--sprintf在vs和gcc下的行为不一致,gcc的printf,sprintf是四舍六入五成双的舍入(银行家舍入),而vs的则是四舍五入,开始研究起浮点数打印字符串。
## 算法的核心就是把小数点后面的每一个数字提取出来,我分别想到2种方法:对于double a
a - (int)a 这样把小数拿出来,乘以精度,得到的数值用round函数后,就是整数,接下来就容易了。可以参考我的实现,在精度不大的情况下没有问题。
可是,在乘以精度的时候就会有问题,因为浮点数不可以精确表示。比如2.34以35位精度表示成2.3400000...000,当乘以超过一定精度的时候,你就会发现,小数点后不是单纯的0而是其他数字了。
a- (int)a 这样把小数拿出来,用fmod(a, 0.1)取余数,可是问题是,0.789不断除0.1后,
double temp = float_val;
while(temp != 0) {
cout<<"int "<<(int)(temp*10)<
temp = fmod(temp, 0.1);
cout<<"float_val "<
cout<<"float_val*10= "<
temp *= 10.0;
}
运行结果:
cin val
0.789
float_val 0.789
int 7
float_val 0.089
float_val*10= 0.89
1
int 8
float_val 0.09
float_val*10= 0.9
1
int 8
float_val 0.1
float_val*10= 1
1
int 9
float_val 0.1
float_val*10= 1
问题就在最后0.09变成0.9后,除0.1再取小数就变0.1了。尝试好几种方法,无法解决。
最后只可以上网查资料,发现这个问题竟然是持续几十年的难题,要是我无意中完美解决了,就像高斯一样了XD
下一篇,开始讲解Dragon4的背景知识。