-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathHS105.SIF
400 lines (307 loc) · 8.64 KB
/
HS105.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
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
***************************
* SET UP THE INITIAL DATA *
***************************
NAME HS105
* Problem :
* *********
* Source: problem 105 in
* W. Hock and K. Schittkowski,
* "Test examples for nonlinear programming codes",
* Lectures Notes in Economics and Mathematical Systems 187, Springer
* Verlag, Heidelberg, 1981.
* SIF input: Nick Gould, August 1991.
* classification OLR2-AY-8-1
* Number of variables
IE N 8
* Other useful parameters
IE 1 1
IE 235 235
* Problem data
AE Y1 95.0
AE Y2 105.0
IE LOW 3
IE UP 6
DO I LOW UP
AE Y(I) 110.0
ND
IE LOW 7
IE UP 10
DO I LOW UP
AE Y(I) 115.0
ND
IE LOW 11
IE UP 25
DO I LOW UP
AE Y(I) 120.0
ND
IE LOW 26
IE UP 40
DO I LOW UP
AE Y(I) 125.0
ND
IE LOW 41
IE UP 55
DO I LOW UP
AE Y(I) 130.0
ND
IE LOW 56
IE UP 68
DO I LOW UP
AE Y(I) 135.0
ND
IE LOW 69
IE UP 89
DO I LOW UP
AE Y(I) 140.0
ND
IE LOW 90
IE UP 101
DO I LOW UP
AE Y(I) 145.0
ND
IE LOW 102
IE UP 118
DO I LOW UP
AE Y(I) 150.0
ND
IE LOW 119
IE UP 122
DO I LOW UP
AE Y(I) 155.0
ND
IE LOW 123
IE UP 142
DO I LOW UP
AE Y(I) 160.0
ND
IE LOW 143
IE UP 150
DO I LOW UP
AE Y(I) 165.0
ND
IE LOW 151
IE UP 167
DO I LOW UP
AE Y(I) 170.0
ND
IE LOW 168
IE UP 175
DO I LOW UP
AE Y(I) 175.0
ND
IE LOW 176
IE UP 181
DO I LOW UP
AE Y(I) 180.0
ND
IE LOW 182
IE UP 187
DO I LOW UP
AE Y(I) 185.0
ND
IE LOW 188
IE UP 194
DO I LOW UP
AE Y(I) 190.0
ND
IE LOW 195
IE UP 198
DO I LOW UP
AE Y(I) 195.0
ND
IE LOW 199
IE UP 201
DO I LOW UP
AE Y(I) 200.0
ND
IE LOW 202
IE UP 204
DO I LOW UP
AE Y(I) 205.0
ND
IE LOW 205
IE UP 212
DO I LOW UP
AE Y(I) 210.0
ND
AE Y213 215.0
IE LOW 214
IE UP 219
DO I LOW UP
AE Y(I) 220.0
ND
IE LOW 220
IE UP 224
DO I LOW UP
AE Y(I) 230.0
ND
AE Y225 235.0
IE LOW 226
IE UP 232
DO I LOW UP
AE Y(I) 240.0
ND
AE Y233 245.0
IE LOW 234
IE UP 235
DO I LOW UP
AE Y(I) 250.0
ND
VARIABLES
DO I 1 N
X X(I)
ND
GROUPS
DO I 1 235
XN OBJ(I)
ND
G C1 X1 -1.0D+0 X2 -1.0D+0
CONSTANTS
HS105 C1 -1.0D+0
BOUNDS
LO HS105 X1 0.001
UP HS105 X1 0.499
LO HS105 X2 0.001
UP HS105 X2 0.499
LO HS105 X3 100.0
UP HS105 X3 180.0
LO HS105 X4 130.0
UP HS105 X4 210.0
LO HS105 X5 170.0
UP HS105 X5 240.0
LO HS105 X6 5.0
UP HS105 X6 25.0
LO HS105 X7 5.0
UP HS105 X7 25.0
LO HS105 X8 5.0
UP HS105 X8 25.0
START POINT
HS105 X1 0.1 X2 0.2
HS105 X3 100.0 X4 125.0
HS105 X5 175.0 X6 11.2
HS105 X7 13.2 X8 15.8
HS105SOL X1 0.4128928 X2 0.4033526
HS105SOL X3 131.2613 X4 164.3135
HS105SOL X5 217.4222 X6 12.28018
HS105SOL X7 15.77170 X8 20.74682
ELEMENT TYPE
* ABI
EV ABI V1 V2
EV ABI V3
EP ABI YI
* CI
EV CI X1 X2
EV CI X5 X8
IV CI V1 V2
IV CI V3
EP CI YI
ELEMENT USES
DO I 1 235
XT A(I) ABI
ZV A(I) V1 X1
ZV A(I) V2 X6
ZV A(I) V3 X3
ZP A(I) YI Y(I)
XT B(I) ABI
ZV B(I) V1 X2
ZV B(I) V2 X7
ZV B(I) V3 X4
ZP B(I) YI Y(I)
XT C(I) CI
ZV C(I) X1 X1
ZV C(I) X2 X2
ZV C(I) X5 X5
ZV C(I) X8 X8
ZP C(I) YI Y(I)
ND
GROUP TYPE
* LOG
GV LOG T
GROUP USES
DO I 1 235
XT OBJ(I) LOG
XE OBJ(I) A(I) B(I)
XE OBJ(I) C(I)
ND
OBJECT BOUND
* Solution
*LO SOLTN 1138.416240
ENDATA
***********************
* SET UP THE FUNCTION *
* AND RANGE ROUTINES *
***********************
ELEMENTS HS105
TEMPORARIES
R R
R D
R E
R DDV2
R DDV3
R DEV2
R DEV3
M EXP
INDIVIDUALS
* Element type ABI
T ABI
A R V1 / V2
A D ( YI - V3 ) / V2
A E EXP( - 5.0D-1 * D * D )
A DDV2 - D / V2
A DEV2 E * ( - D ) * DDV2
A DDV3 - 1.0D+0 / V2
A DEV3 E * ( - D ) * DDV3
F R * E
G V1 E / V2
G V2 ( D * D - 1.0D+0 ) * R * E / V2
G V3 D * R * E / V2
H V1 V2 ( DEV2 - E / V2 ) / V2
H V1 V3 DEV3 / V2
H V2 V2 ( 2.0D+0 * D * DDV2 * E +
H+ ( D * D - 1.0D+0 ) * ( DEV2 -
H+ 2.0D+0 * E / V2 ) ) * R / V2
H V2 V3 ( DDV2 * E + D * DEV2 -
H+ 2.0D+0 * D * E / V2 ) * R / V2
H V3 V3 ( DDV3 * E + D * DEV3 ) * R / V2
* Element type CI
T CI
R V1 X1 -1.0 X1 -1.0
R V2 X8 1.0
R V3 X5 1.0
A R ( 1.0D+0 + V1 ) / V2
A D ( YI - V3 ) / V2
A E EXP( - 5.0D-1 * D * D )
A DDV2 - D / V2
A DEV2 E * ( - D ) * DDV2
A DDV3 - 1.0D+0 / V2
A DEV3 E * ( - D ) * DDV3
F R * E
G V1 E / V2
G V2 ( D * D - 1.0D+0 ) * R * E / V2
G V3 D * R * E / V2
H V1 V2 ( DEV2 - E / V2 ) / V2
H V1 V3 DEV3 / V2
H V2 V2 ( 2.0D+0 * D * DDV2 * E +
H+ ( D * D - 1.0D+0 ) * ( DEV2 -
H+ 2.0D+0 * E / V2 ) ) * R / V2
H V2 V3 ( DDV2 * E + D * DEV2 -
H+ 2.0D+0 * D * E / V2 ) * R / V2
H V3 V3 ( DDV3 * E + D * DEV3 ) * R / V2
ENDATA
***********************
* SET UP THE FUNCTION *
* AND RANGE ROUTINES *
***********************
GROUPS HS105
TEMPORARIES
R P
M LOG
GLOBALS
A P 3.9894228040143270D-01
INDIVIDUALS
* group type log
T LOG
F - LOG( P * T )
G - 1 / T
H 1 / T ** 2
ENDATA