task-reduction-7.c 5.4 KB

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