-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathCOOLHANS.SIF
158 lines (117 loc) · 3.33 KB
/
COOLHANS.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
***************************
* SET UP THE INITIAL DATA *
***************************
NAME COOLHANS
* Problem :
* *********
* A problem arising from the analysis of a Cooley-Hansen economy with
* loglinear approximation. The problem is to solve the matrix equation
*
* A * X * X + B * X + C = 0
*
* where A, B and C are known N times N matrices and X an unknown matrix
* of matching dimension. The instance considered here has N = 3.
* Source:
* S. Ceria, private communication, 1995.
* SIF input: Ph. Toint, Feb 1995.
* classification NQR2-RN-9-9
* order of the matrix equation
IE N 3
* matrix data
RE A1,1 0.0
RE A2,1 0.13725D-6
RE A3,1 0.0
RE A1,2 0.0
RE A2,2 937.62
RE A3,2 0.0
RE A1,3 0.0
RE A2,3 -42.207
RE A3,3 0.0
RE B1,1 0.0060893
RE B2,1 0.13880D-6
RE B3,1 -0.13877d-6
RE B1,2 -44.292
RE B2,2 -1886.0
RE B3,2 42.362
RE B1,3 2.0011
RE B2,3 42.362
RE B3,3 -2.0705
RE C1,1 0.0
RE C2,1 0.0
RE C3,1 0.0
RE C1,2 44.792
RE C2,2 948.21
RE C3,2 -42.684
RE C1,3 0.0
RE C2,3 0.0
RE C3,3 0.0
* useful constants
IE 1 1
IE 2 2
IE 3 3
VARIABLES
DO I 1 N
DO J 1 N
X X(I,J)
ND
GROUPS
DO K 1 N
DO L 1 N
DO M 1 N
ZE G(K,L) X(M,L) B(K,M)
ND
CONSTANTS
DO K 1 N
DO L 1 N
AM -C C(K,L) -1.0
Z COOLHANS G(K,L) -C
ND
BOUNDS
FR COOLHANS 'DEFAULT'
* DO K 1 N
* XX COOLHANS X(1,K) 0.0
* XX COOLHANS X(3,K) 0.0
* OD K
ELEMENT TYPE
EV 2PR XX YY
ELEMENT USES
DO K 1 N
DO L 1 N
DO M 1 N
XT E(K,M,L) 2PR
ZV E(K,M,L) XX X(K,M)
ZV E(K,M,L) YY X(M,L)
ND
GROUP USES
DO L 1 N
DO P 1 N
DO M 1 N
ZE G(1,L) E(P,M,L) A(1,P)
ND
DO L 1 N
DO P 1 N
DO M 1 N
ZE G(2,L) E(P,M,L) A(2,P)
ND
DO L 1 N
DO P 1 N
DO M 1 N
ZE G(3,L) E(P,M,L) A(3,P)
ND
OBJECT BOUND
LO COOLHANS 0.0
* Solution
*LO SOLTN 0.0
ENDATA
***********************
* SET UP THE FUNCTION *
* AND RANGE ROUTINES *
***********************
ELEMENTS COOLHANS
INDIVIDUALS
T 2PR
F XX * YY
G XX YY
G YY XX
H XX YY 1.0
ENDATA