File: f2c_9.f90

package info (click to toggle)
gcc-arm-none-eabi 15%3A14.2.rel1-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 1,099,328 kB
  • sloc: cpp: 3,627,108; ansic: 2,571,498; ada: 834,230; f90: 235,082; makefile: 79,231; asm: 74,984; xml: 51,692; exp: 39,736; sh: 33,298; objc: 15,629; python: 15,069; fortran: 14,429; pascal: 7,003; awk: 5,070; perl: 3,106; ml: 285; lisp: 253; lex: 204; haskell: 135
file content (52 lines) | stat: -rw-r--r-- 979 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
! { 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) STOP 1
  end subroutine check
end program test