c语言实现矩阵乘法的循环展开

矩阵乘法的循环展开是一种优化技术，可以提高计算效率。下面是用 C 语言实现矩阵乘法的循环展开的示例代码：

#include <stdio.h>

#define N 4

void matrix_multiply(int A[][N], int B[][N], int C[][N]) {
    for (int i = 0; i < N; i++) {
        for (int j = 0; j < N; j++) {
            C[i][j] = 0;
            for (int k = 0; k < N; k += 2) {
                C[i][j] += A[i][k] * B[k][j];
                C[i][j] += A[i][k+1] * B[k+1][j];
            }
        }
    }
}

void print_matrix(int matrix[][N]) {
    for (int i = 0; i < N; i++) {
        for (int j = 0; j < N; j++) {
            printf("%d ", matrix[i][j]);
        }
        printf("\n");
    }
}

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

    int B[N][N] = {{1, 2, 3, 4},
                   {5, 6, 7, 8},
                   {9, 10, 11, 12},
                   {13, 14, 15, 16}};

    int C[N][N];

    matrix_multiply(A, B, C);

    printf("Matrix C:\n");
    print_matrix(C);

    return 0;
}

这段代码中，`matrix_multiply` 函数实现了矩阵乘法的循环展开。循环展开是通过将内层循环中的迭代次数减半，将两次乘法操作合并到一次循环中，从而提高计算效率。在这个例子中，内层循环使用 `k += 2` 来减半迭代次数，并在每次迭代中进行两次乘法操作。

注意：上面的代码仅为示例，实际应用中可能还需要进行边界检查、内存分配等操作。