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
|
! { dg-do run }
! { dg-options "-ff2c" }
! PR 34868
function f(a) result(res)
implicit none
real(8), intent(in) :: a(:)
complex(8) :: res
res = cmplx(sum(a),product(a),8)
end function f
function g(a)
implicit none
real(8), intent(in) :: a(:)
complex(8) :: g
g = cmplx(sum(a),product(a),8)
end function g
program test
real(8) :: a(1,5)
complex(8) :: c
integer :: i
interface
complex(8) function f(a)
real(8), intent(in) :: a(:)
end function f
function g(a) result(res)
real(8), intent(in) :: a(:)
complex(8) :: res
end function g
end interface
do i = 1, 5
a(1,i) = sqrt(real(i,kind(a)))
end do
c = f(a(1,:))
call check (real(c), sum(a))
call check (imag(c), product(a))
c = g(a(1,:))
call check (real(c), sum(a))
call check (imag(c), product(a))
contains
subroutine check (a, b)
real(8), intent(in) :: a, b
if (abs(a - b) > 1.e-10_8) call abort
end subroutine check
end program test
|
