Можно ли уменьшить массив с помощью openmp?

Поддерживает ли OpenMP изначально сокращение переменной, представляющей массив?

Это могло бы работать примерно так ...

float* a = (float*) calloc(4*sizeof(float));
omp_set_num_threads(13);
#pragma omp parallel reduction(+:a)
for(i=0;i<4;i++){
   a[i] += 1;  // Thread-local copy of a incremented by something interesting
}
// a now contains [13 13 13 13]

В идеале, было бы что-то подобное для параллельного OMP, и если у вас достаточно большое количество потоков, чтобы это имело смысл, накопление происходило бы через двоичное дерево.