Add *.win.dia.ref in optimization module. 16/12616/2
Adeline CARNIS [Fri, 20 Sep 2013 09:13:36 +0000 (11:13 +0200)]
* test_run('optimization', 'bug_12702')
* test_run('optimization', 'bug_9690')

Change-Id: I362e1bece884a81a0b668007fee067d831e7185a

scilab/modules/optimization/tests/nonreg_tests/bug_12702.win.dia.ref [new file with mode: 0644]
scilab/modules/optimization/tests/nonreg_tests/bug_9690.win.dia.ref [new file with mode: 0644]

diff --git a/scilab/modules/optimization/tests/nonreg_tests/bug_12702.win.dia.ref b/scilab/modules/optimization/tests/nonreg_tests/bug_12702.win.dia.ref
new file mode 100644 (file)
index 0000000..d43d33e
--- /dev/null
@@ -0,0 +1,1405 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2013 - Scilab Enterprises - Paul Bignier
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+//
+// <-- ENGLISH IMPOSED -->
+//
+// <-- Non-regression test for bug 12702 -->
+//
+// <-- Bugzilla URL -->
+// http://bugzilla.scilab.org/show_bug.cgi?id=12702
+//
+// <-- Short Description -->
+// When no extra parameters are needed in the cost function, NDcost does not work:
+//
+function f = rosenbrock(x)
+    f = 100.0 * (x(2)-x(1)^2)^2 + (1-x(1))^2;
+endfunction
+function [f, g, ind] = rosenbrockCost(x, ind)
+    if ((ind == 1) | (ind == 4)) then
+        f = rosenbrock ( x );
+    end
+    if ((ind == 1) | (ind == 4)) then
+        g = derivative ( rosenbrock , x(:) );
+    end
+endfunction
+x0 = [-1.2 1.0];
+lines(0);
+for impval = 1:5
+    [ fopt , xopt ] = optim ( rosenbrockCost , x0 , "gc" , imp=impval);
+end
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   1
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   2
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+ n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
+ n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
+ n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
+ n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
+ n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
+ n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
+ n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
+ n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
+ n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
+ n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
+ n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
+ n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
+ n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
+ n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
+ n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
+ n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
+ n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
+ n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
+ n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
+ n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
+ n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
+ n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
+ n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
+ n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
+ n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
+ n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
+ n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
+ n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
+ n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
+ n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
+ n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
+ n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
+ n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
+ n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
+ n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
+ n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
+ n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
+ n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
+ n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
+ n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
+ n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
+ n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   3
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+ -------------------------------------------------------------------------------
+
+
+ n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
+     nlis0       fpn=-0.200D+01 d2= 0.74D-04  tmin= 0.28D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.194D+01 -0.189D+01
+     nlis0     0.100D+02 -0.147D+02 -0.973D+00
+
+ n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
+     nlis0       fpn=-0.924D+01 d2= 0.66D-02  tmin= 0.30D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.520D+01 -0.138D+01
+
+ n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
+     nlis0       fpn=-0.245D+00 d2= 0.21D-03  tmin= 0.16D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.135D+00 -0.256D-01
+
+ n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
+     nlis0       fpn=-0.581D-02 d2= 0.59D-05  tmin= 0.98D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.435D-02 -0.290D-02
+
+ n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
+     nlis0       fpn=-0.833D-02 d2= 0.26D-04  tmin= 0.61D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.690D-02 -0.549D-02
+
+ n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
+     nlis0       fpn=-0.278D-01 d2= 0.37D-03  tmin= 0.14D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.228D-01 -0.180D-01
+
+ n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
+     nlis0       fpn=-0.828D-01 d2= 0.40D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.713D-01 -0.624D-01
+
+ n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
+     nlis0       fpn=-0.405D+00 d2= 0.11D+00  tmin= 0.79D-15 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.531D+00 -0.721D+00
+     nlis0                                      0.100D+02  0.500D+03  0.222D+03
+     nlis0                                      0.398D+01  0.333D+01  0.774D+01
+     nlis0     0.249D+01 -0.140D+01  0.190D+00
+
+ n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
+     nlis0       fpn=-0.161D+00 d2= 0.57D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.285D+00  0.571D+00
+     nlis0     0.140D+00 -0.107D-01  0.342D-02
+
+ n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
+     nlis0       fpn=-0.384D-01 d2= 0.21D-03  tmin= 0.21D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.376D-01 -0.369D-01
+     nlis0     0.100D+02 -0.318D+00 -0.257D-01
+
+ n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
+     nlis0       fpn=-0.627D+00 d2= 0.71D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.315D+00  0.212D+00
+
+ n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
+     nlis0       fpn=-0.707D+00 d2= 0.27D-02  tmin= 0.52D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.402D+00 -0.117D+00
+
+ n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
+     nlis0       fpn=-0.499D+00 d2= 0.52D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.225D+00  0.427D+00
+
+ n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
+     nlis0       fpn=-0.339D+00 d2= 0.36D-02  tmin= 0.47D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.151D+00  0.252D-01
+
+ n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
+     nlis0       fpn=-0.207D+00 d2= 0.84D-02  tmin= 0.24D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.186D+00 -0.162D+00
+
+ n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
+     nlis0       fpn=-0.650D+00 d2= 0.12D+00  tmin= 0.64D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.169D+01  0.738D+01
+     nlis0     0.374D+00 -0.132D+00  0.166D+00
+
+ n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
+     nlis0       fpn=-0.187D+00 d2= 0.20D-02  tmin= 0.62D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.149D+00 -0.115D+00
+
+ n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
+     nlis0       fpn=-0.311D+00 d2= 0.34D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.102D+00  0.268D+00
+
+ n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
+     nlis0       fpn=-0.421D+00 d2= 0.20D-02  tmin= 0.51D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.181D+00  0.576D-01
+
+ n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
+     nlis0       fpn=-0.126D+00 d2= 0.11D-01  tmin= 0.25D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.110D+00 -0.725D-01
+
+ n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
+     nlis0       fpn=-0.228D+00 d2= 0.83D-01  tmin= 0.95D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.601D+00  0.234D+01
+     nlis0     0.321D+00 -0.373D-01  0.472D-01
+
+ n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
+     nlis0       fpn=-0.554D-01 d2= 0.12D-02  tmin= 0.80D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.477D-01 -0.406D-01
+
+ n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
+     nlis0       fpn=-0.143D+00 d2= 0.30D-01  tmin= 0.17D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.878D-01 -0.226D-01
+
+ n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
+     nlis0       fpn=-0.128D+00 d2= 0.84D-02  tmin= 0.27D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.619D-01 -0.487D-02
+
+ n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
+     nlis0       fpn=-0.100D+00 d2= 0.49D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.155D-01  0.179D+00
+
+ n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
+     nlis0       fpn=-0.120D+00 d2= 0.12D-01  tmin= 0.26D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.251D-01  0.437D-01
+
+ n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
+     nlis0       fpn=-0.208D-01 d2= 0.45D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.183D-01 -0.148D-01
+
+ n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
+     nlis0       fpn=-0.555D-01 d2= 0.58D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.507D-01  0.260D+00
+     nlis0     0.413D+00 -0.132D-01  0.278D-02
+
+ n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
+     nlis0       fpn=-0.918D-02 d2= 0.87D-03  tmin= 0.86D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.813D-02 -0.718D-02
+
+ n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
+     nlis0       fpn=-0.259D-01 d2= 0.13D-01  tmin= 0.22D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.174D-01 -0.109D-01
+
+ n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
+     nlis0       fpn=-0.949D-02 d2= 0.41D-02  tmin= 0.39D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.648D-02 -0.412D-02
+
+ n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
+     nlis0       fpn=-0.127D-01 d2= 0.37D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.236D-01  0.889D-01
+     nlis0     0.308D+00 -0.202D-02  0.129D-02
+
+ n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
+     nlis0       fpn=-0.270D-02 d2= 0.13D-02  tmin= 0.68D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.193D-02 -0.124D-02
+
+ n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
+     nlis0       fpn=-0.106D-02 d2= 0.11D-02  tmin= 0.73D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.680D-03 -0.342D-03
+
+ n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
+     nlis0       fpn=-0.165D-03 d2= 0.22D-03  tmin= 0.17D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.994D-04 -0.378D-04
+
+ n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
+     nlis0       fpn=-0.128D-04 d2= 0.34D-04  tmin= 0.43D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.592D-05  0.983D-06
+
+ n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
+     nlis0       fpn=-0.241D-05 d2= 0.96D-07  tmin= 0.92D-12 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.168D-06  0.208D-05
+
+ n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
+     nlis0       fpn=-0.970D-06 d2= 0.12D-08  tmin= 0.80D-11 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.485D-06 -0.183D-09
+
+ n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
+     nlis0       fpn=-0.143D-11 d2= 0.35D-11  tmin= 0.13D-09 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.734D-12 -0.116D-13
+
+ n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
+     nlis0       fpn=-0.991D-16 d2= 0.25D-15  tmin= 0.16D-07 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.152D-15 -0.112D-20
+
+ n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
+     nlis0       fpn=-0.474D-24 d2= 0.19D-25  tmin= 0.17D-02 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.599D-21 -0.107D-26
+
+ n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
+     nlis0       fpn=-0.248D-29 d2= 0.33D-32  tmin= 0.39D+01 tmax= 0.10D+21
+     nlis0                                      0.387D+01  0.326D-23  0.743D-29
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   4
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+ n1qn3a: descent direction -g: precon =  0.369D-04
+ -------------------------------------------------------------------------------
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
+ n1qn3: line search
+     nlis0       fpn=-0.200D+01 d2= 0.74D-04  tmin= 0.28D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.194D+01 -0.189D+01
+     nlis0     0.100D+02 -0.147D+02 -0.973D+00
+ n1qn3: stopping criterion on g:  0.48674D+00
+ n1qn3: matrix update:
+     Oren-Spedicato factor (not used) =  0.718D-03     diagonal: average value =  0.718D-03
+ n1qn3: descent direction d: angle(-g,d) =   0.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
+ n1qn3: line search
+     nlis0       fpn=-0.924D+01 d2= 0.66D-02  tmin= 0.30D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.520D+01 -0.138D+01
+ n1qn3: stopping criterion on g:  0.73079D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.94908D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.118D+01
+     Oren-Spedicato factor (not used) =  0.846D-03     diagonal: average value =  0.846D-03
+ n1qn3: descent direction d: angle(-g,d) =   0.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
+ n1qn3: line search
+     nlis0       fpn=-0.245D+00 d2= 0.21D-03  tmin= 0.16D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.135D+00 -0.256D-01
+ n1qn3: stopping criterion on g:  0.10270D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.17647D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.111D+01
+     Oren-Spedicato factor (not used) =  0.935D-03     diagonal: average value =  0.942D-03
+ n1qn3: descent direction d: angle(-g,d) =   1.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
+ n1qn3: line search
+     nlis0       fpn=-0.581D-02 d2= 0.59D-05  tmin= 0.98D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.435D-02 -0.290D-02
+ n1qn3: stopping criterion on g:  0.76379D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.16885D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.969D+00
+     Oren-Spedicato factor (not used) =  0.952D-03     diagonal: average value =  0.915D-03
+ n1qn3: descent direction d: angle(-g,d) =  23.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
+ n1qn3: line search
+     nlis0       fpn=-0.833D-02 d2= 0.26D-04  tmin= 0.61D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.690D-02 -0.549D-02
+ n1qn3: stopping criterion on g:  0.11072D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.21032D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.985D+00
+     Oren-Spedicato factor (not used) =  0.954D-03     diagonal: average value =  0.102D-02
+ n1qn3: descent direction d: angle(-g,d) =  56.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
+ n1qn3: line search
+     nlis0       fpn=-0.278D-01 d2= 0.37D-03  tmin= 0.14D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.228D-01 -0.180D-01
+ n1qn3: stopping criterion on g:  0.23226D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.37768D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.113D+01
+     Oren-Spedicato factor (not used) =  0.941D-03     diagonal: average value =  0.127D-02
+ n1qn3: descent direction d: angle(-g,d) =  76.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
+ n1qn3: line search
+     nlis0       fpn=-0.828D-01 d2= 0.40D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.713D-01 -0.624D-01
+ n1qn3: stopping criterion on g:  0.44894D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.32918D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.817D+00
+     Oren-Spedicato factor (not used) =  0.762D-03     diagonal: average value =  0.106D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
+ n1qn3: line search
+     nlis0       fpn=-0.405D+00 d2= 0.11D+00  tmin= 0.79D-15 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.531D+00 -0.721D+00
+     nlis0                                      0.100D+02  0.500D+03  0.222D+03
+     nlis0                                      0.398D+01  0.333D+01  0.774D+01
+     nlis0     0.249D+01 -0.140D+01  0.190D+00
+ n1qn3: stopping criterion on g:  0.60049D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.12764D+02
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.287D+01
+     Oren-Spedicato factor (not used) =  0.248D-02     diagonal: average value =  0.292D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
+ n1qn3: line search
+     nlis0       fpn=-0.161D+00 d2= 0.57D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.285D+00  0.571D+00
+     nlis0     0.140D+00 -0.107D-01  0.342D-02
+ n1qn3: stopping criterion on g:  0.48319D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.41360D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.102D+01
+     Oren-Spedicato factor (not used) =  0.293D-02     diagonal: average value =  0.285D-02
+ n1qn3: descent direction d: angle(-g,d) =  76.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
+ n1qn3: line search
+     nlis0       fpn=-0.384D-01 d2= 0.21D-03  tmin= 0.21D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.376D-01 -0.369D-01
+     nlis0     0.100D+02 -0.318D+00 -0.257D-01
+ n1qn3: stopping criterion on g:  0.41045D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43358D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.193D+02
+     Oren-Spedicato factor (not used) =  0.403D-01     diagonal: average value =  0.645D-01
+ n1qn3: descent direction d: angle(-g,d) =  75.7 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
+ n1qn3: line search
+     nlis0       fpn=-0.627D+00 d2= 0.71D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.315D+00  0.212D+00
+ n1qn3: stopping criterion on g:  0.60708D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.18312D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.259D+00
+     Oren-Spedicato factor (not used) =  0.251D-01     diagonal: average value =  0.248D-01
+ n1qn3: descent direction d: angle(-g,d) =  15.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
+ n1qn3: line search
+     nlis0       fpn=-0.707D+00 d2= 0.27D-02  tmin= 0.52D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.402D+00 -0.117D+00
+ n1qn3: stopping criterion on g:  0.11240D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.19476D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.190D+00
+     Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.468D-02
+ n1qn3: descent direction d: angle(-g,d) =  33.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
+ n1qn3: line search
+     nlis0       fpn=-0.499D+00 d2= 0.52D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.225D+00  0.427D+00
+ n1qn3: stopping criterion on g:  0.43961D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43931D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.364D+01
+     Oren-Spedicato factor (not used) =  0.102D-01     diagonal: average value =  0.311D-01
+ n1qn3: descent direction d: angle(-g,d) =  56.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
+ n1qn3: line search
+     nlis0       fpn=-0.339D+00 d2= 0.36D-02  tmin= 0.47D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.151D+00  0.252D-01
+ n1qn3: stopping criterion on g:  0.12127D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.26394D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.759D+00
+     Oren-Spedicato factor (not used) =  0.522D-02     diagonal: average value =  0.235D-01
+ n1qn3: descent direction d: angle(-g,d) =  36.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
+ n1qn3: line search
+     nlis0       fpn=-0.207D+00 d2= 0.84D-02  tmin= 0.24D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.186D+00 -0.162D+00
+ n1qn3: stopping criterion on g:  0.11958D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.15231D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.493D+01
+     Oren-Spedicato factor (not used) =  0.996D-01     diagonal: average value =  0.118D+00
+ n1qn3: descent direction d: angle(-g,d) =  47.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
+ n1qn3: line search
+     nlis0       fpn=-0.650D+00 d2= 0.12D+00  tmin= 0.64D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.169D+01  0.738D+01
+     nlis0     0.374D+00 -0.132D+00  0.166D+00
+ n1qn3: stopping criterion on g:  0.29691D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.14137D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.211D+00
+     Oren-Spedicato factor (not used) =  0.105D-01     diagonal: average value =  0.372D-01
+ n1qn3: descent direction d: angle(-g,d) =  53.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
+ n1qn3: line search
+     nlis0       fpn=-0.187D+00 d2= 0.20D-02  tmin= 0.62D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.149D+00 -0.115D+00
+ n1qn3: stopping criterion on g:  0.15621D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.34738D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.123D+01
+     Oren-Spedicato factor (not used) =  0.634D-02     diagonal: average value =  0.386D-01
+ n1qn3: descent direction d: angle(-g,d) =  62.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
+ n1qn3: line search
+     nlis0       fpn=-0.311D+00 d2= 0.34D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.102D+00  0.268D+00
+ n1qn3: stopping criterion on g:  0.42709D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.41295D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.285D+00
+     Oren-Spedicato factor (not used) =  0.108D-01     diagonal: average value =  0.841D-02
+ n1qn3: descent direction d: angle(-g,d) =  17.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
+ n1qn3: line search
+     nlis0       fpn=-0.421D+00 d2= 0.20D-02  tmin= 0.51D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.181D+00  0.576D-01
+ n1qn3: stopping criterion on g:  0.94247D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.23849D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.526D+00
+     Oren-Spedicato factor (not used) =  0.345D-02     diagonal: average value =  0.461D-02
+ n1qn3: descent direction d: angle(-g,d) =  57.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
+ n1qn3: line search
+     nlis0       fpn=-0.126D+00 d2= 0.11D-01  tmin= 0.25D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.110D+00 -0.725D-01
+ n1qn3: stopping criterion on g:  0.67714D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.23948D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.105D+01
+     Oren-Spedicato factor (not used) =  0.489D-02     diagonal: average value =  0.547D-02
+ n1qn3: descent direction d: angle(-g,d) =  59.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
+ n1qn3: line search
+     nlis0       fpn=-0.228D+00 d2= 0.83D-01  tmin= 0.95D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.601D+00  0.234D+01
+     nlis0     0.321D+00 -0.373D-01  0.472D-01
+ n1qn3: stopping criterion on g:  0.24287D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.87291D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.807D+00
+     Oren-Spedicato factor (not used) =  0.458D-02     diagonal: average value =  0.431D-02
+ n1qn3: descent direction d: angle(-g,d) =  73.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
+ n1qn3: line search
+     nlis0       fpn=-0.554D-01 d2= 0.12D-02  tmin= 0.80D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.477D-01 -0.406D-01
+ n1qn3: stopping criterion on g:  0.15832D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.36912D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.894D+00
+     Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.392D-02
+ n1qn3: descent direction d: angle(-g,d) =  77.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
+ n1qn3: line search
+     nlis0       fpn=-0.143D+00 d2= 0.30D-01  tmin= 0.17D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.878D-01 -0.226D-01
+ n1qn3: stopping criterion on g:  0.23266D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.51063D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.865D+01
+     Oren-Spedicato factor (not used) =  0.318D-01     diagonal: average value =  0.403D-01
+ n1qn3: descent direction d: angle(-g,d) =  75.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
+ n1qn3: line search
+     nlis0       fpn=-0.128D+00 d2= 0.84D-02  tmin= 0.27D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.619D-01 -0.487D-02
+ n1qn3: stopping criterion on g:  0.10083D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.52582D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.572D-01
+     Oren-Spedicato factor (not used) =  0.208D-02     diagonal: average value =  0.227D-02
+ n1qn3: descent direction d: angle(-g,d) =  78.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
+ n1qn3: line search
+     nlis0       fpn=-0.100D+00 d2= 0.49D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.155D-01  0.179D+00
+ n1qn3: stopping criterion on g:  0.37918D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.24062D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.120D+01
+     Oren-Spedicato factor (not used) =  0.224D-02     diagonal: average value =  0.267D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
+ n1qn3: line search
+     nlis0       fpn=-0.120D+00 d2= 0.12D-01  tmin= 0.26D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.251D-01  0.437D-01
+ n1qn3: stopping criterion on g:  0.43947D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.50600D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.757D+00
+     Oren-Spedicato factor (not used) =  0.170D-02     diagonal: average value =  0.201D-02
+ n1qn3: descent direction d: angle(-g,d) =  72.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
+ n1qn3: line search
+     nlis0       fpn=-0.208D-01 d2= 0.45D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.183D-01 -0.148D-01
+ n1qn3: stopping criterion on g:  0.21889D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.59772D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.170D+01
+     Oren-Spedicato factor (not used) =  0.290D-02     diagonal: average value =  0.351D-02
+ n1qn3: descent direction d: angle(-g,d) =  63.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
+ n1qn3: line search
+     nlis0       fpn=-0.555D-01 d2= 0.58D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.507D-01  0.260D+00
+     nlis0     0.413D+00 -0.132D-01  0.278D-02
+ n1qn3: stopping criterion on g:  0.14651D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.14924D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.980D+00
+     Oren-Spedicato factor (not used) =  0.274D-02     diagonal: average value =  0.355D-02
+ n1qn3: descent direction d: angle(-g,d) =  84.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
+ n1qn3: line search
+     nlis0       fpn=-0.918D-02 d2= 0.87D-03  tmin= 0.86D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.813D-02 -0.718D-02
+ n1qn3: stopping criterion on g:  0.11306D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.29590D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.112D+01
+     Oren-Spedicato factor (not used) =  0.329D-02     diagonal: average value =  0.403D-02
+ n1qn3: descent direction d: angle(-g,d) =  85.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
+ n1qn3: line search
+     nlis0       fpn=-0.259D-01 d2= 0.13D-01  tmin= 0.22D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.174D-01 -0.109D-01
+ n1qn3: stopping criterion on g:  0.64294D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.38281D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.337D+01
+     Oren-Spedicato factor (not used) =  0.116D-01     diagonal: average value =  0.134D-01
+ n1qn3: descent direction d: angle(-g,d) =  84.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
+ n1qn3: line search
+     nlis0       fpn=-0.949D-02 d2= 0.41D-02  tmin= 0.39D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.648D-02 -0.412D-02
+ n1qn3: stopping criterion on g:  0.32255D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.56905D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.218D+00
+     Oren-Spedicato factor (not used) =  0.229D-02     diagonal: average value =  0.316D-02
+ n1qn3: descent direction d: angle(-g,d) =  28.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
+ n1qn3: line search
+     nlis0       fpn=-0.127D-01 d2= 0.37D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.236D-01  0.889D-01
+     nlis0     0.308D+00 -0.202D-02  0.129D-02
+ n1qn3: stopping criterion on g:  0.59635D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.92424D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.995D+00
+     Oren-Spedicato factor (not used) =  0.211D-02     diagonal: average value =  0.320D-02
+ n1qn3: descent direction d: angle(-g,d) =  86.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
+ n1qn3: line search
+     nlis0       fpn=-0.270D-02 d2= 0.13D-02  tmin= 0.68D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.193D-02 -0.124D-02
+ n1qn3: stopping criterion on g:  0.26088D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.61219D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.107D+01
+     Oren-Spedicato factor (not used) =  0.240D-02     diagonal: average value =  0.346D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
+ n1qn3: line search
+     nlis0       fpn=-0.106D-02 d2= 0.11D-02  tmin= 0.73D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.680D-03 -0.342D-03
+ n1qn3: stopping criterion on g:  0.93209D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.93175D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.195D+01
+     Oren-Spedicato factor (not used) =  0.472D-02     diagonal: average value =  0.669D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
+ n1qn3: line search
+     nlis0       fpn=-0.165D-03 d2= 0.22D-03  tmin= 0.17D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.994D-04 -0.378D-04
+ n1qn3: stopping criterion on g:  0.13605D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43415D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.797D+00
+     Oren-Spedicato factor (not used) =  0.369D-02     diagonal: average value =  0.540D-02
+ n1qn3: descent direction d: angle(-g,d) =  86.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
+ n1qn3: line search
+     nlis0       fpn=-0.128D-04 d2= 0.34D-04  tmin= 0.43D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.592D-05  0.983D-06
+ n1qn3: stopping criterion on g:  0.15294D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.39644D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.278D+03
+     Oren-Spedicato factor (not used) =  0.637D+00     diagonal: average value =  0.256D+01
+ n1qn3: descent direction d: angle(-g,d) =  77.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
+ n1qn3: line search
+     nlis0       fpn=-0.241D-05 d2= 0.96D-07  tmin= 0.92D-12 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.168D-06  0.208D-05
+ n1qn3: stopping criterion on g:  0.13388D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.53137D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.794D-03
+     Oren-Spedicato factor (not used) =  0.101D-02     diagonal: average value =  0.159D-02
+ n1qn3: descent direction d: angle(-g,d) =  26.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
+ n1qn3: line search
+     nlis0       fpn=-0.970D-06 d2= 0.12D-08  tmin= 0.80D-11 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.485D-06 -0.183D-09
+ n1qn3: stopping criterion on g:  0.27051D-07
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.11218D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.866D+00
+     Oren-Spedicato factor (not used) =  0.999D-03     diagonal: average value =  0.137D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
+ n1qn3: line search
+     nlis0       fpn=-0.143D-11 d2= 0.35D-11  tmin= 0.13D-09 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.734D-12 -0.116D-13
+ n1qn3: stopping criterion on g:  0.61631D-10
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.53526D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.399D+02
+     Oren-Spedicato factor (not used) =  0.356D-01     diagonal: average value =  0.717D-01
+ n1qn3: descent direction d: angle(-g,d) =  64.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
+ n1qn3: line search
+     nlis0       fpn=-0.991D-16 d2= 0.25D-15  tmin= 0.16D-07 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.152D-15 -0.112D-20
+ n1qn3: stopping criterion on g:  0.92861D-13
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.84647D-02
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.571D+01
+     Oren-Spedicato factor (not used) =  0.482D+00     diagonal: average value =  0.356D+00
+ n1qn3: descent direction d: angle(-g,d) =  80.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
+ n1qn3: line search
+     nlis0       fpn=-0.474D-24 d2= 0.19D-25  tmin= 0.17D-02 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.599D-21 -0.107D-26
+ n1qn3: stopping criterion on g:  0.21365D-15
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.86648D-05
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.352D-02
+     Oren-Spedicato factor (not used) =  0.102D-02     diagonal: average value =  0.198D-02
+ n1qn3: descent direction d: angle(-g,d) =  29.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
+ n1qn3: line search
+     nlis0       fpn=-0.248D-29 d2= 0.33D-32  tmin= 0.39D+01 tmax= 0.10D+21
+     nlis0                                      0.387D+01  0.326D-23  0.743D-29
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   5
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+ n1qn3a: descent direction -g: precon =  0.369D-04
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
+ n1qn3: line search
+     nlis0       fpn=-0.200D+01 d2= 0.74D-04  tmin= 0.28D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.194D+01 -0.189D+01
+     nlis0     0.100D+02 -0.147D+02 -0.973D+00
+ n1qn3: stopping criterion on g:  0.48674D+00
+ n1qn3: matrix update:
+     Oren-Spedicato factor (not used) =  0.718D-03     diagonal: average value =  0.718D-03
+ n1qn3: descent direction d: angle(-g,d) =   0.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
+ n1qn3: line search
+     nlis0       fpn=-0.924D+01 d2= 0.66D-02  tmin= 0.30D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.520D+01 -0.138D+01
+ n1qn3: stopping criterion on g:  0.73079D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.94908D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.118D+01
+     Oren-Spedicato factor (not used) =  0.846D-03     diagonal: average value =  0.846D-03
+ n1qn3: descent direction d: angle(-g,d) =   0.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
+ n1qn3: line search
+     nlis0       fpn=-0.245D+00 d2= 0.21D-03  tmin= 0.16D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.135D+00 -0.256D-01
+ n1qn3: stopping criterion on g:  0.10270D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.17647D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.111D+01
+     Oren-Spedicato factor (not used) =  0.935D-03     diagonal: average value =  0.942D-03
+ n1qn3: descent direction d: angle(-g,d) =   1.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
+ n1qn3: line search
+     nlis0       fpn=-0.581D-02 d2= 0.59D-05  tmin= 0.98D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.435D-02 -0.290D-02
+ n1qn3: stopping criterion on g:  0.76379D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.16885D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.969D+00
+     Oren-Spedicato factor (not used) =  0.952D-03     diagonal: average value =  0.915D-03
+ n1qn3: descent direction d: angle(-g,d) =  23.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
+ n1qn3: line search
+     nlis0       fpn=-0.833D-02 d2= 0.26D-04  tmin= 0.61D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.690D-02 -0.549D-02
+ n1qn3: stopping criterion on g:  0.11072D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.21032D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.985D+00
+     Oren-Spedicato factor (not used) =  0.954D-03     diagonal: average value =  0.102D-02
+ n1qn3: descent direction d: angle(-g,d) =  56.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
+ n1qn3: line search
+     nlis0       fpn=-0.278D-01 d2= 0.37D-03  tmin= 0.14D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.228D-01 -0.180D-01
+ n1qn3: stopping criterion on g:  0.23226D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.37768D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.113D+01
+     Oren-Spedicato factor (not used) =  0.941D-03     diagonal: average value =  0.127D-02
+ n1qn3: descent direction d: angle(-g,d) =  76.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
+ n1qn3: line search
+     nlis0       fpn=-0.828D-01 d2= 0.40D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.713D-01 -0.624D-01
+ n1qn3: stopping criterion on g:  0.44894D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.32918D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.817D+00
+     Oren-Spedicato factor (not used) =  0.762D-03     diagonal: average value =  0.106D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
+ n1qn3: line search
+     nlis0       fpn=-0.405D+00 d2= 0.11D+00  tmin= 0.79D-15 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.531D+00 -0.721D+00
+     nlis0                                      0.100D+02  0.500D+03  0.222D+03
+     nlis0                                      0.398D+01  0.333D+01  0.774D+01
+     nlis0     0.249D+01 -0.140D+01  0.190D+00
+ n1qn3: stopping criterion on g:  0.60049D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.12764D+02
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.287D+01
+     Oren-Spedicato factor (not used) =  0.248D-02     diagonal: average value =  0.292D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
+ n1qn3: line search
+     nlis0       fpn=-0.161D+00 d2= 0.57D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.285D+00  0.571D+00
+     nlis0     0.140D+00 -0.107D-01  0.342D-02
+ n1qn3: stopping criterion on g:  0.48319D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.41360D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.102D+01
+     Oren-Spedicato factor (not used) =  0.293D-02     diagonal: average value =  0.285D-02
+ n1qn3: descent direction d: angle(-g,d) =  76.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
+ n1qn3: line search
+     nlis0       fpn=-0.384D-01 d2= 0.21D-03  tmin= 0.21D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.376D-01 -0.369D-01
+     nlis0     0.100D+02 -0.318D+00 -0.257D-01
+ n1qn3: stopping criterion on g:  0.41045D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43358D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.193D+02
+     Oren-Spedicato factor (not used) =  0.403D-01     diagonal: average value =  0.645D-01
+ n1qn3: descent direction d: angle(-g,d) =  75.7 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
+ n1qn3: line search
+     nlis0       fpn=-0.627D+00 d2= 0.71D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.315D+00  0.212D+00
+ n1qn3: stopping criterion on g:  0.60708D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.18312D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.259D+00
+     Oren-Spedicato factor (not used) =  0.251D-01     diagonal: average value =  0.248D-01
+ n1qn3: descent direction d: angle(-g,d) =  15.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
+ n1qn3: line search
+     nlis0       fpn=-0.707D+00 d2= 0.27D-02  tmin= 0.52D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.402D+00 -0.117D+00
+ n1qn3: stopping criterion on g:  0.11240D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.19476D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.190D+00
+     Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.468D-02
+ n1qn3: descent direction d: angle(-g,d) =  33.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
+ n1qn3: line search
+     nlis0       fpn=-0.499D+00 d2= 0.52D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.225D+00  0.427D+00
+ n1qn3: stopping criterion on g:  0.43961D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43931D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.364D+01
+     Oren-Spedicato factor (not used) =  0.102D-01     diagonal: average value =  0.311D-01
+ n1qn3: descent direction d: angle(-g,d) =  56.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
+ n1qn3: line search
+     nlis0       fpn=-0.339D+00 d2= 0.36D-02  tmin= 0.47D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.151D+00  0.252D-01
+ n1qn3: stopping criterion on g:  0.12127D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.26394D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.759D+00
+     Oren-Spedicato factor (not used) =  0.522D-02     diagonal: average value =  0.235D-01
+ n1qn3: descent direction d: angle(-g,d) =  36.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
+ n1qn3: line search
+     nlis0       fpn=-0.207D+00 d2= 0.84D-02  tmin= 0.24D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.186D+00 -0.162D+00
+ n1qn3: stopping criterion on g:  0.11958D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.15231D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.493D+01
+     Oren-Spedicato factor (not used) =  0.996D-01     diagonal: average value =  0.118D+00
+ n1qn3: descent direction d: angle(-g,d) =  47.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
+ n1qn3: line search
+     nlis0       fpn=-0.650D+00 d2= 0.12D+00  tmin= 0.64D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.169D+01  0.738D+01
+     nlis0     0.374D+00 -0.132D+00  0.166D+00
+ n1qn3: stopping criterion on g:  0.29691D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.14137D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.211D+00
+     Oren-Spedicato factor (not used) =  0.105D-01     diagonal: average value =  0.372D-01
+ n1qn3: descent direction d: angle(-g,d) =  53.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
+ n1qn3: line search
+     nlis0       fpn=-0.187D+00 d2= 0.20D-02  tmin= 0.62D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.149D+00 -0.115D+00
+ n1qn3: stopping criterion on g:  0.15621D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.34738D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.123D+01
+     Oren-Spedicato factor (not used) =  0.634D-02     diagonal: average value =  0.386D-01
+ n1qn3: descent direction d: angle(-g,d) =  62.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
+ n1qn3: line search
+     nlis0       fpn=-0.311D+00 d2= 0.34D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.102D+00  0.268D+00
+ n1qn3: stopping criterion on g:  0.42709D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.41295D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.285D+00
+     Oren-Spedicato factor (not used) =  0.108D-01     diagonal: average value =  0.841D-02
+ n1qn3: descent direction d: angle(-g,d) =  17.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
+ n1qn3: line search
+     nlis0       fpn=-0.421D+00 d2= 0.20D-02  tmin= 0.51D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.181D+00  0.576D-01
+ n1qn3: stopping criterion on g:  0.94247D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.23849D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.526D+00
+     Oren-Spedicato factor (not used) =  0.345D-02     diagonal: average value =  0.461D-02
+ n1qn3: descent direction d: angle(-g,d) =  57.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
+ n1qn3: line search
+     nlis0       fpn=-0.126D+00 d2= 0.11D-01  tmin= 0.25D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.110D+00 -0.725D-01
+ n1qn3: stopping criterion on g:  0.67714D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.23948D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.105D+01
+     Oren-Spedicato factor (not used) =  0.489D-02     diagonal: average value =  0.547D-02
+ n1qn3: descent direction d: angle(-g,d) =  59.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
+ n1qn3: line search
+     nlis0       fpn=-0.228D+00 d2= 0.83D-01  tmin= 0.95D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.601D+00  0.234D+01
+     nlis0     0.321D+00 -0.373D-01  0.472D-01
+ n1qn3: stopping criterion on g:  0.24287D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.87291D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.807D+00
+     Oren-Spedicato factor (not used) =  0.458D-02     diagonal: average value =  0.431D-02
+ n1qn3: descent direction d: angle(-g,d) =  73.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
+ n1qn3: line search
+     nlis0       fpn=-0.554D-01 d2= 0.12D-02  tmin= 0.80D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.477D-01 -0.406D-01
+ n1qn3: stopping criterion on g:  0.15832D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.36912D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.894D+00
+     Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.392D-02
+ n1qn3: descent direction d: angle(-g,d) =  77.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
+ n1qn3: line search
+     nlis0       fpn=-0.143D+00 d2= 0.30D-01  tmin= 0.17D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.878D-01 -0.226D-01
+ n1qn3: stopping criterion on g:  0.23266D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.51063D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.865D+01
+     Oren-Spedicato factor (not used) =  0.318D-01     diagonal: average value =  0.403D-01
+ n1qn3: descent direction d: angle(-g,d) =  75.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
+ n1qn3: line search
+     nlis0       fpn=-0.128D+00 d2= 0.84D-02  tmin= 0.27D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.619D-01 -0.487D-02
+ n1qn3: stopping criterion on g:  0.10083D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.52582D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.572D-01
+     Oren-Spedicato factor (not used) =  0.208D-02     diagonal: average value =  0.227D-02
+ n1qn3: descent direction d: angle(-g,d) =  78.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
+ n1qn3: line search
+     nlis0       fpn=-0.100D+00 d2= 0.49D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.155D-01  0.179D+00
+ n1qn3: stopping criterion on g:  0.37918D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.24062D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.120D+01
+     Oren-Spedicato factor (not used) =  0.224D-02     diagonal: average value =  0.267D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
+ n1qn3: line search
+     nlis0       fpn=-0.120D+00 d2= 0.12D-01  tmin= 0.26D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.251D-01  0.437D-01
+ n1qn3: stopping criterion on g:  0.43947D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.50600D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.757D+00
+     Oren-Spedicato factor (not used) =  0.170D-02     diagonal: average value =  0.201D-02
+ n1qn3: descent direction d: angle(-g,d) =  72.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
+ n1qn3: line search
+     nlis0       fpn=-0.208D-01 d2= 0.45D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.183D-01 -0.148D-01
+ n1qn3: stopping criterion on g:  0.21889D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.59772D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.170D+01
+     Oren-Spedicato factor (not used) =  0.290D-02     diagonal: average value =  0.351D-02
+ n1qn3: descent direction d: angle(-g,d) =  63.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
+ n1qn3: line search
+     nlis0       fpn=-0.555D-01 d2= 0.58D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.507D-01  0.260D+00
+     nlis0     0.413D+00 -0.132D-01  0.278D-02
+ n1qn3: stopping criterion on g:  0.14651D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.14924D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.980D+00
+     Oren-Spedicato factor (not used) =  0.274D-02     diagonal: average value =  0.355D-02
+ n1qn3: descent direction d: angle(-g,d) =  84.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
+ n1qn3: line search
+     nlis0       fpn=-0.918D-02 d2= 0.87D-03  tmin= 0.86D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.813D-02 -0.718D-02
+ n1qn3: stopping criterion on g:  0.11306D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.29590D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.112D+01
+     Oren-Spedicato factor (not used) =  0.329D-02     diagonal: average value =  0.403D-02
+ n1qn3: descent direction d: angle(-g,d) =  85.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
+ n1qn3: line search
+     nlis0       fpn=-0.259D-01 d2= 0.13D-01  tmin= 0.22D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.174D-01 -0.109D-01
+ n1qn3: stopping criterion on g:  0.64294D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.38281D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.337D+01
+     Oren-Spedicato factor (not used) =  0.116D-01     diagonal: average value =  0.134D-01
+ n1qn3: descent direction d: angle(-g,d) =  84.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
+ n1qn3: line search
+     nlis0       fpn=-0.949D-02 d2= 0.41D-02  tmin= 0.39D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.648D-02 -0.412D-02
+ n1qn3: stopping criterion on g:  0.32255D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.56905D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.218D+00
+     Oren-Spedicato factor (not used) =  0.229D-02     diagonal: average value =  0.316D-02
+ n1qn3: descent direction d: angle(-g,d) =  28.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
+ n1qn3: line search
+     nlis0       fpn=-0.127D-01 d2= 0.37D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.236D-01  0.889D-01
+     nlis0     0.308D+00 -0.202D-02  0.129D-02
+ n1qn3: stopping criterion on g:  0.59635D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.92424D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.995D+00
+     Oren-Spedicato factor (not used) =  0.211D-02     diagonal: average value =  0.320D-02
+ n1qn3: descent direction d: angle(-g,d) =  86.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
+ n1qn3: line search
+     nlis0       fpn=-0.270D-02 d2= 0.13D-02  tmin= 0.68D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.193D-02 -0.124D-02
+ n1qn3: stopping criterion on g:  0.26088D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.61219D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.107D+01
+     Oren-Spedicato factor (not used) =  0.240D-02     diagonal: average value =  0.346D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
+ n1qn3: line search
+     nlis0       fpn=-0.106D-02 d2= 0.11D-02  tmin= 0.73D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.680D-03 -0.342D-03
+ n1qn3: stopping criterion on g:  0.93209D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.93175D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.195D+01
+     Oren-Spedicato factor (not used) =  0.472D-02     diagonal: average value =  0.669D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
+ n1qn3: line search
+     nlis0       fpn=-0.165D-03 d2= 0.22D-03  tmin= 0.17D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.994D-04 -0.378D-04
+ n1qn3: stopping criterion on g:  0.13605D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43415D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.797D+00
+     Oren-Spedicato factor (not used) =  0.369D-02     diagonal: average value =  0.540D-02
+ n1qn3: descent direction d: angle(-g,d) =  86.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
+ n1qn3: line search
+     nlis0       fpn=-0.128D-04 d2= 0.34D-04  tmin= 0.43D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.592D-05  0.983D-06
+ n1qn3: stopping criterion on g:  0.15294D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.39644D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.278D+03
+     Oren-Spedicato factor (not used) =  0.637D+00     diagonal: average value =  0.256D+01
+ n1qn3: descent direction d: angle(-g,d) =  77.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
+ n1qn3: line search
+     nlis0       fpn=-0.241D-05 d2= 0.96D-07  tmin= 0.92D-12 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.168D-06  0.208D-05
+ n1qn3: stopping criterion on g:  0.13388D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.53137D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.794D-03
+     Oren-Spedicato factor (not used) =  0.101D-02     diagonal: average value =  0.159D-02
+ n1qn3: descent direction d: angle(-g,d) =  26.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
+ n1qn3: line search
+     nlis0       fpn=-0.970D-06 d2= 0.12D-08  tmin= 0.80D-11 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.485D-06 -0.183D-09
+ n1qn3: stopping criterion on g:  0.27051D-07
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.11218D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.866D+00
+     Oren-Spedicato factor (not used) =  0.999D-03     diagonal: average value =  0.137D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
+ n1qn3: line search
+     nlis0       fpn=-0.143D-11 d2= 0.35D-11  tmin= 0.13D-09 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.734D-12 -0.116D-13
+ n1qn3: stopping criterion on g:  0.61631D-10
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.53526D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.399D+02
+     Oren-Spedicato factor (not used) =  0.356D-01     diagonal: average value =  0.717D-01
+ n1qn3: descent direction d: angle(-g,d) =  64.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
+ n1qn3: line search
+     nlis0       fpn=-0.991D-16 d2= 0.25D-15  tmin= 0.16D-07 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.152D-15 -0.112D-20
+ n1qn3: stopping criterion on g:  0.92861D-13
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.84647D-02
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.571D+01
+     Oren-Spedicato factor (not used) =  0.482D+00     diagonal: average value =  0.356D+00
+ n1qn3: descent direction d: angle(-g,d) =  80.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
+ n1qn3: line search
+     nlis0       fpn=-0.474D-24 d2= 0.19D-25  tmin= 0.17D-02 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.599D-21 -0.107D-26
+ n1qn3: stopping criterion on g:  0.21365D-15
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.86648D-05
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.352D-02
+     Oren-Spedicato factor (not used) =  0.102D-02     diagonal: average value =  0.198D-02
+ n1qn3: descent direction d: angle(-g,d) =  29.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
+ n1qn3: line search
+     nlis0       fpn=-0.248D-29 d2= 0.33D-32  tmin= 0.39D+01 tmax= 0.10D+21
+     nlis0                                      0.387D+01  0.326D-23  0.743D-29
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+
diff --git a/scilab/modules/optimization/tests/nonreg_tests/bug_9690.win.dia.ref b/scilab/modules/optimization/tests/nonreg_tests/bug_9690.win.dia.ref
new file mode 100644 (file)
index 0000000..3740471
--- /dev/null
@@ -0,0 +1,1407 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2013 - Scilab Enterprises - Paul Bignier
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+//
+// <-- CLI SHELL MODE -->
+//
+// <-- ENGLISH IMPOSED -->
+//
+// <-- Non-regression test for bug 9690 -->
+//
+// <-- Bugzilla URL -->
+// http://bugzilla.scilab.org/show_bug.cgi?id=9690
+//
+// <-- Short Description -->
+// optim(): option "imp"=5 could crash Scilab
+//
+function f = rosenbrock(x)
+    f = 100.0 * (x(2)-x(1)^2)^2 + (1-x(1))^2;
+endfunction
+function [f, g, ind] = rosenbrockCost(x, ind)
+    if ((ind == 1) | (ind == 4)) then
+        f = rosenbrock ( x );
+    end
+    if ((ind == 1) | (ind == 4)) then
+        g = derivative ( rosenbrock , x(:) );
+    end
+endfunction
+x0 = [-1.2 1.0];
+lines(0);
+for impval = 1:5
+    [ fopt , xopt ] = optim ( rosenbrockCost , x0 , "gc" , imp=impval);
+end
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   1
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   2
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+ n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
+ n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
+ n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
+ n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
+ n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
+ n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
+ n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
+ n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
+ n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
+ n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
+ n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
+ n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
+ n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
+ n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
+ n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
+ n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
+ n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
+ n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
+ n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
+ n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
+ n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
+ n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
+ n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
+ n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
+ n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
+ n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
+ n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
+ n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
+ n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
+ n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
+ n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
+ n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
+ n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
+ n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
+ n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
+ n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
+ n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
+ n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
+ n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
+ n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
+ n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
+ n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   3
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+ -------------------------------------------------------------------------------
+
+
+ n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
+     nlis0       fpn=-0.200D+01 d2= 0.74D-04  tmin= 0.28D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.194D+01 -0.189D+01
+     nlis0     0.100D+02 -0.147D+02 -0.973D+00
+
+ n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
+     nlis0       fpn=-0.924D+01 d2= 0.66D-02  tmin= 0.30D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.520D+01 -0.138D+01
+
+ n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
+     nlis0       fpn=-0.245D+00 d2= 0.21D-03  tmin= 0.16D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.135D+00 -0.256D-01
+
+ n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
+     nlis0       fpn=-0.581D-02 d2= 0.59D-05  tmin= 0.98D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.435D-02 -0.290D-02
+
+ n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
+     nlis0       fpn=-0.833D-02 d2= 0.26D-04  tmin= 0.61D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.690D-02 -0.549D-02
+
+ n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
+     nlis0       fpn=-0.278D-01 d2= 0.37D-03  tmin= 0.14D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.228D-01 -0.180D-01
+
+ n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
+     nlis0       fpn=-0.828D-01 d2= 0.40D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.713D-01 -0.624D-01
+
+ n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
+     nlis0       fpn=-0.405D+00 d2= 0.11D+00  tmin= 0.79D-15 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.531D+00 -0.721D+00
+     nlis0                                      0.100D+02  0.500D+03  0.222D+03
+     nlis0                                      0.398D+01  0.333D+01  0.774D+01
+     nlis0     0.249D+01 -0.140D+01  0.190D+00
+
+ n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
+     nlis0       fpn=-0.161D+00 d2= 0.57D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.285D+00  0.571D+00
+     nlis0     0.140D+00 -0.107D-01  0.342D-02
+
+ n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
+     nlis0       fpn=-0.384D-01 d2= 0.21D-03  tmin= 0.21D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.376D-01 -0.369D-01
+     nlis0     0.100D+02 -0.318D+00 -0.257D-01
+
+ n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
+     nlis0       fpn=-0.627D+00 d2= 0.71D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.315D+00  0.212D+00
+
+ n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
+     nlis0       fpn=-0.707D+00 d2= 0.27D-02  tmin= 0.52D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.402D+00 -0.117D+00
+
+ n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
+     nlis0       fpn=-0.499D+00 d2= 0.52D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.225D+00  0.427D+00
+
+ n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
+     nlis0       fpn=-0.339D+00 d2= 0.36D-02  tmin= 0.47D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.151D+00  0.252D-01
+
+ n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
+     nlis0       fpn=-0.207D+00 d2= 0.84D-02  tmin= 0.24D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.186D+00 -0.162D+00
+
+ n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
+     nlis0       fpn=-0.650D+00 d2= 0.12D+00  tmin= 0.64D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.169D+01  0.738D+01
+     nlis0     0.374D+00 -0.132D+00  0.166D+00
+
+ n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
+     nlis0       fpn=-0.187D+00 d2= 0.20D-02  tmin= 0.62D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.149D+00 -0.115D+00
+
+ n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
+     nlis0       fpn=-0.311D+00 d2= 0.34D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.102D+00  0.268D+00
+
+ n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
+     nlis0       fpn=-0.421D+00 d2= 0.20D-02  tmin= 0.51D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.181D+00  0.576D-01
+
+ n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
+     nlis0       fpn=-0.126D+00 d2= 0.11D-01  tmin= 0.25D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.110D+00 -0.725D-01
+
+ n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
+     nlis0       fpn=-0.228D+00 d2= 0.83D-01  tmin= 0.95D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.601D+00  0.234D+01
+     nlis0     0.321D+00 -0.373D-01  0.472D-01
+
+ n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
+     nlis0       fpn=-0.554D-01 d2= 0.12D-02  tmin= 0.80D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.477D-01 -0.406D-01
+
+ n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
+     nlis0       fpn=-0.143D+00 d2= 0.30D-01  tmin= 0.17D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.878D-01 -0.226D-01
+
+ n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
+     nlis0       fpn=-0.128D+00 d2= 0.84D-02  tmin= 0.27D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.619D-01 -0.487D-02
+
+ n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
+     nlis0       fpn=-0.100D+00 d2= 0.49D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.155D-01  0.179D+00
+
+ n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
+     nlis0       fpn=-0.120D+00 d2= 0.12D-01  tmin= 0.26D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.251D-01  0.437D-01
+
+ n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
+     nlis0       fpn=-0.208D-01 d2= 0.45D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.183D-01 -0.148D-01
+
+ n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
+     nlis0       fpn=-0.555D-01 d2= 0.58D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.507D-01  0.260D+00
+     nlis0     0.413D+00 -0.132D-01  0.278D-02
+
+ n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
+     nlis0       fpn=-0.918D-02 d2= 0.87D-03  tmin= 0.86D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.813D-02 -0.718D-02
+
+ n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
+     nlis0       fpn=-0.259D-01 d2= 0.13D-01  tmin= 0.22D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.174D-01 -0.109D-01
+
+ n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
+     nlis0       fpn=-0.949D-02 d2= 0.41D-02  tmin= 0.39D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.648D-02 -0.412D-02
+
+ n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
+     nlis0       fpn=-0.127D-01 d2= 0.37D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.236D-01  0.889D-01
+     nlis0     0.308D+00 -0.202D-02  0.129D-02
+
+ n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
+     nlis0       fpn=-0.270D-02 d2= 0.13D-02  tmin= 0.68D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.193D-02 -0.124D-02
+
+ n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
+     nlis0       fpn=-0.106D-02 d2= 0.11D-02  tmin= 0.73D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.680D-03 -0.342D-03
+
+ n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
+     nlis0       fpn=-0.165D-03 d2= 0.22D-03  tmin= 0.17D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.994D-04 -0.378D-04
+
+ n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
+     nlis0       fpn=-0.128D-04 d2= 0.34D-04  tmin= 0.43D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.592D-05  0.983D-06
+
+ n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
+     nlis0       fpn=-0.241D-05 d2= 0.96D-07  tmin= 0.92D-12 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.168D-06  0.208D-05
+
+ n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
+     nlis0       fpn=-0.970D-06 d2= 0.12D-08  tmin= 0.80D-11 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.485D-06 -0.183D-09
+
+ n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
+     nlis0       fpn=-0.143D-11 d2= 0.35D-11  tmin= 0.13D-09 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.734D-12 -0.116D-13
+
+ n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
+     nlis0       fpn=-0.991D-16 d2= 0.25D-15  tmin= 0.16D-07 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.152D-15 -0.112D-20
+
+ n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
+     nlis0       fpn=-0.474D-24 d2= 0.19D-25  tmin= 0.17D-02 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.599D-21 -0.107D-26
+
+ n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
+     nlis0       fpn=-0.248D-29 d2= 0.33D-32  tmin= 0.39D+01 tmax= 0.10D+21
+     nlis0                                      0.387D+01  0.326D-23  0.743D-29
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   4
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+ n1qn3a: descent direction -g: precon =  0.369D-04
+ -------------------------------------------------------------------------------
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
+ n1qn3: line search
+     nlis0       fpn=-0.200D+01 d2= 0.74D-04  tmin= 0.28D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.194D+01 -0.189D+01
+     nlis0     0.100D+02 -0.147D+02 -0.973D+00
+ n1qn3: stopping criterion on g:  0.48674D+00
+ n1qn3: matrix update:
+     Oren-Spedicato factor (not used) =  0.718D-03     diagonal: average value =  0.718D-03
+ n1qn3: descent direction d: angle(-g,d) =   0.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
+ n1qn3: line search
+     nlis0       fpn=-0.924D+01 d2= 0.66D-02  tmin= 0.30D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.520D+01 -0.138D+01
+ n1qn3: stopping criterion on g:  0.73079D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.94908D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.118D+01
+     Oren-Spedicato factor (not used) =  0.846D-03     diagonal: average value =  0.846D-03
+ n1qn3: descent direction d: angle(-g,d) =   0.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
+ n1qn3: line search
+     nlis0       fpn=-0.245D+00 d2= 0.21D-03  tmin= 0.16D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.135D+00 -0.256D-01
+ n1qn3: stopping criterion on g:  0.10270D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.17647D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.111D+01
+     Oren-Spedicato factor (not used) =  0.935D-03     diagonal: average value =  0.942D-03
+ n1qn3: descent direction d: angle(-g,d) =   1.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
+ n1qn3: line search
+     nlis0       fpn=-0.581D-02 d2= 0.59D-05  tmin= 0.98D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.435D-02 -0.290D-02
+ n1qn3: stopping criterion on g:  0.76379D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.16885D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.969D+00
+     Oren-Spedicato factor (not used) =  0.952D-03     diagonal: average value =  0.915D-03
+ n1qn3: descent direction d: angle(-g,d) =  23.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
+ n1qn3: line search
+     nlis0       fpn=-0.833D-02 d2= 0.26D-04  tmin= 0.61D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.690D-02 -0.549D-02
+ n1qn3: stopping criterion on g:  0.11072D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.21032D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.985D+00
+     Oren-Spedicato factor (not used) =  0.954D-03     diagonal: average value =  0.102D-02
+ n1qn3: descent direction d: angle(-g,d) =  56.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
+ n1qn3: line search
+     nlis0       fpn=-0.278D-01 d2= 0.37D-03  tmin= 0.14D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.228D-01 -0.180D-01
+ n1qn3: stopping criterion on g:  0.23226D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.37768D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.113D+01
+     Oren-Spedicato factor (not used) =  0.941D-03     diagonal: average value =  0.127D-02
+ n1qn3: descent direction d: angle(-g,d) =  76.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
+ n1qn3: line search
+     nlis0       fpn=-0.828D-01 d2= 0.40D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.713D-01 -0.624D-01
+ n1qn3: stopping criterion on g:  0.44894D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.32918D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.817D+00
+     Oren-Spedicato factor (not used) =  0.762D-03     diagonal: average value =  0.106D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
+ n1qn3: line search
+     nlis0       fpn=-0.405D+00 d2= 0.11D+00  tmin= 0.79D-15 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.531D+00 -0.721D+00
+     nlis0                                      0.100D+02  0.500D+03  0.222D+03
+     nlis0                                      0.398D+01  0.333D+01  0.774D+01
+     nlis0     0.249D+01 -0.140D+01  0.190D+00
+ n1qn3: stopping criterion on g:  0.60049D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.12764D+02
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.287D+01
+     Oren-Spedicato factor (not used) =  0.248D-02     diagonal: average value =  0.292D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
+ n1qn3: line search
+     nlis0       fpn=-0.161D+00 d2= 0.57D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.285D+00  0.571D+00
+     nlis0     0.140D+00 -0.107D-01  0.342D-02
+ n1qn3: stopping criterion on g:  0.48319D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.41360D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.102D+01
+     Oren-Spedicato factor (not used) =  0.293D-02     diagonal: average value =  0.285D-02
+ n1qn3: descent direction d: angle(-g,d) =  76.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
+ n1qn3: line search
+     nlis0       fpn=-0.384D-01 d2= 0.21D-03  tmin= 0.21D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.376D-01 -0.369D-01
+     nlis0     0.100D+02 -0.318D+00 -0.257D-01
+ n1qn3: stopping criterion on g:  0.41045D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43358D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.193D+02
+     Oren-Spedicato factor (not used) =  0.403D-01     diagonal: average value =  0.645D-01
+ n1qn3: descent direction d: angle(-g,d) =  75.7 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
+ n1qn3: line search
+     nlis0       fpn=-0.627D+00 d2= 0.71D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.315D+00  0.212D+00
+ n1qn3: stopping criterion on g:  0.60708D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.18312D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.259D+00
+     Oren-Spedicato factor (not used) =  0.251D-01     diagonal: average value =  0.248D-01
+ n1qn3: descent direction d: angle(-g,d) =  15.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
+ n1qn3: line search
+     nlis0       fpn=-0.707D+00 d2= 0.27D-02  tmin= 0.52D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.402D+00 -0.117D+00
+ n1qn3: stopping criterion on g:  0.11240D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.19476D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.190D+00
+     Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.468D-02
+ n1qn3: descent direction d: angle(-g,d) =  33.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
+ n1qn3: line search
+     nlis0       fpn=-0.499D+00 d2= 0.52D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.225D+00  0.427D+00
+ n1qn3: stopping criterion on g:  0.43961D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43931D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.364D+01
+     Oren-Spedicato factor (not used) =  0.102D-01     diagonal: average value =  0.311D-01
+ n1qn3: descent direction d: angle(-g,d) =  56.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
+ n1qn3: line search
+     nlis0       fpn=-0.339D+00 d2= 0.36D-02  tmin= 0.47D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.151D+00  0.252D-01
+ n1qn3: stopping criterion on g:  0.12127D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.26394D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.759D+00
+     Oren-Spedicato factor (not used) =  0.522D-02     diagonal: average value =  0.235D-01
+ n1qn3: descent direction d: angle(-g,d) =  36.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
+ n1qn3: line search
+     nlis0       fpn=-0.207D+00 d2= 0.84D-02  tmin= 0.24D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.186D+00 -0.162D+00
+ n1qn3: stopping criterion on g:  0.11958D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.15231D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.493D+01
+     Oren-Spedicato factor (not used) =  0.996D-01     diagonal: average value =  0.118D+00
+ n1qn3: descent direction d: angle(-g,d) =  47.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
+ n1qn3: line search
+     nlis0       fpn=-0.650D+00 d2= 0.12D+00  tmin= 0.64D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.169D+01  0.738D+01
+     nlis0     0.374D+00 -0.132D+00  0.166D+00
+ n1qn3: stopping criterion on g:  0.29691D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.14137D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.211D+00
+     Oren-Spedicato factor (not used) =  0.105D-01     diagonal: average value =  0.372D-01
+ n1qn3: descent direction d: angle(-g,d) =  53.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
+ n1qn3: line search
+     nlis0       fpn=-0.187D+00 d2= 0.20D-02  tmin= 0.62D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.149D+00 -0.115D+00
+ n1qn3: stopping criterion on g:  0.15621D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.34738D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.123D+01
+     Oren-Spedicato factor (not used) =  0.634D-02     diagonal: average value =  0.386D-01
+ n1qn3: descent direction d: angle(-g,d) =  62.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
+ n1qn3: line search
+     nlis0       fpn=-0.311D+00 d2= 0.34D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.102D+00  0.268D+00
+ n1qn3: stopping criterion on g:  0.42709D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.41295D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.285D+00
+     Oren-Spedicato factor (not used) =  0.108D-01     diagonal: average value =  0.841D-02
+ n1qn3: descent direction d: angle(-g,d) =  17.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
+ n1qn3: line search
+     nlis0       fpn=-0.421D+00 d2= 0.20D-02  tmin= 0.51D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.181D+00  0.576D-01
+ n1qn3: stopping criterion on g:  0.94247D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.23849D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.526D+00
+     Oren-Spedicato factor (not used) =  0.345D-02     diagonal: average value =  0.461D-02
+ n1qn3: descent direction d: angle(-g,d) =  57.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
+ n1qn3: line search
+     nlis0       fpn=-0.126D+00 d2= 0.11D-01  tmin= 0.25D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.110D+00 -0.725D-01
+ n1qn3: stopping criterion on g:  0.67714D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.23948D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.105D+01
+     Oren-Spedicato factor (not used) =  0.489D-02     diagonal: average value =  0.547D-02
+ n1qn3: descent direction d: angle(-g,d) =  59.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
+ n1qn3: line search
+     nlis0       fpn=-0.228D+00 d2= 0.83D-01  tmin= 0.95D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.601D+00  0.234D+01
+     nlis0     0.321D+00 -0.373D-01  0.472D-01
+ n1qn3: stopping criterion on g:  0.24287D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.87291D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.807D+00
+     Oren-Spedicato factor (not used) =  0.458D-02     diagonal: average value =  0.431D-02
+ n1qn3: descent direction d: angle(-g,d) =  73.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
+ n1qn3: line search
+     nlis0       fpn=-0.554D-01 d2= 0.12D-02  tmin= 0.80D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.477D-01 -0.406D-01
+ n1qn3: stopping criterion on g:  0.15832D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.36912D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.894D+00
+     Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.392D-02
+ n1qn3: descent direction d: angle(-g,d) =  77.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
+ n1qn3: line search
+     nlis0       fpn=-0.143D+00 d2= 0.30D-01  tmin= 0.17D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.878D-01 -0.226D-01
+ n1qn3: stopping criterion on g:  0.23266D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.51063D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.865D+01
+     Oren-Spedicato factor (not used) =  0.318D-01     diagonal: average value =  0.403D-01
+ n1qn3: descent direction d: angle(-g,d) =  75.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
+ n1qn3: line search
+     nlis0       fpn=-0.128D+00 d2= 0.84D-02  tmin= 0.27D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.619D-01 -0.487D-02
+ n1qn3: stopping criterion on g:  0.10083D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.52582D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.572D-01
+     Oren-Spedicato factor (not used) =  0.208D-02     diagonal: average value =  0.227D-02
+ n1qn3: descent direction d: angle(-g,d) =  78.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
+ n1qn3: line search
+     nlis0       fpn=-0.100D+00 d2= 0.49D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.155D-01  0.179D+00
+ n1qn3: stopping criterion on g:  0.37918D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.24062D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.120D+01
+     Oren-Spedicato factor (not used) =  0.224D-02     diagonal: average value =  0.267D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
+ n1qn3: line search
+     nlis0       fpn=-0.120D+00 d2= 0.12D-01  tmin= 0.26D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.251D-01  0.437D-01
+ n1qn3: stopping criterion on g:  0.43947D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.50600D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.757D+00
+     Oren-Spedicato factor (not used) =  0.170D-02     diagonal: average value =  0.201D-02
+ n1qn3: descent direction d: angle(-g,d) =  72.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
+ n1qn3: line search
+     nlis0       fpn=-0.208D-01 d2= 0.45D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.183D-01 -0.148D-01
+ n1qn3: stopping criterion on g:  0.21889D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.59772D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.170D+01
+     Oren-Spedicato factor (not used) =  0.290D-02     diagonal: average value =  0.351D-02
+ n1qn3: descent direction d: angle(-g,d) =  63.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
+ n1qn3: line search
+     nlis0       fpn=-0.555D-01 d2= 0.58D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.507D-01  0.260D+00
+     nlis0     0.413D+00 -0.132D-01  0.278D-02
+ n1qn3: stopping criterion on g:  0.14651D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.14924D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.980D+00
+     Oren-Spedicato factor (not used) =  0.274D-02     diagonal: average value =  0.355D-02
+ n1qn3: descent direction d: angle(-g,d) =  84.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
+ n1qn3: line search
+     nlis0       fpn=-0.918D-02 d2= 0.87D-03  tmin= 0.86D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.813D-02 -0.718D-02
+ n1qn3: stopping criterion on g:  0.11306D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.29590D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.112D+01
+     Oren-Spedicato factor (not used) =  0.329D-02     diagonal: average value =  0.403D-02
+ n1qn3: descent direction d: angle(-g,d) =  85.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
+ n1qn3: line search
+     nlis0       fpn=-0.259D-01 d2= 0.13D-01  tmin= 0.22D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.174D-01 -0.109D-01
+ n1qn3: stopping criterion on g:  0.64294D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.38281D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.337D+01
+     Oren-Spedicato factor (not used) =  0.116D-01     diagonal: average value =  0.134D-01
+ n1qn3: descent direction d: angle(-g,d) =  84.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
+ n1qn3: line search
+     nlis0       fpn=-0.949D-02 d2= 0.41D-02  tmin= 0.39D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.648D-02 -0.412D-02
+ n1qn3: stopping criterion on g:  0.32255D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.56905D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.218D+00
+     Oren-Spedicato factor (not used) =  0.229D-02     diagonal: average value =  0.316D-02
+ n1qn3: descent direction d: angle(-g,d) =  28.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
+ n1qn3: line search
+     nlis0       fpn=-0.127D-01 d2= 0.37D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.236D-01  0.889D-01
+     nlis0     0.308D+00 -0.202D-02  0.129D-02
+ n1qn3: stopping criterion on g:  0.59635D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.92424D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.995D+00
+     Oren-Spedicato factor (not used) =  0.211D-02     diagonal: average value =  0.320D-02
+ n1qn3: descent direction d: angle(-g,d) =  86.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
+ n1qn3: line search
+     nlis0       fpn=-0.270D-02 d2= 0.13D-02  tmin= 0.68D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.193D-02 -0.124D-02
+ n1qn3: stopping criterion on g:  0.26088D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.61219D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.107D+01
+     Oren-Spedicato factor (not used) =  0.240D-02     diagonal: average value =  0.346D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
+ n1qn3: line search
+     nlis0       fpn=-0.106D-02 d2= 0.11D-02  tmin= 0.73D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.680D-03 -0.342D-03
+ n1qn3: stopping criterion on g:  0.93209D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.93175D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.195D+01
+     Oren-Spedicato factor (not used) =  0.472D-02     diagonal: average value =  0.669D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
+ n1qn3: line search
+     nlis0       fpn=-0.165D-03 d2= 0.22D-03  tmin= 0.17D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.994D-04 -0.378D-04
+ n1qn3: stopping criterion on g:  0.13605D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43415D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.797D+00
+     Oren-Spedicato factor (not used) =  0.369D-02     diagonal: average value =  0.540D-02
+ n1qn3: descent direction d: angle(-g,d) =  86.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
+ n1qn3: line search
+     nlis0       fpn=-0.128D-04 d2= 0.34D-04  tmin= 0.43D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.592D-05  0.983D-06
+ n1qn3: stopping criterion on g:  0.15294D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.39644D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.278D+03
+     Oren-Spedicato factor (not used) =  0.637D+00     diagonal: average value =  0.256D+01
+ n1qn3: descent direction d: angle(-g,d) =  77.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
+ n1qn3: line search
+     nlis0       fpn=-0.241D-05 d2= 0.96D-07  tmin= 0.92D-12 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.168D-06  0.208D-05
+ n1qn3: stopping criterion on g:  0.13388D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.53137D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.794D-03
+     Oren-Spedicato factor (not used) =  0.101D-02     diagonal: average value =  0.159D-02
+ n1qn3: descent direction d: angle(-g,d) =  26.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
+ n1qn3: line search
+     nlis0       fpn=-0.970D-06 d2= 0.12D-08  tmin= 0.80D-11 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.485D-06 -0.183D-09
+ n1qn3: stopping criterion on g:  0.27051D-07
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.11218D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.866D+00
+     Oren-Spedicato factor (not used) =  0.999D-03     diagonal: average value =  0.137D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
+ n1qn3: line search
+     nlis0       fpn=-0.143D-11 d2= 0.35D-11  tmin= 0.13D-09 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.734D-12 -0.116D-13
+ n1qn3: stopping criterion on g:  0.61631D-10
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.53526D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.399D+02
+     Oren-Spedicato factor (not used) =  0.356D-01     diagonal: average value =  0.717D-01
+ n1qn3: descent direction d: angle(-g,d) =  64.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
+ n1qn3: line search
+     nlis0       fpn=-0.991D-16 d2= 0.25D-15  tmin= 0.16D-07 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.152D-15 -0.112D-20
+ n1qn3: stopping criterion on g:  0.92861D-13
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.84647D-02
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.571D+01
+     Oren-Spedicato factor (not used) =  0.482D+00     diagonal: average value =  0.356D+00
+ n1qn3: descent direction d: angle(-g,d) =  80.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
+ n1qn3: line search
+     nlis0       fpn=-0.474D-24 d2= 0.19D-25  tmin= 0.17D-02 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.599D-21 -0.107D-26
+ n1qn3: stopping criterion on g:  0.21365D-15
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.86648D-05
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.352D-02
+     Oren-Spedicato factor (not used) =  0.102D-02     diagonal: average value =  0.198D-02
+ n1qn3: descent direction d: angle(-g,d) =  29.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
+ n1qn3: line search
+     nlis0       fpn=-0.248D-29 d2= 0.33D-32  tmin= 0.39D+01 tmax= 0.10D+21
+     nlis0                                      0.387D+01  0.326D-23  0.743D-29
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+
+ n1qn3: entry point
+     dimension of the problem (n):     2
+     absolute precision on x (dxmin): 0.22D-15
+     expected decrease for f (df1): 0.10D+01
+     relative precision on g (epsg): 0.95D-18
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   5
+     allocated memory (nrz) :     58
+     used memory :                58
+     number of updates :          10
+
+     f         =  0.24200000D+02
+     norm of g =  0.23286769D+03
+ n1qn3a: descent direction -g: precon =  0.369D-04
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   1, simul   1, f= 0.24200000D+02, h'(0)=-0.20000D+01
+ n1qn3: line search
+     nlis0       fpn=-0.200D+01 d2= 0.74D-04  tmin= 0.28D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.194D+01 -0.189D+01
+     nlis0     0.100D+02 -0.147D+02 -0.973D+00
+ n1qn3: stopping criterion on g:  0.48674D+00
+ n1qn3: matrix update:
+     Oren-Spedicato factor (not used) =  0.718D-03     diagonal: average value =  0.718D-03
+ n1qn3: descent direction d: angle(-g,d) =   0.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   2, simul   3, f= 0.94705115D+01, h'(0)=-0.92378D+01
+ n1qn3: line search
+     nlis0       fpn=-0.924D+01 d2= 0.66D-02  tmin= 0.30D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.520D+01 -0.138D+01
+ n1qn3: stopping criterion on g:  0.73079D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.94908D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.118D+01
+     Oren-Spedicato factor (not used) =  0.846D-03     diagonal: average value =  0.846D-03
+ n1qn3: descent direction d: angle(-g,d) =   0.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   3, simul   4, f= 0.42678391D+01, h'(0)=-0.24478D+00
+ n1qn3: line search
+     nlis0       fpn=-0.245D+00 d2= 0.21D-03  tmin= 0.16D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.135D+00 -0.256D-01
+ n1qn3: stopping criterion on g:  0.10270D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.17647D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.111D+01
+     Oren-Spedicato factor (not used) =  0.935D-03     diagonal: average value =  0.942D-03
+ n1qn3: descent direction d: angle(-g,d) =   1.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   4, simul   5, f= 0.41332370D+01, h'(0)=-0.58072D-02
+ n1qn3: line search
+     nlis0       fpn=-0.581D-02 d2= 0.59D-05  tmin= 0.98D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.435D-02 -0.290D-02
+ n1qn3: stopping criterion on g:  0.76379D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.16885D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.969D+00
+     Oren-Spedicato factor (not used) =  0.952D-03     diagonal: average value =  0.915D-03
+ n1qn3: descent direction d: angle(-g,d) =  23.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   5, simul   6, f= 0.41288855D+01, h'(0)=-0.83281D-02
+ n1qn3: line search
+     nlis0       fpn=-0.833D-02 d2= 0.26D-04  tmin= 0.61D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.690D-02 -0.549D-02
+ n1qn3: stopping criterion on g:  0.11072D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.21032D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.985D+00
+     Oren-Spedicato factor (not used) =  0.954D-03     diagonal: average value =  0.102D-02
+ n1qn3: descent direction d: angle(-g,d) =  56.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   6, simul   7, f= 0.41219818D+01, h'(0)=-0.27767D-01
+ n1qn3: line search
+     nlis0       fpn=-0.278D-01 d2= 0.37D-03  tmin= 0.14D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.228D-01 -0.180D-01
+ n1qn3: stopping criterion on g:  0.23226D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.37768D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.113D+01
+     Oren-Spedicato factor (not used) =  0.941D-03     diagonal: average value =  0.127D-02
+ n1qn3: descent direction d: angle(-g,d) =  76.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   7, simul   8, f= 0.40992137D+01, h'(0)=-0.82808D-01
+ n1qn3: line search
+     nlis0       fpn=-0.828D-01 d2= 0.40D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.713D-01 -0.624D-01
+ n1qn3: stopping criterion on g:  0.44894D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.32918D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.817D+00
+     Oren-Spedicato factor (not used) =  0.762D-03     diagonal: average value =  0.106D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   8, simul   9, f= 0.40279192D+01, h'(0)=-0.40499D+00
+ n1qn3: line search
+     nlis0       fpn=-0.405D+00 d2= 0.11D+00  tmin= 0.79D-15 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.531D+00 -0.721D+00
+     nlis0                                      0.100D+02  0.500D+03  0.222D+03
+     nlis0                                      0.398D+01  0.333D+01  0.774D+01
+     nlis0     0.249D+01 -0.140D+01  0.190D+00
+ n1qn3: stopping criterion on g:  0.60049D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.12764D+02
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.287D+01
+     Oren-Spedicato factor (not used) =  0.248D-02     diagonal: average value =  0.292D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter   9, simul  13, f= 0.26255752D+01, h'(0)=-0.16053D+00
+ n1qn3: line search
+     nlis0       fpn=-0.161D+00 d2= 0.57D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.285D+00  0.571D+00
+     nlis0     0.140D+00 -0.107D-01  0.342D-02
+ n1qn3: stopping criterion on g:  0.48319D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.41360D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.102D+01
+     Oren-Spedicato factor (not used) =  0.293D-02     diagonal: average value =  0.285D-02
+ n1qn3: descent direction d: angle(-g,d) =  76.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  10, simul  15, f= 0.26148282D+01, h'(0)=-0.38378D-01
+ n1qn3: line search
+     nlis0       fpn=-0.384D-01 d2= 0.21D-03  tmin= 0.21D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.376D-01 -0.369D-01
+     nlis0     0.100D+02 -0.318D+00 -0.257D-01
+ n1qn3: stopping criterion on g:  0.41045D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43358D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.193D+02
+     Oren-Spedicato factor (not used) =  0.403D-01     diagonal: average value =  0.645D-01
+ n1qn3: descent direction d: angle(-g,d) =  75.7 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  11, simul  17, f= 0.22964498D+01, h'(0)=-0.62657D+00
+ n1qn3: line search
+     nlis0       fpn=-0.627D+00 d2= 0.71D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.315D+00  0.212D+00
+ n1qn3: stopping criterion on g:  0.60708D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.18312D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.259D+00
+     Oren-Spedicato factor (not used) =  0.251D-01     diagonal: average value =  0.248D-01
+ n1qn3: descent direction d: angle(-g,d) =  15.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  12, simul  18, f= 0.19818830D+01, h'(0)=-0.70680D+00
+ n1qn3: line search
+     nlis0       fpn=-0.707D+00 d2= 0.27D-02  tmin= 0.52D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.402D+00 -0.117D+00
+ n1qn3: stopping criterion on g:  0.11240D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.19476D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.190D+00
+     Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.468D-02
+ n1qn3: descent direction d: angle(-g,d) =  33.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  13, simul  19, f= 0.15797116D+01, h'(0)=-0.49943D+00
+ n1qn3: line search
+     nlis0       fpn=-0.499D+00 d2= 0.52D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.225D+00  0.427D+00
+ n1qn3: stopping criterion on g:  0.43961D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43931D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.364D+01
+     Oren-Spedicato factor (not used) =  0.102D-01     diagonal: average value =  0.311D-01
+ n1qn3: descent direction d: angle(-g,d) =  56.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  14, simul  20, f= 0.13542804D+01, h'(0)=-0.33867D+00
+ n1qn3: line search
+     nlis0       fpn=-0.339D+00 d2= 0.36D-02  tmin= 0.47D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.151D+00  0.252D-01
+ n1qn3: stopping criterion on g:  0.12127D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.26394D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.759D+00
+     Oren-Spedicato factor (not used) =  0.522D-02     diagonal: average value =  0.235D-01
+ n1qn3: descent direction d: angle(-g,d) =  36.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  15, simul  21, f= 0.12033041D+01, h'(0)=-0.20678D+00
+ n1qn3: line search
+     nlis0       fpn=-0.207D+00 d2= 0.84D-02  tmin= 0.24D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.186D+00 -0.162D+00
+ n1qn3: stopping criterion on g:  0.11958D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.15231D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.493D+01
+     Oren-Spedicato factor (not used) =  0.996D-01     diagonal: average value =  0.118D+00
+ n1qn3: descent direction d: angle(-g,d) =  47.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  16, simul  22, f= 0.10169710D+01, h'(0)=-0.64958D+00
+ n1qn3: line search
+     nlis0       fpn=-0.650D+00 d2= 0.12D+00  tmin= 0.64D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.169D+01  0.738D+01
+     nlis0     0.374D+00 -0.132D+00  0.166D+00
+ n1qn3: stopping criterion on g:  0.29691D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.14137D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.211D+00
+     Oren-Spedicato factor (not used) =  0.105D-01     diagonal: average value =  0.372D-01
+ n1qn3: descent direction d: angle(-g,d) =  53.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  17, simul  24, f= 0.88479943D+00, h'(0)=-0.18657D+00
+ n1qn3: line search
+     nlis0       fpn=-0.187D+00 d2= 0.20D-02  tmin= 0.62D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.149D+00 -0.115D+00
+ n1qn3: stopping criterion on g:  0.15621D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.34738D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.123D+01
+     Oren-Spedicato factor (not used) =  0.634D-02     diagonal: average value =  0.386D-01
+ n1qn3: descent direction d: angle(-g,d) =  62.6 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  18, simul  25, f= 0.73594993D+00, h'(0)=-0.31085D+00
+ n1qn3: line search
+     nlis0       fpn=-0.311D+00 d2= 0.34D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.102D+00  0.268D+00
+ n1qn3: stopping criterion on g:  0.42709D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.41295D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.285D+00
+     Oren-Spedicato factor (not used) =  0.108D-01     diagonal: average value =  0.841D-02
+ n1qn3: descent direction d: angle(-g,d) =  17.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  19, simul  26, f= 0.63444382D+00, h'(0)=-0.42075D+00
+ n1qn3: line search
+     nlis0       fpn=-0.421D+00 d2= 0.20D-02  tmin= 0.51D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.181D+00  0.576D-01
+ n1qn3: stopping criterion on g:  0.94247D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.23849D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.526D+00
+     Oren-Spedicato factor (not used) =  0.345D-02     diagonal: average value =  0.461D-02
+ n1qn3: descent direction d: angle(-g,d) =  57.2 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  20, simul  27, f= 0.45317916D+00, h'(0)=-0.12594D+00
+ n1qn3: line search
+     nlis0       fpn=-0.126D+00 d2= 0.11D-01  tmin= 0.25D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.110D+00 -0.725D-01
+ n1qn3: stopping criterion on g:  0.67714D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.23948D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.105D+01
+     Oren-Spedicato factor (not used) =  0.489D-02     diagonal: average value =  0.547D-02
+ n1qn3: descent direction d: angle(-g,d) =  59.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  21, simul  28, f= 0.34356278D+00, h'(0)=-0.22838D+00
+ n1qn3: line search
+     nlis0       fpn=-0.228D+00 d2= 0.83D-01  tmin= 0.95D-15 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.601D+00  0.234D+01
+     nlis0     0.321D+00 -0.373D-01  0.472D-01
+ n1qn3: stopping criterion on g:  0.24287D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.87291D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.807D+00
+     Oren-Spedicato factor (not used) =  0.458D-02     diagonal: average value =  0.431D-02
+ n1qn3: descent direction d: angle(-g,d) =  73.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  22, simul  30, f= 0.30624976D+00, h'(0)=-0.55446D-01
+ n1qn3: line search
+     nlis0       fpn=-0.554D-01 d2= 0.12D-02  tmin= 0.80D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.477D-01 -0.406D-01
+ n1qn3: stopping criterion on g:  0.15832D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.36912D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.894D+00
+     Oren-Spedicato factor (not used) =  0.380D-02     diagonal: average value =  0.392D-02
+ n1qn3: descent direction d: angle(-g,d) =  77.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  23, simul  31, f= 0.25851794D+00, h'(0)=-0.14344D+00
+ n1qn3: line search
+     nlis0       fpn=-0.143D+00 d2= 0.30D-01  tmin= 0.17D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.878D-01 -0.226D-01
+ n1qn3: stopping criterion on g:  0.23266D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.51063D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.865D+01
+     Oren-Spedicato factor (not used) =  0.318D-01     diagonal: average value =  0.403D-01
+ n1qn3: descent direction d: angle(-g,d) =  75.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  24, simul  32, f= 0.17073490D+00, h'(0)=-0.12793D+00
+ n1qn3: line search
+     nlis0       fpn=-0.128D+00 d2= 0.84D-02  tmin= 0.27D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.619D-01 -0.487D-02
+ n1qn3: stopping criterion on g:  0.10083D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.52582D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.572D-01
+     Oren-Spedicato factor (not used) =  0.208D-02     diagonal: average value =  0.227D-02
+ n1qn3: descent direction d: angle(-g,d) =  78.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  25, simul  33, f= 0.10886882D+00, h'(0)=-0.10034D+00
+ n1qn3: line search
+     nlis0       fpn=-0.100D+00 d2= 0.49D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.155D-01  0.179D+00
+ n1qn3: stopping criterion on g:  0.37918D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.24062D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.120D+01
+     Oren-Spedicato factor (not used) =  0.224D-02     diagonal: average value =  0.267D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  26, simul  34, f= 0.93394212D-01, h'(0)=-0.12026D+00
+ n1qn3: line search
+     nlis0       fpn=-0.120D+00 d2= 0.12D-01  tmin= 0.26D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.251D-01  0.437D-01
+ n1qn3: stopping criterion on g:  0.43947D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.50600D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.757D+00
+     Oren-Spedicato factor (not used) =  0.170D-02     diagonal: average value =  0.201D-02
+ n1qn3: descent direction d: angle(-g,d) =  72.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  27, simul  35, f= 0.68247217D-01, h'(0)=-0.20762D-01
+ n1qn3: line search
+     nlis0       fpn=-0.208D-01 d2= 0.45D-02  tmin= 0.41D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.183D-01 -0.148D-01
+ n1qn3: stopping criterion on g:  0.21889D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.59772D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.170D+01
+     Oren-Spedicato factor (not used) =  0.290D-02     diagonal: average value =  0.351D-02
+ n1qn3: descent direction d: angle(-g,d) =  63.1 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  28, simul  36, f= 0.49926693D-01, h'(0)=-0.55493D-01
+ n1qn3: line search
+     nlis0       fpn=-0.555D-01 d2= 0.58D-01  tmin= 0.11D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.507D-01  0.260D+00
+     nlis0     0.413D+00 -0.132D-01  0.278D-02
+ n1qn3: stopping criterion on g:  0.14651D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.14924D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.980D+00
+     Oren-Spedicato factor (not used) =  0.274D-02     diagonal: average value =  0.355D-02
+ n1qn3: descent direction d: angle(-g,d) =  84.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  29, simul  38, f= 0.36768990D-01, h'(0)=-0.91754D-02
+ n1qn3: line search
+     nlis0       fpn=-0.918D-02 d2= 0.87D-03  tmin= 0.86D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.813D-02 -0.718D-02
+ n1qn3: stopping criterion on g:  0.11306D-01
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.29590D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.112D+01
+     Oren-Spedicato factor (not used) =  0.329D-02     diagonal: average value =  0.403D-02
+ n1qn3: descent direction d: angle(-g,d) =  85.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  30, simul  39, f= 0.28635823D-01, h'(0)=-0.25872D-01
+ n1qn3: line search
+     nlis0       fpn=-0.259D-01 d2= 0.13D-01  tmin= 0.22D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.174D-01 -0.109D-01
+ n1qn3: stopping criterion on g:  0.64294D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.38281D+01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.337D+01
+     Oren-Spedicato factor (not used) =  0.116D-01     diagonal: average value =  0.134D-01
+ n1qn3: descent direction d: angle(-g,d) =  84.3 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  31, simul  40, f= 0.11210318D-01, h'(0)=-0.94861D-02
+ n1qn3: line search
+     nlis0       fpn=-0.949D-02 d2= 0.41D-02  tmin= 0.39D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.648D-02 -0.412D-02
+ n1qn3: stopping criterion on g:  0.32255D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.56905D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.218D+00
+     Oren-Spedicato factor (not used) =  0.229D-02     diagonal: average value =  0.316D-02
+ n1qn3: descent direction d: angle(-g,d) =  28.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  32, simul  41, f= 0.47312381D-02, h'(0)=-0.12658D-01
+ n1qn3: line search
+     nlis0       fpn=-0.127D-01 d2= 0.37D-01  tmin= 0.13D-14 tmax= 0.10D+21
+     nlis0                                      0.100D+01  0.236D-01  0.889D-01
+     nlis0     0.308D+00 -0.202D-02  0.129D-02
+ n1qn3: stopping criterion on g:  0.59635D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.92424D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.995D+00
+     Oren-Spedicato factor (not used) =  0.211D-02     diagonal: average value =  0.320D-02
+ n1qn3: descent direction d: angle(-g,d) =  86.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  33, simul  43, f= 0.27115837D-02, h'(0)=-0.27047D-02
+ n1qn3: line search
+     nlis0       fpn=-0.270D-02 d2= 0.13D-02  tmin= 0.68D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.193D-02 -0.124D-02
+ n1qn3: stopping criterion on g:  0.26088D-02
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.61219D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.107D+01
+     Oren-Spedicato factor (not used) =  0.240D-02     diagonal: average value =  0.346D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  34, simul  44, f= 0.78628297D-03, h'(0)=-0.10616D-02
+ n1qn3: line search
+     nlis0       fpn=-0.106D-02 d2= 0.11D-02  tmin= 0.73D-14 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.680D-03 -0.342D-03
+ n1qn3: stopping criterion on g:  0.93209D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.93175D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.195D+01
+     Oren-Spedicato factor (not used) =  0.472D-02     diagonal: average value =  0.669D-02
+ n1qn3: descent direction d: angle(-g,d) =  87.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  35, simul  45, f= 0.10596353D-03, h'(0)=-0.16466D-03
+ n1qn3: line search
+     nlis0       fpn=-0.165D-03 d2= 0.22D-03  tmin= 0.17D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.994D-04 -0.378D-04
+ n1qn3: stopping criterion on g:  0.13605D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.43415D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.797D+00
+     Oren-Spedicato factor (not used) =  0.369D-02     diagonal: average value =  0.540D-02
+ n1qn3: descent direction d: angle(-g,d) =  86.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  36, simul  46, f= 0.65777606D-05, h'(0)=-0.12835D-04
+ n1qn3: line search
+     nlis0       fpn=-0.128D-04 d2= 0.34D-04  tmin= 0.43D-13 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.592D-05  0.983D-06
+ n1qn3: stopping criterion on g:  0.15294D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.39644D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.278D+03
+     Oren-Spedicato factor (not used) =  0.637D+00     diagonal: average value =  0.256D+01
+ n1qn3: descent direction d: angle(-g,d) =  77.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  37, simul  47, f= 0.65293408D-06, h'(0)=-0.24132D-05
+ n1qn3: line search
+     nlis0       fpn=-0.241D-05 d2= 0.96D-07  tmin= 0.92D-12 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.168D-06  0.208D-05
+ n1qn3: stopping criterion on g:  0.13388D-03
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.53137D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.794D-03
+     Oren-Spedicato factor (not used) =  0.101D-02     diagonal: average value =  0.159D-02
+ n1qn3: descent direction d: angle(-g,d) =  26.4 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  38, simul  48, f= 0.48527633D-06, h'(0)=-0.97036D-06
+ n1qn3: line search
+     nlis0       fpn=-0.970D-06 d2= 0.12D-08  tmin= 0.80D-11 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.485D-06 -0.183D-09
+ n1qn3: stopping criterion on g:  0.27051D-07
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.11218D+00
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.866D+00
+     Oren-Spedicato factor (not used) =  0.999D-03     diagonal: average value =  0.137D-02
+ n1qn3: descent direction d: angle(-g,d) =  83.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  39, simul  49, f= 0.73395789D-12, h'(0)=-0.14314D-11
+ n1qn3: line search
+     nlis0       fpn=-0.143D-11 d2= 0.35D-11  tmin= 0.13D-09 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.734D-12 -0.116D-13
+ n1qn3: stopping criterion on g:  0.61631D-10
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.53526D-01
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.399D+02
+     Oren-Spedicato factor (not used) =  0.356D-01     diagonal: average value =  0.717D-01
+ n1qn3: descent direction d: angle(-g,d) =  64.0 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  40, simul  50, f= 0.20622903D-15, h'(0)=-0.99086D-16
+ n1qn3: line search
+     nlis0       fpn=-0.991D-16 d2= 0.25D-15  tmin= 0.16D-07 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.152D-15 -0.112D-20
+ n1qn3: stopping criterion on g:  0.92861D-13
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.84647D-02
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.571D+01
+     Oren-Spedicato factor (not used) =  0.482D+00     diagonal: average value =  0.356D+00
+ n1qn3: descent direction d: angle(-g,d) =  80.8 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  41, simul  51, f= 0.53784084D-16, h'(0)=-0.47400D-24
+ n1qn3: line search
+     nlis0       fpn=-0.474D-24 d2= 0.19D-25  tmin= 0.17D-02 tmax= 0.10D+21
+     nlis0     0.100D+01 -0.599D-21 -0.107D-26
+ n1qn3: stopping criterion on g:  0.21365D-15
+ n1qn3: convergence rate, s(k)/s(k-1) =  0.86648D-05
+ n1qn3: matrix update:
+     fitting the ellipsoid: factor  0.352D-02
+     Oren-Spedicato factor (not used) =  0.102D-02     diagonal: average value =  0.198D-02
+ n1qn3: descent direction d: angle(-g,d) =  29.9 degrees
+ -------------------------------------------------------------------------------
+
+ n1qn3: iter  42, simul  52, f= 0.53783484D-16, h'(0)=-0.24769D-29
+ n1qn3: line search
+     nlis0       fpn=-0.248D-29 d2= 0.33D-32  tmin= 0.39D+01 tmax= 0.10D+21
+     nlis0                                      0.387D+01  0.326D-23  0.743D-29
+     nlis0   fin sur tmin        pas            fonctions     derivees
+     nlis0                  0.00000000D+00    0.53783484D-16 -0.248D-29
+     nlis0                  0.38722999D+01    0.53783488D-16  0.743D-29
+ -------------------------------------------------------------------------------
+ n1qn3: output mode is  6
+     number of iterations:   42
+     number of simulations:     53
+     realized relative precision on g:  0.21D-15
+     norm of x =  0.14142135D+01
+     f         =  0.53783484D-16
+     norm of x =  0.14924691D-12
+Optimization stops because too small variations for x.
+