如何正确地比较两个double类型的数据

今天写代码发现一个很有意思的问题。

代码如下：

double sum = 0;
for(i = 11213; i < 11325; i++) {
sum += v1[i];
}
printf("sum = %f\n", sum);
double sum2 = 66.144421;
if(sum == sum2) {
printf("sum = 66.144421 .\n");
} else {
printf("sum != 66.144421 .\n");
}
double a = 66.144421, b = 66.144421;
if (a == b) {
printf("a = b = 66.144421\n");
}
exit(0);

运行结果如下：

sum打印出来的值和sum2相等，但是与sum2进行比较的时候得到的结果却是不等的。这种现象应该跟double类型表示的精度有关。
所以当我们比较两个double值是否相等时，应该采用相减的方式。即：

double sum = 0;
for(i = 11213; i < 11325; i++) {
sum += v1[i];
}
printf("sum = %f\n", sum);
double sum2 = 66.144421;
if(sum - sum2 < 1e-6) {
printf("sum = 66.144421 .\n");
} else {
printf("sum != 66.144421 .\n");
}
double a = 66.144421, b = 66.144421;
if (a == b) {
printf("a = b = 66.144421\n");
}
exit(0);

在写代码的时候要特别注意这一点，如果程序因为这种原因发生问题，排查起来很困难。