OpenMP缩减,变量不是私有的?

我有一个这样的数组(0,0在左下角):

0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0
0 0 1 0 1 0 1 0 0
0 0 1 1 1 1 1 1 1
1 0 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1

我的目标是得到更高线的指数 0 . 为此,我制定了以下代码(工作正常):

max=0;
for (i=0 ; i<width ; ++i) {
  for (j=max ; j<height ; ++j) {
    if (array[i*height+j]!=0) {
      max=j;
    }
  }
}

对于第二个循环i初始化 j MAX,因为全局最大值不能小于局部最大值。这样我可以减少测试次数。

我试图把它与之并行化。 OpenMp . 我的代码是:

max=0;
#pragma omp parallel for  default(none)                 \
                          shared(spec, width, height)   \
                          collapse(2)                   \
                          reduction(max:max)
for (i=0 ; i<width ; ++i) {
  for (j=max ; j<height ; ++j) {
    if (array[i*height+j]!=0) {
      max=j;
    }
  }
}

这会导致分割错误。为了让它起作用,我改变了 j=max 到 j=0 . 所以问题似乎来自 max 变量。

我不明白为什么,因为随着减少,这个变量在每个线程之间应该是私有的(或lastprivate)。为什么它会崩溃?我怎样才能用OpenMP来使用我的“优化”呢?