замена #pragma OMP критической (C++)

Question

Я использую OpenMP на вложенном цикле, который работает как этот

#pragma omp parallel shared(vector1) private(i,j) 
{ 
#pragma omp for schedule(dynamic) 
    for (i = 0; i < vector1.size(); ++i){ 

     //some code here 

     for (j = 0; j < vector1.size(); ++j){ 

      //some other code goes here 
     #pragma omp critical 
     A+=B; 
     } 
    C +=A; 
    } 
} 

Проблема здесь состоит в том, что мой код делает много вычислений в A+=B части код. Поэтому, делая это критическим, я не достигаю ускорения, которое я бы хотел. (На самом деле, как представляется, некоторые накладные расходы, так как моя программа занимает больше времени, чтобы выполнить ее, она записывается последовательно).

Я попытался с помощью

#pragma omp reduction private(B) reduction(+:A) 
    A+=B 

это ускоряет время выполнения, однако, кажется, что он не заботится о гоночных условиях, как п critical, так как я не получаю те же результаты А.

Есть ли альтернатива этому, я могу попробовать?

Answer 1

Если вы не хотите, чтобы ваш класс Vector3 был потокобезопасным или переписал ваши операции для использования с std::atomic<Vector3>, оба из которых по-прежнему будут иметь недостатки в производительности (хотя это и не так серьезно, как использование критической секции) , вы можете имитировать behaviour of OpenMP reduction:

#pragma omp parallel // no need to declare variables declared outside/inside as shared/private 
{ 

    Vector3 A{}, LocalC{}; // both thread-private 

    #pragma omp for schedule(dynamic) 
    for (i = 0; i < vector1.size(); ++i){ 

     //some code here 

     for (j = 0; j < vector1.size(); ++j){ 

      //some other code goes here 

      A += B; // does not need a barrier 
     } 
     LocalC += A; // does not need a barrier 
    } 

    #pragma omp critical 
    C += LocalC; 
} 

NB, что это предполагает, что у вас нет доступа A для чтения в ваших «некоторый код» комментарии, но не все равно следует, если вы когда-нибудь думали о использовании reduction пункт.