Я использую две разные версии сокращения в openmp и получаю совершенно разные результаты. Что из следующего не так?
omp_set_num_threads(t);
long long unsigned int d = 0;
#pragma omp parallel for default(none) shared(some_stuff) reduction(+:d)
for (int i=start; i< n; i++)
{
d += calc(i,some_stuff);
}
cout << d << endl;
и вторая версия такова:
omp_set_num_threads(t);
//reduction array
long long unsigned int* d = new long long unsigned int[t];
for(int i = 0; i < t; i++)
d[i] = 0;
#pragma omp parallel for default(none) shared(somestuff, d)
for (int i=start; i< n; i++)
{
long long unsigned dd = calc(i, somestuff);
d[omp_get_thread_num()] += dd;
}
long long unsigned int res = 0;
for(int i = 0; i < omp_get_num_threads(); i++){
res += d[i];
}
delete[] d;
cout << res << endl;
Второй код неверен. omp_get_num_threads()
возвращается 1
при вызове вне параллельной области, и, следовательно, ваш код не сводит все значения в конечный результат. Поскольку вы явно фиксируете количество потоков, которые будут t
вместо этого вы должны использовать:
for(int i = 0; i < t; i++){
res += d[i];
}
В качестве альтернативы, вы можете использовать omp_get_max_threads()
,