-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathEXPLIN.SIF
107 lines (70 loc) · 1.91 KB
/
EXPLIN.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
***************************
* SET UP THE INITIAL DATA *
***************************
NAME EXPLIN
* Problem :
* *********
* A problem involving exponential and linear terms.
* SIF input: Ph. Toint, 1992.
* minor correction by Ph. Shott, Jan 1995.
* classification OBR2-AN-V-V
*IE N 12 $-PARAMETER
*IE M 6 $-PARAMETER
*IE N 120 $-PARAMETER original value
*IE M 10 $-PARAMETER original value
IE N 1200 $-PARAMETER
IE M 100 $-PARAMETER
* Define constants
IE 0 0
IE 1 1
IE 2 2
VARIABLES
DO I 1 N
X X(I)
ND
GROUPS
DO I 1 N
RI RI I
RM C RI - 10.0
ZN OBJ X(I) C
ND
CONSTANTS
BOUNDS
XU EXPLIN 'DEFAULT' 10.0
START POINT
ELEMENT TYPE
EV EXP X Y
ELEMENT USES
XT 'DEFAULT' EXP
DO I 1 M
I+ I+1 I 1
ZV E(I) X X(I)
ZV E(I) Y X(I+1)
OD I
GROUP USES
DO I 1 M
XE OBJ E(I)
OD I
OBJECT BOUND
LO EXPLIN 0.0
* Solution
ENDATA
***********************
* SET UP THE FUNCTION *
* AND RANGE ROUTINES *
***********************
ELEMENTS EXPLIN
TEMPORARIES
M EXP
R F
INDIVIDUALS
* Exponential of the product of X and Y divided by 10
T EXP
A F EXP( 0.1 * X * Y )
F F
G X 0.1 * Y * F
G Y 0.1 * X * F
H X X ( 0.1 * Y ) ** 2 * F
H Y Y ( 0.1 * X ) ** 2 * F
H Y X ( 0.1 + 0.01 * X * Y ) * F
ENDATA