-
Notifications
You must be signed in to change notification settings - Fork 16
/
Copy pathserial_firstprivate.F90
138 lines (126 loc) · 2.93 KB
/
serial_firstprivate.F90
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
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
#ifndef T1
!T1:serial,firstprivate,V:2.6-2.7
LOGICAL FUNCTION test1()
IMPLICIT NONE
INCLUDE "acc_testsuite.Fh"
INTEGER:: errors
REAL(8),DIMENSION(10, LOOPCOUNT):: a, b, d
REAL(8),DIMENSION(10):: c, c_copy
INTEGER:: x, y
errors = 0
SEEDDIM(1) = 1
# ifdef SEED
SEEDDIM(1) = SEED
# endif
CALL RANDOM_SEED(PUT=SEEDDIM)
CALL RANDOM_NUMBER(a)
CALL RANDOM_NUMBER(b)
CALL RANDOM_NUMBER(c)
c_copy = c
d = 0
!$acc data copyin(a(1:10, 1:LOOPCOUNT), b(1:10, 1:LOOPCOUNT)) copy(d(1:10, 1:LOOPCOUNT))
!$acc serial firstprivate(c(1:10))
!$acc loop gang
DO y = 1, LOOPCOUNT
!$acc loop worker
DO x = 1, 10
d(x, y) = a(x, y) + b(x, y) + c(x)
END DO
END DO
!$acc end serial
!$acc end data
DO y = 1, LOOPCOUNT
DO x = 1, 10
IF (abs(d(x, y) - (a(x, y) + b(x, y) + c_copy(x))) .gt. PRECISION) THEN
errors = errors + 1
END IF
END DO
END DO
IF (errors .eq. 0) THEN
test1 = .FALSE.
ELSE
test1 = .TRUE.
END IF
END
#endif
#ifndef T2
!T2:serial,firstprivate,V:2.6-2.7
LOGICAL FUNCTION test2()
IMPLICIT NONE
INCLUDE "acc_testsuite.Fh"
INTEGER:: errors
REAL(8),DIMENSION(10, LOOPCOUNT):: a, b, d
REAL(8),DIMENSION(10):: c, c_copy
INTEGER:: x, y
errors = 0
SEEDDIM(1) = 1
# ifdef SEED
SEEDDIM(1) = SEED
# endif
CALL RANDOM_SEED(PUT=SEEDDIM)
CALL RANDOM_NUMBER(a)
CALL RANDOM_NUMBER(b)
d = 0
c = 0
!$acc data copyin(a(1:10, 1:LOOPCOUNT), b(1:10, 1:LOOPCOUNT)) copy(d(1:10, 1:LOOPCOUNT))
!$acc serial firstprivate(c(1:10))
!$acc loop gang independent
DO y = 1, LOOPCOUNT
!$acc loop worker independent
DO x = 1, 10
c(x) = a(x, y) - b(x, y)
END DO
!$acc loop worker independent
DO x = 1, 10
d(x, y) = a(x, y) + b(x, y) + c(x)
END DO
END DO
!$acc end serial
!$acc end data
DO y = 1, LOOPCOUNT
DO x = 1, 10
IF (abs(d(x, y) - (2 * a(x, y))) .gt. PRECISION) THEN
errors = errors + 1
END IF
END DO
END DO
IF (errors .eq. 0) THEN
test2 = .FALSE.
ELSE
test2 = .TRUE.
END IF
END
#endif
PROGRAM main
IMPLICIT NONE
INTEGER :: failcode, testrun
LOGICAL :: failed
INCLUDE "acc_testsuite.Fh"
#ifndef T1
LOGICAL :: test1
#endif
#ifndef T2
LOGICAL :: test2
#endif
failed = .FALSE.
failcode = 0
#ifndef T1
DO testrun = 1, NUM_TEST_CALLS
failed = failed .or. test1()
END DO
IF (failed) THEN
failcode = failcode + 2 ** 0
failed = .FALSE.
END IF
#endif
#ifndef T2
DO testrun = 1, NUM_TEST_CALLS
failed = failed .or. test2()
END DO
IF (failed) THEN
failcode = failcode + 2 ** 1
failed = .FALSE.
END IF
#endif
CALL EXIT (failcode)
END PROGRAM