8 main (int argc
, char **argv
[])
10 int n
= argc
< 5 ? 12 : 31, i
, m
, l
;
12 unsigned short b
[n
/ 2 - 1];
15 for (i
= 0; i
< n
+ 3; i
++)
17 for (i
= 0; i
< n
/ 2 - 1; i
++)
19 for (i
= 0; i
< n
* 2 + 1; i
++)
23 #pragma omp parallel default (shared) num_threads (4) \
24 firstprivate (a, m) private (b, i) reduction (+:l)
26 for (i
= 0; i
< m
+ 3; i
++)
29 for (i
= 0; i
< m
* 2 + 1; i
++)
30 if (c
[i
] != ((i
<< 24) | i
))
33 memset (a
, omp_get_thread_num (), m
+ 3);
34 for (i
= 0; i
< m
/ 2 - 1; i
++)
38 for (i
= 0; i
< m
* 2 + 1; i
++)
42 if (a
[0] != omp_get_thread_num ())
44 for (i
= 1; i
< m
+ 3; i
++)
47 for (i
= 0; i
< m
/ 2 - 1; i
++)
50 for (i
= 0; i
< m
* 2 + 1; i
++)
56 for (i
= 0; i
< n
* 2 + 1; i
++)