File: pass_array_by_data_10.f90

package info (click to toggle)
lfortran 0.45.0-1
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 46,332 kB
  • sloc: cpp: 137,068; f90: 51,260; python: 6,444; ansic: 4,277; yacc: 2,285; fortran: 806; sh: 524; makefile: 30; javascript: 15
file content (55 lines) | stat: -rw-r--r-- 1,288 bytes parent folder | download
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
module pass_array_by_data_10_pintrf_module
public :: OBJCON

contains

    subroutine OBJCON(x, out)
        implicit none
        real, intent(in) :: x(:)
        real, intent(out) :: out
    end subroutine OBJCON

end module pass_array_by_data_10_pintrf_module

subroutine cobyla(x, out)
    implicit none

    real, intent(inout) :: x(:)
    real, intent(out) :: out
    call evaluate(calcfc_internal, out)
    contains

    subroutine calcfc_internal(x_internal, out)
        implicit none
        real, intent(in) :: x_internal(:)
        real, intent(out) :: out

        out = sum(x_internal**2)
    end subroutine calcfc_internal

    subroutine evaluate(calcfc, out)
        use, non_intrinsic :: pass_array_by_data_10_pintrf_module, only : OBJCON
        implicit none
        procedure(OBJCON) :: calcfc
        real, intent(out) :: out
        call calcfc(x, out)
    end subroutine evaluate

end subroutine

program pass_array_by_data_10
    implicit none
    interface
        subroutine cobyla(x, out)
            real, intent(inout) :: x(:)
            real, intent(out) :: out
        end subroutine cobyla
    end interface

    real :: x(2) = [1.0, 2.0]
    real :: out
    call cobyla(x, out)
    print *, out
    if ( abs(out - 5.0 ) > 1e-8 ) error stop
end program