C语言数据科学项目矩阵乘法精度误差问题排查求助

• 严重内存分配错误（核心原因）：ddg函数中二级指针分配逻辑错误，原代码res = (double**)calloc(n*n, sizeof(double*));错误分配了n*n个一级指针，实际仅需分配n个一级指针。你仅对前n个指针分配了实际的double数组空间，剩余指针为未初始化的野指针，矩阵乘法时若访问到对应位置会读取内存垃圾值，这就是本应全零的位置出现异常小数值的核心诱因。
• 浮点数判断逻辑错误：mat_mult中使用assert(res[i][j] == 0)判断结果矩阵是否初始化完成，浮点数受精度限制不能直接用==做等值判断，若初始化残留极小非零值会导致断言误判，也可能引入初始误差。
• 打印精度与实际存储精度不匹配：当前打印矩阵仅输出8位小数，你看到的输入矩阵是截断后的值，实际内存中double类型有15~17位有效数字，用截断后的值手动计算的预期结果本身就和全精度计算的结果存在偏差。
• 浮点数累积误差：当前mat_mult使用ijk循环顺序，乘加操作的累积误差会逐步放大，小矩阵场景下误差虽小但会传递到后续计算。
• 冗余操作隐患：lnorm中调用mat_mult(copyMatrix(res,n), ddg,tmp,n)时，copyMatrix返回的临时矩阵未释放会造成内存泄漏，若copyMatrix实现存在逻辑错误也会导致数值拷贝异常。

1. 修复ddg函数的内存分配错误，将res = (double**)calloc(n*n, sizeof(double*));修改为res = (double**)calloc(n, sizeof(double*));，匹配二级矩阵的分配逻辑。
2. 修正浮点数判断逻辑，将assert(res[i][j] == 0)修改为assert(fabs(res[i][j]) < 1e-12)，同时确保每次调用mat_mult前结果矩阵已被正确初始化为全0。
3. 调整矩阵打印函数的输出精度，使用%.15g格式符输出全精度数值，确认输入矩阵的真实存储值后再核对预期结果，避免截断误差造成的误判。
4. 优化矩阵乘法的精度：
   • 调整循环顺序为kij，提升缓存命中率的同时减少浮点数累积误差
   • 开启编译器FMA（融合乘加）优化选项（如GCC加-mfma参数），减少单步乘加的误差
   • 必要时将累加中间变量改为long double类型，进一步降低累积误差
5. 优化lnorm中的内存操作，去掉不必要的copyMatrix调用，直接传入res作为矩阵乘法的输入参数，避免临时对象引入的异常和内存泄漏。
6. 删除mat_mult中res[i][j] += mat1[i][k] * mat2[k][j];;末尾多余的分号，避免后续代码维护时出现误判。

内容的提问来源于stack exchange，提问作者Itai Zemah