| ! { dg-do run } |
| |
| program udr7 |
| implicit none |
| interface |
| elemental subroutine omp_priv (x, y, z) |
| real, intent (in) :: x |
| real, intent (inout) :: y |
| real, intent (in) :: z |
| end subroutine omp_priv |
| elemental real function omp_orig (x) |
| real, intent (in) :: x |
| end function omp_orig |
| end interface |
| !$omp declare reduction (omp_priv : real : & |
| !$omp & omp_priv (omp_orig (omp_in), omp_out, 1.0)) & |
| !$omp & initializer (omp_out (omp_priv, omp_in (omp_orig))) |
| real :: x (2:4, 1:1, -2:0) |
| integer :: i |
| x = 0 |
| !$omp parallel do reduction (omp_priv : x) |
| do i = 1, 64 |
| x = x + i |
| end do |
| if (any (x /= 2080.0)) stop 1 |
| contains |
| elemental subroutine omp_out (x, y) |
| real, intent (out) :: x |
| real, intent (in) :: y |
| x = y - 4.0 |
| end subroutine omp_out |
| elemental real function omp_in (x) |
| real, intent (in) :: x |
| omp_in = x + 4.0 |
| end function omp_in |
| end program udr7 |
| elemental subroutine omp_priv (x, y, z) |
| real, intent (in) :: x |
| real, intent (inout) :: y |
| real, intent (in) :: z |
| y = y + (x - 4.0) + (z - 1.0) |
| end subroutine omp_priv |
| elemental real function omp_orig (x) |
| real, intent (in) :: x |
| omp_orig = x + 4.0 |
| end function omp_orig |