c语言 4x4矩阵乘法

下面是两个4x4的矩阵相乘的C语言代码：

#include <stdio.h>

void matrix_multiply(float A[][4], float B[][4], float C[][4]) {
    for(int i=0; i<4; i++) {
        for(int j=0; j<4; j++) {
            float sum = 0;
            for(int k=0; k<4; k++) {
                sum += A[i][k] * B[k][j];
            }
            C[i][j] = sum;
        }
    }
}

int main() {
    float A[4][4] = {{1, 2, 3, 4},
                     {5, 6, 7, 8},
                     {9, 10, 11, 12},
                     {13, 14, 15, 16}};
    float B[4][4] = {{1, 0, 0, 0},
                     {0, 1, 0, 0},
                     {0, 0, 1, 0},
                     {0, 0, 0, 1}};
    float C[4][4];

    matrix_multiply(A, B, C);

    for(int i=0; i<4; i++) {
        for(int j=0; j<4; j++) {
            printf("%f ", C[i][j]);
        }
        printf("\n");
    }

    return 0;
}

其中 `matrix_multiply` 函数用于计算两个4x4矩阵的乘积，并将结果存储在C矩阵中。在主函数中，我们先定义两个矩阵A和B，然后调用 `matrix_multiply` 函数计算它们的乘积，最后打印输出结果矩阵C。