在C语言中，以二维数组作为参数的函数必须在参数列表中说明列数，所以C程序员通常只使用一维数组，然后用代码中显示地将二维下标转换为--维下标。修改本书中的OpenMP矩阵–向量乘法程序，采用一维数组表示矩阵。

以下是使用一维数组表示矩阵的OpenMP矩阵-向量乘法程序：

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

#define N 1000

int main() {
    int i, j, tid;
    double start_time, end_time;
    double a[N*N], x[N], y[N];
    omp_set_num_threads(4); // 设置线程数为4

    // 初始化矩阵 a 和向量 x
    for (i = 0; i < N; i++) {
        x[i] = i * 1.0;
        for (j = 0; j < N; j++) {
            a[i*N+j] = i * j * 1.0;
        }
    }

    // 计时开始
    start_time = omp_get_wtime();

    // 矩阵-向量乘法
    #pragma omp parallel for private(j, tid)
    for (i = 0; i < N; i++) {
        tid = omp_get_thread_num();
        y[i] = 0.0;
        for (j = 0; j < N; j++) {
            y[i] += a[i*N+j] * x[j];
        }
        printf("Thread %d: i=%d, y[%d]=%f\n", tid, i, i, y[i]);
    }

    // 计时结束
    end_time = omp_get_wtime();

    // 输出计算结果和计算时间
    printf("Calculation time: %f seconds\n", end_time - start_time);
    printf("y[0]=%f, y[%d]=%f\n", y[0], N-1, y[N-1]);

    return 0;
}

在这个代码中，我们将矩阵 a 表示为一维数组 a[N*N]，其中 a[i*N+j] 表示矩阵的第 i 行和第 j 列的元素。这样，我们就可以使用一维数组表示矩阵，并且不需要在函数参数列表中指定列数。在矩阵-向量乘法中，我们仍然使用了 OpenMP 的并行 for 循环指令 #pragma omp parallel for，但是我们需要注意到在循环中使用了一维数组表示矩阵元素，即 a[i*N+j]，而不是二维数组表示矩阵元素，即 a[i][j]。