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
|
! { dg-do run }
!
! PR119540 comment2: REDUCE was getting the shape wrong. This testcase also
! verifies that the longest possible name for the OPERATION wrapper function
! is catered for.
!
! Contributed by Harald Anlauf <anlauf@gcc.gnu.org>
!
program p2345678901234567890123456789012345678901234567890123456789_123
implicit none
integer, parameter :: n = 3
integer, parameter :: vec(n) = [2, 5, 10]
integer, parameter :: mat(n,2) = reshape([vec,2*vec],[n,2])
integer :: mat_shape(2), reduce_shape(1), r
integer, dimension(:), allocatable :: res1
mat_shape = shape (mat)
reduce_shape = shape (reduce (mat, add, 1), 1)
if (reduce_shape(1) /= mat_shape(2)) stop 1
reduce_shape = shape (reduce (mat, add, 1), 1)
if (reduce_shape(1) /= mat_shape(2)) stop 2
res1 = reduce (mat, add, 1)
if (any (res1 /= [17, 34])) stop 3
res1 = reduce (mat, add, 2)
if (any (res1 /= [6, 15, 30])) stop 4
r = reduce (vec, &
o2345678901234567890123456789012345678901234567890123456789_123)
if (r /= 17) stop 5
deallocate (res1)
contains
pure function add(i,j) result(sum_ij)
integer, intent(in) :: i, j
integer :: sum_ij
sum_ij = i + j
end function add
pure function o2345678901234567890123456789012345678901234567890123456789_123 (i, j) &
result (sum)
integer, intent(in) :: i, j
integer :: sum
sum = i + j
end function
end
|
