1 char z[10] = { 0 };
2
3 template <int N>
4 __attribute__((noinline, noclone)) void
foo(int (* & x)[3][N],int * y,long (& w)[1][N])5 foo (int (*&x)[3][N], int *y, long (&w)[1][N])
6 {
7 unsigned long long a[9] = {};
8 short b[5] = {};
9 #pragma omp parallel for reduction(+:x[0:N][:][0:N], z[:4]) \
10 reduction(*:y[:3]) reduction(|:a[:4]) \
11 reduction(&:w[0:][:N]) reduction(max:b)
12 for (int i = 0; i < 128; i++)
13 {
14 x[i / 64][i % 3][(i / 4) & 1] += i;
15 if ((i & 15) == 1)
16 y[0] *= 3;
17 if ((i & 31) == N)
18 y[1] *= 7;
19 if ((i & 63) == 3)
20 y[N] *= 17;
21 z[i / 32] += (i & 3);
22 if (i < 4)
23 z[i] += i;
24 a[i / 32] |= 1ULL << (i & 30);
25 w[0][i & 1] &= ~(1L << (i / 17 * 3));
26 if ((i % 79) > b[0])
27 b[0] = i % 79;
28 if ((i % 13) > b[1])
29 b[1] = i % 13;
30 if ((i % 23) > b[N])
31 b[N] = i % 23;
32 if ((i % 85) > b[3])
33 b[3] = i % 85;
34 if ((i % 192) > b[4])
35 b[4] = i % 192;
36 }
37 for (int i = 0; i < 9; i++)
38 if (a[i] != (i < 4 ? 0x55555555ULL : 0))
39 __builtin_abort ();
40 if (b[0] != 78 || b[1] != 12 || b[N] != 22 || b[3] != 84 || b[4] != 127)
41 __builtin_abort ();
42 }
43
44 int a3[4][3][2];
45 int (*p3)[3][2] = &a3[1];
46 int y3[5] = { 0, 1, 1, 1, 0 };
47 long w3[1][2] = { ~0L, ~0L };
48 short bb[5];
49
50 template <int N>
51 struct S
52 {
53 int (*&x)[3][N];
54 int *y;
55 long (&w)[1][N];
56 char z[10];
57 short (&b)[5];
58 unsigned long long a[9];
SS59 S() : x(p3), y(y3+1), w(w3), z(), a(), b(bb) {}
60 __attribute__((noinline, noclone)) void foo ();
61 };
62
63 template <int N>
64 void
foo()65 S<N>::foo ()
66 {
67 #pragma omp parallel for reduction(+:x[0:N][:][0:N], z[:4]) \
68 reduction(*:y[:3]) reduction(|:a[:4]) \
69 reduction(&:w[0:][:N]) reduction(max:b)
70 for (int i = 0; i < 128; i++)
71 {
72 x[i / 64][i % 3][(i / 4) & 1] += i;
73 if ((i & 15) == 1)
74 y[0] *= 3;
75 if ((i & 31) == N)
76 y[1] *= 7;
77 if ((i & 63) == 3)
78 y[N] *= 17;
79 z[i / 32] += (i & 3);
80 if (i < 4)
81 z[i] += i;
82 a[i / 32] |= 1ULL << (i & 30);
83 w[0][i & 1] &= ~(1L << (i / 17 * 3));
84 if ((i % 79) > b[0])
85 b[0] = i % 79;
86 if ((i % 13) > b[1])
87 b[1] = i % 13;
88 if ((i % 23) > b[N])
89 b[N] = i % 23;
90 if ((i % 85) > b[3])
91 b[3] = i % 85;
92 if ((i % 192) > b[4])
93 b[4] = i % 192;
94 }
95 }
96
97 int
main()98 main ()
99 {
100 int a[4][3][2] = {};
101 static int a2[4][3][2] = {{{ 0, 0 }, { 0, 0 }, { 0, 0 }},
102 {{ 312, 381 }, { 295, 356 }, { 337, 335 }},
103 {{ 1041, 975 }, { 1016, 1085 }, { 935, 1060 }},
104 {{ 0, 0 }, { 0, 0 }, { 0, 0 }}};
105 int (*p)[3][2] = &a[1];
106 int y[5] = { 0, 1, 1, 1, 0 };
107 int y2[5] = { 0, 6561, 2401, 289, 0 };
108 char z2[10] = { 48, 49, 50, 51, 0, 0, 0, 0, 0, 0 };
109 long w[1][2] = { ~0L, ~0L };
110 foo<2> (p, y + 1, w);
111 if (__builtin_memcmp (a, a2, sizeof (a))
112 || __builtin_memcmp (y, y2, sizeof (y))
113 || __builtin_memcmp (z, z2, sizeof (z))
114 || w[0][0] != ~0x249249L
115 || w[0][1] != ~0x249249L)
116 __builtin_abort ();
117 S<2> s;
118 s.foo ();
119 for (int i = 0; i < 9; i++)
120 if (s.a[i] != (i < 4 ? 0x55555555ULL : 0))
121 __builtin_abort ();
122 if (__builtin_memcmp (a3, a2, sizeof (a3))
123 || __builtin_memcmp (y3, y2, sizeof (y3))
124 || __builtin_memcmp (s.z, z2, sizeof (s.z))
125 || w3[0][0] != ~0x249249L
126 || w3[0][1] != ~0x249249L)
127 __builtin_abort ();
128 if (bb[0] != 78 || bb[1] != 12 || bb[2] != 22 || bb[3] != 84 || bb[4] != 127)
129 __builtin_abort ();
130 }
131