pthread矩阵向量乘法

下面是一个使用pthread进行矩阵向量乘法的示例代码：

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

#define MAX_SIZE 1000

int A[MAX_SIZE][MAX_SIZE];
int x[MAX_SIZE];
int y[MAX_SIZE];
int row, col;

struct thread_data {
    int start;
    int end;
};

void *multiply(void *args) {
    struct thread_data *data = (struct thread_data *) args;
    for (int i = data->start; i <= data->end; i++) {
        int sum = 0;
        for (int j = 0; j < col; j++) {
            sum += A[i][j] * x[j];
        }
        y[i] = sum;
    }
    pthread_exit(NULL);
}

int main() {
    printf("Enter the number of rows and columns of the matrix (row, col): ");
    scanf("%d %d", &row, &col);

    printf("Enter the elements of the matrix: \n");
    for (int i = 0; i < row; i++) {
        for (int j = 0; j < col; j++) {
            scanf("%d", &A[i][j]);
        }
    }

    printf("Enter the elements of the vector: \n");
    for (int i = 0; i < col; i++) {
        scanf("%d", &x[i]);
    }

    int num_threads;
    printf("Enter the number of threads: ");
    scanf("%d", &num_threads);

    pthread_t threads[num_threads];
    struct thread_data data[num_threads];
    int chunk_size = row / num_threads;

    for (int i = 0; i < num_threads; i++) {
        data[i].start = i * chunk_size;
        data[i].end = i == num_threads - 1 ? row - 1 : (i + 1) * chunk_size - 1;
        pthread_create(&threads[i], NULL, multiply, (void *) &data[i]);
    }

    for (int i = 0; i < num_threads; i++) {
        pthread_join(threads[i], NULL);
    }

    printf("Output vector: \n");
    for (int i = 0; i < row; i++) {
        printf("%d\n", y[i]);
    }

    return 0;
}

该程序首先要求用户输入矩阵和向量的大小及其元素。然后，它会要求用户输入线程的数量，创建线程并将它们分配到矩阵行的子集中。每个线程计算其分配的行的乘积，并将结果存储在输出向量中。最后，主线程打印输出向量。

在该程序中，`struct thread_data` 结构体用于传递线程的起始行和终止行。该结构体作为参数传递给 `multiply` 函数。函数计算其分配的行的乘积，并将结果存储在输出向量中。当线程完成计算时，它将退出。

为了创建线程，程序使用 `pthread_create` 函数。该函数接受四个参数：指向线程标识符的指针，线程属性的指针（通常设置为 `NULL`），指向线程函数的指针和函数参数的指针。在本例中，线程函数是 `multiply`，并且函数参数是指向 `struct thread_data` 结构体的指针。

为了等待所有线程完成，程序使用 `pthread_join` 函数。该函数接受两个参数：要等待的线程的标识符和指向线程的返回值的指针（通常设置为 `NULL`）。

最后，程序打印输出向量的内容。