C语言通用未知大小矩阵乘法函数的跨编译器实现方案

你的代码使用了C99标准的变长数组（VLA），这也是它能在支持C99的编译器（比如GCC、Clang）上运行，但在旧版MSVC上报错的原因——旧版MSVC不支持C99的VLA特性，它要求数组维度必须是编译期常量。

要实现兼容所有C编译器的矩阵乘法函数，推荐以下两种方案：

方案一：用一维数组模拟二维矩阵（最通用） #

直接将二维矩阵当作一维数组传递，手动计算元素的一维索引，完全避开VLA的使用，所有C编译器都能兼容。

void matrix_multiply(int rows1, int cols1, int cols2, const float* matrix1, const float* matrix2, float* result)
{
    for (int i = 0; i < rows1; i++)
    {
        for (int j = 0; j < cols2; j++)
        {
            // 计算结果矩阵当前元素的一维索引
            int res_idx = i * cols2 + j;
            result[res_idx] = 0.0f;
            for (int k = 0; k < cols1; k++)
            {
                // 分别计算两个输入矩阵对应元素的一维索引
                int mat1_idx = i * cols1 + k;
                int mat2_idx = k * cols2 + j;
                result[res_idx] += matrix1[mat1_idx] * matrix2[mat2_idx];
            }
        }
    }
}

调用时直接传入二维数组的数组名即可（C语言中数组会自动退化为指向首元素的指针），比如：

float mat1[3][4] = {...};
float mat2[4][2] = {...};
float result[3][2];
matrix_multiply(3, 4, 2, (float*)mat1, (float*)mat2, (float*)result);

方案二：用宏封装二维索引（保留二维语法） #

如果想让代码看起来更贴近二维矩阵的写法，可以用宏封装索引计算逻辑，本质还是基于一维数组的实现：

// 宏：获取一维数组模拟的二维矩阵元素
#define MAT_GET(mat, cols, i, j) (mat[(i) * (cols) + (j)])

void matrix_multiply(int rows1, int cols1, int cols2, const float* matrix1, const float* matrix2, float* result)
{
    for (int i = 0; i < rows1; i++)
    {
        for (int j = 0; j < cols2; j++)
        {
            MAT_GET(result, cols2, i, j) = 0.0f;
            for (int k = 0; k < cols1; k++)
            {
                MAT_GET(result, cols2, i, j) += MAT_GET(matrix1, cols1, i, k) * MAT_GET(matrix2, cols2, k, j);
            }
        }
    }
}

这种写法的调用方式和方案一完全一致，但代码可读性更接近原生二维数组的写法。

额外说明：针对新版MSVC的兼容 #

如果你只需要兼容2019及以后的MSVC，可以开启编译选项/std:c99或/std:c11，这样你的原始VLA代码就能正常编译运行。但如果要兼容所有编译器（包括旧版MSVC），还是优先选择前两种方案。

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