1 #include <omp.h>
2 #include <stdlib.h>
3 
4 int
main(void)5 main (void)
6 {
7   int i = 0, j = 0, k = ~0, l;
8   double d = 1.0;
9 #pragma omp parallel num_threads(4)
10   {
11 #pragma omp single
12     {
13       i = 16;
14       k ^= (1 << 16);
15       d += 32.0;
16     }
17 
18 #pragma omp for reduction(+:i) reduction(*:d) reduction(&:k) nowait
19     for (l = 0; l < 4; l++)
20       {
21 	if (omp_get_num_threads () == 4 && (i != 0 || d != 1.0 || k != ~0))
22 #pragma omp atomic
23 	  j |= 1;
24 
25 	if (l == omp_get_thread_num ())
26 	  {
27 	    i = omp_get_thread_num ();
28 	    d = i + 1;
29 	    k = ~(1 << (2 * i));
30 	  }
31       }
32 
33     if (omp_get_num_threads () == 4)
34       {
35 #pragma omp barrier
36 	if (i != (16 + 0 + 1 + 2 + 3))
37 #pragma omp atomic
38 	  j |= 2;
39 	if (d != (33.0 * 1.0 * 2.0 * 3.0 * 4.0))
40 #pragma omp atomic
41 	  j |= 4;
42 	if (k != (~0 ^ 0x55 ^ (1 << 16)))
43 #pragma omp atomic
44 	  j |= 8;
45       }
46   }
47 
48   if (j)
49     abort ();
50   return 0;
51 }
52