Mercurial > hg > CbC > CbC_gcc

comparison libgomp/testsuite/libgomp.fortran/reduction1.f90 @ 0:a06113de4d67

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
first commit
author kent <kent@cr.ie.u-ryukyu.ac.jp>
date Fri, 17 Jul 2009 14:47:48 +0900
parents
children 84e7813d76e9
comparison
equal deleted inserted replaced
-1:000000000000 0:a06113de4d67
1 ! { dg-do run }
2 !$ use omp_lib
3
4 integer :: i, ia (6), n, cnt
5 real :: r, ra (4)
6 double precision :: d, da (5)
7 complex :: c, ca (3)
8 logical :: v
9
10 i = 1
11 ia = 2
12 r = 3
13 ra = 4
14 d = 5.5
15 da = 6.5
16 c = cmplx (7.5, 1.5)
17 ca = cmplx (8.5, -3.0)
18 v = .false.
19 cnt = -1
20
21 !$omp parallel num_threads (3) private (n) reduction (.or.:v) &
22 !$omp & reduction (+:i, ia, r, ra, d, da, c, ca)
23 !$ if (i .ne. 0 .or. any (ia .ne. 0)) v = .true.
24 !$ if (r .ne. 0 .or. any (ra .ne. 0)) v = .true.
25 !$ if (d .ne. 0 .or. any (da .ne. 0)) v = .true.
26 !$ if (c .ne. cmplx (0) .or. any (ca .ne. cmplx (0))) v = .true.
27 n = omp_get_thread_num ()
28 if (n .eq. 0) then
29 cnt = omp_get_num_threads ()
30 i = 4
31 ia(3:5) = -2
32 r = 5
33 ra(1:2) = 6.5
34 d = -2.5
35 da(2:4) = 8.5
36 c = cmplx (2.5, -3.5)
37 ca(1) = cmplx (4.5, 5)
38 else if (n .eq. 1) then
39 i = 2
40 ia(4:6) = 5
41 r = 1
42 ra(2:4) = -1.5
43 d = 8.5
44 da(1:3) = 2.5
45 c = cmplx (0.5, -3)
46 ca(2:3) = cmplx (-1, 6)
47 else
48 i = 1
49 ia = 1
50 r = -1
51 ra = -1
52 d = 1
53 da = -1
54 c = 1
55 ca = cmplx (-1, 0)
56 end if
57 !$omp end parallel
58 if (v) call abort
59 if (cnt .eq. 3) then
60 if (i .ne. 8 .or. any (ia .ne. (/3, 3, 1, 6, 6, 8/))) call abort
61 if (r .ne. 8 .or. any (ra .ne. (/9.5, 8.0, 1.5, 1.5/))) call abort
62 if (d .ne. 12.5 .or. any (da .ne. (/8.0, 16.5, 16.5, 14.0, 5.5/))) call abort
63 if (c .ne. cmplx (11.5, -5)) call abort
64 if (ca(1) .ne. cmplx (12, 2)) call abort
65 if (ca(2) .ne. cmplx (6.5, 3) .or. ca(2) .ne. ca(3)) call abort
66 end if
67
68 i = 1
69 ia = 2
70 r = 3
71 ra = 4
72 d = 5.5
73 da = 6.5
74 c = cmplx (7.5, 1.5)
75 ca = cmplx (8.5, -3.0)
76 v = .false.
77 cnt = -1
78
79 !$omp parallel num_threads (3) private (n) reduction (.or.:v) &
80 !$omp & reduction (-:i, ia, r, ra, d, da, c, ca)
81 !$ if (i .ne. 0 .or. any (ia .ne. 0)) v = .true.
82 !$ if (r .ne. 0 .or. any (ra .ne. 0)) v = .true.
83 !$ if (d .ne. 0 .or. any (da .ne. 0)) v = .true.
84 !$ if (c .ne. cmplx (0) .or. any (ca .ne. cmplx (0))) v = .true.
85 n = omp_get_thread_num ()
86 if (n .eq. 0) then
87 cnt = omp_get_num_threads ()
88 i = 4
89 ia(3:5) = -2
90 r = 5
91 ra(1:2) = 6.5
92 d = -2.5
93 da(2:4) = 8.5
94 c = cmplx (2.5, -3.5)
95 ca(1) = cmplx (4.5, 5)
96 else if (n .eq. 1) then
97 i = 2
98 ia(4:6) = 5
99 r = 1
100 ra(2:4) = -1.5
101 d = 8.5
102 da(1:3) = 2.5
103 c = cmplx (0.5, -3)
104 ca(2:3) = cmplx (-1, 6)
105 else
106 i = 1
107 ia = 1
108 r = -1
109 ra = -1
110 d = 1
111 da = -1
112 c = 1
113 ca = cmplx (-1, 0)
114 end if
115 !$omp end parallel
116 if (v) call abort
117 if (cnt .eq. 3) then
118 if (i .ne. 8 .or. any (ia .ne. (/3, 3, 1, 6, 6, 8/))) call abort
119 if (r .ne. 8 .or. any (ra .ne. (/9.5, 8.0, 1.5, 1.5/))) call abort
120 if (d .ne. 12.5 .or. any (da .ne. (/8.0, 16.5, 16.5, 14.0, 5.5/))) call abort
121 if (c .ne. cmplx (11.5, -5)) call abort
122 if (ca(1) .ne. cmplx (12, 2)) call abort
123 if (ca(2) .ne. cmplx (6.5, 3) .or. ca(2) .ne. ca(3)) call abort
124 end if
125
126 i = 1
127 ia = 2
128 r = 4
129 ra = 8
130 d = 16
131 da = 32
132 c = 2
133 ca = cmplx (0, 2)
134 v = .false.
135 cnt = -1
136
137 !$omp parallel num_threads (3) private (n) reduction (.or.:v) &
138 !$omp & reduction (*:i, ia, r, ra, d, da, c, ca)
139 !$ if (i .ne. 1 .or. any (ia .ne. 1)) v = .true.
140 !$ if (r .ne. 1 .or. any (ra .ne. 1)) v = .true.
141 !$ if (d .ne. 1 .or. any (da .ne. 1)) v = .true.
142 !$ if (c .ne. cmplx (1) .or. any (ca .ne. cmplx (1))) v = .true.
143 n = omp_get_thread_num ()
144 if (n .eq. 0) then
145 cnt = omp_get_num_threads ()
146 i = 3
147 ia(3:5) = 2
148 r = 0.5
149 ra(1:2) = 2
150 d = -1
151 da(2:4) = -2
152 c = 2.5
153 ca(1) = cmplx (-5, 0)
154 else if (n .eq. 1) then
155 i = 2
156 ia(4:6) = -2
157 r = 8
158 ra(2:4) = -0.5
159 da(1:3) = -1
160 c = -3
161 ca(2:3) = cmplx (0, -1)
162 else
163 ia = 2
164 r = 0.5
165 ra = 0.25
166 d = 2.5
167 da = -1
168 c = cmplx (0, -1)
169 ca = cmplx (-1, 0)
170 end if
171 !$omp end parallel
172 if (v) call abort
173 if (cnt .eq. 3) then
174 if (i .ne. 6 .or. any (ia .ne. (/4, 4, 8, -16, -16, -8/))) call abort
175 if (r .ne. 8 .or. any (ra .ne. (/4., -2., -1., -1./))) call abort
176 if (d .ne. -40 .or. any (da .ne. (/32., -64., -64., 64., -32./))) call abort
177 if (c .ne. cmplx (0, 15)) call abort
178 if (ca(1) .ne. cmplx (0, 10)) call abort
179 if (ca(2) .ne. cmplx (-2, 0) .or. ca(2) .ne. ca(3)) call abort
180 end if
181 end