PHP超大/极小浮点数算术运算返回NAN/INF如何处理

问题原因 #

• PHP原生64位双精度浮点的最大可表示值约为1.7976931348623E+308，示例中-8.3802985809867E+217与4.8047258326981E+215相乘结果约为4.0265E+432，远超出上限触发溢出，返回INF或NAN。
• 常规%运算符仅支持整数运算，超出整数范围就会出现计算错误；fmod仅适配普通范围浮点运算，数值溢出后直接失效；原生bc_mod仅支持正整数的模运算，不支持浮点、负数、科学计数法格式的输入。

具体实现方案 #

方案1：使用Decimal扩展（推荐，适配线性方程组矩阵求解场景） #

该扩展为高精度十进制运算专用扩展，支持任意精度的正负浮点、科学计数法输入，原生覆盖加减乘除、模运算、比较等所有常用算术操作，不会出现精度丢失或溢出问题，完全满足矩阵求解的数值计算需求。
运算示例如下：

// 直接传入科学计数法格式的字符串初始化高精度数
$a = new Decimal\Decimal("-8.3802985809867E+217");
$b = new Decimal\Decimal("4.8047258326981E+215");

// 乘法运算，精度可自定义，默认保留28位有效数字
$result = $a->mul($b);
// 输出科学计数法格式的结果：-4.0265082638315717136541407644E+433
echo $result->toScientificString();

// 高精度模运算示例，支持负数、浮点输入
$modResult = $a->mod($b);

方案2：基于BC Math扩展二次封装（无需额外安装第三方扩展） #

如果环境不支持安装额外扩展，可以对原生BC Math扩展进行封装适配：

1. 先实现科学计数法转普通数字字符串的转换逻辑，解决BC Math无法识别科学计数法输入的问题
2. 运算前统一记录所有操作数的小数位数，运算完成后做精度对齐
3. 自行实现浮点、负数的模运算逻辑：先将数值放大为整数完成模运算，再缩小对应倍数还原结果
   示例封装乘法函数参考：

function bcmul_scientific(string $a, string $b, int $precision = 50): string
{
    $a_normal = normalize_scientific($a); // 自定义函数，实现科学计数法转普通数字字符串
    $b_normal = normalize_scientific($b);
    return bcmul($a_normal, $b_normal, $precision);
}

矩阵求解场景适配注意事项 #

• 所有矩阵元素统一使用字符串类型存储，不要转换为PHP原生浮点类型，避免运算前就丢失精度
• 矩阵行变换、消元等所有操作全部使用高精度运算API实现，不要混用PHP原生算术运算符

内容的提问来源于stack exchange，提问作者Tomas Cuellar Escalona