taskloop-reduction-2.c 5.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212
  1. #ifdef __cplusplus
  2. extern "C"
  3. #endif
  4. void abort (void);
  5. int a[2];
  6. long long int b[7] = { 9, 11, 1, 1, 1, 13, 15 };
  7. int e[3] = { 5, 0, 5 };
  8. int f[5] = { 6, 7, 0, 0, 9 };
  9. int g[4] = { 1, 0, 0, 2 };
  10. int h[3] = { 0, 1, 4 };
  11. int k[4][2] = { { 5, 6 }, { 0, 0 }, { 0, 0 }, { 7, 8 } };
  12. long long *s;
  13. long long (*t)[2];
  14. void
  15. foo (int n, int *c, long long int *d, int m[3], int *r, int o[4], int *p, int q[4][2])
  16. {
  17. int i;
  18. #pragma omp taskloop in_reduction (+: a, c[:2]) in_reduction (*: b[2 * n:3 * n], d[0:2]) \
  19. in_reduction (+: o[n:n*2], m[1], k[1:2][:], p[0], f[2:2]) \
  20. in_reduction (+: q[1:2][:], g[n:n*2], e[1], h[0], r[2:2]) \
  21. in_reduction (*: s[1:2], t[2:2][:]) nogroup
  22. for (i = 0; i < 2; i++)
  23. {
  24. a[0] += 7;
  25. a[1] += 17;
  26. b[2] *= 2;
  27. b[4] *= 2;
  28. c[0] += 6;
  29. d[1] *= 2;
  30. e[1] += 19;
  31. f[2] += 21;
  32. f[3] += 23;
  33. g[1] += 25;
  34. g[2] += 27;
  35. h[0] += 29;
  36. k[1][0] += 31;
  37. k[2][1] += 33;
  38. m[1] += 19;
  39. r[2] += 21;
  40. r[3] += 23;
  41. o[1] += 25;
  42. o[2] += 27;
  43. p[0] += 29;
  44. q[1][0] += 31;
  45. q[2][1] += 33;
  46. s[1] *= 2;
  47. t[2][0] *= 2;
  48. t[3][1] *= 2;
  49. }
  50. }
  51. void
  52. test (int n)
  53. {
  54. int c[2] = { 0, 0 };
  55. int p[3] = { 0, 1, 4 };
  56. int q[4][2] = { { 5, 6 }, { 0, 0 }, { 0, 0 }, { 7, 8 } };
  57. long long ss[4] = { 5, 1, 1, 6 };
  58. long long tt[5][2] = { { 9, 10 }, { 11, 12 }, { 1, 1 }, { 1, 1 }, { 13, 14 } };
  59. s = ss;
  60. t = tt;
  61. #pragma omp parallel
  62. #pragma omp single
  63. {
  64. long long int d[] = { 1, 1 };
  65. int m[3] = { 5, 0, 5 };
  66. int r[5] = { 6, 7, 0, 0, 9 };
  67. int o[4] = { 1, 0, 0, 2 };
  68. int i;
  69. #pragma omp taskloop reduction (+: a, c) reduction (default, *: b[2 * n:3 * n], d) \
  70. reduction (+: e[1], f[2:2], g[n:n*2], h[0], k[1:2][0:2]) \
  71. reduction (default, +: o[n:n*2], m[1], q[1:2][:], p[0], r[2:2]) \
  72. reduction (*: t[2:2][:], s[1:n + 1])
  73. for (i = 0; i < 4; i++)
  74. {
  75. int j;
  76. a[0] += 2;
  77. a[1] += 3;
  78. b[2] *= 2;
  79. f[3] += 8;
  80. g[1] += 9;
  81. g[2] += 10;
  82. h[0] += 11;
  83. k[1][1] += 13;
  84. k[2][1] += 15;
  85. m[1] += 16;
  86. r[2] += 8;
  87. s[1] *= 2;
  88. t[2][1] *= 2;
  89. t[3][1] *= 2;
  90. #pragma omp taskloop in_reduction (+: a, c[:2]) \
  91. in_reduction (*: b[2 * n:3 * n], d[n - 1:n + 1]) \
  92. in_reduction (+: e[1], f[2:2], g[n:n*2], h[0], k[1:2][:2]) \
  93. in_reduction (+: m[1], r[2:2], o[n:n*2], p[0], q[1:2][:2]) \
  94. in_reduction (*: s[n:2], t[2:2][:]) nogroup
  95. for (j = 0; j < 2; j++)
  96. {
  97. m[1] += 6;
  98. r[2] += 7;
  99. q[1][0] += 17;
  100. q[2][0] += 19;
  101. a[0] += 4;
  102. a[1] += 5;
  103. b[3] *= 2;
  104. b[4] *= 2;
  105. f[3] += 18;
  106. g[1] += 29;
  107. g[2] += 18;
  108. h[0] += 19;
  109. s[2] *= 2;
  110. t[2][0] *= 2;
  111. t[3][0] *= 2;
  112. foo (n, c, d, m, r, o, p, q);
  113. r[3] += 18;
  114. o[1] += 29;
  115. o[2] += 18;
  116. p[0] += 19;
  117. c[0] += 4;
  118. c[1] += 5;
  119. d[0] *= 2;
  120. e[1] += 6;
  121. f[2] += 7;
  122. k[1][0] += 17;
  123. k[2][0] += 19;
  124. }
  125. r[3] += 8;
  126. o[1] += 9;
  127. o[2] += 10;
  128. p[0] += 11;
  129. q[1][1] += 13;
  130. q[2][1] += 15;
  131. b[3] *= 2;
  132. c[0] += 4;
  133. c[1] += 9;
  134. d[0] *= 2;
  135. e[1] += 16;
  136. f[2] += 8;
  137. }
  138. if (d[0] != 1LL << (8 + 4)
  139. || d[1] != 1LL << 16
  140. || m[0] != 5
  141. || m[1] != 19 * 16 + 6 * 8 + 16 * 4
  142. || m[2] != 5
  143. || r[0] != 6
  144. || r[1] != 7
  145. || r[2] != 21 * 16 + 7 * 8 + 8 * 4
  146. || r[3] != 23 * 16 + 18 * 8 + 8 * 4
  147. || r[4] != 9
  148. || o[0] != 1
  149. || o[1] != 25 * 16 + 29 * 8 + 9 * 4
  150. || o[2] != 27 * 16 + 18 * 8 + 10 * 4
  151. || o[3] != 2)
  152. abort ();
  153. }
  154. if (a[0] != 7 * 16 + 4 * 8 + 2 * 4
  155. || a[1] != 17 * 16 + 5 * 8 + 3 * 4
  156. || b[0] != 9 || b[1] != 11
  157. || b[2] != 1LL << (16 + 4)
  158. || b[3] != 1LL << (8 + 4)
  159. || b[4] != 1LL << (16 + 8)
  160. || b[5] != 13 || b[6] != 15
  161. || c[0] != 6 * 16 + 4 * 8 + 4 * 4
  162. || c[1] != 5 * 8 + 9 * 4
  163. || e[0] != 5
  164. || e[1] != 19 * 16 + 6 * 8 + 16 * 4
  165. || e[2] != 5
  166. || f[0] != 6
  167. || f[1] != 7
  168. || f[2] != 21 * 16 + 7 * 8 + 8 * 4
  169. || f[3] != 23 * 16 + 18 * 8 + 8 * 4
  170. || f[4] != 9
  171. || g[0] != 1
  172. || g[1] != 25 * 16 + 29 * 8 + 9 * 4
  173. || g[2] != 27 * 16 + 18 * 8 + 10 * 4
  174. || g[3] != 2
  175. || h[0] != 29 * 16 + 19 * 8 + 11 * 4
  176. || h[1] != 1 || h[2] != 4
  177. || k[0][0] != 5 || k[0][1] != 6
  178. || k[1][0] != 31 * 16 + 17 * 8
  179. || k[1][1] != 13 * 4
  180. || k[2][0] != 19 * 8
  181. || k[2][1] != 33 * 16 + 15 * 4
  182. || k[3][0] != 7 || k[3][1] != 8
  183. || p[0] != 29 * 16 + 19 * 8 + 11 * 4
  184. || p[1] != 1 || p[2] != 4
  185. || q[0][0] != 5 || q[0][1] != 6
  186. || q[1][0] != 31 * 16 + 17 * 8
  187. || q[1][1] != 13 * 4
  188. || q[2][0] != 19 * 8
  189. || q[2][1] != 33 * 16 + 15 * 4
  190. || q[3][0] != 7 || q[3][1] != 8
  191. || ss[0] != 5
  192. || ss[1] != 1LL << (16 + 4)
  193. || ss[2] != 1LL << 8
  194. || ss[3] != 6
  195. || tt[0][0] != 9 || tt[0][1] != 10 || tt[1][0] != 11 || tt[1][1] != 12
  196. || tt[2][0] != 1LL << (16 + 8)
  197. || tt[2][1] != 1LL << 4
  198. || tt[3][0] != 1LL << 8
  199. || tt[3][1] != 1LL << (16 + 4)
  200. || tt[4][0] != 13 || tt[4][1] != 14)
  201. abort ();
  202. }
  203. int
  204. main ()
  205. {
  206. test (1);
  207. return 0;
  208. }