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
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
! { dg-do run }
!
! PR fortran/51972
!
! Contributed by Damian Rouson
!
module surrogate_module
type ,abstract :: surrogate
end type
end module
module strategy_module
use surrogate_module
type :: strategy
end type
end module
module integrand_module
use surrogate_module
use strategy_module
implicit none
type ,abstract, extends(surrogate) :: integrand
class(strategy), allocatable :: quadrature
end type
end module integrand_module
module lorenz_module
use strategy_module
use integrand_module
implicit none
type ,extends(integrand) :: lorenz
real, dimension(:), allocatable :: state
contains
procedure ,public :: assign => assign_lorenz
end type
contains
type(lorenz) function constructor(initial_state, this_strategy)
real ,dimension(:) ,intent(in) :: initial_state
class(strategy) ,intent(in) :: this_strategy
constructor%state=initial_state
allocate (constructor%quadrature, source=this_strategy)
end function
subroutine assign_lorenz(lhs,rhs)
class(lorenz) ,intent(inout) :: lhs
class(integrand) ,intent(in) :: rhs
select type(rhs)
class is (lorenz)
allocate (lhs%quadrature, source=rhs%quadrature)
lhs%state=rhs%state
end select
end subroutine
end module lorenz_module
module runge_kutta_2nd_module
use surrogate_module,only : surrogate
use strategy_module ,only : strategy
use integrand_module,only : integrand
implicit none
type, extends(strategy) ,public :: runge_kutta_2nd
contains
procedure, nopass :: integrate
end type
contains
subroutine integrate(this)
class(surrogate) ,intent(inout) :: this
class(integrand) ,allocatable :: this_half
select type (this)
class is (integrand)
allocate (this_half, source=this)
end select
end subroutine
end module
program main
use lorenz_module
use runge_kutta_2nd_module ,only : runge_kutta_2nd, integrate
implicit none
type(runge_kutta_2nd) :: timed_lorenz_integrator
type(lorenz) :: attractor
attractor = constructor( [1., 1., 1.] , timed_lorenz_integrator)
call integrate(attractor)
end program main
|
