晟松的代码学习日志之fsum.c

4:fsum.c
这个代码是为了测试精度的问题，当一个极大地数加一个较小的数，最后保存的到底精不精确。
代码为：

#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;
}
/*
@ubuntu:/mnt/hgfs/share/csapp_code$ ./a.out 1e20 -1e20 3.14
1e+20 + -1e+20 = 0
1e+20 + -1e+20 + 3.14 = 3.14
@ubuntu:/mnt/hgfs/share/csapp_code$ ./a.out -1e20 3.14  
-1e+20 + 3.14 = -1e+20
@ubuntu:/mnt/hgfs/share/csapp_code$ ./a.out -1e20 3.14  1e20
-1e+20 + 3.14 = -1e+20
-1e+20 + 3.14 + 1e+20 = 0

*/

从最后的测试结果可以看出当一个极大的数加一个较小的数时，最后会忽略掉那个较小的数，这是因为浮点数表示的小数点后的数的精度不够，当超过了这个精度能表示的范围的时候，就会发生截断，所以就会发生大数吃小数的情况。
其中的三个数的加法只是运算的顺序不同，所以在运算的时候的数的准确性就发生了改变，最后的结果也跟着改变。