1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
! { dg-do run }
!
MODULE reduction_test
CONTAINS
SUBROUTINE reduction_kernel(x_min,x_max,y_min,y_max,arr,sum)
IMPLICIT NONE
INTEGER :: x_min,x_max,y_min,y_max
REAL(KIND=8), DIMENSION(x_min-2:x_max+2,y_min-2:y_max+2) :: arr
REAL(KIND=8) :: sum
INTEGER :: j,k
!$ACC DATA PRESENT(arr) COPY(sum)
!$ACC PARALLEL LOOP REDUCTION(+ : sum)
DO k=y_min,y_max
DO j=x_min,x_max
sum=sum+arr(j,k)
ENDDO
ENDDO
!$ACC END PARALLEL LOOP
!$ACC END DATA
END SUBROUTINE reduction_kernel
END MODULE reduction_test
program main
use reduction_test
integer :: x_min,x_max,y_min,y_max
real(kind=8), dimension(1:10,1:10) :: arr
real(kind=8) :: sum
x_min = 1
x_max = 2
y_min = 1
y_max = 2
arr(:,:) = 1.0
sum = 0.0
!$acc data copy(arr)
call reduction_kernel(x_min,x_max,y_min,y_max,arr,sum)
!$acc end data
if (sum .ne. 4.0) STOP 1
end program
|