vla5.f90 7.9 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205
  1. ! { dg-do run }
  2. call test
  3. contains
  4. subroutine check (x, y, l)
  5. integer :: x, y
  6. logical :: l
  7. l = l .or. x .ne. y
  8. end subroutine check
  9. subroutine foo (c, d, e, f, g, h, i, j, k, n)
  10. use omp_lib
  11. interface
  12. subroutine GOMP_barrier () bind(c, name="GOMP_barrier")
  13. end subroutine
  14. end interface
  15. integer :: n
  16. character (len = *) :: c
  17. character (len = n) :: d
  18. integer, dimension (2, 3:5, n) :: e
  19. integer, dimension (2, 3:n, n) :: f
  20. character (len = *), dimension (5, 3:n) :: g
  21. character (len = n), dimension (5, 3:n) :: h
  22. real, dimension (:, :, :) :: i
  23. double precision, dimension (3:, 5:, 7:) :: j
  24. integer, dimension (:, :, :) :: k
  25. logical :: l
  26. integer :: p, q, r
  27. character (len = n) :: s
  28. integer, dimension (2, 3:5, n) :: t
  29. integer, dimension (2, 3:n, n) :: u
  30. character (len = n), dimension (5, 3:n) :: v
  31. character (len = 2 * n + 24) :: w
  32. integer :: x, z, z2
  33. character (len = 1) :: y
  34. s = 'PQRSTUV'
  35. forall (p = 1:2, q = 3:5, r = 1:7) t(p, q, r) = -10 + p - q + 2 * r
  36. forall (p = 1:2, q = 3:7, r = 1:7) u(p, q, r) = 30 - p + q - 2 * r
  37. forall (p = 1:5, q = 3:7, p + q .le. 8) v(p, q) = '_+|/Oo_'
  38. forall (p = 1:5, q = 3:7, p + q .gt. 8) v(p, q) = '///|||!'
  39. l = .false.
  40. call omp_set_dynamic (.false.)
  41. call omp_set_num_threads (6)
  42. !$omp parallel do default (none) lastprivate (c, d, e, f, g, h, i, j, k) &
  43. !$omp & lastprivate (s, t, u, v) reduction (.or.:l) num_threads (6) &
  44. !$omp private (p, q, r, w, x, y) schedule (static) shared (z2)
  45. do 110 z = 0, omp_get_num_threads () - 1
  46. if (omp_get_thread_num () .eq. 0) z2 = omp_get_num_threads ()
  47. x = omp_get_thread_num ()
  48. w = ''
  49. if (x .eq. 0) w = 'thread0thr_number_0THREAD0THR_NUMBER_0'
  50. if (x .eq. 1) w = 'thread1thr_number_1THREAD1THR_NUMBER_1'
  51. if (x .eq. 2) w = 'thread2thr_number_2THREAD2THR_NUMBER_2'
  52. if (x .eq. 3) w = 'thread3thr_number_3THREAD3THR_NUMBER_3'
  53. if (x .eq. 4) w = 'thread4thr_number_4THREAD4THR_NUMBER_4'
  54. if (x .eq. 5) w = 'thread5thr_number_5THREAD5THR_NUMBER_5'
  55. c = w(8:19)
  56. d = w(1:7)
  57. forall (p = 1:2, q = 3:5, r = 1:7) e(p, q, r) = 5 * x + p + q + 2 * r
  58. forall (p = 1:2, q = 3:7, r = 1:7) f(p, q, r) = 25 * x + p + q + 2 * r
  59. forall (p = 1:5, q = 3:7, p + q .le. 8) g(p, q) = w(8:19)
  60. forall (p = 1:5, q = 3:7, p + q .gt. 8) g(p, q) = w(27:38)
  61. forall (p = 1:5, q = 3:7, p + q .le. 8) h(p, q) = w(1:7)
  62. forall (p = 1:5, q = 3:7, p + q .gt. 8) h(p, q) = w(20:26)
  63. forall (p = 3:5, q = 2:6, r = 1:7) i(p - 2, q - 1, r) = (7.5 + x) * p * q * r
  64. forall (p = 3:5, q = 2:6, r = 1:7) j(p, q + 3, r + 6) = (9.5 + x) * p * q * r
  65. forall (p = 1:5, q = 7:7, r = 4:6) k(p, q - 6, r - 3) = 19 + x + p + q + 3 * r
  66. s = w(20:26)
  67. forall (p = 1:2, q = 3:5, r = 1:7) t(p, q, r) = -10 + x + p - q + 2 * r
  68. forall (p = 1:2, q = 3:7, r = 1:7) u(p, q, r) = 30 - x - p + q - 2 * r
  69. forall (p = 1:5, q = 3:7, p + q .le. 8) v(p, q) = w(1:7)
  70. forall (p = 1:5, q = 3:7, p + q .gt. 8) v(p, q) = w(20:26)
  71. call GOMP_barrier
  72. y = ''
  73. if (x .eq. 0) y = '0'
  74. if (x .eq. 1) y = '1'
  75. if (x .eq. 2) y = '2'
  76. if (x .eq. 3) y = '3'
  77. if (x .eq. 4) y = '4'
  78. if (x .eq. 5) y = '5'
  79. l = l .or. w(7:7) .ne. y
  80. l = l .or. w(19:19) .ne. y
  81. l = l .or. w(26:26) .ne. y
  82. l = l .or. w(38:38) .ne. y
  83. l = l .or. c .ne. w(8:19)
  84. l = l .or. d .ne. w(1:7)
  85. l = l .or. s .ne. w(20:26)
  86. do 103, p = 1, 2
  87. do 103, q = 3, 7
  88. do 103, r = 1, 7
  89. if (q .lt. 6) l = l .or. e(p, q, r) .ne. 5 * x + p + q + 2 * r
  90. l = l .or. f(p, q, r) .ne. 25 * x + p + q + 2 * r
  91. if (r .lt. 6 .and. q + r .le. 8) l = l .or. g(r, q) .ne. w(8:19)
  92. if (r .lt. 6 .and. q + r .gt. 8) l = l .or. g(r, q) .ne. w(27:38)
  93. if (r .lt. 6 .and. q + r .le. 8) l = l .or. h(r, q) .ne. w(1:7)
  94. if (r .lt. 6 .and. q + r .gt. 8) l = l .or. h(r, q) .ne. w(20:26)
  95. if (q .lt. 6) l = l .or. t(p, q, r) .ne. -10 + x + p - q + 2 * r
  96. l = l .or. u(p, q, r) .ne. 30 - x - p + q - 2 * r
  97. if (r .lt. 6 .and. q + r .le. 8) l = l .or. v(r, q) .ne. w(1:7)
  98. if (r .lt. 6 .and. q + r .gt. 8) l = l .or. v(r, q) .ne. w(20:26)
  99. 103 continue
  100. do 104, p = 3, 5
  101. do 104, q = 2, 6
  102. do 104, r = 1, 7
  103. l = l .or. i(p - 2, q - 1, r) .ne. (7.5 + x) * p * q * r
  104. l = l .or. j(p, q + 3, r + 6) .ne. (9.5 + x) * p * q * r
  105. 104 continue
  106. do 105, p = 1, 5
  107. do 105, q = 4, 6
  108. l = l .or. k(p, 1, q - 3) .ne. 19 + x + p + 7 + 3 * q
  109. 105 continue
  110. call check (size (e, 1), 2, l)
  111. call check (size (e, 2), 3, l)
  112. call check (size (e, 3), 7, l)
  113. call check (size (e), 42, l)
  114. call check (size (f, 1), 2, l)
  115. call check (size (f, 2), 5, l)
  116. call check (size (f, 3), 7, l)
  117. call check (size (f), 70, l)
  118. call check (size (g, 1), 5, l)
  119. call check (size (g, 2), 5, l)
  120. call check (size (g), 25, l)
  121. call check (size (h, 1), 5, l)
  122. call check (size (h, 2), 5, l)
  123. call check (size (h), 25, l)
  124. call check (size (i, 1), 3, l)
  125. call check (size (i, 2), 5, l)
  126. call check (size (i, 3), 7, l)
  127. call check (size (i), 105, l)
  128. call check (size (j, 1), 4, l)
  129. call check (size (j, 2), 5, l)
  130. call check (size (j, 3), 7, l)
  131. call check (size (j), 140, l)
  132. call check (size (k, 1), 5, l)
  133. call check (size (k, 2), 1, l)
  134. call check (size (k, 3), 3, l)
  135. call check (size (k), 15, l)
  136. 110 continue
  137. !$omp end parallel do
  138. if (l) call abort
  139. if (z2 == 6) then
  140. x = 5
  141. w = 'thread5thr_number_5THREAD5THR_NUMBER_5'
  142. y = '5'
  143. l = l .or. w(7:7) .ne. y
  144. l = l .or. w(19:19) .ne. y
  145. l = l .or. w(26:26) .ne. y
  146. l = l .or. w(38:38) .ne. y
  147. l = l .or. c .ne. w(8:19)
  148. l = l .or. d .ne. w(1:7)
  149. l = l .or. s .ne. w(20:26)
  150. do 113, p = 1, 2
  151. do 113, q = 3, 7
  152. do 113, r = 1, 7
  153. if (q .lt. 6) l = l .or. e(p, q, r) .ne. 5 * x + p + q + 2 * r
  154. l = l .or. f(p, q, r) .ne. 25 * x + p + q + 2 * r
  155. if (r .lt. 6 .and. q + r .le. 8) l = l .or. g(r, q) .ne. w(8:19)
  156. if (r .lt. 6 .and. q + r .gt. 8) l = l .or. g(r, q) .ne. w(27:38)
  157. if (r .lt. 6 .and. q + r .le. 8) l = l .or. h(r, q) .ne. w(1:7)
  158. if (r .lt. 6 .and. q + r .gt. 8) l = l .or. h(r, q) .ne. w(20:26)
  159. if (q .lt. 6) l = l .or. t(p, q, r) .ne. -10 + x + p - q + 2 * r
  160. l = l .or. u(p, q, r) .ne. 30 - x - p + q - 2 * r
  161. if (r .lt. 6 .and. q + r .le. 8) l = l .or. v(r, q) .ne. w(1:7)
  162. if (r .lt. 6 .and. q + r .gt. 8) l = l .or. v(r, q) .ne. w(20:26)
  163. 113 continue
  164. do 114, p = 3, 5
  165. do 114, q = 2, 6
  166. do 114, r = 1, 7
  167. l = l .or. i(p - 2, q - 1, r) .ne. (7.5 + x) * p * q * r
  168. l = l .or. j(p, q + 3, r + 6) .ne. (9.5 + x) * p * q * r
  169. 114 continue
  170. do 115, p = 1, 5
  171. do 115, q = 4, 6
  172. l = l .or. k(p, 1, q - 3) .ne. 19 + x + p + 7 + 3 * q
  173. 115 continue
  174. if (l) call abort
  175. end if
  176. end subroutine foo
  177. subroutine test
  178. character (len = 12) :: c
  179. character (len = 7) :: d
  180. integer, dimension (2, 3:5, 7) :: e
  181. integer, dimension (2, 3:7, 7) :: f
  182. character (len = 12), dimension (5, 3:7) :: g
  183. character (len = 7), dimension (5, 3:7) :: h
  184. real, dimension (3:5, 2:6, 1:7) :: i
  185. double precision, dimension (3:6, 2:6, 1:7) :: j
  186. integer, dimension (1:5, 7:7, 4:6) :: k
  187. integer :: p, q, r
  188. c = 'abcdefghijkl'
  189. d = 'ABCDEFG'
  190. forall (p = 1:2, q = 3:5, r = 1:7) e(p, q, r) = 5 + p + q + 2 * r
  191. forall (p = 1:2, q = 3:7, r = 1:7) f(p, q, r) = 25 + p + q + 2 * r
  192. forall (p = 1:5, q = 3:7, p + q .le. 8) g(p, q) = '0123456789AB'
  193. forall (p = 1:5, q = 3:7, p + q .gt. 8) g(p, q) = '9876543210ZY'
  194. forall (p = 1:5, q = 3:7, p + q .le. 8) h(p, q) = '0123456'
  195. forall (p = 1:5, q = 3:7, p + q .gt. 8) h(p, q) = '9876543'
  196. forall (p = 3:5, q = 2:6, r = 1:7) i(p, q, r) = 7.5 * p * q * r
  197. forall (p = 3:6, q = 2:6, r = 1:7) j(p, q, r) = 9.5 * p * q * r
  198. forall (p = 1:5, q = 7:7, r = 4:6) k(p, q, r) = 19 + p + q + 3 * r
  199. call foo (c, d, e, f, g, h, i, j, k, 7)
  200. end subroutine test
  201. end