使用Pthreads進行動態矩陣乘法

Question

我是一名初學者，主要使用線程編程和C，我試圖找出如何使用Pthreads進行簡單的矩陣乘法。我想爲每列創建一個線程，並將結果放入結果矩陣中。我試圖動態地做到這一點，這意味着用戶可以使用輸入作爲大小來創建兩個n×n矩陣。

我的代碼，現在，不包括填充所述矩陣和讀取大小n爲執行以下操作：

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

typedef struct Matrix { 
int line, col, size; 
double (*MA)[]; 
double (*MB)[]; 
double (*MC)[]; 
} Matrix; 


void *multiply(void *arg) { 
Matrix* work = (Matrix*) arg; 
int s, z; 
s = work->col; 
z = work->line; 
work->MC[0][0] = 0.0.//can't use MC, MB, MA here!! 
return 0; 
} 

int main() { 
Matrix* m; 
//read size and set it to int size (miissing here, does work) 

double MA[size][size], MB[size][size], MC[size][size]; 
int i, j; 
//filling the matrices (missing here, does work) 

pthread_t threads[size]; 

for (i = 0; i < size; i++) { 
    m = malloc(sizeof(Matrix*)); 
    m->size = size; 
    m->col = i; 
    pthread_create(&threads[i], NULL, multiply, m); 

} 

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

} 

的問題是，我不能使用既不MA，MB，也不NC（：=結果）在乘法與代碼中顯示的一樣。 即使我在主要方法中聲明瞭所有這三個元素，我也只是得到錯誤「數組的非特定邊界無效使用」。

我瞭解這裏有什麼不對嗎，或者我該如何解決這個問題？我試圖修改我講座的一個例子，在那裏爲每個元素創建一個線程。 在此先感謝！

Answer 1

就有關該錯誤：

work->MC[0][0] = 0.0.//can't use MC, MB, MA here!! 

MC被宣佈爲double (*MC)[]，並嘗試使用它作爲一個二維數組像你已經宣佈它double MC[N]{M]。您可以使用兩個（或更多）維數組，如果且僅當第一維是固定的，或者如果您逐行分配它。

所以，你的程序可能是：

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

typedef struct Matrix { 
    int line, col, size; 
    double MA[][]; 
    double MB[][]; 
    double MC[][]; 
} Matrix; 

void *multiply(void *arg) { 
    Matrix* work = (Matrix*) arg; 
    int s, z; 
    s = work->col; 
    z = work->line; 
    work->MC[0][0] = 0.0 
    return 0; 
} 

int main() { 
    Matrix* m; 
    //read size and set it to int size (miissing here, does work) 

    double MA[][], MB[][], MC[][]; 
    int i, j; 
    pthread_t threads[size]; 

    MA = (double **) malloc(size * sizeof(double *)); 
    MB = (double **) malloc(size * sizeof(double *)); 
    MC = (double **) malloc(size * sizeof(double *)); 
    for(int i=0;i<size;++i){ 
     MA[i] = (double *) malloc(size * sizeof(double)); 
     MB[i] = (double *) malloc(size * sizeof(double)); 
     MC[i] = (double *) malloc(size * sizeof(double)); 
    } 

    for (i = 0; i < size; i++) { 
     m = malloc(sizeof(Matrix*)); 
     m->MA = MA; 
     m->MB = MB; 
     m->MC = MC; 
     m->size = size; 
     m->col = i; 
     pthread_create(&threads[i], NULL, multiply, m); 
    } 

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

但你必須要小心的是，線程可以訪問到的數據同時，所以你應該使用一些鎖，如果不同的線程可以使用和更改相同的值。