593 / 531 / 76
Регистрация: 22.03.2011
Сообщений: 1,585
1

Распараллелить скалярное произведение векторов

08.11.2012, 23:22. Показов 8610. Ответов 4

Задание - распараллелить скалярное произведение векторов.
вектора vec1 и vec2 заполняются не рандомом для однозначности результата.

в общем все работает,но препод впихнул симуляцию задержки 1го потока (в нашем случае который выполняет 10ую итерацию), и все накрывается. не знаю как это исправить (место обозначено (!!!)).
пробывал поставить барьеры - не работает.
как решить проблему? нельзя убирать слип.
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
const int N = 20;
 
int res = 0;
short vec1[N], vec2[N];
 
srand(time(0));
 
for (int i=0; i<N; i++) {
    vec1[i] = 1;
    vec2[i] = 2;
    printf("%3d %3d\n",vec1[i], vec2[i]);
}
 
::omp_set_num_threads(omp_get_num_procs()); // ?
 
int chunk = 1;
 
int tmp2, tmp;
 
#pragma omp parallel shared(N, res) private(tmp, tmp2)
 {
     #pragma omp for schedule(dynamic, chunk) nowait
     for (int i=0; i<N; i++) {
 
         tmp = vec1[i] * vec2[i];
         tmp2 = res + tmp;
         if (i == 10)  {          // !!!
             Sleep(1000);
         }
 
         res = tmp2;
     }
 
 }
 
 
 printf("\n\nres: %d\n", res);
спс.
1
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
08.11.2012, 23:22
Ответы с готовыми решениями:

Скалярное произведение векторов
Объясните, пожалуйста, значение вот этой строки: int scalar(int V1, int V2, int n)

Скалярное произведение векторов
Здравствуйте! Помогите пожалуйста найти ошибки #include &lt;iostream&gt; using namespace std; int...

Найти скалярное произведение векторов
Помогите пожалуйста, у самой никак не получается, надо написать программу на С++. Даны координаты...

Скалярное произведение двух векторов
Даны два вектора (одномерных массива), содержащих n вещественных элементов. Найти скалярное ...

4
Эксперт С++
2373 / 1657 / 279
Регистрация: 29.05.2011
Сообщений: 3,387
11.11.2012, 02:20 2
А редукцию использовать тоже нельзя?
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
#pragma omp parallel shared(res) private(tmp)
 {
     #pragma omp for schedule(dynamic, chunk) reduction(+:res) nowait
     for (int i=0; i<N; i++) {
 
         tmp = vec1[i] * vec2[i];
         if (i == 10)  {          // !!!
             Sleep(1000);
         }
 
         res += tmp;
     }
 }
1
593 / 531 / 76
Регистрация: 22.03.2011
Сообщений: 1,585
18.11.2012, 00:39  [ТС] 3
grizlik78, да, я думаю. конкретно в OpenMP с редукцией не работал, и вообще слабо разбираюсь в сабже... мне интересно как это решает проблему торможения потока. ...
0
Эксперт С++
2373 / 1657 / 279
Регистрация: 29.05.2011
Сообщений: 3,387
18.11.2012, 00:59 4
Строго говоря, даже если убрать задержку в исходном варианте, то и тогда программа может выдавать неправильный результат. И в моём случае она выдаёт то 40, то 38, то 36... Это происходит из-за конкурирующего изменения одной переменной. Один из потоков может при вычислении tmp прочитать значение res раньше другого потока, а изменить значение res позже, затирая тем самым результат второго потока. А при наличии задержки это происходит практически гарантировано. Редукция — стандартный механизм для таких случаев, когда каждый поток вычисляет независимо свою частичную сумму, а объединением частичных сумм в полную занимается уже код OpenMP.
1
593 / 531 / 76
Регистрация: 22.03.2011
Сообщений: 1,585
22.11.2012, 21:39  [ТС] 5
grizlik78, Спасибо !
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
22.11.2012, 21:39
Помогаю со студенческими работами здесь

[Геометрия]Скалярное произведение векторов
Добрый день. Есть вот такая задача: Я написал программу, но не уверен верно ли написал. ...

Найти скалярное произведение векторов
нужно написать две программы 1) нахождение скалярного произведения векторов ...

Найти скалярное произведение векторов
Даны два вектора x, y есть Rn. Найти скалярное произведение векторов.

Найти скалярное и векторное произведение векторов
Вот программка для операций над векторами, но почему-то не хочет работать, посмотрите плиз: ...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
5
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2022, CyberForum.ru