-
Notifications
You must be signed in to change notification settings - Fork 0
/
BDVALUE.SIF
131 lines (92 loc) · 2.62 KB
/
BDVALUE.SIF
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
***************************
* SET UP THE INITIAL DATA *
***************************
NAME BDVALUE
* Problem :
* *********
* The Boundary Value problem.
* This is a nonlinear equations problems.
* Source: problem 28 in
* J.J. More', B.S. Garbow and K.E. Hillstrom,
* "Testing Unconstrained Optimization Software",
* ACM Transactions on Mathematical Software, vol. 7(1), pp. 17-41, 1981.
* SIF input: Ph. Toint, Dec 1989.
* classification NOR2-MN-V-V
* NDP is the number of discretization points (variable).
* The number of free variables is NDP-2
*IE NDP 12 $-PARAMETER original value
*IE NDP 52 $-PARAMETER
*IE NDP 102 $-PARAMETER
*IE NDP 502 $-PARAMETER
*IE NDP 1002 $-PARAMETER
IE NDP 5002 $-PARAMETER
* Define useful parameters
IE 1 1
IE 2 2
RE 1.0 1.0
IA NDP-1 NDP -1
RI RNDP-1 NDP-1
R/ H 1.0 RNDP-1
R* H2 H H
RM HALFH2 H2 0.5D0
VARIABLES
DO I 1 NDP
X X(I)
ND
GROUPS
DO I 2 NDP-1
IA I-1 I -1
IA I+1 I 1
XE G(I) X(I-1) -1.0 X(I) 2.0
XE G(I) X(I+1) -1.0
ND
BOUNDS
FR BDVALUE 'DEFAULT'
XX BDVALUE X(1) 0.0
XX BDVALUE X(NDP) 0.0
START POINT
X BDVALUE X(1) 0.0
DO I 2 NDP-1
IA I-1 I -1
RI RI-1 I-1
R* IH RI-1 H
RA IH-1 IH -1.0
R* TI IH IH-1
Z BDVALUE X(I) TI
ND
X BDVALUE X(NDP) 0.0
ELEMENT TYPE
EV WCUBE V
EP WCUBE B
ELEMENT USES
T 'DEFAULT' WCUBE
DO I 2 NDP-1
RI REALI I
RA RI+1 REALI 1.0
R* IH REALI H
RA IH+1 IH 1.0
ZV E(I) V X(I)
ZP E(I) B IH+1
ND
GROUP USES
DO I 2 NDP-1
ZE G(I) E(I) HALFH2
ND
OBJECT BOUND
* Solution
*LO SOLTN 0.0
ENDATA
***********************
* SET UP THE FUNCTION *
* AND RANGE ROUTINES *
***********************
ELEMENTS BDVALUE
TEMPORARIES
R VPLUSB
INDIVIDUALS
T WCUBE
A VPLUSB V + B
F VPLUSB**3
G V 3.0 * VPLUSB**2
H V V 6.0 * VPLUSB
ENDATA