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
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
|
!----------------------------------------------------------------------
!----------------------------------------------------------------------
! pvl : Setting functional parameter values (for Bratu's equation)
!----------------------------------------------------------------------
!----------------------------------------------------------------------
SUBROUTINE FUNC(NDIM,U,ICP,PAR,IJAC,F,DFDU,DFDP)
! ---------- ----
IMPLICIT NONE
INTEGER, INTENT(IN) :: NDIM, IJAC, ICP(*)
DOUBLE PRECISION, INTENT(IN) :: U(NDIM), PAR(*)
DOUBLE PRECISION, INTENT(OUT) :: F(NDIM)
DOUBLE PRECISION, INTENT(INOUT) :: DFDU(NDIM,*), DFDP(NDIM,*)
DOUBLE PRECISION E
E=EXP(U(1))
F(1)=U(2)
F(2)=-PAR(1)*E
END SUBROUTINE FUNC
SUBROUTINE STPNT(NDIM,U,PAR,T)
! ---------- -----
IMPLICIT NONE
INTEGER, INTENT(IN) :: NDIM
DOUBLE PRECISION, INTENT(INOUT) :: U(NDIM), PAR(*)
DOUBLE PRECISION, INTENT(IN) :: T
PAR(1)=0
U(1)=0.0
U(2)=0.0
END SUBROUTINE STPNT
SUBROUTINE BCND(NDIM,PAR,ICP,NBC,U0,U1,FB,IJAC,DBC)
! ---------- ----
IMPLICIT NONE
INTEGER, INTENT(IN) :: NDIM, ICP(*), NBC, IJAC
DOUBLE PRECISION, INTENT(IN) :: PAR(*), U0(NDIM), U1(NDIM)
DOUBLE PRECISION, INTENT(OUT) :: FB(NBC)
DOUBLE PRECISION, INTENT(INOUT) :: DBC(NBC,*)
FB(1)=U0(1)
FB(2)=U1(1)
END SUBROUTINE BCND
DOUBLE PRECISION FUNCTION GETU2(U,NDX,NTST,NCOL)
! ------ --------- -------- -----
INTEGER, INTENT(IN) :: NDX,NCOL,NTST
DOUBLE PRECISION, INTENT(IN) :: U(NDX,0:NCOL*NTST)
GETU2 = U(2,0)
END FUNCTION GETU2
SUBROUTINE PVLS(NDIM,U,PAR)
! ---------- ----
IMPLICIT NONE
INTEGER, INTENT(IN) :: NDIM
DOUBLE PRECISION, INTENT(IN) :: U(NDIM)
DOUBLE PRECISION, INTENT(INOUT) :: PAR(*)
DOUBLE PRECISION, EXTERNAL :: GETP,GETU2
INTEGER NDX,NCOL,NTST
!----------------------------------------------------------------------
! NOTE :
! Parameters set in this subroutine should be considered as ``solution
! measures'' and be used for output purposes only.
!
! They should never be used as `true'' continuation parameters.
!
! They may, however, be added as ``over-specified parameters'' in the
! parameter list associated with the AUTO-Constant NICP, in order to
! print their values on the screen and in the ``p.xxx file.
!
! They may also appear in the list associated with AUTO-constant NUZR.
!
!----------------------------------------------------------------------
! For algebraic problems the argument U is, as usual, the state vector.
! For differential equations the argument U represents the approximate
! solution on the entire interval [0,1]. In this case its values can
! be accessed indirectly by calls to GETP, as illustrated below, or
! by obtaining NDIM, NCOL, NTST via GETP and then dimensioning U as
! U(NDIM,0:NCOL*NTST) in a seperate subroutine that is called by PVLS.
!----------------------------------------------------------------------
! Set PAR(2) equal to the L2-norm of U(1)
PAR(2)=GETP('NRM',1,U)
! Set PAR(3) equal to the minimum of U(2)
PAR(3)=GETP('MIN',2,U)
! Set PAR(4) equal to the value of U(2) at the left boundary.
PAR(4)=GETP('BV0',2,U)
! Set PAR(5) equal to the value of U(2) at the left boundary using
! another method
NDX=NINT(GETP('NDX',0,U))
NTST=NINT(GETP('NTST',0,U))
NCOL=NINT(GETP('NCOL',0,U))
PAR(5)=GETU2(U,NDX,NTST,NCOL)
!----------------------------------------------------------------------
! The first argument of GETP may be one of the following:
! 'NRM' (L2-norm), 'MAX' (maximum),
! 'INT' (integral), 'BV0 (left boundary value),
! 'MIN' (minimum), 'BV1' (right boundary value).
! 'MNT' (t value for minimum)
! 'MXT' (t value for maximum)
! 'NDIM', 'NDX' (effective (active) number of dimensions)
! 'NTST' (NTST from constant file)
! 'NCOL' (NCOL from constant file)
! 'NBC' (active NBC)
! 'NINT' (active NINT)
! 'DTM' (delta t for all t values, I=1...NTST)
! 'WINT' (integration weights used for interpolation, I=0...NCOL)
!
! Also available are
! 'STP' (Pseudo-arclength step size used).
! 'FLD' (`Fold function', which vanishes at folds).
! 'BIF' (`Bifurcation function', which vanishes at singular points).
! 'HBF' (`Hopf function'; which vanishes at Hopf points).
! 'SPB' ( Function which vanishes at secondary periodic bifurcations).
! 'EIG' ( Eigenvalues/multipliers, I=1...2*NDIM, alternates real/imag parts).
! 'STA' ( Number of stable eigenvalues/multipliers).
!----------------------------------------------------------------------
END SUBROUTINE PVLS
SUBROUTINE ICND
END SUBROUTINE ICND
SUBROUTINE FOPT
END SUBROUTINE FOPT
|
