|
145
|
1 /* { dg-require-effective-target size32plus } */
|
|
|
2
|
|
|
3 extern void abort (void);
|
|
|
4 float r = 1.0f, a[1024], b[1024];
|
|
|
5
|
|
|
6 __attribute__((noipa)) void
|
|
|
7 foo (float *a, float *b)
|
|
|
8 {
|
|
|
9 #pragma omp for reduction (inscan, *:r)
|
|
|
10 for (int i = 0; i < 1024; i++)
|
|
|
11 {
|
|
|
12 r *= a[i];
|
|
|
13 #pragma omp scan inclusive(r)
|
|
|
14 b[i] = r;
|
|
|
15 }
|
|
|
16 }
|
|
|
17
|
|
|
18 __attribute__((noipa)) float
|
|
|
19 bar (void)
|
|
|
20 {
|
|
|
21 float s = -__builtin_inff ();
|
|
|
22 #pragma omp parallel for reduction (inscan, max:s)
|
|
|
23 for (int i = 0; i < 1024; i++)
|
|
|
24 {
|
|
|
25 s = s > a[i] ? s : a[i];
|
|
|
26 #pragma omp scan inclusive(s)
|
|
|
27 b[i] = s;
|
|
|
28 }
|
|
|
29 return s;
|
|
|
30 }
|
|
|
31
|
|
|
32 int
|
|
|
33 main ()
|
|
|
34 {
|
|
|
35 float s = 1.0f;
|
|
|
36 for (int i = 0; i < 1024; ++i)
|
|
|
37 {
|
|
|
38 if (i < 80)
|
|
|
39 a[i] = (i & 1) ? 0.25f : 0.5f;
|
|
|
40 else if (i < 200)
|
|
|
41 a[i] = (i % 3) == 0 ? 2.0f : (i % 3) == 1 ? 4.0f : 1.0f;
|
|
|
42 else if (i < 280)
|
|
|
43 a[i] = (i & 1) ? 0.25f : 0.5f;
|
|
|
44 else if (i < 380)
|
|
|
45 a[i] = (i % 3) == 0 ? 2.0f : (i % 3) == 1 ? 4.0f : 1.0f;
|
|
|
46 else
|
|
|
47 switch (i % 6)
|
|
|
48 {
|
|
|
49 case 0: a[i] = 0.25f; break;
|
|
|
50 case 1: a[i] = 2.0f; break;
|
|
|
51 case 2: a[i] = -1.0f; break;
|
|
|
52 case 3: a[i] = -4.0f; break;
|
|
|
53 case 4: a[i] = 0.5f; break;
|
|
|
54 case 5: a[i] = 1.0f; break;
|
|
|
55 default: a[i] = 0.0f; break;
|
|
|
56 }
|
|
|
57 b[i] = -19.0f;
|
|
|
58 asm ("" : "+g" (i));
|
|
|
59 }
|
|
|
60 #pragma omp parallel
|
|
|
61 foo (a, b);
|
|
|
62 if (r * 16384.0f != 0.125f)
|
|
|
63 abort ();
|
|
|
64 float m = -175.25f;
|
|
|
65 for (int i = 0; i < 1024; ++i)
|
|
|
66 {
|
|
|
67 s *= a[i];
|
|
|
68 if (b[i] != s)
|
|
|
69 abort ();
|
|
|
70 else
|
|
|
71 {
|
|
|
72 a[i] = m - ((i % 3) == 1 ? 2.0f : (i % 3) == 2 ? 4.0f : 0.0f);
|
|
|
73 b[i] = -231.75f;
|
|
|
74 m += 0.75f;
|
|
|
75 }
|
|
|
76 }
|
|
|
77 if (bar () != 592.0f)
|
|
|
78 abort ();
|
|
|
79 s = -__builtin_inff ();
|
|
|
80 for (int i = 0; i < 1024; ++i)
|
|
|
81 {
|
|
|
82 if (s < a[i])
|
|
|
83 s = a[i];
|
|
|
84 if (b[i] != s)
|
|
|
85 abort ();
|
|
|
86 }
|
|
|
87 return 0;
|
|
|
88 }
|
