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
|
!----------------------------------------------------------------------
!----------------------------------------------------------------------
! cir : Homoclinic Bifurcation in an Electronic Circuit
! (the same equations as in demo tor)
!----------------------------------------------------------------------
!----------------------------------------------------------------------
SUBROUTINE FUNC(NDIM,U,ICP,PAR,IJAC,F,DFDU,DFDP)
! ---------- ----
IMPLICIT NONE
INTEGER, INTENT(IN) :: NDIM, ICP(*), IJAC
DOUBLE PRECISION, INTENT(IN) :: U(NDIM), PAR(*)
DOUBLE PRECISION, INTENT(OUT) :: F(NDIM)
DOUBLE PRECISION, INTENT(INOUT) :: DFDU(NDIM,NDIM), DFDP(NDIM,*)
DOUBLE PRECISION P
P=PAR(1)
F(1)= U(2)
F(2)= U(3)
F(3)= U(4)
F(4)= -P*U(3)-U(1)+U(1)**3
IF(IJAC.EQ.0)RETURN
DFDU(1,1)=0.0d0
DFDU(1,2)=1.0d0
DFDU(1,3)=0.0d0
DFDU(1,4)=0.0d0
DFDU(2,1)=0.0d0
DFDU(2,2)=0.0d0
DFDU(2,3)=1.0d0
DFDU(2,4)=0.0d0
DFDU(3,1)=0.0d0
DFDU(3,2)=0.0d0
DFDU(3,3)=0.0d0
DFDU(3,4)=1.0d0
DFDU(4,1)=-1.0+3.0d0*U(1)**2
DFDU(4,2)=0.0d0
DFDU(4,3)=-P
DFDU(4,4)=0.0d0
IF(IJAC.EQ.1)RETURN
! *Parameter derivatives
DFDP(1,1)=0.0d0
DFDP(2,1)=0.0d0
DFDP(3,1)=0.0d0
DFDP(4,1)=-U(3)
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
!----------------------------------------------------------------------
! Problem parameters (only PAR(1-9) are available to the user) :
PAR(1)=1.6 ! P
!----------------------------------------------------------------------
! If IEQUIB >0 put initial equilibrium in PAR(11+i), i=1,...,NDIM :
PAR(12) = 0.0
PAR(13) = 0.0
PAR(14) = 0.0
PAR(15) = 0.0
END SUBROUTINE STPNT
SUBROUTINE PVLS(NDIM,U,PAR)
! ---------- ----
IMPLICIT NONE
INTEGER, INTENT(IN) :: NDIM
DOUBLE PRECISION, INTENT(IN) :: U(NDIM)
DOUBLE PRECISION, INTENT(INOUT) :: PAR(*)
! Homoclinic bifurcations COMMON block needed here :
COMMON /BLHOM/ ITWIST,ISTART,IEQUIB,NFIXED,NPSI,NUNSTAB,NSTAB,NREV
INTEGER ITWIST,ISTART,IEQUIB,NFIXED,NPSI,NUNSTAB,NSTAB,NREV
INTEGER I
! If IEQUIB =0 put analytic equilibrium in PAR(11+i), i=1..NDIM
IF(IEQUIB.EQ.0)THEN
DO I=1,NDIM
PAR(11+I)= 0.0
ENDDO
ENDIF
END SUBROUTINE PVLS
SUBROUTINE BCND
END SUBROUTINE BCND
SUBROUTINE ICND
END SUBROUTINE ICND
SUBROUTINE FOPT
END SUBROUTINE FOPT
|
