blob: 2b8005dac150e59ae12939e0dba4b13d22ca0e8a [file] [log] [blame]
! { dg-do run }
! double precision reductions
program reduction_3
implicit none
integer, parameter :: n = 10, vl = 2
integer :: i
double precision, parameter :: e = .001
double precision :: vresult, result
logical :: lresult, lvresult
double precision, dimension (n) :: array
do i = 1, n
array(i) = i
end do
result = 0
vresult = 0
! '+' reductions
!$acc parallel vector_length(vl) num_gangs(1)
!$acc loop reduction(+:result)
do i = 1, n
result = result + array(i)
end do
!$acc end parallel
! Verify the results
do i = 1, n
vresult = vresult + array(i)
end do
if (abs (result - vresult) .ge. e) call abort
result = 1
vresult = 1
! '*' reductions
!$acc parallel vector_length(vl) num_gangs(1)
!$acc loop reduction(*:result)
do i = 1, n
result = result * array(i)
end do
!$acc end parallel
! Verify the results
do i = 1, n
vresult = vresult * array(i)
end do
if (result.ne.vresult) call abort
result = 0
vresult = 0
! 'max' reductions
!$acc parallel vector_length(vl) num_gangs(1)
!$acc loop reduction(max:result)
do i = 1, n
result = max (result, array(i))
end do
!$acc end parallel
! Verify the results
do i = 1, n
vresult = max (vresult, array(i))
end do
if (result.ne.vresult) call abort
result = 1
vresult = 1
! 'min' reductions
!$acc parallel vector_length(vl) num_gangs(1)
!$acc loop reduction(min:result)
do i = 1, n
result = min (result, array(i))
end do
!$acc end parallel
! Verify the results
do i = 1, n
vresult = min (vresult, array(i))
end do
if (result.ne.vresult) call abort
result = 1
vresult = 1
! '.and.' reductions
!$acc parallel vector_length(vl) num_gangs(1)
!$acc loop reduction(.and.:lresult)
do i = 1, n
lresult = lresult .and. (array(i) .ge. 5)
end do
!$acc end parallel
! Verify the results
do i = 1, n
lvresult = lvresult .and. (array(i) .ge. 5)
end do
if (result.ne.vresult) call abort
lresult = .false.
lvresult = .false.
! '.or.' reductions
!$acc parallel vector_length(vl) num_gangs(1)
!$acc loop reduction(.or.:lresult)
do i = 1, n
lresult = lresult .or. (array(i) .ge. 5)
end do
!$acc end parallel
! Verify the results
do i = 1, n
lvresult = lvresult .or. (array(i) .ge. 5)
end do
if (result.ne.vresult) call abort
lresult = .false.
lvresult = .false.
! '.eqv.' reductions
!$acc parallel vector_length(vl) num_gangs(1)
!$acc loop reduction(.eqv.:lresult)
do i = 1, n
lresult = lresult .eqv. (array(i) .ge. 5)
end do
!$acc end parallel
! Verify the results
do i = 1, n
lvresult = lvresult .eqv. (array(i) .ge. 5)
end do
if (result.ne.vresult) call abort
lresult = .false.
lvresult = .false.
! '.neqv.' reductions
!$acc parallel vector_length(vl) num_gangs(1)
!$acc loop reduction(.neqv.:lresult)
do i = 1, n
lresult = lresult .neqv. (array(i) .ge. 5)
end do
!$acc end parallel
! Verify the results
do i = 1, n
lvresult = lvresult .neqv. (array(i) .ge. 5)
end do
if (result.ne.vresult) call abort
end program reduction_3