OpenMP Segfault

Question

我想將OpenMP並行化添加到工作代碼中（僅適用於單個for循環），但是我無法擺脫分段錯誤。該問題起因於這一行：

pos += sprintf(com + pos, "%d ", i); 

com是一個字符陣列，我試圖將其定義爲char com[255]或char *com = malloc(255*sizeof(char))，內部和for循環之前。當我在循環之前定義com時，我添加了private(com)到#pragma omp parallel for指令。我也嘗試初始化它，並使用firstprivate。 （pos是一個整數，初始化爲0）

當我不加-fopenmp一切工作正常，但-fopenmp它給段錯誤。我錯過了什麼？

Answer 1

段故障來自多個線程同時更新的pos值，因此將其設置爲某個值，輪流com + pos成指向超出或com所分配的存儲器之前的指針。到parallelise這樣一個循環的正確方法是串聯在私人字符串值，然後串聯私人字符串有序方式：

char com[255]; 
int pos = 0; 

#pragma omp parallel 
{ 
    char mycom[255]; 
    int mypos = 0; 

    #pragma omp for schedule(static) nowait 
    for (int i = 0; i < N; i++) 
     mypos += sprintf(mycom + mypos, "%d ", i); 

    // Concatenate the strings in an ordered fashion 
    #pragma omp for schedule(static) ordered 
    for (int i = 0; i < omp_get_num_threads(); i++) 
    { 
     #pragma omp ordered 
     pos += sprintf(com + pos, "%s", mycom); 
    } 
} 

的ordered結構可確保正確的同步，以便一個不需要critical。使用schedule(static)對於確保每個線程處理迭代空間的單個連續部分非常重要。