OpenMP 段错误

Question

我正在尝试将 OpenMP 并行化添加到工作代码中（仅添加到单个 for 循环中），但是我无法摆脱分段错误。问题出在这一行：

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

com 是一个字符数组，我尝试在 for 循环内部和之前将其定义为 char com[255] 或 char *com = malloc(255*sizeof(char))。当我在循环之前定义 com 时，我将 private(com) 添加到 #pragma omp parallel for 指令。我也尝试初始化它并使用 firstprivate。 (pos为整数，初始化为0)

当我不添加 -fopenmp 时，一切正常，但使用 -fopenmp 时会出现段错误。我错过了什么？

Answer 1

分段错误来自多个线程同时更新 pos 的值，因此将其设置为某个值，将 com + pos 变成指向分配内存之外或之前的指针com。并行化此类循环的正确方法是连接私有字符串中的值，然后以有序的方式连接私有字符串：

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) 的使用很重要，以保证每个线程处理迭代的单个连续部分 space.

OpenMP 段错误

OpenMP Segfault

c

malloc

char

openmp