e942d1ce32473ec2473797828b3c9e7ecb37da4b
[scilab.git] / scilab / modules / optimization / tests / nonreg_tests / bug_4638.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2009 - DIGITEO - Allan CORNET
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 // <-- Non-regression test for bug 4638 -->
8 //
9 // <-- Bugzilla URL -->
10 // http://bugzilla.scilab.org/show_bug.cgi?id=4638
11 //
12 // <-- Short Description -->
13 // this example with optim freezes scilab
14 // =============================================================================
15 function A = build_matrix_A()
16 A = spzeros(36,36);
17 A(    1,    1) =         12.356155 ;
18 A(    1,    2) =         1.8561553 ;
19 A(    1,    5) =       - 5.25 ;
20 A(    1,    7) =       - 1.8561553 ;
21 A(    1,    8) =       - 1.8561553 ;
22 A(    2,    1) =         1.8561553 ;
23 A(    2,    2) =         7.1061553 ;
24 A(    2,    4) =       - 5.25 ;
25 A(    2,    7) =       - 1.8561553 ;
26 A(    2,    8) =       - 1.8561553 ;
27 A(    3,    3) =         8.9623106 ;
28 A(    3,    5) =       - 1.8561553 ;
29 A(    3,    6) =         1.8561553 ;
30 A(    3,    7) =       - 5.25 ;
31 A(    4,    2) =       - 5.25 ;
32 A(    4,    4) =         8.9623106 ;
33 A(    4,    5) =         1.8561553 ;
34 A(    4,    6) =       - 1.8561553 ;
35 A(    5,    1) =       - 5.25 ;
36 A(    5,    3) =       - 1.8561553 ;
37 A(    5,    4) =         1.8561553 ;
38 A(    5,    5) =         14.212311 ;
39 A(    5,    9) =       - 5.25 ;
40 A(    5,   11) =       - 1.8561553 ;
41 A(    5,   12) =       - 1.8561553 ;
42 A(    6,    3) =         1.8561553 ;
43 A(    6,    4) =       - 1.8561553 ;
44 A(    6,    6) =         8.9623106 ;
45 A(    6,    8) =       - 5.25 ;
46 A(    6,   11) =       - 1.8561553 ;
47 A(    6,   12) =       - 1.8561553 ;
48 A(    7,    1) =       - 1.8561553 ;
49 A(    7,    2) =       - 1.8561553 ;
50 A(    7,    3) =       - 5.25 ;
51 A(    7,    7) =         14.212311 ;
52 A(    7,    9) =       - 1.8561553 ;
53 A(    7,   10) =         1.8561553 ;
54 A(    7,   11) =       - 5.25 ;
55 A(    8,    1) =       - 1.8561553 ;
56 A(    8,    2) =       - 1.8561553 ;
57 A(    8,    6) =       - 5.25 ;
58 A(    8,    8) =         8.9623106 ;
59 A(    8,    9) =         1.8561553 ;
60 A(    8,   10) =       - 1.8561553 ;
61 A(    9,    5) =       - 5.25 ;
62 A(    9,    7) =       - 1.8561553 ;
63 A(    9,    8) =         1.8561553 ;
64 A(    9,    9) =         14.212311 ;
65 A(    9,   13) =       - 5.25 ;
66 A(    9,   15) =       - 1.8561553 ;
67 A(    9,   16) =       - 1.8561553 ;
68 A(   10,    7) =         1.8561553 ;
69 A(   10,    8) =       - 1.8561553 ;
70 A(   10,   10) =         8.9623106 ;
71 A(   10,   12) =       - 5.25 ;
72 A(   10,   15) =       - 1.8561553 ;
73 A(   10,   16) =       - 1.8561553 ;
74 A(   11,    5) =       - 1.8561553 ;
75 A(   11,    6) =       - 1.8561553 ;
76 A(   11,    7) =       - 5.25 ;
77 A(   11,   11) =         14.212311 ;
78 A(   11,   13) =       - 1.8561553 ;
79 A(   11,   14) =         1.8561553 ;
80 A(   11,   15) =       - 5.25 ;
81 A(   12,    5) =       - 1.8561553 ;
82 A(   12,    6) =       - 1.8561553 ;
83 A(   12,   10) =       - 5.25 ;
84 A(   12,   12) =         8.9623106 ;
85 A(   12,   13) =         1.8561553 ;
86 A(   12,   14) =       - 1.8561553 ;
87 A(   13,    9) =       - 5.25 ;
88 A(   13,   11) =       - 1.8561553 ;
89 A(   13,   12) =         1.8561553 ;
90 A(   13,   13) =         14.212311 ;
91 A(   13,   17) =       - 5.25 ;
92 A(   13,   19) =       - 1.8561553 ;
93 A(   13,   20) =       - 1.8561553 ;
94 A(   14,   11) =         1.8561553 ;
95 A(   14,   12) =       - 1.8561553 ;
96 A(   14,   14) =         8.9623106 ;
97 A(   14,   16) =       - 5.25 ;
98 A(   14,   19) =       - 1.8561553 ;
99 A(   14,   20) =       - 1.8561553 ;
100 A(   15,    9) =       - 1.8561553 ;
101 A(   15,   10) =       - 1.8561553 ;
102 A(   15,   11) =       - 5.25 ;
103 A(   15,   15) =         14.212311 ;
104 A(   15,   17) =       - 1.8561553 ;
105 A(   15,   18) =         1.8561553 ;
106 A(   15,   19) =       - 5.25 ;
107 A(   16,    9) =       - 1.8561553 ;
108 A(   16,   10) =       - 1.8561553 ;
109 A(   16,   14) =       - 5.25 ;
110 A(   16,   16) =         8.9623106 ;
111 A(   16,   17) =         1.8561553 ;
112 A(   16,   18) =       - 1.8561553 ;
113 A(   17,   13) =       - 5.25 ;
114 A(   17,   15) =       - 1.8561553 ;
115 A(   17,   16) =         1.8561553 ;
116 A(   17,   17) =         14.212311 ;
117 A(   17,   21) =       - 5.25 ;
118 A(   17,   23) =       - 1.8561553 ;
119 A(   17,   24) =       - 1.8561553 ;
120 A(   18,   15) =         1.8561553 ;
121 A(   18,   16) =       - 1.8561553 ;
122 A(   18,   18) =         8.9623106 ;
123 A(   18,   20) =       - 5.25 ;
124 A(   18,   23) =       - 1.8561553 ;
125 A(   18,   24) =       - 1.8561553 ;
126 A(   19,   13) =       - 1.8561553 ;
127 A(   19,   14) =       - 1.8561553 ;
128 A(   19,   15) =       - 5.25 ;
129 A(   19,   19) =         14.212311 ;
130 A(   19,   21) =       - 1.8561553 ;
131 A(   19,   22) =         1.8561553 ;
132 A(   19,   23) =       - 5.25 ;
133 A(   20,   13) =       - 1.8561553 ;
134 A(   20,   14) =       - 1.8561553 ;
135 A(   20,   18) =       - 5.25 ;
136 A(   20,   20) =         8.9623106 ;
137 A(   20,   21) =         1.8561553 ;
138 A(   20,   22) =       - 1.8561553 ;
139 A(   21,   17) =       - 5.25 ;
140 A(   21,   19) =       - 1.8561553 ;
141 A(   21,   20) =         1.8561553 ;
142 A(   21,   21) =         14.212311 ;
143 A(   21,   25) =       - 5.25 ;
144 A(   21,   27) =       - 1.8561553 ;
145 A(   21,   28) =       - 1.8561553 ;
146 A(   22,   19) =         1.8561553 ;
147 A(   22,   20) =       - 1.8561553 ;
148 A(   22,   22) =         8.9623106 ;
149 A(   22,   24) =       - 5.25 ;
150 A(   22,   27) =       - 1.8561553 ;
151 A(   22,   28) =       - 1.8561553 ;
152 A(   23,   17) =       - 1.8561553 ;
153 A(   23,   18) =       - 1.8561553 ;
154 A(   23,   19) =       - 5.25 ;
155 A(   23,   23) =         14.212311 ;
156 A(   23,   25) =       - 1.8561553 ;
157 A(   23,   26) =         1.8561553 ;
158 A(   23,   27) =       - 5.25 ;
159 A(   24,   17) =       - 1.8561553 ;
160 A(   24,   18) =       - 1.8561553 ;
161 A(   24,   22) =       - 5.25 ;
162 A(   24,   24) =         8.9623106 ;
163 A(   24,   25) =         1.8561553 ;
164 A(   24,   26) =       - 1.8561553 ;
165 A(   25,   21) =       - 5.25 ;
166 A(   25,   23) =       - 1.8561553 ;
167 A(   25,   24) =         1.8561553 ;
168 A(   25,   25) =         14.212311 ;
169 A(   25,   29) =       - 5.25 ;
170 A(   25,   31) =       - 1.8561553 ;
171 A(   25,   32) =       - 1.8561553 ;
172 A(   26,   23) =         1.8561553 ;
173 A(   26,   24) =       - 1.8561553 ;
174 A(   26,   26) =         8.9623106 ;
175 A(   26,   28) =       - 5.25 ;
176 A(   26,   31) =       - 1.8561553 ;
177 A(   26,   32) =       - 1.8561553 ;
178 A(   27,   21) =       - 1.8561553 ;
179 A(   27,   22) =       - 1.8561553 ;
180 A(   27,   23) =       - 5.25 ;
181 A(   27,   27) =         14.212311 ;
182 A(   27,   29) =       - 1.8561553 ;
183 A(   27,   30) =         1.8561553 ;
184 A(   27,   31) =       - 5.25 ;
185 A(   28,   21) =       - 1.8561553 ;
186 A(   28,   22) =       - 1.8561553 ;
187 A(   28,   26) =       - 5.25 ;
188 A(   28,   28) =         8.9623106 ;
189 A(   28,   29) =         1.8561553 ;
190 A(   28,   30) =       - 1.8561553 ;
191 A(   29,   25) =       - 5.25 ;
192 A(   29,   27) =       - 1.8561553 ;
193 A(   29,   28) =         1.8561553 ;
194 A(   29,   29) =         14.212311 ;
195 A(   29,   33) =       - 5.25 ;
196 A(   29,   35) =       - 1.8561553 ;
197 A(   29,   36) =       - 1.8561553 ;
198 A(   30,   27) =         1.8561553 ;
199 A(   30,   28) =       - 1.8561553 ;
200 A(   30,   30) =         8.9623106 ;
201 A(   30,   32) =       - 5.25 ;
202 A(   30,   35) =       - 1.8561553 ;
203 A(   30,   36) =       - 1.8561553 ;
204 A(   31,   25) =       - 1.8561553 ;
205 A(   31,   26) =       - 1.8561553 ;
206 A(   31,   27) =       - 5.25 ;
207 A(   31,   31) =         14.212311 ;
208 A(   31,   33) =       - 1.8561553 ;
209 A(   31,   34) =         1.8561553 ;
210 A(   31,   35) =       - 5.25 ;
211 A(   32,   25) =       - 1.8561553 ;
212 A(   32,   26) =       - 1.8561553 ;
213 A(   32,   30) =       - 5.25 ;
214 A(   32,   32) =         8.9623106 ;
215 A(   32,   33) =         1.8561553 ;
216 A(   32,   34) =       - 1.8561553 ;
217 A(   33,   29) =       - 5.25 ;
218 A(   33,   31) =       - 1.8561553 ;
219 A(   33,   32) =         1.8561553 ;
220 A(   33,   33) =         12.356155 ;
221 A(   33,   34) =       - 1.8561553 ;
222 A(   34,   31) =         1.8561553 ;
223 A(   34,   32) =       - 1.8561553 ;
224 A(   34,   33) =       - 1.8561553 ;
225 A(   34,   34) =         7.1061553 ;
226 A(   34,   36) =       - 5.25 ;
227 A(   35,   29) =       - 1.8561553 ;
228 A(   35,   30) =       - 1.8561553 ;
229 A(   35,   31) =       - 5.25 ;
230 A(   35,   35) =         8.9623106 ;
231 A(   36,   29) =       - 1.8561553 ;
232 A(   36,   30) =       - 1.8561553 ;
233 A(   36,   34) =       - 5.25 ;
234 A(   36,   36) =         8.9623106 ;
235 A = A * 1.0D+08;
236 endfunction
237 function [y, dy, ind] = optim_test(x,ind)
238   printf('hello from optim_test\n');
239   y = sum(x.^2);
240   dy = 2*x;
241   A = build_matrix_A();
242   Ainv = inv(A);
243 endfunction
244 [f_opt, x_opt] = optim(optim_test,ones(10,1));
245 hellofromoptim_test
246
247 hellofromoptim_test
248
249 hellofromoptim_test
250
251 // =============================================================================