浮点运算

原代码：

#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;
}

结果：guest-ep7ljC@ubuntu:~/Desktop$ ./a 1e20 -1e20 3.14
1e+20 + -1e+20 = 0
1e+20 + -1e+20 + 3.14 = 3.14
guest-ep7ljC@ubuntu:~/Desktop$ ./a 1e30 1e20
1e+30 + 1e+20 = 1e+30
guest-ep7ljC@ubuntu:~/Desktop$ ./a 1e20 1e19
1e+20 + 1e+19 = 1.1e+20
guest-ep7ljC@ubuntu:~/Desktop$ ./a 1e30 3.54 -1e30
1e+30 + 3.54 = 1e+30
1e+30 + 3.54 + -1e+30 = 0
guest-ep7ljC@ubuntu:~/Desktop$
我的感想：一直认为“a+b-a=b“是必然正确的，然而在电脑里面却不全是这样的，就如上面的结果一样。让我觉得有趣的是，电脑似乎也像我一样喜欢取近似值，如果让我算1e30+3.54，我也会忽略掉3.54，因为1e30相比于3.54相差太大了，因而在舍入后，3.54的值丢失了。