упражнение openMP omp_bug2.c

это упражнение с сайта OpenMP:
https://computing.llnl.gov/tutorials/openMP/exercise.html

#include "stdafx.h"#include <omp.h>
#include <stdio.h>
#include <stdlib.h>

int _tmain(int argc, _TCHAR* argv[])
{
int nthreads, i, tid;
float total;

/*** Spawn parallel region ***/
#pragma omp parallel private(i, tid) // i changed this line
{
/* Obtain thread number */
tid = omp_get_thread_num();
/* Only master thread does this */
if (tid == 0) {
nthreads = omp_get_num_threads();
printf("Number of threads = %d\n", nthreads);
}
printf("Thread %d is starting...\n",tid);

#pragma omp barrier

/* do some work */
total = 0.0;
#pragma omp for schedule(dynamic,10)
for (i=0; i<1000000; i++)
total = total + i*1.0;

printf ("Thread %d is done! Total= %e\n",tid,total);

}
}

выход для этого

Number of threads = 4
Thread 0 is starting...
Thread 3 is starting...
Thread 2 is starting...
Thread 1 is starting...
Thread 0 is done! Total= 0.000000e+000
Thread 3 is done! Total= 0.000000e+000
Thread 2 is done! Total= 0.000000e+000
Thread 1 is done! Total= 0.000000e+000

Это означает, что у нас есть проблема с переменной «Всего»

это помощь на сайте введите описание изображения здесь

Вот мое решение: как вы думаете, это правильный способ сделать это?

#include "stdafx.h"#include <omp.h>
#include <stdio.h>
#include <stdlib.h>

int _tmain(int argc, _TCHAR* argv[])
{
int nthreads, i, tid;
float total;

/*** Spawn parallel region ***/
#pragma omp parallel private(total,tid)
{
/* Obtain thread number */
tid = omp_get_thread_num();
total= 0.0;
/* Only master thread does this */
if (tid == 0) {
nthreads = omp_get_num_threads();
printf("Number of threads = %d\n", nthreads);
}
printf("Thread %d is starting...\n",tid);#pragma omp parallel for schedule(static,10)\
private(i)\
reduction(+:total)
for (i=0; i<1000000; i++)
total = total + i*1.0;

printf ("Thread %d is done! Total= %e\n",tid,total);

} /*** End of parallel region ***/
}

Вот мой новый вывод:

Number of threads = 4
Thread 0 is starting...
Thread 1 is starting...
Thread 0 is done! Total= 4.999404e+011
Thread 2 is starting...
Thread 1 is done! Total= 4.999404e+011
Thread 2 is done! Total= 4.999404e+011
Thread 3 is starting...
Thread 3 is done! Total= 4.999404e+011

0

Решение

Да ты конечно хочешь total быть частной переменной потока. Одна вещь, которую вы, вероятно, сделали бы в реальном примере, это уменьшить частный поток totals до единого глобального итога в конце (и только тогда один поток может напечатать результат). Одним из способов сделать это является простой

#pragma omp atomic
global_total += total

в конце (есть лучшие способы, хотя и с использованием сокращений).

PS: счетчики петель для omp for по умолчанию являются частными, поэтому вам не нужно явно указывать это.

1

Другие решения

Других решений пока нет …

По вопросам рекламы ammmcru@yandex.ru
Adblock
detector