OpenMP: имея полный «для» петли в каждую нить

Question

У меня есть этот код:

#pragma omp parallel 
{ 
    #pragma omp single 
    { 
    for (int i=0; i<given_number; ++i) myBuffer_1[i] = myObject_1->myFunction(); 
    } 

    #pragma omp single 
    { 
    for (int i=0; i<given_number; ++i) myBuffer_2[i] = myObject_2->myFunction(); 
    } 
} 

// and so on... up to 5 or 6 of myObject_x 

// Then I sum up the buffers and do something with them 
float result; 
for (int i=0; i<given_number; ++i) 
    result = myBuffer_1[i] + myBuffer_2[i]; 

// do something with result 

Если я запускаю этот код, я получаю то, что я ожидал, но загрузка процессора выглядит довольно высокой. Вместо этого, если я запускаю его нормально без OpenMP, я получаю те же результаты, но использование ЦП намного ниже, несмотря на то, что работает в одном потоке.

Я не хочу указывать количество потоков, я хочу, чтобы программа выбирала максимальное количество потоков в соответствии с возможностями ЦП, но я хочу, чтобы каждый цикл for полностью запускался в своем потоке. Как я могу это сделать?

Кроме того, мое ожидание состоит в том, что цикл for для myBuffer_1 запускает поток, другой цикл for запускает другой поток, а остальное выполняется в потоке «мастер». Это верно?

Answer 1

1. #pragma omp single имеет неявный барьер в конце концов, вам нужно использовать #pragma omp single nowait, если вы хотите, чтобы два single блока могут работать одновременно.

2. Однако, для вашего требования, используя section может быть лучшей идеей

   #pragma omp parallel 
   { 
       #pragma omp sections 
       { 
        #pragma omp section 
        { 
         for (int i=0; i<given_number; ++i) myBuffer_1[i] = myObject_1->myFunction(); 
        } 
        #pragma omp section 
        { 
         for (int i=0; i<given_number; ++i) myBuffer_2[i] = myObject_2->myFunction(); 
        } 
       } 
   
   }