c - openmp 没有正确地减少

Question

我想使用 openmp 缩减运算符添加循环进行的所有迭代。

这是我的代码:

#include <omp.h>
#include <stdio.h>
#include <stdlib.h>
#define CHUNKSIZE   2
#define N       100

int main (int argc, char *argv[]) 
{
int nthreads, tid, i, chunk;
float a[N], b[N], c[N];

/* Some initializations */
for (i=0; i < N; i++)
  a[i] = b[i] = i * 1.0;
chunk = CHUNKSIZE;


int x=0;

#pragma omp parallel shared(a,b,c,nthreads,chunk) private(i,tid) reduction(+ : x)
  {
  tid = omp_get_thread_num();
  if (tid == 0)
    {
    nthreads = omp_get_num_threads();
    printf("Number of threads = %d\n", nthreads);
    }
  printf("Thread %d starting...\n",tid);

  #pragma omp for schedule(static,chunk)
  for (i=0; i<N; i++)
    {
    c[i] = a[i] + b[i];
    printf("Thread %d: c[%d]= %f\n",tid,i,c[i]);
    x++;
    }

  }  /* end of parallel section */


printf("Value of x is %d" , x);
}

问题是 x 的最终值是 100，而不是 200。我不明白为什么我没有得到 200 的预期值。有人可以帮我吗？

Answer 1

100 是预期的结果:

#pragma omp for schedule(static,chunk)
for (i=0; i<N; i++) {
  ...
}

每个线程将从 [0..N[ 间隔中获取一组双元素 block ，并且只为其分配的值递增 x。因此 x++ 将在所有线程中执行的总次数为 N。

即假设三个线程，第一个将运行 i=0,1,6,7,12,13,... 的循环体，第二个线程将运行 i=2,3,8， 9,... 和第三个 i=4,5,10,11,...。