Neden OpenMP atomik ve kritik doğru sonucu vermiyor?

Question

Ben atomic sınamak için aşağıdaki Fortran kod yazdım

ve critical

program test 
    implicit none 
    integer::i 
    integer::a(10),b(10),atmp(10),btmp(10) 
    a=[1,2,3,4,5,6,7,8,9,10] 
    b=[12,32,54,77,32,19,34,1,75,45] 
    atmp=a 
    btmp=b 
    write(*,'(1X,10I4)') a+b 
    print*,'------------------' 
    !$omp parallel 
    !$omp do 
    do i=1,10 
     B(I) = B(I)+A(I) 
    end do 
    !$omp end do 
    !$omp single 
    write(*,'(1X,10I4)') b 
    !$omp end single 

    a=atmp 
    b=btmp 
    !$omp do 
    do i=1,10 
     !$omp critical 
     B(I) = B(I)+A(I) 
     !$omp end critical 
    end do 
    !$omp end do 
    !$omp single 
    write(*,'(1X,10I4)') b 
    !$omp end single 

    a=atmp 
    b=btmp 
    !$omp do 
    do i=1,10 
     !$omp atomic 
     B(I) = B(I)+A(I) 
     !$omp end atomic 
    end do 
    !$omp end do 
    !$omp single 
    write(*,'(1X,10I4)') b 
    !$omp end single 

    !$omp end parallel 
end program 

çıktı O atomic ve critical yanlıştır o sonucunu demektir

olduğunu. Bu garip, onlara yarış koşullarını engelleyebileceklerini düşündüm. Ancak, senkronizasyon olmadan ilk döngü doğru cevap verir, burada yarış yok mu? Kodumdaki sorun nedir?

Answer 1

Kodunuzdaki sorun yarış durumu

!$omp parallel 

... 
    a=atmp 
    b=btmp 
... 
    !$omp end parallel 

bütün konuları bu işlemi yapmak ve onlar çatışma olduğunu. Bu satırların etrafında omp single olmasını istiyorsunuz. Her iş parçacığı farklı bir dizi öğesi çalışır çünkü

Sen

!$omp do 
do i=1,10 
    B(I) = B(I)+A(I) 
end do 
!$omp end do 

herhangi atomic veya critical gerekmez. OpenMP şartname adresinin örnekte

konu

!$OMP PARALLEL DO SHARED(X, Y, INDEX, N) 
DO I=1,N 
    !$OMP  ATOMIC UPDATE 
    X(INDEX(I)) = X(INDEX(I)) + WORK1(I) 

yılında dizi veya işlev INDEX(I) farklı I ile iki farklı diş için aynı değeri döndürebilir ve bu potansiyel yarış durumu korumak gerektiğidir.