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
|
extern void abort (void);
int a[16], b[16], c[16], d[5][2];
__attribute__((noinline, noclone)) void
foo (int x, int y)
{
int i;
#pragma omp for schedule (static, 1) reduction (+:a[:3])
for (i = 0; i < 64; i++)
{
a[0] += i;
a[1] += 2 * i;
a[2] += 3 * i;
}
#pragma omp for schedule (guided) reduction (+:b[4:3])
for (i = 0; i < 64; i++)
{
b[4] += i;
b[5] += 2 * i;
b[6] += 3 * i;
}
#pragma omp for schedule (static) reduction (+:c[x:4])
for (i = 0; i < 64; i++)
{
c[9] += i;
c[10] += 2 * i;
c[11] += 3 * i;
c[12] += 4 * i;
}
#pragma omp for reduction (+:d[x - 8:2][y:])
for (i = 0; i < 64; i++)
{
d[1][0] += i;
d[1][1] += 2 * i;
d[2][0] += 3 * i;
d[2][1] += 4 * i;
}
}
int
main ()
{
int i;
#pragma omp parallel
foo (9, 0);
for (i = 0; i < 16; i++)
if (a[i] != (i < 3 ? 64 * 63 / 2 * (i + 1) : 0)
|| b[i] != ((i >= 4 && i < 7) ? 64 * 63 / 2 * (i - 3) : 0)
|| c[i] != ((i >= 9 && i < 13) ? 64 * 63 / 2 * (i - 8) : 0))
abort ();
for (i = 0; i < 5; i++)
if (d[i][0] != ((i && i <= 2) ? 64 * 63 / 2 * (2 * i - 1) : 0)
|| d[i][1] != ((i && i <= 2) ? 64 * 63 / 2 * (2 * i) : 0))
abort ();
return 0;
}
|
