数值计算是计算机科学的核心领域之一，C语言提供了丰富的数值类型和运算方法，但同时也存在一些陷阱。本文详细介绍C语言数值计算的要点。
浮点数精度问题
浮点数采用二进制表示，很多十进制小数无法精确表示：
/* 浮点数精度问题 */
float f = 0.1f;
double d = 0.1;
/* 0.1 在二进制中是无限循环小数 */
printf("float: %.20f\n", f); /* 0.10000000149011611938 */
printf("double: %.20f\n", d); /* 0.10000000000000000555 */
/* 避免直接比较浮点数是否相等 */
double a = 0.1 + 0.2;
double b = 0.3;
/* a == b 的结果为 false，因为 0.1+0.2 != 0.3 */
/* 正确方式：使用误差范围比较 */
#define EPSILON 1e-9
if (fabs(a - b) < EPSILON) {
/* 认为相等 */
}整数溢出
整数运算可能发生溢出，有符号和无符号整数的溢出行为不同：
/* 有符号整数溢出是未定义行为 */
int i = 2147483647;
i++; /* 未定义行为：可能崩溃、可能静默回绕 */
/* 无符号整数溢出是定义行为（模 2^n） */
unsigned int u = 4294967295;
u++; /* u 变为 0，定义行为 */
/* 检测加法溢出 */
int safe_add(int a, int b) {
/* 检查是否会溢出 */
if (b > 0 && a > INT_MAX - b) {
/* 处理溢出 */
}
return a + b;
}数值稳定性
有些算法在数学上正确，但在计算机上可能不稳定：
/* 不稳定的求根公式 */
/* x^2 - 4x + 3.9999999 = 0 的两个根 */
double a = 1, b = -4, c = 3.9999999;
double r1 = (-b + sqrt(b*b - 4*a*c)) / (2*a); /* 可能丢失精度 */
double r2 = (-b - sqrt(b*b - 4*a*c)) / (2*a);
/* 稳定的求根方式：使用韦达定理 */
double r1_stable = (2*c) / (-b - sqrt(b*b - 4*a*c));
double r2_stable = (2*c) / (-b + sqrt(b*b - 4*a*c));
/* 大数吃小数问题 */
double big = 1e20;
double small = 1e-20;
double result = big + small; /* 结果还是 1e20，丢失了 small */类型选择
根据需要选择合适的数值类型：
/* 整数类型选择 */
int8_t /* -128 到 127 */
uint8_t /* 0 到 255 */
int16_t /* -32768 到 32767 */
int32_t /* 约 ±21 亿 */
int64_t /* 约 ±9.2 quintillion */
/* 浮点数类型选择 */
float /* 32 位，约 7 位十进制精度 */
double /* 64 位，约 15 位十进制精度 */
long double /* 80 位或 128 位，更高精度（但跨平台性差） */数值计算最佳实践
- 避免直接比较浮点数相等，使用误差范围
- 有符号整数运算前检查溢出
- 注意大数吃小数的问题
- 选择合适的数值类型
- 对于高精度计算，考虑使用定点数或专用数学库
安全建议
在安全敏感的代码中，整数溢出可能导致严重的安全漏洞。使用编译器提供的溢出检查选项（如 GCC 的 -ftrapv）可以帮助发现潜在问题。
总结
- 浮点数精度有限，避免直接比较
- 有符号整数溢出是未定义行为
- 注意数值稳定性和精度损失
- 选择合适的数据类型