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