1*404b540aSrobert #include <omp.h>
2*404b540aSrobert #include <stdlib.h>
3*404b540aSrobert 
4*404b540aSrobert int
main(void)5*404b540aSrobert main (void)
6*404b540aSrobert {
7*404b540aSrobert   int i = 0, j = 0, k = ~0;
8*404b540aSrobert   double d = 1.0;
9*404b540aSrobert #pragma omp parallel num_threads(4) reduction(+:i) reduction(*:d) reduction(&:k)
10*404b540aSrobert   {
11*404b540aSrobert     if (i != 0 || d != 1.0 || k != ~0)
12*404b540aSrobert #pragma omp atomic
13*404b540aSrobert       j |= 1;
14*404b540aSrobert 
15*404b540aSrobert     if (omp_get_num_threads () != 4)
16*404b540aSrobert #pragma omp atomic
17*404b540aSrobert       j |= 2;
18*404b540aSrobert 
19*404b540aSrobert     i = omp_get_thread_num ();
20*404b540aSrobert     d = i + 1;
21*404b540aSrobert     k = ~(1 << (2 * i));
22*404b540aSrobert   }
23*404b540aSrobert 
24*404b540aSrobert   if (j & 1)
25*404b540aSrobert     abort ();
26*404b540aSrobert   if ((j & 2) == 0)
27*404b540aSrobert     {
28*404b540aSrobert       if (i != (0 + 1 + 2 + 3))
29*404b540aSrobert 	abort ();
30*404b540aSrobert       if (d != (1.0 * 2.0 * 3.0 * 4.0))
31*404b540aSrobert 	abort ();
32*404b540aSrobert       if (k != (~0 ^ 0x55))
33*404b540aSrobert 	abort ();
34*404b540aSrobert     }
35*404b540aSrobert   return 0;
36*404b540aSrobert }
37