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