OpenMP. Нетривиальный случай

@prozaeck · Регистрация: 26.01.2019

Author24 — интернет-сервис помощи студентам

Аналог реальной программы. Как переместить создание потоков (#pragma omp parallel) в main_, но оставить funcSerial однопоточным (последовательным)? funcSerial я распараллеливать пробовал - толку мало. В реальной программе там идет очень интенсивный обмен с памятью, и более чем в два раза она не ускоряется на любом количестве ядер.

C++

double funcPar()
{
    double sum[] = { 0, 0, 0, 0 };
#pragma omp parallel
    {
#pragma omp for
        for (int i = 0; i < 4; i++)  //loop body should be parallel (multy-threaded)
        {
            auto a = (double)i;
            sum[i] += sin(a)*cos(a); // real code is more complicated and works much longer
        }
    }//omp parallel
    return sum[0] + sum[1] + sum[2] + sum[3];
}
double funcSerial() //function should remain serial (single-threaded)
{
    int K = 5000;
    double sum = 0.0;
    for (int k = 0; k < K; k++)
    {
        sum += funcPar();
    }
    return sum;
}
void main_() // //function body should remain serial (single-threaded)
{
    int J = 1000;
    double sum = 0.0;
    for (int j = 0; j < J; j++)
    {
        sum += funcSerial();
    }
    printf("\n\n sum = %f \n", sum);
}

Добавлено через 28 минут
А можно мне рыбу вместо удочек? Как то не получилось у меня с рыбалкой на этот раз. Пробовал я #pragma omp single, но кроме зависаний ничего не получил.

@mat_for_c · 27.01.2019, 11:44

C++

#pragma omp parallel for reduction(+: sum)
for (int j = 0; j < J; j++)
{
   sum += funcSerial();
}

@prozaeck · 27.01.2019, 14:37 **[ТС]**

Это не верно. Такой вариант я делал. Он приводит к параллельной работе funcSerial, а поскольку в ней чтение и запись больших объемов данных, то шина памяти перегружается уже на трех потоках. Ускорение максимум в два раза. Можно сформулировать задачу так: нужно #pragma omp parallel перенести в main, но #pragma omp for оставить в funcPar, дописав в main и funcSerial некий код ( возможно основанный на использовании #pragma omp serial ) который оставит функции main и funcSerial однопоточными.

@mat_for_c · 27.01.2019, 16:28

Сообщение от prozaeck

но #pragma omp for оставить в funcPar

а зачем вам это?
ваша функция funcPar() эквивалентна следующему коду:

C++

double funcPar()
{
   return (1 + sin(2.) + sin(4.) + sin(6.)) / 2.;
}

а если соптимизировать и остальные циклы, то последовательный вариант будет значительно быстрее параллельной реализации с циклами

@prozaeck · 27.01.2019, 16:41 **[ТС]**

Вы вообще читаете, что я писал в описании? Написано же - АНАЛОГ реальной программы. А также "// real code is more complicated and works much longer"

@prozaeck 0 / 0 / 0 Регистрация: 26.01.2019 Сообщений: 3
	27.01.2019, 14:37 [ТС]	3
	Это не верно. Такой вариант я делал. Он приводит к параллельной работе funcSerial, а поскольку в ней чтение и запись больших объемов данных, то шина памяти перегружается уже на трех потоках. Ускорение максимум в два раза. Можно сформулировать задачу так: нужно #pragma omp parallel перенести в main, но #pragma omp for оставить в funcPar, дописав в main и funcSerial некий код ( возможно основанный на использовании #pragma omp serial ) который оставит функции main и funcSerial однопоточными. 0

@prozaeck 0 / 0 / 0 Регистрация: 26.01.2019 Сообщений: 3
	27.01.2019, 16:41 [ТС]	5
	Вы вообще читаете, что я писал в описании? Написано же - АНАЛОГ реальной программы. А также "// real code is more complicated and works much longer" 0