C语言向函数传递二维数组报错及矩阵乘法实现问题求助

1. 代码逻辑错误 #

你在输入矩阵b的代码行写错了参数：

input(row_b, col_b, a); // 错误，应该传入b的地址

这行代码把用户输入的矩阵b的值写到了矩阵a的内存空间里，矩阵b全程没有被赋值，处于未初始化的垃圾值状态，运行自然得不到正确结果。当你用预定义的固定大小矩阵时没有走这段输入逻辑，所以不会触发这个错误。

2. GDB调试报错原因 #

你用了C99标准的变长数组（VLA）作为函数参数，GCC 6.3.0本身对VLA的编译支持是完善的，但配套的GDB版本对函数参数形式的VLA调试支持存在缺陷：GDB无法在运行时获取变长数组的列跨度信息，无法计算指针偏移，所以输出cannot perform pointer math on incomplete type try casting的报错。

1. 修正输入逻辑的参数错误，将输入矩阵b的代码改为：

input(row_b, col_b, b);

如果需要验证输入是否正确，可以打开你注释掉的output函数调用，打印输入后的矩阵值。

2. 解决GDB调试VLA的问题，两种可选方案：

• 调试时手动强转指针类型，比如要查看a数组第i行第k列的值，在GDB中输入：

  p ((double (*)[col_a])a)[i][k]
  

• 放弃栈上变长数组的实现，改用堆内存动态分配的二维数组，编译期类型明确，GDB可以直接识别。

1. 循环顺序优化（缓存友好）
   当前循环顺序是i-j-k，访问b矩阵时是按列遍历，内存步长等于b矩阵的列数，缓存命中率极低。将循环顺序调整为i-k-j，a、b矩阵都按行连续访问，缓存命中率可以提升数倍：

// 注意res矩阵需要提前初始化为0
for (int i = 0; i < row_a; i++) {
    for (int k = 0; k < col_a; k++) {
        double a_ik = a[i][k]; // 寄存器缓存a[i][k]的值，避免重复读取
        for (int j = 0; j < col_b; j++) {
            res[i][j] += a_ik * b[k][j];
        }
    }
}

2. 编译选项优化
   GCC 6.3.0可以添加以下编译参数提升性能：
   -O2 -march=native -ftree-vectorize，开启二级优化、适配本地CPU指令集、自动向量化循环运算。

3. 避免栈溢出风险
   当前所有矩阵都分配在栈空间，栈空间默认只有几MB，矩阵稍大就会触发栈溢出。建议改用malloc分配堆内存存储矩阵，支持更大规模的矩阵运算。

4. 错误提示完善
   当前矩阵维度不匹配时直接静默返回，用户无法感知错误原因，建议添加错误打印：

if (col_a != row_b) {
    fprintf(stderr, "矩阵维度不匹配，a的列数必须等于b的行数\n");
    return;
}

5. 大规模矩阵可选分块优化
   如果矩阵大小超过CPU缓存容量，可以采用分块（阻塞）乘法，将矩阵拆分为多个小块运算，进一步提升缓存利用率。

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