C语言写出使用OpenMP计算大矩阵相乘并找出最佳线程数的完整程序

下面是一个使用OpenMP计算大矩阵相乘并找出最佳线程数的C语言程序：

#include <stdio.h>
#include <stdlib.h>
#include <omp.h>

#define N 1000

int main() {
    int i, j, k, num_threads, best_threads;
    double start_time, end_time, best_time;
    double a[N][N], b[N][N], c[N][N];

    // 初始化矩阵a和b
    for (i = 0; i < N; i++) {
        for (j = 0; j < N; j++) {
            a[i][j] = i + j;
            b[i][j] = i - j;
        }
    }

    best_time = 1e9; // 初始最佳时间设置为一个大数

    // 计算最佳线程数
    for (num_threads = 1; num_threads <= 8; num_threads++) {
        omp_set_num_threads(num_threads); // 设置线程数
        start_time = omp_get_wtime(); // 记录开始时间

        // 矩阵相乘
        #pragma omp parallel for private(j, k) shared(a, b, c)
        for (i = 0; i < N; i++) {
            for (j = 0; j < N; j++) {
                c[i][j] = 0;
                for (k = 0; k < N; k++) {
                    c[i][j] += a[i][k] * b[k][j];
                }
            }
        }

        end_time = omp_get_wtime(); // 记录结束时间
        printf("%d threads: %f seconds\n", num_threads, end_time - start_time);

        // 更新最佳时间和最佳线程数
        if (end_time - start_time < best_time) {
            best_time = end_time - start_time;
            best_threads = num_threads;
        }
    }

    printf("Best time: %f seconds with %d threads.\n", best_time, best_threads);

    return 0;
}

这个程序使用双重循环计算矩阵相乘，并使用OpenMP进行并行化。在计算过程中，通过不断改变线程数来寻找最佳线程数，并记录最佳时间和最佳线程数。最后输出最佳时间和最佳线程数。