111
|
1 ! { dg-do run }
|
|
2
|
|
3 module m
|
|
4 interface operator(.add.)
|
|
5 module procedure do_add
|
|
6 end interface
|
|
7 type dt
|
|
8 real :: r = 0.0
|
|
9 end type
|
|
10 contains
|
|
11 function do_add(x, y)
|
|
12 type (dt), intent (in) :: x, y
|
|
13 type (dt) :: do_add
|
|
14 do_add%r = x%r + y%r
|
|
15 end function
|
|
16 subroutine dp_add(x, y)
|
|
17 double precision :: x, y
|
|
18 x = x + y
|
|
19 end subroutine
|
|
20 subroutine dp_init(x)
|
|
21 double precision :: x
|
|
22 x = 0.0
|
|
23 end subroutine
|
|
24 end module
|
|
25
|
|
26 program udr5
|
|
27 use m, only : operator(.add.), dt, dp_add, dp_init
|
|
28 type(dt) :: xdt, one
|
|
29 real :: r
|
|
30 integer (kind = 4) :: i4
|
|
31 integer (kind = 8) :: i8
|
|
32 real (kind = 4) :: r4
|
|
33 double precision :: dp
|
|
34 !$omp declare reduction(.add.:dt:omp_out=omp_out.add.omp_in)
|
|
35 !$omp declare reduction(foo:integer(4),integer(kind=8),real (kind = 4) &
|
|
36 !$omp & :omp_out = omp_out + omp_in) initializer (omp_priv = 0)
|
|
37 !$omp declare reduction(foo:double precision:dp_add (omp_out, omp_in)) &
|
|
38 !$omp & initializer (dp_init (omp_priv))
|
|
39
|
|
40 one%r = 1.0
|
|
41 r = 0.0
|
|
42 i4 = 0
|
|
43 i8 = 0
|
|
44 r4 = 0.0
|
|
45 call dp_init (dp)
|
|
46 !$omp parallel reduction(.add.: xdt) reduction(+: r) &
|
|
47 !$omp & reduction(foo: i4, i8, r4, dp)
|
|
48 xdt = xdt.add.one
|
|
49 r = r + 1.0
|
|
50 i4 = i4 + 1
|
|
51 i8 = i8 + 1
|
|
52 r4 = r4 + 1.0
|
|
53 call dp_add (dp, 1.0d0)
|
|
54 !$omp end parallel
|
131
|
55 if (xdt%r .ne. r) STOP 1
|
|
56 if (i4.ne.r.or.i8.ne.r.or.r4.ne.r.or.dp.ne.r) STOP 2
|
111
|
57 end program udr5
|