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[ и только приращение x для присвоенных ему значений. Таким образом, общее число выполнений x++ во всех потоках составит N.

т.е. при условии, что три потока будут запускать тело цикла для i=0,1,6,7,12,13,..., второй поток для i=2,3,8,9,... и третий для i=4,5,10,11,....