File: interface_18.f90

package info (click to toggle)
lfortran 0.59.0-3
  • links: PTS, VCS
  • area: main
  • in suites: sid
  • size: 56,736 kB
  • sloc: cpp: 168,052; f90: 74,272; python: 17,537; ansic: 7,705; yacc: 2,345; sh: 1,334; fortran: 895; makefile: 37; javascript: 15
file content (72 lines) | stat: -rw-r--r-- 1,459 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
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
module module_interface_18
    implicit none

    abstract interface
        subroutine OBJ(x, f)
            implicit none
            real(4), intent(in) :: x(:)
            real(4), intent(out) :: f
        end subroutine OBJ
    end interface

    interface evaluate
        module procedure evaluatef
    end interface evaluate

    contains

    function moderatex(x) result(y)
        implicit none
        real, intent(in) :: x(:)
        real :: y(size(x))
        y = 2
    end function moderatex

    subroutine evaluatef(calfun, x, f)
        implicit none

        procedure(OBJ) :: calfun  
        real, intent(in) :: x(:)
        real, intent(out) :: f
        call calfun(moderatex(x), f)
    end subroutine evaluatef

    subroutine initxf(calfun, xpt)
        implicit none

        procedure(OBJ) :: calfun
        real, intent(out) :: xpt(:, :)

        real :: f
        real :: x(size(xpt, 1))

        x = 1
        f = 2
        call evaluate(calfun, x, f)
    end subroutine initxf

end module module_interface_18

program interface_18
    use module_interface_18
    implicit none

    real :: xpt(1, 5)
    
    call initxf(calfun, xpt)

contains

    subroutine calfun(x, f)
        implicit none
        real, intent(in) :: x(:)
        real, intent(out) :: f

        print * , "x = " , x
        print * , "f = " , f
        
        if (x(1) /= 2.0) error stop
        if (f /= 2.0) error stop
    end subroutine

end program