73 lines
1.9 KiB
Fortran
73 lines
1.9 KiB
Fortran
! { dg-do run }
|
|
!$ use omp_lib
|
|
call test_workshare
|
|
|
|
contains
|
|
subroutine test_workshare
|
|
integer :: i, j, k, l, m
|
|
double precision, dimension (64) :: d, e
|
|
integer, dimension (10) :: f, g
|
|
integer, dimension (16, 16) :: a, b, c
|
|
integer, dimension (16) :: n
|
|
d(:) = 1
|
|
e = 7
|
|
f = 10
|
|
l = 256
|
|
m = 512
|
|
g(1:3) = -1
|
|
g(4:6) = 0
|
|
g(7:8) = 5
|
|
g(9:10) = 10
|
|
forall (i = 1:16, j = 1:16) a (i, j) = i * 16 + j
|
|
forall (j = 1:16) n (j) = j
|
|
!$omp parallel num_threads (4) private (j, k)
|
|
!$omp barrier
|
|
!$omp workshare
|
|
i = 6
|
|
e(:) = d(:)
|
|
where (g .lt. 0)
|
|
f = 100
|
|
elsewhere (g .eq. 0)
|
|
f = 200 + f
|
|
elsewhere
|
|
where (g .gt. 6) f = f + sum (g)
|
|
f = 300 + f
|
|
end where
|
|
where (f .gt. 210) g = 0
|
|
!$omp end workshare nowait
|
|
!$omp workshare
|
|
forall (j = 1:16, k = 1:16) b (k, j) = a (j, k)
|
|
forall (k = 1:16) c (k, 1:16) = a (1:16, k)
|
|
forall (j = 2:16, n (17 - j) / 4 * 4 .ne. n (17 - j))
|
|
n (j) = n (j - 1) * n (j)
|
|
end forall
|
|
!$omp endworkshare
|
|
!$omp workshare
|
|
!$omp atomic
|
|
i = i + 8 + 6
|
|
!$omp critical
|
|
!$omp critical (critical_foox)
|
|
l = 128
|
|
!$omp end critical (critical_foox)
|
|
!$omp endcritical
|
|
!$omp parallel num_threads (2)
|
|
!$ if (omp_get_thread_num () .eq. 0) m = omp_get_num_threads ()
|
|
!$omp atomic
|
|
l = 1 + l
|
|
!$omp end parallel
|
|
!$omp end workshare
|
|
!$omp end parallel
|
|
|
|
if (any (f .ne. (/100, 100, 100, 210, 210, 210, 310, 310, 337, 337/))) &
|
|
& call abort
|
|
if (any (g .ne. (/-1, -1, -1, 0, 0, 0, 0, 0, 0, 0/))) call abort
|
|
if (i .ne. 20) call abort
|
|
!$ if (l .ne. 128 + m) call abort
|
|
if (any (d .ne. 1 .or. e .ne. 1)) call abort
|
|
if (any (b .ne. transpose (a))) call abort
|
|
if (any (c .ne. b)) call abort
|
|
if (any (n .ne. (/1, 2, 6, 12, 5, 30, 42, 56, 9, 90, &
|
|
& 110, 132, 13, 182, 210, 240/))) call abort
|
|
end subroutine test_workshare
|
|
end
|