-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathHS56.SIF
137 lines (97 loc) · 2.56 KB
/
HS56.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
***************************
* SET UP THE INITIAL DATA *
***************************
NAME HS56
* Problem :
* *********
* Source: problem 56 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: A.R. Conn, April 1990
* classification OOR2-AN-7-4
* some useful parameters, including N, the number of variables.
IE N 7
IE 1 1
IE 2 2
IE 3 3
IE 4 4
AE PAR(1) 4.2
AE PAR(2) 4.2
AE PAR(3) 4.2
AE PAR(4) 7.2
VARIABLES
DO I 1 N
X X(I)
ND
GROUPS
N OBJ
XE CON(1) X1 1.0
XE CON(2) X2 1.0
XE CON(3) X3 1.0
XE CON(4) X1 1.0 X2 2.0
XE CON(4) X3 2.0
BOUNDS
FR HS56 'DEFAULT'
START POINT
* 0.05097397 = arcsin( sqrt( 1/4.2 ) )
* 0.09851108 = arcsin( sqrt( 5/7.2 ) )
* Note there is a typo in Hock & Schittkowski
* (incorrect leading decimal point).
XV HS56 'DEFAULT' 1.0
HS56 X4 0.50973968
HS56 X5 0.50973968
HS56 X6 0.50973968
HS56 X7 0.98511078
ELEMENT TYPE
EV 3PROD V1 V2
EV 3PROD V3
EV PSNSQ V1
EP PSNSQ P
ELEMENT USES
T E1 3PROD
V E1 V1 X1
V E1 V2 X2
V E1 V3 X3
DO I 1 4
IA J I 1
IA K I 3
XT E(J) PSNSQ
ZV E(J) V1 X(K)
ZP E(J) P PAR(I)
ND
GROUP USES
E OBJ E1 -1.0
DO I 1 4
IA J I 1
XE CON(I) E(J)
ND
OBJECT BOUND
* Solution
*LO SOLTN -3.456
ENDATA
***********************
* SET UP THE FUNCTION *
* AND RANGE ROUTINES *
***********************
ELEMENTS HS56
TEMPORARIES
R SUM
M SIN
M COS
INDIVIDUALS
T 3PROD
F V1 * V2 * V3
G V1 V2 * V3
G V2 V1 * V3
G V3 V1 * V2
H V1 V2 V3
H V1 V3 V2
H V2 V3 V1
T PSNSQ
A SUM V1 + V1
F -P * SIN(V1)**2
G V1 -P * SIN(SUM)
H V1 V1 -2.0 * P * COS(SUM)
ENDATA