Всем сяп! Подскажите пж, не понимаю, что за бред. Есть функция

Question

Всем сяп! Подскажите пж, не понимаю, что за бред. Есть функция

умножения матрицы (вектор векторов) на вектор, в котором при паралельном выполнении ловлю сегфолт. Если запускаю последовательно - все ок. В чем проблема? Выход из индексов не вижу

vector<double> matrix_vec_mul(vector<vector<double>> m, vector<double> b) {
vector<double> res;
res.reserve(b.size());

#pragma omp parallel for shared(b, m, res) default(none)
for (int i = 0; i < b.size(); ++i) {
double row_res = 0;
for (int j = 0; j < b.size(); ++j) {
row_res += m.at(i).at(j) * b.at(j);
}
res.insert(res.begin() + i, row_res);
}

return res;
}

#c #embedded #programming #russian

0

15.06.2022

13 ответов

21 просмотр

Max Azatian Автор вопроса

Ilya Zviagin
А ты точно хочешь копировать матрицы при вызове фу...

это второй вопрос, его уже тоже увидел

0

15.06.2022

Ilya Zviagin

У тебя с постановкой задачи нет тут проблем ?

0

15.06.2022

Ilya Zviagin

Так а сегфолт же поймал - можно сразу в дебагере глянуть, что там и как...

0

15.06.2022

Max Azatian Автор вопроса

Ilya Zviagin
Так а сегфолт же поймал - можно сразу в дебагере ...

В дебаг версии CLion запускаю - только out of bounds, а где - хрен пойми

0

15.06.2022

Ilya Zviagin

res тут не может быть shared никак...

0

15.06.2022

Ilya Zviagin

Ты тут будешь параллельно вставлять в res. Это недопустимо.

0

15.06.2022

Max Azatian Автор вопроса

Ilya Zviagin
Ты тут будешь параллельно вставлять в res. Это нед...

Так по разным индексам ж

0

15.06.2022

Ilya Zviagin

Max Azatian
Так по разным индексам ж

Где ?

0

15.06.2022

Max Azatian Автор вопроса

Ilya Zviagin
Где ?

Res, в конце Мы ж I-loop параллелим

0

15.06.2022

Ilya Zviagin

res.insert(res.begin() + i, row_res); эту операцию нельзя делать параллельно

0

15.06.2022

Max Azatian Автор вопроса

Ilya Zviagin
res.insert(res.begin() + i, row_res); эту операц...

А res[i] = row_res; ?

0

15.06.2022

Ilya Zviagin

Max Azatian
Res, в конце Мы ж I-loop параллелим

#pragma omp parallel for shared(b, m, res) default(none) for (int i = 0; i < b.size(); ++i) { double row_res = 0; for (int j = 0; j < b.size(); ++j) { row_res += m.at(i).at(j) * b.at(j); } res.insert(res.begin() + i, row_res); } Вот это параллелицо

0

15.06.2022

Ilya Zviagin · Accepted Answer

Ilya Zviagin

А ты точно хочешь копировать матрицы при вызове фукнции ?

0

15.06.2022

170 похожих чатов

Всем сяп! Подскажите пж, не понимаю, что за бред. Есть функция

13 ответов

Похожие вопросы