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
    26 мая в 10:30
  • @megorit нужно через openMP
    – Freakgen
    26 мая в 13:59
  • Всякие такие штуки больше годятся для распределенных вычислений, т.е. по сети.
    – megorit
    26 мая в 14:05

1 ответ 1

0

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

Ваш ответ

By clicking “Отправить ответ”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

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