xref: /openbsd-src/gnu/gcc/libgomp/testsuite/libgomp.c++/reduction-2.C (revision 404b540a9034ac75a6199ad1a32d1bbc7a0d4210) 
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, l;
8*404b540aSrobert   double d = 1.0;
9*404b540aSrobert #pragma omp parallel num_threads(4)
10*404b540aSrobert   {
11*404b540aSrobert #pragma omp single
12*404b540aSrobert     {
13*404b540aSrobert       i = 16;
14*404b540aSrobert       k ^= (1 << 16);
15*404b540aSrobert       d += 32.0;
16*404b540aSrobert     }
17*404b540aSrobert 
18*404b540aSrobert #pragma omp for reduction(+:i) reduction(*:d) reduction(&:k)
19*404b540aSrobert     for (l = 0; l < 4; l++)
20*404b540aSrobert       {
21*404b540aSrobert 	if (omp_get_num_threads () == 4 && (i != 0 || d != 1.0 || k != ~0))
22*404b540aSrobert #pragma omp atomic
23*404b540aSrobert 	  j |= 1;
24*404b540aSrobert 
25*404b540aSrobert 	if (l == omp_get_thread_num ())
26*404b540aSrobert 	  {
27*404b540aSrobert 	    i = omp_get_thread_num ();
28*404b540aSrobert 	    d = i + 1;
29*404b540aSrobert 	    k = ~(1 << (2 * i));
30*404b540aSrobert 	  }
31*404b540aSrobert       }
32*404b540aSrobert 
33*404b540aSrobert     if (omp_get_num_threads () == 4)
34*404b540aSrobert       {
35*404b540aSrobert 	if (i != (16 + 0 + 1 + 2 + 3))
36*404b540aSrobert #pragma omp atomic
37*404b540aSrobert 	  j |= 2;
38*404b540aSrobert 	if (d != (33.0 * 1.0 * 2.0 * 3.0 * 4.0))
39*404b540aSrobert #pragma omp atomic
40*404b540aSrobert 	  j |= 4;
41*404b540aSrobert 	if (k != (~0 ^ 0x55 ^ (1 << 16)))
42*404b540aSrobert #pragma omp atomic
43*404b540aSrobert 	  j |= 8;
44*404b540aSrobert       }
45*404b540aSrobert   }
46*404b540aSrobert 
47*404b540aSrobert   if (j)
48*404b540aSrobert     abort ();
49*404b540aSrobert   return 0;
50*404b540aSrobert }
51