0

Функция принимает на вход матрицу, и решает СЛАУ. Но время выполнения без omp и с omp не отличается (т.е если комментировать строку с #pragma parallel). В свойствах проекта поддержка openMP включена, если ставлю num_threads = 1, то код замедляется (т.е параллельный алгоритм работает). Время при n = 3000, 4.5 секунды на обоих алгоритмах, а на параллельном должно быть в 2 раза быстрее. В чем может быть ошибка?

void gaussMethod(vector<vector<double>> a, vector<double> b, int n, vector<double> y)
{
    int i, j, k;
    double ratio;
    for (i = 0; i < n - 1; i++) 
    {
#pragma omp parallel for private(j,k,ratio) num_threads(12)
        for (j = i+1; j < n; j++) 
        {
            ratio = a[j][i] / a[i][i];
            for (k = i; k <= n; k++) 
                a[j][k] -= ratio * a[i][k];
            b[j] -= ratio * b[i];
        }
    }
    for (int k = n - 1; k >= 0; k--) 
    {
        double d = 0;
        for (int j = k + 1; j < n; j++)
            d += a[k][j] * y[j];
        y[k] = (b[k] - d) / a[k][k];
    }
}
3
  • std::thread не?
    – megorit
    Commented 26 мая 2023 в 10:30
  • @megorit нужно через openMP
    – Freakgen
    Commented 26 мая 2023 в 13:59
  • Всякие такие штуки больше годятся для распределенных вычислений, т.е. по сети.
    – megorit
    Commented 26 мая 2023 в 14:05

1 ответ 1

0

Оказывается всё дело было в конфигурации. Если запускать в Release, то эффекта вроде как и нет. А если в Debug, то есть ускорении в разы.

Ваш ответ

Нажимая «Отправить ответ», вы соглашаетесь с условиями пользования и подтверждаете, что прочитали политику конфиденциальности.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.