为什么 OpenMP atomic 和 critical 没有给出正确的结果？

Question

我编写了以下 Fortran 代码来测试 atomic 和 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

输出为

表示atomic和critical的结果是错误的。这很奇怪，我认为添加它们可以避免赛车状况。然而，没有同步的第一个循环给出了正确的答案，这里没有比赛吗？我的代码有什么问题？

Answer 1

您的代码中的问题是竞争条件

   !$omp parallel

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

所有线程都执行该操作并且它们发生冲突。您需要 omp single 围绕这些行。

您不需要 atomic 或 critical

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

因为每个线程对不同的数组元素进行操作。

在您的 OpenMP 规范示例中，问题在于

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

数组或函数 INDEX(I) 可以 return 具有不同 I 的两个不同线程的相同值，您必须保护这种潜在的竞争条件。

为什么 OpenMP atomic 和 critical 没有给出正确的结果？

Why OpenMP atomic and critical don't give the right result?

fortran

openmp