Add *.win.dia.ref in optimization module.
[scilab.git] / scilab / modules / optimization / tests / nonreg_tests / bug_9690.win.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2013 - Scilab Enterprises - Paul Bignier
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 //
8 // <-- CLI SHELL MODE -->
9 //
10 // <-- ENGLISH IMPOSED -->
11 //
12 // <-- Non-regression test for bug 9690 -->
13 //
14 // <-- Bugzilla URL -->
15 // http://bugzilla.scilab.org/show_bug.cgi?id=9690
16 //
17 // <-- Short Description -->
18 // optim(): option "imp"=5 could crash Scilab
19 //
20 function f = rosenbrock(x)
21     f = 100.0 * (x(2)-x(1)^2)^2 + (1-x(1))^2;
22 endfunction
23 function [f, g, ind] = rosenbrockCost(x, ind)
24     if ((ind == 1) | (ind == 4)) then
25         f = rosenbrock ( x );
26     end
27     if ((ind == 1) | (ind == 4)) then
28         g = derivative ( rosenbrock , x(:) );
29     end
30 endfunction
31 x0 = [-1.2 1.0];
32 lines(0);
33 for impval = 1:5
34     [ fopt , xopt ] = optim ( rosenbrockCost , x0 , "gc" , imp=impval);
35 end
36  n1qn3: entry point
37      dimension of the problem (n):     2
38      absolute precision on x (dxmin): 0.22D-15
39      expected decrease for f (df1): 0.10D+01
40      relative precision on g (epsg): 0.95D-18
41      maximal number of iterations (niter):   100
42      maximal number of simulations (nsim):   100
43      printing level (impres):   1
44      allocated memory (nrz) :     58
45      used memory :                58
46      number of updates :          10
47
48      f         =  0.24200000D+02
49      norm of g =  0.23286769D+03
50      nlis0   fin sur tmin        pas            fonctions     derivees
51      nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
52      nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
53  -------------------------------------------------------------------------------
54  n1qn3: output mode is  6
55      number of iterations:   42
56      number of simulations:     53
57      realized relative precision on g:  0.21D-15
58      norm of x =  0.14142135D+01
59      f         =  0.53783484D-16
60      norm of x =  0.14924691D-12
61 Optimization stops because too small variations for x.
62
63  n1qn3: entry point
64      dimension of the problem (n):     2
65      absolute precision on x (dxmin): 0.22D-15
66      expected decrease for f (df1): 0.10D+01
67      relative precision on g (epsg): 0.95D-18
68      maximal number of iterations (niter):   100
69      maximal number of simulations (nsim):   100
70      printing level (impres):   2
71      allocated memory (nrz) :     58
72      used memory :                58
73      number of updates :          10
74
75      f         =  0.24200000D+02
76      norm of g =  0.23286769D+03
77  n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
78  n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
79  n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
80  n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
81  n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
82  n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
83  n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
84  n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
85  n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
86  n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
87  n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
88  n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
89  n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
90  n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
91  n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
92  n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
93  n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
94  n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
95  n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
96  n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
97  n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
98  n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
99  n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
100  n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
101  n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
102  n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
103  n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
104  n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
105  n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
106  n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
107  n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
108  n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
109  n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
110  n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
111  n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
112  n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
113  n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
114  n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
115  n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
116  n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
117  n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
118  n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
119      nlis0   fin sur tmin        pas            fonctions     derivees
120      nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
121      nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
122  -------------------------------------------------------------------------------
123  n1qn3: output mode is  6
124      number of iterations:   42
125      number of simulations:     53
126      realized relative precision on g:  0.21D-15
127      norm of x =  0.14142135D+01
128      f         =  0.53783484D-16
129      norm of x =  0.14924691D-12
130 Optimization stops because too small variations for x.
131
132  n1qn3: entry point
133      dimension of the problem (n):     2
134      absolute precision on x (dxmin): 0.22D-15
135      expected decrease for f (df1): 0.10D+01
136      relative precision on g (epsg): 0.95D-18
137      maximal number of iterations (niter):   100
138      maximal number of simulations (nsim):   100
139      printing level (impres):   3
140      allocated memory (nrz) :     58
141      used memory :                58
142      number of updates :          10
143
144      f         =  0.24200000D+02
145      norm of g =  0.23286769D+03
146  -------------------------------------------------------------------------------
147
148
149  n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
150      nlis0       fpn=-0.200D+01 d2= 0.74D-04  tmin= 0.28D-13 tmax= 0.10D+21
151      nlis0     0.100D+01 -0.194D+01 -0.189D+01
152      nlis0     0.100D+02 -0.147D+02 -0.973D+00
153
154  n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
155      nlis0       fpn=-0.924D+01 d2= 0.66D-02  tmin= 0.30D-14 tmax= 0.10D+21
156      nlis0     0.100D+01 -0.520D+01 -0.138D+01
157
158  n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
159      nlis0       fpn=-0.245D+00 d2= 0.21D-03  tmin= 0.16D-13 tmax= 0.10D+21
160      nlis0     0.100D+01 -0.135D+00 -0.256D-01
161
162  n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
163      nlis0       fpn=-0.581D-02 d2= 0.59D-05  tmin= 0.98D-13 tmax= 0.10D+21
164      nlis0     0.100D+01 -0.435D-02 -0.290D-02
165
166  n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
167      nlis0       fpn=-0.833D-02 d2= 0.26D-04  tmin= 0.61D-13 tmax= 0.10D+21
168      nlis0     0.100D+01 -0.690D-02 -0.549D-02
169
170  n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
171      nlis0       fpn=-0.278D-01 d2= 0.37D-03  tmin= 0.14D-13 tmax= 0.10D+21
172      nlis0     0.100D+01 -0.228D-01 -0.180D-01
173
174  n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
175      nlis0       fpn=-0.828D-01 d2= 0.40D-02  tmin= 0.41D-14 tmax= 0.10D+21
176      nlis0     0.100D+01 -0.713D-01 -0.624D-01
177
178  n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
179      nlis0       fpn=-0.405D+00 d2= 0.11D+00  tmin= 0.79D-15 tmax= 0.10D+21
180      nlis0     0.100D+01 -0.531D+00 -0.721D+00
181      nlis0                                      0.100D+02  0.500D+03  0.222D+03
182      nlis0                                      0.398D+01  0.333D+01  0.774D+01
183      nlis0     0.249D+01 -0.140D+01  0.190D+00
184
185  n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
186      nlis0       fpn=-0.161D+00 d2= 0.57D-01  tmin= 0.11D-14 tmax= 0.10D+21
187      nlis0                                      0.100D+01  0.285D+00  0.571D+00
188      nlis0     0.140D+00 -0.107D-01  0.342D-02
189
190  n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
191      nlis0       fpn=-0.384D-01 d2= 0.21D-03  tmin= 0.21D-13 tmax= 0.10D+21
192      nlis0     0.100D+01 -0.376D-01 -0.369D-01
193      nlis0     0.100D+02 -0.318D+00 -0.257D-01
194
195  n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
196      nlis0       fpn=-0.627D+00 d2= 0.71D-01  tmin= 0.11D-14 tmax= 0.10D+21
197      nlis0     0.100D+01 -0.315D+00  0.212D+00
198
199  n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
200      nlis0       fpn=-0.707D+00 d2= 0.27D-02  tmin= 0.52D-14 tmax= 0.10D+21
201      nlis0     0.100D+01 -0.402D+00 -0.117D+00
202
203  n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
204      nlis0       fpn=-0.499D+00 d2= 0.52D-01  tmin= 0.11D-14 tmax= 0.10D+21
205      nlis0     0.100D+01 -0.225D+00  0.427D+00
206
207  n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
208      nlis0       fpn=-0.339D+00 d2= 0.36D-02  tmin= 0.47D-14 tmax= 0.10D+21
209      nlis0     0.100D+01 -0.151D+00  0.252D-01
210
211  n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
212      nlis0       fpn=-0.207D+00 d2= 0.84D-02  tmin= 0.24D-14 tmax= 0.10D+21
213      nlis0     0.100D+01 -0.186D+00 -0.162D+00
214
215  n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
216      nlis0       fpn=-0.650D+00 d2= 0.12D+00  tmin= 0.64D-15 tmax= 0.10D+21
217      nlis0                                      0.100D+01  0.169D+01  0.738D+01
218      nlis0     0.374D+00 -0.132D+00  0.166D+00
219
220  n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
221      nlis0       fpn=-0.187D+00 d2= 0.20D-02  tmin= 0.62D-14 tmax= 0.10D+21
222      nlis0     0.100D+01 -0.149D+00 -0.115D+00
223
224  n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
225      nlis0       fpn=-0.311D+00 d2= 0.34D-01  tmin= 0.13D-14 tmax= 0.10D+21
226      nlis0     0.100D+01 -0.102D+00  0.268D+00
227
228  n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
229      nlis0       fpn=-0.421D+00 d2= 0.20D-02  tmin= 0.51D-14 tmax= 0.10D+21
230      nlis0     0.100D+01 -0.181D+00  0.576D-01
231
232  n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
233      nlis0       fpn=-0.126D+00 d2= 0.11D-01  tmin= 0.25D-14 tmax= 0.10D+21
234      nlis0     0.100D+01 -0.110D+00 -0.725D-01
235
236  n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
237      nlis0       fpn=-0.228D+00 d2= 0.83D-01  tmin= 0.95D-15 tmax= 0.10D+21
238      nlis0                                      0.100D+01  0.601D+00  0.234D+01
239      nlis0     0.321D+00 -0.373D-01  0.472D-01
240
241  n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
242      nlis0       fpn=-0.554D-01 d2= 0.12D-02  tmin= 0.80D-14 tmax= 0.10D+21
243      nlis0     0.100D+01 -0.477D-01 -0.406D-01
244
245  n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
246      nlis0       fpn=-0.143D+00 d2= 0.30D-01  tmin= 0.17D-14 tmax= 0.10D+21
247      nlis0     0.100D+01 -0.878D-01 -0.226D-01
248
249  n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
250      nlis0       fpn=-0.128D+00 d2= 0.84D-02  tmin= 0.27D-14 tmax= 0.10D+21
251      nlis0     0.100D+01 -0.619D-01 -0.487D-02
252
253  n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
254      nlis0       fpn=-0.100D+00 d2= 0.49D-01  tmin= 0.13D-14 tmax= 0.10D+21
255      nlis0     0.100D+01 -0.155D-01  0.179D+00
256
257  n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
258      nlis0       fpn=-0.120D+00 d2= 0.12D-01  tmin= 0.26D-14 tmax= 0.10D+21
259      nlis0     0.100D+01 -0.251D-01  0.437D-01
260
261  n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
262      nlis0       fpn=-0.208D-01 d2= 0.45D-02  tmin= 0.41D-14 tmax= 0.10D+21
263      nlis0     0.100D+01 -0.183D-01 -0.148D-01
264
265  n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
266      nlis0       fpn=-0.555D-01 d2= 0.58D-01  tmin= 0.11D-14 tmax= 0.10D+21
267      nlis0                                      0.100D+01  0.507D-01  0.260D+00
268      nlis0     0.413D+00 -0.132D-01  0.278D-02
269
270  n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
271      nlis0       fpn=-0.918D-02 d2= 0.87D-03  tmin= 0.86D-14 tmax= 0.10D+21
272      nlis0     0.100D+01 -0.813D-02 -0.718D-02
273
274  n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
275      nlis0       fpn=-0.259D-01 d2= 0.13D-01  tmin= 0.22D-14 tmax= 0.10D+21
276      nlis0     0.100D+01 -0.174D-01 -0.109D-01
277
278  n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
279      nlis0       fpn=-0.949D-02 d2= 0.41D-02  tmin= 0.39D-14 tmax= 0.10D+21
280      nlis0     0.100D+01 -0.648D-02 -0.412D-02
281
282  n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
283      nlis0       fpn=-0.127D-01 d2= 0.37D-01  tmin= 0.13D-14 tmax= 0.10D+21
284      nlis0                                      0.100D+01  0.236D-01  0.889D-01
285      nlis0     0.308D+00 -0.202D-02  0.129D-02
286
287  n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
288      nlis0       fpn=-0.270D-02 d2= 0.13D-02  tmin= 0.68D-14 tmax= 0.10D+21
289      nlis0     0.100D+01 -0.193D-02 -0.124D-02
290
291  n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
292      nlis0       fpn=-0.106D-02 d2= 0.11D-02  tmin= 0.73D-14 tmax= 0.10D+21
293      nlis0     0.100D+01 -0.680D-03 -0.342D-03
294
295  n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
296      nlis0       fpn=-0.165D-03 d2= 0.22D-03  tmin= 0.17D-13 tmax= 0.10D+21
297      nlis0     0.100D+01 -0.994D-04 -0.378D-04
298
299  n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
300      nlis0       fpn=-0.128D-04 d2= 0.34D-04  tmin= 0.43D-13 tmax= 0.10D+21
301      nlis0     0.100D+01 -0.592D-05  0.983D-06
302
303  n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
304      nlis0       fpn=-0.241D-05 d2= 0.96D-07  tmin= 0.92D-12 tmax= 0.10D+21
305      nlis0     0.100D+01 -0.168D-06  0.208D-05
306
307  n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
308      nlis0       fpn=-0.970D-06 d2= 0.12D-08  tmin= 0.80D-11 tmax= 0.10D+21
309      nlis0     0.100D+01 -0.485D-06 -0.183D-09
310
311  n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
312      nlis0       fpn=-0.143D-11 d2= 0.35D-11  tmin= 0.13D-09 tmax= 0.10D+21
313      nlis0     0.100D+01 -0.734D-12 -0.116D-13
314
315  n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
316      nlis0       fpn=-0.991D-16 d2= 0.25D-15  tmin= 0.16D-07 tmax= 0.10D+21
317      nlis0     0.100D+01 -0.152D-15 -0.112D-20
318
319  n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
320      nlis0       fpn=-0.474D-24 d2= 0.19D-25  tmin= 0.17D-02 tmax= 0.10D+21
321      nlis0     0.100D+01 -0.599D-21 -0.107D-26
322
323  n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
324      nlis0       fpn=-0.248D-29 d2= 0.33D-32  tmin= 0.39D+01 tmax= 0.10D+21
325      nlis0                                      0.387D+01  0.326D-23  0.743D-29
326      nlis0   fin sur tmin        pas            fonctions     derivees
327      nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
328      nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
329  -------------------------------------------------------------------------------
330  n1qn3: output mode is  6
331      number of iterations:   42
332      number of simulations:     53
333      realized relative precision on g:  0.21D-15
334      norm of x =  0.14142135D+01
335      f         =  0.53783484D-16
336      norm of x =  0.14924691D-12
337 Optimization stops because too small variations for x.
338
339  n1qn3: entry point
340      dimension of the problem (n):     2
341      absolute precision on x (dxmin): 0.22D-15
342      expected decrease for f (df1): 0.10D+01
343      relative precision on g (epsg): 0.95D-18
344      maximal number of iterations (niter):   100
345      maximal number of simulations (nsim):   100
346      printing level (impres):   4
347      allocated memory (nrz) :     58
348      used memory :                58
349      number of updates :          10
350
351      f         =  0.24200000D+02
352      norm of g =  0.23286769D+03
353  n1qn3a: descent direction -g: precon =  0.369D-04
354  -------------------------------------------------------------------------------
355  -------------------------------------------------------------------------------
356
357  n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
358  n1qn3: line search
359      nlis0       fpn=-0.200D+01 d2= 0.74D-04  tmin= 0.28D-13 tmax= 0.10D+21
360      nlis0     0.100D+01 -0.194D+01 -0.189D+01
361      nlis0     0.100D+02 -0.147D+02 -0.973D+00
362  n1qn3: stopping criterion on g:  0.48674D+00
363  n1qn3: matrix update:
364      Oren-Spedicato factor (not used) =  0.718D-03     diagonal: average value =  0.718D-03
365  n1qn3: descent direction d: angle(-g,d) =   0.9 degrees
366  -------------------------------------------------------------------------------
367
368  n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
369  n1qn3: line search
370      nlis0       fpn=-0.924D+01 d2= 0.66D-02  tmin= 0.30D-14 tmax= 0.10D+21
371      nlis0     0.100D+01 -0.520D+01 -0.138D+01
372  n1qn3: stopping criterion on g:  0.73079D-01
373  n1qn3: convergence rate, s(k)/s(k-1) =  0.94908D+00
374  n1qn3: matrix update:
375      fitting the ellipsoid: factor  0.118D+01
376      Oren-Spedicato factor (not used) =  0.846D-03     diagonal: average value =  0.846D-03
377  n1qn3: descent direction d: angle(-g,d) =   0.3 degrees
378  -------------------------------------------------------------------------------
379
380  n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
381  n1qn3: line search
382      nlis0       fpn=-0.245D+00 d2= 0.21D-03  tmin= 0.16D-13 tmax= 0.10D+21
383      nlis0     0.100D+01 -0.135D+00 -0.256D-01
384  n1qn3: stopping criterion on g:  0.10270D-01
385  n1qn3: convergence rate, s(k)/s(k-1) =  0.17647D+00
386  n1qn3: matrix update:
387      fitting the ellipsoid: factor  0.111D+01
388      Oren-Spedicato factor (not used) =  0.935D-03     diagonal: average value =  0.942D-03
389  n1qn3: descent direction d: angle(-g,d) =   1.1 degrees
390  -------------------------------------------------------------------------------
391
392  n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
393  n1qn3: line search
394      nlis0       fpn=-0.581D-02 d2= 0.59D-05  tmin= 0.98D-13 tmax= 0.10D+21
395      nlis0     0.100D+01 -0.435D-02 -0.290D-02
396  n1qn3: stopping criterion on g:  0.76379D-02
397  n1qn3: convergence rate, s(k)/s(k-1) =  0.16885D+00
398  n1qn3: matrix update:
399      fitting the ellipsoid: factor  0.969D+00
400      Oren-Spedicato factor (not used) =  0.952D-03     diagonal: average value =  0.915D-03
401  n1qn3: descent direction d: angle(-g,d) =  23.6 degrees
402  -------------------------------------------------------------------------------
403
404  n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
405  n1qn3: line search
406      nlis0       fpn=-0.833D-02 d2= 0.26D-04  tmin= 0.61D-13 tmax= 0.10D+21
407      nlis0     0.100D+01 -0.690D-02 -0.549D-02
408  n1qn3: stopping criterion on g:  0.11072D-01
409  n1qn3: convergence rate, s(k)/s(k-1) =  0.21032D+01
410  n1qn3: matrix update:
411      fitting the ellipsoid: factor  0.985D+00
412      Oren-Spedicato factor (not used) =  0.954D-03     diagonal: average value =  0.102D-02
413  n1qn3: descent direction d: angle(-g,d) =  56.1 degrees
414  -------------------------------------------------------------------------------
415
416  n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
417  n1qn3: line search
418      nlis0       fpn=-0.278D-01 d2= 0.37D-03  tmin= 0.14D-13 tmax= 0.10D+21
419      nlis0     0.100D+01 -0.228D-01 -0.180D-01
420  n1qn3: stopping criterion on g:  0.23226D-01
421  n1qn3: convergence rate, s(k)/s(k-1) =  0.37768D+01
422  n1qn3: matrix update:
423      fitting the ellipsoid: factor  0.113D+01
424      Oren-Spedicato factor (not used) =  0.941D-03     diagonal: average value =  0.127D-02
425  n1qn3: descent direction d: angle(-g,d) =  76.0 degrees
426  -------------------------------------------------------------------------------
427
428  n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
429  n1qn3: line search
430      nlis0       fpn=-0.828D-01 d2= 0.40D-02  tmin= 0.41D-14 tmax= 0.10D+21
431      nlis0     0.100D+01 -0.713D-01 -0.624D-01
432  n1qn3: stopping criterion on g:  0.44894D-01
433  n1qn3: convergence rate, s(k)/s(k-1) =  0.32918D+01
434  n1qn3: matrix update:
435      fitting the ellipsoid: factor  0.817D+00
436      Oren-Spedicato factor (not used) =  0.762D-03     diagonal: average value =  0.106D-02
437  n1qn3: descent direction d: angle(-g,d) =  83.2 degrees
438  -------------------------------------------------------------------------------
439
440  n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
441  n1qn3: line search
442      nlis0       fpn=-0.405D+00 d2= 0.11D+00  tmin= 0.79D-15 tmax= 0.10D+21
443      nlis0     0.100D+01 -0.531D+00 -0.721D+00
444      nlis0                                      0.100D+02  0.500D+03  0.222D+03
445      nlis0                                      0.398D+01  0.333D+01  0.774D+01
446      nlis0     0.249D+01 -0.140D+01  0.190D+00
447  n1qn3: stopping criterion on g:  0.60049D-01
448  n1qn3: convergence rate, s(k)/s(k-1) =  0.12764D+02
449  n1qn3: matrix update:
450      fitting the ellipsoid: factor  0.287D+01
451      Oren-Spedicato factor (not used) =  0.248D-02     diagonal: average value =  0.292D-02
452  n1qn3: descent direction d: angle(-g,d) =  87.2 degrees
453  -------------------------------------------------------------------------------
454
455  n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
456  n1qn3: line search
457      nlis0       fpn=-0.161D+00 d2= 0.57D-01  tmin= 0.11D-14 tmax= 0.10D+21
458      nlis0                                      0.100D+01  0.285D+00  0.571D+00
459      nlis0     0.140D+00 -0.107D-01  0.342D-02
460  n1qn3: stopping criterion on g:  0.48319D-01
461  n1qn3: convergence rate, s(k)/s(k-1) =  0.41360D-01
462  n1qn3: matrix update:
463      fitting the ellipsoid: factor  0.102D+01
464      Oren-Spedicato factor (not used) =  0.293D-02     diagonal: average value =  0.285D-02
465  n1qn3: descent direction d: angle(-g,d) =  76.4 degrees
466  -------------------------------------------------------------------------------
467
468  n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
469  n1qn3: line search
470      nlis0       fpn=-0.384D-01 d2= 0.21D-03  tmin= 0.21D-13 tmax= 0.10D+21
471      nlis0     0.100D+01 -0.376D-01 -0.369D-01
472      nlis0     0.100D+02 -0.318D+00 -0.257D-01
473  n1qn3: stopping criterion on g:  0.41045D-01
474  n1qn3: convergence rate, s(k)/s(k-1) =  0.43358D+01
475  n1qn3: matrix update:
476      fitting the ellipsoid: factor  0.193D+02
477      Oren-Spedicato factor (not used) =  0.403D-01     diagonal: average value =  0.645D-01
478  n1qn3: descent direction d: angle(-g,d) =  75.7 degrees
479  -------------------------------------------------------------------------------
480
481  n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
482  n1qn3: line search
483      nlis0       fpn=-0.627D+00 d2= 0.71D-01  tmin= 0.11D-14 tmax= 0.10D+21
484      nlis0     0.100D+01 -0.315D+00  0.212D+00
485  n1qn3: stopping criterion on g:  0.60708D-01
486  n1qn3: convergence rate, s(k)/s(k-1) =  0.18312D+01
487  n1qn3: matrix update:
488      fitting the ellipsoid: factor  0.259D+00
489      Oren-Spedicato factor (not used) =  0.251D-01     diagonal: average value =  0.248D-01
490  n1qn3: descent direction d: angle(-g,d) =  15.3 degrees
491  -------------------------------------------------------------------------------
492
493  n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
494  n1qn3: line search
495      nlis0       fpn=-0.707D+00 d2= 0.27D-02  tmin= 0.52D-14 tmax= 0.10D+21
496      nlis0     0.100D+01 -0.402D+00 -0.117D+00
497  n1qn3: stopping criterion on g:  0.11240D-01
498  n1qn3: convergence rate, s(k)/s(k-1) =  0.19476D+00
499  n1qn3: matrix update:
500      fitting the ellipsoid: factor  0.190D+00
501      Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.468D-02
502  n1qn3: descent direction d: angle(-g,d) =  33.1 degrees
503  -------------------------------------------------------------------------------
504
505  n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
506  n1qn3: line search
507      nlis0       fpn=-0.499D+00 d2= 0.52D-01  tmin= 0.11D-14 tmax= 0.10D+21
508      nlis0     0.100D+01 -0.225D+00  0.427D+00
509  n1qn3: stopping criterion on g:  0.43961D-01
510  n1qn3: convergence rate, s(k)/s(k-1) =  0.43931D+01
511  n1qn3: matrix update:
512      fitting the ellipsoid: factor  0.364D+01
513      Oren-Spedicato factor (not used) =  0.102D-01     diagonal: average value =  0.311D-01
514  n1qn3: descent direction d: angle(-g,d) =  56.6 degrees
515  -------------------------------------------------------------------------------
516
517  n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
518  n1qn3: line search
519      nlis0       fpn=-0.339D+00 d2= 0.36D-02  tmin= 0.47D-14 tmax= 0.10D+21
520      nlis0     0.100D+01 -0.151D+00  0.252D-01
521  n1qn3: stopping criterion on g:  0.12127D-01
522  n1qn3: convergence rate, s(k)/s(k-1) =  0.26394D+00
523  n1qn3: matrix update:
524      fitting the ellipsoid: factor  0.759D+00
525      Oren-Spedicato factor (not used) =  0.522D-02     diagonal: average value =  0.235D-01
526  n1qn3: descent direction d: angle(-g,d) =  36.9 degrees
527  -------------------------------------------------------------------------------
528
529  n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
530  n1qn3: line search
531      nlis0       fpn=-0.207D+00 d2= 0.84D-02  tmin= 0.24D-14 tmax= 0.10D+21
532      nlis0     0.100D+01 -0.186D+00 -0.162D+00
533  n1qn3: stopping criterion on g:  0.11958D-01
534  n1qn3: convergence rate, s(k)/s(k-1) =  0.15231D+01
535  n1qn3: matrix update:
536      fitting the ellipsoid: factor  0.493D+01
537      Oren-Spedicato factor (not used) =  0.996D-01     diagonal: average value =  0.118D+00
538  n1qn3: descent direction d: angle(-g,d) =  47.6 degrees
539  -------------------------------------------------------------------------------
540
541  n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
542  n1qn3: line search
543      nlis0       fpn=-0.650D+00 d2= 0.12D+00  tmin= 0.64D-15 tmax= 0.10D+21
544      nlis0                                      0.100D+01  0.169D+01  0.738D+01
545      nlis0     0.374D+00 -0.132D+00  0.166D+00
546  n1qn3: stopping criterion on g:  0.29691D-01
547  n1qn3: convergence rate, s(k)/s(k-1) =  0.14137D+01
548  n1qn3: matrix update:
549      fitting the ellipsoid: factor  0.211D+00
550      Oren-Spedicato factor (not used) =  0.105D-01     diagonal: average value =  0.372D-01
551  n1qn3: descent direction d: angle(-g,d) =  53.1 degrees
552  -------------------------------------------------------------------------------
553
554  n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
555  n1qn3: line search
556      nlis0       fpn=-0.187D+00 d2= 0.20D-02  tmin= 0.62D-14 tmax= 0.10D+21
557      nlis0     0.100D+01 -0.149D+00 -0.115D+00
558  n1qn3: stopping criterion on g:  0.15621D-01
559  n1qn3: convergence rate, s(k)/s(k-1) =  0.34738D+00
560  n1qn3: matrix update:
561      fitting the ellipsoid: factor  0.123D+01
562      Oren-Spedicato factor (not used) =  0.634D-02     diagonal: average value =  0.386D-01
563  n1qn3: descent direction d: angle(-g,d) =  62.6 degrees
564  -------------------------------------------------------------------------------
565
566  n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
567  n1qn3: line search
568      nlis0       fpn=-0.311D+00 d2= 0.34D-01  tmin= 0.13D-14 tmax= 0.10D+21
569      nlis0     0.100D+01 -0.102D+00  0.268D+00
570  n1qn3: stopping criterion on g:  0.42709D-01
571  n1qn3: convergence rate, s(k)/s(k-1) =  0.41295D+01
572  n1qn3: matrix update:
573      fitting the ellipsoid: factor  0.285D+00
574      Oren-Spedicato factor (not used) =  0.108D-01     diagonal: average value =  0.841D-02
575  n1qn3: descent direction d: angle(-g,d) =  17.2 degrees
576  -------------------------------------------------------------------------------
577
578  n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
579  n1qn3: line search
580      nlis0       fpn=-0.421D+00 d2= 0.20D-02  tmin= 0.51D-14 tmax= 0.10D+21
581      nlis0     0.100D+01 -0.181D+00  0.576D-01
582  n1qn3: stopping criterion on g:  0.94247D-02
583  n1qn3: convergence rate, s(k)/s(k-1) =  0.23849D+00
584  n1qn3: matrix update:
585      fitting the ellipsoid: factor  0.526D+00
586      Oren-Spedicato factor (not used) =  0.345D-02     diagonal: average value =  0.461D-02
587  n1qn3: descent direction d: angle(-g,d) =  57.2 degrees
588  -------------------------------------------------------------------------------
589
590  n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
591  n1qn3: line search
592      nlis0       fpn=-0.126D+00 d2= 0.11D-01  tmin= 0.25D-14 tmax= 0.10D+21
593      nlis0     0.100D+01 -0.110D+00 -0.725D-01
594  n1qn3: stopping criterion on g:  0.67714D-02
595  n1qn3: convergence rate, s(k)/s(k-1) =  0.23948D+01
596  n1qn3: matrix update:
597      fitting the ellipsoid: factor  0.105D+01
598      Oren-Spedicato factor (not used) =  0.489D-02     diagonal: average value =  0.547D-02
599  n1qn3: descent direction d: angle(-g,d) =  59.9 degrees
600  -------------------------------------------------------------------------------
601
602  n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
603  n1qn3: line search
604      nlis0       fpn=-0.228D+00 d2= 0.83D-01  tmin= 0.95D-15 tmax= 0.10D+21
605      nlis0                                      0.100D+01  0.601D+00  0.234D+01
606      nlis0     0.321D+00 -0.373D-01  0.472D-01
607  n1qn3: stopping criterion on g:  0.24287D-01
608  n1qn3: convergence rate, s(k)/s(k-1) =  0.87291D+00
609  n1qn3: matrix update:
610      fitting the ellipsoid: factor  0.807D+00
611      Oren-Spedicato factor (not used) =  0.458D-02     diagonal: average value =  0.431D-02
612  n1qn3: descent direction d: angle(-g,d) =  73.3 degrees
613  -------------------------------------------------------------------------------
614
615  n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
616  n1qn3: line search
617      nlis0       fpn=-0.554D-01 d2= 0.12D-02  tmin= 0.80D-14 tmax= 0.10D+21
618      nlis0     0.100D+01 -0.477D-01 -0.406D-01
619  n1qn3: stopping criterion on g:  0.15832D-01
620  n1qn3: convergence rate, s(k)/s(k-1) =  0.36912D+00
621  n1qn3: matrix update:
622      fitting the ellipsoid: factor  0.894D+00
623      Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.392D-02
624  n1qn3: descent direction d: angle(-g,d) =  77.1 degrees
625  -------------------------------------------------------------------------------
626
627  n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
628  n1qn3: line search
629      nlis0       fpn=-0.143D+00 d2= 0.30D-01  tmin= 0.17D-14 tmax= 0.10D+21
630      nlis0     0.100D+01 -0.878D-01 -0.226D-01
631  n1qn3: stopping criterion on g:  0.23266D-01
632  n1qn3: convergence rate, s(k)/s(k-1) =  0.51063D+01
633  n1qn3: matrix update:
634      fitting the ellipsoid: factor  0.865D+01
635      Oren-Spedicato factor (not used) =  0.318D-01     diagonal: average value =  0.403D-01
636  n1qn3: descent direction d: angle(-g,d) =  75.1 degrees
637  -------------------------------------------------------------------------------
638
639  n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
640  n1qn3: line search
641      nlis0       fpn=-0.128D+00 d2= 0.84D-02  tmin= 0.27D-14 tmax= 0.10D+21
642      nlis0     0.100D+01 -0.619D-01 -0.487D-02
643  n1qn3: stopping criterion on g:  0.10083D-01
644  n1qn3: convergence rate, s(k)/s(k-1) =  0.52582D+00
645  n1qn3: matrix update:
646      fitting the ellipsoid: factor  0.572D-01
647      Oren-Spedicato factor (not used) =  0.208D-02     diagonal: average value =  0.227D-02
648  n1qn3: descent direction d: angle(-g,d) =  78.8 degrees
649  -------------------------------------------------------------------------------
650
651  n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
652  n1qn3: line search
653      nlis0       fpn=-0.100D+00 d2= 0.49D-01  tmin= 0.13D-14 tmax= 0.10D+21
654      nlis0     0.100D+01 -0.155D-01  0.179D+00
655  n1qn3: stopping criterion on g:  0.37918D-01
656  n1qn3: convergence rate, s(k)/s(k-1) =  0.24062D+01
657  n1qn3: matrix update:
658      fitting the ellipsoid: factor  0.120D+01
659      Oren-Spedicato factor (not used) =  0.224D-02     diagonal: average value =  0.267D-02
660  n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
661  -------------------------------------------------------------------------------
662
663  n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
664  n1qn3: line search
665      nlis0       fpn=-0.120D+00 d2= 0.12D-01  tmin= 0.26D-14 tmax= 0.10D+21
666      nlis0     0.100D+01 -0.251D-01  0.437D-01
667  n1qn3: stopping criterion on g:  0.43947D-02
668  n1qn3: convergence rate, s(k)/s(k-1) =  0.50600D+00
669  n1qn3: matrix update:
670      fitting the ellipsoid: factor  0.757D+00
671      Oren-Spedicato factor (not used) =  0.170D-02     diagonal: average value =  0.201D-02
672  n1qn3: descent direction d: angle(-g,d) =  72.3 degrees
673  -------------------------------------------------------------------------------
674
675  n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
676  n1qn3: line search
677      nlis0       fpn=-0.208D-01 d2= 0.45D-02  tmin= 0.41D-14 tmax= 0.10D+21
678      nlis0     0.100D+01 -0.183D-01 -0.148D-01
679  n1qn3: stopping criterion on g:  0.21889D-02
680  n1qn3: convergence rate, s(k)/s(k-1) =  0.59772D+00
681  n1qn3: matrix update:
682      fitting the ellipsoid: factor  0.170D+01
683      Oren-Spedicato factor (not used) =  0.290D-02     diagonal: average value =  0.351D-02
684  n1qn3: descent direction d: angle(-g,d) =  63.1 degrees
685  -------------------------------------------------------------------------------
686
687  n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
688  n1qn3: line search
689      nlis0       fpn=-0.555D-01 d2= 0.58D-01  tmin= 0.11D-14 tmax= 0.10D+21
690      nlis0                                      0.100D+01  0.507D-01  0.260D+00
691      nlis0     0.413D+00 -0.132D-01  0.278D-02
692  n1qn3: stopping criterion on g:  0.14651D-01
693  n1qn3: convergence rate, s(k)/s(k-1) =  0.14924D+01
694  n1qn3: matrix update:
695      fitting the ellipsoid: factor  0.980D+00
696      Oren-Spedicato factor (not used) =  0.274D-02     diagonal: average value =  0.355D-02
697  n1qn3: descent direction d: angle(-g,d) =  84.8 degrees
698  -------------------------------------------------------------------------------
699
700  n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
701  n1qn3: line search
702      nlis0       fpn=-0.918D-02 d2= 0.87D-03  tmin= 0.86D-14 tmax= 0.10D+21
703      nlis0     0.100D+01 -0.813D-02 -0.718D-02
704  n1qn3: stopping criterion on g:  0.11306D-01
705  n1qn3: convergence rate, s(k)/s(k-1) =  0.29590D+00
706  n1qn3: matrix update:
707      fitting the ellipsoid: factor  0.112D+01
708      Oren-Spedicato factor (not used) =  0.329D-02     diagonal: average value =  0.403D-02
709  n1qn3: descent direction d: angle(-g,d) =  85.0 degrees
710  -------------------------------------------------------------------------------
711
712  n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
713  n1qn3: line search
714      nlis0       fpn=-0.259D-01 d2= 0.13D-01  tmin= 0.22D-14 tmax= 0.10D+21
715      nlis0     0.100D+01 -0.174D-01 -0.109D-01
716  n1qn3: stopping criterion on g:  0.64294D-02
717  n1qn3: convergence rate, s(k)/s(k-1) =  0.38281D+01
718  n1qn3: matrix update:
719      fitting the ellipsoid: factor  0.337D+01
720      Oren-Spedicato factor (not used) =  0.116D-01     diagonal: average value =  0.134D-01
721  n1qn3: descent direction d: angle(-g,d) =  84.3 degrees
722  -------------------------------------------------------------------------------
723
724  n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
725  n1qn3: line search
726      nlis0       fpn=-0.949D-02 d2= 0.41D-02  tmin= 0.39D-14 tmax= 0.10D+21
727      nlis0     0.100D+01 -0.648D-02 -0.412D-02
728  n1qn3: stopping criterion on g:  0.32255D-03
729  n1qn3: convergence rate, s(k)/s(k-1) =  0.56905D+00
730  n1qn3: matrix update:
731      fitting the ellipsoid: factor  0.218D+00
732      Oren-Spedicato factor (not used) =  0.229D-02     diagonal: average value =  0.316D-02
733  n1qn3: descent direction d: angle(-g,d) =  28.9 degrees
734  -------------------------------------------------------------------------------
735
736  n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
737  n1qn3: line search
738      nlis0       fpn=-0.127D-01 d2= 0.37D-01  tmin= 0.13D-14 tmax= 0.10D+21
739      nlis0                                      0.100D+01  0.236D-01  0.889D-01
740      nlis0     0.308D+00 -0.202D-02  0.129D-02
741  n1qn3: stopping criterion on g:  0.59635D-02
742  n1qn3: convergence rate, s(k)/s(k-1) =  0.92424D+00
743  n1qn3: matrix update:
744      fitting the ellipsoid: factor  0.995D+00
745      Oren-Spedicato factor (not used) =  0.211D-02     diagonal: average value =  0.320D-02
746  n1qn3: descent direction d: angle(-g,d) =  86.9 degrees
747  -------------------------------------------------------------------------------
748
749  n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
750  n1qn3: line search
751      nlis0       fpn=-0.270D-02 d2= 0.13D-02  tmin= 0.68D-14 tmax= 0.10D+21
752      nlis0     0.100D+01 -0.193D-02 -0.124D-02
753  n1qn3: stopping criterion on g:  0.26088D-02
754  n1qn3: convergence rate, s(k)/s(k-1) =  0.61219D+00
755  n1qn3: matrix update:
756      fitting the ellipsoid: factor  0.107D+01
757      Oren-Spedicato factor (not used) =  0.240D-02     diagonal: average value =  0.346D-02
758  n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
759  -------------------------------------------------------------------------------
760
761  n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
762  n1qn3: line search
763      nlis0       fpn=-0.106D-02 d2= 0.11D-02  tmin= 0.73D-14 tmax= 0.10D+21
764      nlis0     0.100D+01 -0.680D-03 -0.342D-03
765  n1qn3: stopping criterion on g:  0.93209D-03
766  n1qn3: convergence rate, s(k)/s(k-1) =  0.93175D+00
767  n1qn3: matrix update:
768      fitting the ellipsoid: factor  0.195D+01
769      Oren-Spedicato factor (not used) =  0.472D-02     diagonal: average value =  0.669D-02
770  n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
771  -------------------------------------------------------------------------------
772
773  n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
774  n1qn3: line search
775      nlis0       fpn=-0.165D-03 d2= 0.22D-03  tmin= 0.17D-13 tmax= 0.10D+21
776      nlis0     0.100D+01 -0.994D-04 -0.378D-04
777  n1qn3: stopping criterion on g:  0.13605D-03
778  n1qn3: convergence rate, s(k)/s(k-1) =  0.43415D+00
779  n1qn3: matrix update:
780      fitting the ellipsoid: factor  0.797D+00
781      Oren-Spedicato factor (not used) =  0.369D-02     diagonal: average value =  0.540D-02
782  n1qn3: descent direction d: angle(-g,d) =  86.0 degrees
783  -------------------------------------------------------------------------------
784
785  n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
786  n1qn3: line search
787      nlis0       fpn=-0.128D-04 d2= 0.34D-04  tmin= 0.43D-13 tmax= 0.10D+21
788      nlis0     0.100D+01 -0.592D-05  0.983D-06
789  n1qn3: stopping criterion on g:  0.15294D-03
790  n1qn3: convergence rate, s(k)/s(k-1) =  0.39644D+00
791  n1qn3: matrix update:
792      fitting the ellipsoid: factor  0.278D+03
793      Oren-Spedicato factor (not used) =  0.637D+00     diagonal: average value =  0.256D+01
794  n1qn3: descent direction d: angle(-g,d) =  77.4 degrees
795  -------------------------------------------------------------------------------
796
797  n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
798  n1qn3: line search
799      nlis0       fpn=-0.241D-05 d2= 0.96D-07  tmin= 0.92D-12 tmax= 0.10D+21
800      nlis0     0.100D+01 -0.168D-06  0.208D-05
801  n1qn3: stopping criterion on g:  0.13388D-03
802  n1qn3: convergence rate, s(k)/s(k-1) =  0.53137D-01
803  n1qn3: matrix update:
804      fitting the ellipsoid: factor  0.794D-03
805      Oren-Spedicato factor (not used) =  0.101D-02     diagonal: average value =  0.159D-02
806  n1qn3: descent direction d: angle(-g,d) =  26.4 degrees
807  -------------------------------------------------------------------------------
808
809  n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
810  n1qn3: line search
811      nlis0       fpn=-0.970D-06 d2= 0.12D-08  tmin= 0.80D-11 tmax= 0.10D+21
812      nlis0     0.100D+01 -0.485D-06 -0.183D-09
813  n1qn3: stopping criterion on g:  0.27051D-07
814  n1qn3: convergence rate, s(k)/s(k-1) =  0.11218D+00
815  n1qn3: matrix update:
816      fitting the ellipsoid: factor  0.866D+00
817      Oren-Spedicato factor (not used) =  0.999D-03     diagonal: average value =  0.137D-02
818  n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
819  -------------------------------------------------------------------------------
820
821  n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
822  n1qn3: line search
823      nlis0       fpn=-0.143D-11 d2= 0.35D-11  tmin= 0.13D-09 tmax= 0.10D+21
824      nlis0     0.100D+01 -0.734D-12 -0.116D-13
825  n1qn3: stopping criterion on g:  0.61631D-10
826  n1qn3: convergence rate, s(k)/s(k-1) =  0.53526D-01
827  n1qn3: matrix update:
828      fitting the ellipsoid: factor  0.399D+02
829      Oren-Spedicato factor (not used) =  0.356D-01     diagonal: average value =  0.717D-01
830  n1qn3: descent direction d: angle(-g,d) =  64.0 degrees
831  -------------------------------------------------------------------------------
832
833  n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
834  n1qn3: line search
835      nlis0       fpn=-0.991D-16 d2= 0.25D-15  tmin= 0.16D-07 tmax= 0.10D+21
836      nlis0     0.100D+01 -0.152D-15 -0.112D-20
837  n1qn3: stopping criterion on g:  0.92861D-13
838  n1qn3: convergence rate, s(k)/s(k-1) =  0.84647D-02
839  n1qn3: matrix update:
840      fitting the ellipsoid: factor  0.571D+01
841      Oren-Spedicato factor (not used) =  0.482D+00     diagonal: average value =  0.356D+00
842  n1qn3: descent direction d: angle(-g,d) =  80.8 degrees
843  -------------------------------------------------------------------------------
844
845  n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
846  n1qn3: line search
847      nlis0       fpn=-0.474D-24 d2= 0.19D-25  tmin= 0.17D-02 tmax= 0.10D+21
848      nlis0     0.100D+01 -0.599D-21 -0.107D-26
849  n1qn3: stopping criterion on g:  0.21365D-15
850  n1qn3: convergence rate, s(k)/s(k-1) =  0.86648D-05
851  n1qn3: matrix update:
852      fitting the ellipsoid: factor  0.352D-02
853      Oren-Spedicato factor (not used) =  0.102D-02     diagonal: average value =  0.198D-02
854  n1qn3: descent direction d: angle(-g,d) =  29.9 degrees
855  -------------------------------------------------------------------------------
856
857  n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
858  n1qn3: line search
859      nlis0       fpn=-0.248D-29 d2= 0.33D-32  tmin= 0.39D+01 tmax= 0.10D+21
860      nlis0                                      0.387D+01  0.326D-23  0.743D-29
861      nlis0   fin sur tmin        pas            fonctions     derivees
862      nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
863      nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
864  -------------------------------------------------------------------------------
865  n1qn3: output mode is  6
866      number of iterations:   42
867      number of simulations:     53
868      realized relative precision on g:  0.21D-15
869      norm of x =  0.14142135D+01
870      f         =  0.53783484D-16
871      norm of x =  0.14924691D-12
872 Optimization stops because too small variations for x.
873
874  n1qn3: entry point
875      dimension of the problem (n):     2
876      absolute precision on x (dxmin): 0.22D-15
877      expected decrease for f (df1): 0.10D+01
878      relative precision on g (epsg): 0.95D-18
879      maximal number of iterations (niter):   100
880      maximal number of simulations (nsim):   100
881      printing level (impres):   5
882      allocated memory (nrz) :     58
883      used memory :                58
884      number of updates :          10
885
886      f         =  0.24200000D+02
887      norm of g =  0.23286769D+03
888  n1qn3a: descent direction -g: precon =  0.369D-04
889  -------------------------------------------------------------------------------
890
891  n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
892  n1qn3: line search
893      nlis0       fpn=-0.200D+01 d2= 0.74D-04  tmin= 0.28D-13 tmax= 0.10D+21
894      nlis0     0.100D+01 -0.194D+01 -0.189D+01
895      nlis0     0.100D+02 -0.147D+02 -0.973D+00
896  n1qn3: stopping criterion on g:  0.48674D+00
897  n1qn3: matrix update:
898      Oren-Spedicato factor (not used) =  0.718D-03     diagonal: average value =  0.718D-03
899  n1qn3: descent direction d: angle(-g,d) =   0.9 degrees
900  -------------------------------------------------------------------------------
901
902  n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
903  n1qn3: line search
904      nlis0       fpn=-0.924D+01 d2= 0.66D-02  tmin= 0.30D-14 tmax= 0.10D+21
905      nlis0     0.100D+01 -0.520D+01 -0.138D+01
906  n1qn3: stopping criterion on g:  0.73079D-01
907  n1qn3: convergence rate, s(k)/s(k-1) =  0.94908D+00
908  n1qn3: matrix update:
909      fitting the ellipsoid: factor  0.118D+01
910      Oren-Spedicato factor (not used) =  0.846D-03     diagonal: average value =  0.846D-03
911  n1qn3: descent direction d: angle(-g,d) =   0.3 degrees
912  -------------------------------------------------------------------------------
913
914  n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
915  n1qn3: line search
916      nlis0       fpn=-0.245D+00 d2= 0.21D-03  tmin= 0.16D-13 tmax= 0.10D+21
917      nlis0     0.100D+01 -0.135D+00 -0.256D-01
918  n1qn3: stopping criterion on g:  0.10270D-01
919  n1qn3: convergence rate, s(k)/s(k-1) =  0.17647D+00
920  n1qn3: matrix update:
921      fitting the ellipsoid: factor  0.111D+01
922      Oren-Spedicato factor (not used) =  0.935D-03     diagonal: average value =  0.942D-03
923  n1qn3: descent direction d: angle(-g,d) =   1.1 degrees
924  -------------------------------------------------------------------------------
925
926  n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
927  n1qn3: line search
928      nlis0       fpn=-0.581D-02 d2= 0.59D-05  tmin= 0.98D-13 tmax= 0.10D+21
929      nlis0     0.100D+01 -0.435D-02 -0.290D-02
930  n1qn3: stopping criterion on g:  0.76379D-02
931  n1qn3: convergence rate, s(k)/s(k-1) =  0.16885D+00
932  n1qn3: matrix update:
933      fitting the ellipsoid: factor  0.969D+00
934      Oren-Spedicato factor (not used) =  0.952D-03     diagonal: average value =  0.915D-03
935  n1qn3: descent direction d: angle(-g,d) =  23.6 degrees
936  -------------------------------------------------------------------------------
937
938  n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
939  n1qn3: line search
940      nlis0       fpn=-0.833D-02 d2= 0.26D-04  tmin= 0.61D-13 tmax= 0.10D+21
941      nlis0     0.100D+01 -0.690D-02 -0.549D-02
942  n1qn3: stopping criterion on g:  0.11072D-01
943  n1qn3: convergence rate, s(k)/s(k-1) =  0.21032D+01
944  n1qn3: matrix update:
945      fitting the ellipsoid: factor  0.985D+00
946      Oren-Spedicato factor (not used) =  0.954D-03     diagonal: average value =  0.102D-02
947  n1qn3: descent direction d: angle(-g,d) =  56.1 degrees
948  -------------------------------------------------------------------------------
949
950  n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
951  n1qn3: line search
952      nlis0       fpn=-0.278D-01 d2= 0.37D-03  tmin= 0.14D-13 tmax= 0.10D+21
953      nlis0     0.100D+01 -0.228D-01 -0.180D-01
954  n1qn3: stopping criterion on g:  0.23226D-01
955  n1qn3: convergence rate, s(k)/s(k-1) =  0.37768D+01
956  n1qn3: matrix update:
957      fitting the ellipsoid: factor  0.113D+01
958      Oren-Spedicato factor (not used) =  0.941D-03     diagonal: average value =  0.127D-02
959  n1qn3: descent direction d: angle(-g,d) =  76.0 degrees
960  -------------------------------------------------------------------------------
961
962  n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
963  n1qn3: line search
964      nlis0       fpn=-0.828D-01 d2= 0.40D-02  tmin= 0.41D-14 tmax= 0.10D+21
965      nlis0     0.100D+01 -0.713D-01 -0.624D-01
966  n1qn3: stopping criterion on g:  0.44894D-01
967  n1qn3: convergence rate, s(k)/s(k-1) =  0.32918D+01
968  n1qn3: matrix update:
969      fitting the ellipsoid: factor  0.817D+00
970      Oren-Spedicato factor (not used) =  0.762D-03     diagonal: average value =  0.106D-02
971  n1qn3: descent direction d: angle(-g,d) =  83.2 degrees
972  -------------------------------------------------------------------------------
973
974  n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
975  n1qn3: line search
976      nlis0       fpn=-0.405D+00 d2= 0.11D+00  tmin= 0.79D-15 tmax= 0.10D+21
977      nlis0     0.100D+01 -0.531D+00 -0.721D+00
978      nlis0                                      0.100D+02  0.500D+03  0.222D+03
979      nlis0                                      0.398D+01  0.333D+01  0.774D+01
980      nlis0     0.249D+01 -0.140D+01  0.190D+00
981  n1qn3: stopping criterion on g:  0.60049D-01
982  n1qn3: convergence rate, s(k)/s(k-1) =  0.12764D+02
983  n1qn3: matrix update:
984      fitting the ellipsoid: factor  0.287D+01
985      Oren-Spedicato factor (not used) =  0.248D-02     diagonal: average value =  0.292D-02
986  n1qn3: descent direction d: angle(-g,d) =  87.2 degrees
987  -------------------------------------------------------------------------------
988
989  n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
990  n1qn3: line search
991      nlis0       fpn=-0.161D+00 d2= 0.57D-01  tmin= 0.11D-14 tmax= 0.10D+21
992      nlis0                                      0.100D+01  0.285D+00  0.571D+00
993      nlis0     0.140D+00 -0.107D-01  0.342D-02
994  n1qn3: stopping criterion on g:  0.48319D-01
995  n1qn3: convergence rate, s(k)/s(k-1) =  0.41360D-01
996  n1qn3: matrix update:
997      fitting the ellipsoid: factor  0.102D+01
998      Oren-Spedicato factor (not used) =  0.293D-02     diagonal: average value =  0.285D-02
999  n1qn3: descent direction d: angle(-g,d) =  76.4 degrees
1000  -------------------------------------------------------------------------------
1001
1002  n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
1003  n1qn3: line search
1004      nlis0       fpn=-0.384D-01 d2= 0.21D-03  tmin= 0.21D-13 tmax= 0.10D+21
1005      nlis0     0.100D+01 -0.376D-01 -0.369D-01
1006      nlis0     0.100D+02 -0.318D+00 -0.257D-01
1007  n1qn3: stopping criterion on g:  0.41045D-01
1008  n1qn3: convergence rate, s(k)/s(k-1) =  0.43358D+01
1009  n1qn3: matrix update:
1010      fitting the ellipsoid: factor  0.193D+02
1011      Oren-Spedicato factor (not used) =  0.403D-01     diagonal: average value =  0.645D-01
1012  n1qn3: descent direction d: angle(-g,d) =  75.7 degrees
1013  -------------------------------------------------------------------------------
1014
1015  n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
1016  n1qn3: line search
1017      nlis0       fpn=-0.627D+00 d2= 0.71D-01  tmin= 0.11D-14 tmax= 0.10D+21
1018      nlis0     0.100D+01 -0.315D+00  0.212D+00
1019  n1qn3: stopping criterion on g:  0.60708D-01
1020  n1qn3: convergence rate, s(k)/s(k-1) =  0.18312D+01
1021  n1qn3: matrix update:
1022      fitting the ellipsoid: factor  0.259D+00
1023      Oren-Spedicato factor (not used) =  0.251D-01     diagonal: average value =  0.248D-01
1024  n1qn3: descent direction d: angle(-g,d) =  15.3 degrees
1025  -------------------------------------------------------------------------------
1026
1027  n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
1028  n1qn3: line search
1029      nlis0       fpn=-0.707D+00 d2= 0.27D-02  tmin= 0.52D-14 tmax= 0.10D+21
1030      nlis0     0.100D+01 -0.402D+00 -0.117D+00
1031  n1qn3: stopping criterion on g:  0.11240D-01
1032  n1qn3: convergence rate, s(k)/s(k-1) =  0.19476D+00
1033  n1qn3: matrix update:
1034      fitting the ellipsoid: factor  0.190D+00
1035      Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.468D-02
1036  n1qn3: descent direction d: angle(-g,d) =  33.1 degrees
1037  -------------------------------------------------------------------------------
1038
1039  n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
1040  n1qn3: line search
1041      nlis0       fpn=-0.499D+00 d2= 0.52D-01  tmin= 0.11D-14 tmax= 0.10D+21
1042      nlis0     0.100D+01 -0.225D+00  0.427D+00
1043  n1qn3: stopping criterion on g:  0.43961D-01
1044  n1qn3: convergence rate, s(k)/s(k-1) =  0.43931D+01
1045  n1qn3: matrix update:
1046      fitting the ellipsoid: factor  0.364D+01
1047      Oren-Spedicato factor (not used) =  0.102D-01     diagonal: average value =  0.311D-01
1048  n1qn3: descent direction d: angle(-g,d) =  56.6 degrees
1049  -------------------------------------------------------------------------------
1050
1051  n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
1052  n1qn3: line search
1053      nlis0       fpn=-0.339D+00 d2= 0.36D-02  tmin= 0.47D-14 tmax= 0.10D+21
1054      nlis0     0.100D+01 -0.151D+00  0.252D-01
1055  n1qn3: stopping criterion on g:  0.12127D-01
1056  n1qn3: convergence rate, s(k)/s(k-1) =  0.26394D+00
1057  n1qn3: matrix update:
1058      fitting the ellipsoid: factor  0.759D+00
1059      Oren-Spedicato factor (not used) =  0.522D-02     diagonal: average value =  0.235D-01
1060  n1qn3: descent direction d: angle(-g,d) =  36.9 degrees
1061  -------------------------------------------------------------------------------
1062
1063  n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
1064  n1qn3: line search
1065      nlis0       fpn=-0.207D+00 d2= 0.84D-02  tmin= 0.24D-14 tmax= 0.10D+21
1066      nlis0     0.100D+01 -0.186D+00 -0.162D+00
1067  n1qn3: stopping criterion on g:  0.11958D-01
1068  n1qn3: convergence rate, s(k)/s(k-1) =  0.15231D+01
1069  n1qn3: matrix update:
1070      fitting the ellipsoid: factor  0.493D+01
1071      Oren-Spedicato factor (not used) =  0.996D-01     diagonal: average value =  0.118D+00
1072  n1qn3: descent direction d: angle(-g,d) =  47.6 degrees
1073  -------------------------------------------------------------------------------
1074
1075  n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
1076  n1qn3: line search
1077      nlis0       fpn=-0.650D+00 d2= 0.12D+00  tmin= 0.64D-15 tmax= 0.10D+21
1078      nlis0                                      0.100D+01  0.169D+01  0.738D+01
1079      nlis0     0.374D+00 -0.132D+00  0.166D+00
1080  n1qn3: stopping criterion on g:  0.29691D-01
1081  n1qn3: convergence rate, s(k)/s(k-1) =  0.14137D+01
1082  n1qn3: matrix update:
1083      fitting the ellipsoid: factor  0.211D+00
1084      Oren-Spedicato factor (not used) =  0.105D-01     diagonal: average value =  0.372D-01
1085  n1qn3: descent direction d: angle(-g,d) =  53.1 degrees
1086  -------------------------------------------------------------------------------
1087
1088  n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
1089  n1qn3: line search
1090      nlis0       fpn=-0.187D+00 d2= 0.20D-02  tmin= 0.62D-14 tmax= 0.10D+21
1091      nlis0     0.100D+01 -0.149D+00 -0.115D+00
1092  n1qn3: stopping criterion on g:  0.15621D-01
1093  n1qn3: convergence rate, s(k)/s(k-1) =  0.34738D+00
1094  n1qn3: matrix update:
1095      fitting the ellipsoid: factor  0.123D+01
1096      Oren-Spedicato factor (not used) =  0.634D-02     diagonal: average value =  0.386D-01
1097  n1qn3: descent direction d: angle(-g,d) =  62.6 degrees
1098  -------------------------------------------------------------------------------
1099
1100  n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
1101  n1qn3: line search
1102      nlis0       fpn=-0.311D+00 d2= 0.34D-01  tmin= 0.13D-14 tmax= 0.10D+21
1103      nlis0     0.100D+01 -0.102D+00  0.268D+00
1104  n1qn3: stopping criterion on g:  0.42709D-01
1105  n1qn3: convergence rate, s(k)/s(k-1) =  0.41295D+01
1106  n1qn3: matrix update:
1107      fitting the ellipsoid: factor  0.285D+00
1108      Oren-Spedicato factor (not used) =  0.108D-01     diagonal: average value =  0.841D-02
1109  n1qn3: descent direction d: angle(-g,d) =  17.2 degrees
1110  -------------------------------------------------------------------------------
1111
1112  n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
1113  n1qn3: line search
1114      nlis0       fpn=-0.421D+00 d2= 0.20D-02  tmin= 0.51D-14 tmax= 0.10D+21
1115      nlis0     0.100D+01 -0.181D+00  0.576D-01
1116  n1qn3: stopping criterion on g:  0.94247D-02
1117  n1qn3: convergence rate, s(k)/s(k-1) =  0.23849D+00
1118  n1qn3: matrix update:
1119      fitting the ellipsoid: factor  0.526D+00
1120      Oren-Spedicato factor (not used) =  0.345D-02     diagonal: average value =  0.461D-02
1121  n1qn3: descent direction d: angle(-g,d) =  57.2 degrees
1122  -------------------------------------------------------------------------------
1123
1124  n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
1125  n1qn3: line search
1126      nlis0       fpn=-0.126D+00 d2= 0.11D-01  tmin= 0.25D-14 tmax= 0.10D+21
1127      nlis0     0.100D+01 -0.110D+00 -0.725D-01
1128  n1qn3: stopping criterion on g:  0.67714D-02
1129  n1qn3: convergence rate, s(k)/s(k-1) =  0.23948D+01
1130  n1qn3: matrix update:
1131      fitting the ellipsoid: factor  0.105D+01
1132      Oren-Spedicato factor (not used) =  0.489D-02     diagonal: average value =  0.547D-02
1133  n1qn3: descent direction d: angle(-g,d) =  59.9 degrees
1134  -------------------------------------------------------------------------------
1135
1136  n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
1137  n1qn3: line search
1138      nlis0       fpn=-0.228D+00 d2= 0.83D-01  tmin= 0.95D-15 tmax= 0.10D+21
1139      nlis0                                      0.100D+01  0.601D+00  0.234D+01
1140      nlis0     0.321D+00 -0.373D-01  0.472D-01
1141  n1qn3: stopping criterion on g:  0.24287D-01
1142  n1qn3: convergence rate, s(k)/s(k-1) =  0.87291D+00
1143  n1qn3: matrix update:
1144      fitting the ellipsoid: factor  0.807D+00
1145      Oren-Spedicato factor (not used) =  0.458D-02     diagonal: average value =  0.431D-02
1146  n1qn3: descent direction d: angle(-g,d) =  73.3 degrees
1147  -------------------------------------------------------------------------------
1148
1149  n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
1150  n1qn3: line search
1151      nlis0       fpn=-0.554D-01 d2= 0.12D-02  tmin= 0.80D-14 tmax= 0.10D+21
1152      nlis0     0.100D+01 -0.477D-01 -0.406D-01
1153  n1qn3: stopping criterion on g:  0.15832D-01
1154  n1qn3: convergence rate, s(k)/s(k-1) =  0.36912D+00
1155  n1qn3: matrix update:
1156      fitting the ellipsoid: factor  0.894D+00
1157      Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.392D-02
1158  n1qn3: descent direction d: angle(-g,d) =  77.1 degrees
1159  -------------------------------------------------------------------------------
1160
1161  n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
1162  n1qn3: line search
1163      nlis0       fpn=-0.143D+00 d2= 0.30D-01  tmin= 0.17D-14 tmax= 0.10D+21
1164      nlis0     0.100D+01 -0.878D-01 -0.226D-01
1165  n1qn3: stopping criterion on g:  0.23266D-01
1166  n1qn3: convergence rate, s(k)/s(k-1) =  0.51063D+01
1167  n1qn3: matrix update:
1168      fitting the ellipsoid: factor  0.865D+01
1169      Oren-Spedicato factor (not used) =  0.318D-01     diagonal: average value =  0.403D-01
1170  n1qn3: descent direction d: angle(-g,d) =  75.1 degrees
1171  -------------------------------------------------------------------------------
1172
1173  n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
1174  n1qn3: line search
1175      nlis0       fpn=-0.128D+00 d2= 0.84D-02  tmin= 0.27D-14 tmax= 0.10D+21
1176      nlis0     0.100D+01 -0.619D-01 -0.487D-02
1177  n1qn3: stopping criterion on g:  0.10083D-01
1178  n1qn3: convergence rate, s(k)/s(k-1) =  0.52582D+00
1179  n1qn3: matrix update:
1180      fitting the ellipsoid: factor  0.572D-01
1181      Oren-Spedicato factor (not used) =  0.208D-02     diagonal: average value =  0.227D-02
1182  n1qn3: descent direction d: angle(-g,d) =  78.8 degrees
1183  -------------------------------------------------------------------------------
1184
1185  n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
1186  n1qn3: line search
1187      nlis0       fpn=-0.100D+00 d2= 0.49D-01  tmin= 0.13D-14 tmax= 0.10D+21
1188      nlis0     0.100D+01 -0.155D-01  0.179D+00
1189  n1qn3: stopping criterion on g:  0.37918D-01
1190  n1qn3: convergence rate, s(k)/s(k-1) =  0.24062D+01
1191  n1qn3: matrix update:
1192      fitting the ellipsoid: factor  0.120D+01
1193      Oren-Spedicato factor (not used) =  0.224D-02     diagonal: average value =  0.267D-02
1194  n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
1195  -------------------------------------------------------------------------------
1196
1197  n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
1198  n1qn3: line search
1199      nlis0       fpn=-0.120D+00 d2= 0.12D-01  tmin= 0.26D-14 tmax= 0.10D+21
1200      nlis0     0.100D+01 -0.251D-01  0.437D-01
1201  n1qn3: stopping criterion on g:  0.43947D-02
1202  n1qn3: convergence rate, s(k)/s(k-1) =  0.50600D+00
1203  n1qn3: matrix update:
1204      fitting the ellipsoid: factor  0.757D+00
1205      Oren-Spedicato factor (not used) =  0.170D-02     diagonal: average value =  0.201D-02
1206  n1qn3: descent direction d: angle(-g,d) =  72.3 degrees
1207  -------------------------------------------------------------------------------
1208
1209  n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
1210  n1qn3: line search
1211      nlis0       fpn=-0.208D-01 d2= 0.45D-02  tmin= 0.41D-14 tmax= 0.10D+21
1212      nlis0     0.100D+01 -0.183D-01 -0.148D-01
1213  n1qn3: stopping criterion on g:  0.21889D-02
1214  n1qn3: convergence rate, s(k)/s(k-1) =  0.59772D+00
1215  n1qn3: matrix update:
1216      fitting the ellipsoid: factor  0.170D+01
1217      Oren-Spedicato factor (not used) =  0.290D-02     diagonal: average value =  0.351D-02
1218  n1qn3: descent direction d: angle(-g,d) =  63.1 degrees
1219  -------------------------------------------------------------------------------
1220
1221  n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
1222  n1qn3: line search
1223      nlis0       fpn=-0.555D-01 d2= 0.58D-01  tmin= 0.11D-14 tmax= 0.10D+21
1224      nlis0                                      0.100D+01  0.507D-01  0.260D+00
1225      nlis0     0.413D+00 -0.132D-01  0.278D-02
1226  n1qn3: stopping criterion on g:  0.14651D-01
1227  n1qn3: convergence rate, s(k)/s(k-1) =  0.14924D+01
1228  n1qn3: matrix update:
1229      fitting the ellipsoid: factor  0.980D+00
1230      Oren-Spedicato factor (not used) =  0.274D-02     diagonal: average value =  0.355D-02
1231  n1qn3: descent direction d: angle(-g,d) =  84.8 degrees
1232  -------------------------------------------------------------------------------
1233
1234  n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
1235  n1qn3: line search
1236      nlis0       fpn=-0.918D-02 d2= 0.87D-03  tmin= 0.86D-14 tmax= 0.10D+21
1237      nlis0     0.100D+01 -0.813D-02 -0.718D-02
1238  n1qn3: stopping criterion on g:  0.11306D-01
1239  n1qn3: convergence rate, s(k)/s(k-1) =  0.29590D+00
1240  n1qn3: matrix update:
1241      fitting the ellipsoid: factor  0.112D+01
1242      Oren-Spedicato factor (not used) =  0.329D-02     diagonal: average value =  0.403D-02
1243  n1qn3: descent direction d: angle(-g,d) =  85.0 degrees
1244  -------------------------------------------------------------------------------
1245
1246  n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
1247  n1qn3: line search
1248      nlis0       fpn=-0.259D-01 d2= 0.13D-01  tmin= 0.22D-14 tmax= 0.10D+21
1249      nlis0     0.100D+01 -0.174D-01 -0.109D-01
1250  n1qn3: stopping criterion on g:  0.64294D-02
1251  n1qn3: convergence rate, s(k)/s(k-1) =  0.38281D+01
1252  n1qn3: matrix update:
1253      fitting the ellipsoid: factor  0.337D+01
1254      Oren-Spedicato factor (not used) =  0.116D-01     diagonal: average value =  0.134D-01
1255  n1qn3: descent direction d: angle(-g,d) =  84.3 degrees
1256  -------------------------------------------------------------------------------
1257
1258  n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
1259  n1qn3: line search
1260      nlis0       fpn=-0.949D-02 d2= 0.41D-02  tmin= 0.39D-14 tmax= 0.10D+21
1261      nlis0     0.100D+01 -0.648D-02 -0.412D-02
1262  n1qn3: stopping criterion on g:  0.32255D-03
1263  n1qn3: convergence rate, s(k)/s(k-1) =  0.56905D+00
1264  n1qn3: matrix update:
1265      fitting the ellipsoid: factor  0.218D+00
1266      Oren-Spedicato factor (not used) =  0.229D-02     diagonal: average value =  0.316D-02
1267  n1qn3: descent direction d: angle(-g,d) =  28.9 degrees
1268  -------------------------------------------------------------------------------
1269
1270  n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
1271  n1qn3: line search
1272      nlis0       fpn=-0.127D-01 d2= 0.37D-01  tmin= 0.13D-14 tmax= 0.10D+21
1273      nlis0                                      0.100D+01  0.236D-01  0.889D-01
1274      nlis0     0.308D+00 -0.202D-02  0.129D-02
1275  n1qn3: stopping criterion on g:  0.59635D-02
1276  n1qn3: convergence rate, s(k)/s(k-1) =  0.92424D+00
1277  n1qn3: matrix update:
1278      fitting the ellipsoid: factor  0.995D+00
1279      Oren-Spedicato factor (not used) =  0.211D-02     diagonal: average value =  0.320D-02
1280  n1qn3: descent direction d: angle(-g,d) =  86.9 degrees
1281  -------------------------------------------------------------------------------
1282
1283  n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
1284  n1qn3: line search
1285      nlis0       fpn=-0.270D-02 d2= 0.13D-02  tmin= 0.68D-14 tmax= 0.10D+21
1286      nlis0     0.100D+01 -0.193D-02 -0.124D-02
1287  n1qn3: stopping criterion on g:  0.26088D-02
1288  n1qn3: convergence rate, s(k)/s(k-1) =  0.61219D+00
1289  n1qn3: matrix update:
1290      fitting the ellipsoid: factor  0.107D+01
1291      Oren-Spedicato factor (not used) =  0.240D-02     diagonal: average value =  0.346D-02
1292  n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
1293  -------------------------------------------------------------------------------
1294
1295  n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
1296  n1qn3: line search
1297      nlis0       fpn=-0.106D-02 d2= 0.11D-02  tmin= 0.73D-14 tmax= 0.10D+21
1298      nlis0     0.100D+01 -0.680D-03 -0.342D-03
1299  n1qn3: stopping criterion on g:  0.93209D-03
1300  n1qn3: convergence rate, s(k)/s(k-1) =  0.93175D+00
1301  n1qn3: matrix update:
1302      fitting the ellipsoid: factor  0.195D+01
1303      Oren-Spedicato factor (not used) =  0.472D-02     diagonal: average value =  0.669D-02
1304  n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
1305  -------------------------------------------------------------------------------
1306
1307  n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
1308  n1qn3: line search
1309      nlis0       fpn=-0.165D-03 d2= 0.22D-03  tmin= 0.17D-13 tmax= 0.10D+21
1310      nlis0     0.100D+01 -0.994D-04 -0.378D-04
1311  n1qn3: stopping criterion on g:  0.13605D-03
1312  n1qn3: convergence rate, s(k)/s(k-1) =  0.43415D+00
1313  n1qn3: matrix update:
1314      fitting the ellipsoid: factor  0.797D+00
1315      Oren-Spedicato factor (not used) =  0.369D-02     diagonal: average value =  0.540D-02
1316  n1qn3: descent direction d: angle(-g,d) =  86.0 degrees
1317  -------------------------------------------------------------------------------
1318
1319  n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
1320  n1qn3: line search
1321      nlis0       fpn=-0.128D-04 d2= 0.34D-04  tmin= 0.43D-13 tmax= 0.10D+21
1322      nlis0     0.100D+01 -0.592D-05  0.983D-06
1323  n1qn3: stopping criterion on g:  0.15294D-03
1324  n1qn3: convergence rate, s(k)/s(k-1) =  0.39644D+00
1325  n1qn3: matrix update:
1326      fitting the ellipsoid: factor  0.278D+03
1327      Oren-Spedicato factor (not used) =  0.637D+00     diagonal: average value =  0.256D+01
1328  n1qn3: descent direction d: angle(-g,d) =  77.4 degrees
1329  -------------------------------------------------------------------------------
1330
1331  n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
1332  n1qn3: line search
1333      nlis0       fpn=-0.241D-05 d2= 0.96D-07  tmin= 0.92D-12 tmax= 0.10D+21
1334      nlis0     0.100D+01 -0.168D-06  0.208D-05
1335  n1qn3: stopping criterion on g:  0.13388D-03
1336  n1qn3: convergence rate, s(k)/s(k-1) =  0.53137D-01
1337  n1qn3: matrix update:
1338      fitting the ellipsoid: factor  0.794D-03
1339      Oren-Spedicato factor (not used) =  0.101D-02     diagonal: average value =  0.159D-02
1340  n1qn3: descent direction d: angle(-g,d) =  26.4 degrees
1341  -------------------------------------------------------------------------------
1342
1343  n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
1344  n1qn3: line search
1345      nlis0       fpn=-0.970D-06 d2= 0.12D-08  tmin= 0.80D-11 tmax= 0.10D+21
1346      nlis0     0.100D+01 -0.485D-06 -0.183D-09
1347  n1qn3: stopping criterion on g:  0.27051D-07
1348  n1qn3: convergence rate, s(k)/s(k-1) =  0.11218D+00
1349  n1qn3: matrix update:
1350      fitting the ellipsoid: factor  0.866D+00
1351      Oren-Spedicato factor (not used) =  0.999D-03     diagonal: average value =  0.137D-02
1352  n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
1353  -------------------------------------------------------------------------------
1354
1355  n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
1356  n1qn3: line search
1357      nlis0       fpn=-0.143D-11 d2= 0.35D-11  tmin= 0.13D-09 tmax= 0.10D+21
1358      nlis0     0.100D+01 -0.734D-12 -0.116D-13
1359  n1qn3: stopping criterion on g:  0.61631D-10
1360  n1qn3: convergence rate, s(k)/s(k-1) =  0.53526D-01
1361  n1qn3: matrix update:
1362      fitting the ellipsoid: factor  0.399D+02
1363      Oren-Spedicato factor (not used) =  0.356D-01     diagonal: average value =  0.717D-01
1364  n1qn3: descent direction d: angle(-g,d) =  64.0 degrees
1365  -------------------------------------------------------------------------------
1366
1367  n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
1368  n1qn3: line search
1369      nlis0       fpn=-0.991D-16 d2= 0.25D-15  tmin= 0.16D-07 tmax= 0.10D+21
1370      nlis0     0.100D+01 -0.152D-15 -0.112D-20
1371  n1qn3: stopping criterion on g:  0.92861D-13
1372  n1qn3: convergence rate, s(k)/s(k-1) =  0.84647D-02
1373  n1qn3: matrix update:
1374      fitting the ellipsoid: factor  0.571D+01
1375      Oren-Spedicato factor (not used) =  0.482D+00     diagonal: average value =  0.356D+00
1376  n1qn3: descent direction d: angle(-g,d) =  80.8 degrees
1377  -------------------------------------------------------------------------------
1378
1379  n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
1380  n1qn3: line search
1381      nlis0       fpn=-0.474D-24 d2= 0.19D-25  tmin= 0.17D-02 tmax= 0.10D+21
1382      nlis0     0.100D+01 -0.599D-21 -0.107D-26
1383  n1qn3: stopping criterion on g:  0.21365D-15
1384  n1qn3: convergence rate, s(k)/s(k-1) =  0.86648D-05
1385  n1qn3: matrix update:
1386      fitting the ellipsoid: factor  0.352D-02
1387      Oren-Spedicato factor (not used) =  0.102D-02     diagonal: average value =  0.198D-02
1388  n1qn3: descent direction d: angle(-g,d) =  29.9 degrees
1389  -------------------------------------------------------------------------------
1390
1391  n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
1392  n1qn3: line search
1393      nlis0       fpn=-0.248D-29 d2= 0.33D-32  tmin= 0.39D+01 tmax= 0.10D+21
1394      nlis0                                      0.387D+01  0.326D-23  0.743D-29
1395      nlis0   fin sur tmin        pas            fonctions     derivees
1396      nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
1397      nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
1398  -------------------------------------------------------------------------------
1399  n1qn3: output mode is  6
1400      number of iterations:   42
1401      number of simulations:     53
1402      realized relative precision on g:  0.21D-15
1403      norm of x =  0.14142135D+01
1404      f         =  0.53783484D-16
1405      norm of x =  0.14924691D-12
1406 Optimization stops because too small variations for x.
1407