File: reduction-4.f90

package info (click to toggle)
gcc-avr 1%3A5.4.0%2BAtmel3.6.1-2
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 589,832 kB
  • sloc: ansic: 2,775,567; ada: 756,757; cpp: 723,977; f90: 117,673; asm: 66,896; makefile: 62,755; xml: 44,466; sh: 29,549; exp: 23,315; objc: 15,216; fortran: 10,901; pascal: 4,185; python: 4,093; perl: 2,969; awk: 2,811; ml: 2,385; cs: 879; yacc: 316; lex: 198; haskell: 112; lisp: 8
file content (54 lines) | stat: -rw-r--r-- 943 bytes parent folder | download | duplicates (3) 
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
51
52
53
54
 
! { dg-do run }

! complex reductions

program reduction_4
  implicit none

  integer, parameter    :: n = 10, vl = 32
  integer               :: i
  complex               :: vresult, result
  complex, 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 (result .ne. vresult) call abort

  result = 1
  vresult = 1

!  ! '*' reductions
!
!  !$acc parallel vector_length(vl)
!  !$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
end program reduction_4