-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathSCURLY10.SIF
126 lines (84 loc) · 2.46 KB
/
SCURLY10.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
***************************
* SET UP THE INITIAL DATA *
***************************
NAME SCURLY10
* Problem :
* --------
* A banded function with semi-bandwidth 10 and
* negative curvature near the starting point.
* NB: scaled version of CURLY10
* Source: Nick Gould
* SIF input: Nick Gould, September 1997.
* classification SUR2-AN-V-0
* Number of variables
*IE N 10 $-PARAMETER
*IE N 100 $-PARAMETER
*IE N 1000 $-PARAMETER original value
IE N 10000 $-PARAMETER
*IE N 100000 $-PARAMETER
*IE N 1000000 $-PARAMETER
* Semi-bandwidth
IE K 10
* ratio of smallest to largest scale factors will be exp(scal)
RE SCAL 12.0
* other parameter definitions
IE 1 1
I- N-K N K
IA N-K+1 N-K 1
RI RN N
RA RN-1 RN -1
RA RN+1 RN 1
VARIABLES
DO I 1 N
IA I-1 I -1
RI RI-1 I-1
R/ RAT RI-1 RN-1
R* ARG RAT SCAL
A( S(I) EXP ARG
Z X(I)
ND
GROUPS
DO I 1 N-K
I+ I+K I K
DO J I I+K
ZN Q(I) X(J) S(J)
ND
DO I N-K+1 N
DO J I N
ZN Q(I) X(J) S(J)
ND
BOUNDS
FR SCURLY10 'DEFAULT'
START POINT
* start with X(I) = 0.0001/N+1.
DO I 1 N
RI RI I
R/ T RI RN+1
RM T T 0.0001
AM T T S(I)
ZV SCURLY10 X(I) T
ND
GROUP TYPE
GV P4 GVAR
GROUP USES
XT 'DEFAULT' P4
OBJECT BOUND
* Solution
*ZL SOLTN -1.003163D+5 $ (n=1000)
ENDATA
*********************
* SET UP THE GROUPS *
* ROUTINE *
*********************
GROUPS SCURLY10
TEMPORARIES
R APB
INDIVIDUALS
T P4
A APB 2.0D+1
F GVAR * ( GVAR * ( GVAR ** 2 - APB )
F+ - 1.0D-1 )
G 2.0D+0 * GVAR * ( 2.0D+0 * GVAR ** 2
G+ - APB ) - 1.0D-1
H 1.2D+1 * GVAR ** 2 - 2.0D+0 * APB
ENDATA