1! PR middle-end/66199 2! { dg-do run } 3! { dg-options "-O2" } 4 5 integer :: u(1024), v(1024), w(1024), a, b, c, d, e, a1, b1, a2, b2, d1, d2 6 a = 1 7 b = 1024 8 d = 75 9 !$omp parallel do simd default(none) firstprivate (a, b) shared(u, v, w) 10 do d = a, b 11 u(d) = v(d) + w(d) 12 end do 13 if (d .ne. 1025) STOP 1 14 c = 17 15 d = 75 16 !$omp parallel do simd default(none) firstprivate (a, b) shared(u, v, w) & 17 !$omp& linear(d) linear(c:5) lastprivate(e) 18 do d = a, b 19 u(d) = v(d) + w(d) 20 c = c + 5 21 e = c 22 end do 23 if (d .ne. 1025 .or. c .ne. (17 + 5 * 1024)) STOP 2 24 if (e .ne. (17 + 5 * 1024)) STOP 3 25 a1 = 0 26 a2 = 0 27 b1 = 31 28 b2 = 31 29 d1 = 7 30 d2 = 9 31 !$omp parallel do simd default(none) firstprivate (a1, b1, a2, b2) & 32 !$omp& shared(u, v, w) lastprivate(d1, d2) collapse(2) 33 do d1 = a1, b1 34 do d2 = a2, b2 35 u(d1 * 32 + d2 + 1) = v(d1 * 32 + d2 + 1) + w(d1 * 32 + d2 + 1) 36 end do 37 end do 38 if (d1 .ne. 32 .or. d2 .ne. 32) STOP 4 39 d1 = 7 40 d2 = 9 41 !$omp parallel do simd default(none) firstprivate (a1, b1, a2, b2) & 42 !$omp& shared(u, v, w) collapse(2) 43 do d1 = a1, b1 44 do d2 = a2, b2 45 u(d1 * 32 + d2 + 1) = v(d1 * 32 + d2 + 1) + w(d1 * 32 + d2 + 1) 46 end do 47 end do 48 if (d1 .ne. 32 .or. d2 .ne. 32) STOP 5 49end 50