1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
|
! { dg-do run }
type dt
integer :: g
integer, allocatable :: h(:)
end type
!$omp declare reduction (baz : dt : bar (omp_out, omp_in)) &
!$omp & initializer (foo (omp_priv, omp_orig))
integer :: r
type (dt), allocatable :: a(:)
allocate (a(7:8))
a(:)%g = 0
a(7)%h = (/ 0, 0, 0 /)
r = 0
!$omp parallel reduction(+:r) reduction (baz:a)
if (.not.allocated (a)) call abort
if (lbound (a, 1) /= 7 .or. ubound (a, 1) /= 8) call abort
if (.not.allocated (a(7)%h)) call abort
if (allocated (a(8)%h)) call abort
if (lbound (a(7)%h, 1) /= 1 .or. ubound (a(7)%h, 1) /= 3) call abort
a(:)%g = a(:)%g + 2
a(7)%h = a(7)%h + 3
r = r + 1
!$omp end parallel
if (.not.allocated (a)) call abort
if (lbound (a, 1) /= 7 .or. ubound (a, 1) /= 8) call abort
if (.not.allocated (a(7)%h)) call abort
if (allocated (a(8)%h)) call abort
if (lbound (a(7)%h, 1) /= 1 .or. ubound (a(7)%h, 1) /= 3) call abort
if (any (a(:)%g /= 2 * r) .or. any (a(7)%h(:) /= 3 * r)) call abort
contains
subroutine foo (x, y)
type (dt), allocatable :: x(:), y(:)
if (allocated (x) .neqv. allocated (y)) call abort
if (lbound (x, 1) /= lbound (y, 1)) call abort
if (ubound (x, 1) /= ubound (y, 1)) call abort
if (allocated (x(7)%h) .neqv. allocated (y(7)%h)) call abort
if (allocated (x(8)%h) .neqv. allocated (y(8)%h)) call abort
if (lbound (x(7)%h, 1) /= lbound (y(7)%h, 1)) call abort
if (ubound (x(7)%h, 1) /= ubound (y(7)%h, 1)) call abort
x(7)%g = 0
x(7)%h = 0
x(8)%g = 0
end subroutine
subroutine bar (x, y)
type (dt), allocatable :: x(:), y(:)
x(:)%g = x(:)%g + y(:)%g
x(7)%h(:) = x(7)%h(:) + y(7)%h(:)
end subroutine
end
|
