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