fsum代码
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define BUFSIZE 256
int main(int argc, char *argv[]) {
char prefix[BUFSIZE];
char next[BUFSIZE];
int i;
float sum = 0.0;
for (i = 1; i < argc; i++) {
float x = atof(argv[i]);
sum += x;
if (i == 1) {
sprintf(prefix, "%.4g", x);
} else {
sprintf(next, " + %.4g", x);
strcat(prefix, next);
printf("%s = %.4g\n", prefix, sum);
}
}
return 0;
}
运行结果
代码解释
根据运行结果我们发现第二次和第三次都出现了大数吞小数的状况,为什么会出现这样的错误呢?我们明白浮点数的表示方法是符号位加阶码加尾数,当进行浮点数的加法时需要先对阶,而当两个加数的阶码相差过大,对阶时较小的那个数的尾数可能就只剩下0了,这时也就出现了小数被吞了的情况。而且第一次和第三次的结果对比我们发现浮点数的加法运算不满足交换律。
结论
在进行浮点数的加法运算时,不满足加法的交换律;对于单精度浮点数,当两者阶码相差25及以上时,就不能进行加法运算了,否则会出现错误!