Starting preprocessing of the model file ...
Substitution of Expectation operator: added 2 auxiliary variables and equations.
Substitution of endo leads >= 2: added 2 auxiliary variables and equations.
Substitution of endo lags >= 2: added 13 auxiliary variables and equations.
Found 28 equation(s).
Evaluating expressions...done
Computing static model derivatives (order 1).
Computing dynamic model derivatives (order 1).
Processing outputs ...
done
Preprocessing completed.
STEADY-STATE RESULTS:
pi_p 0
pi_w 0
x_tilde 0
pi_p_bar 0
r 0
r_bar 0
mu_p 0
g 0
s 0
U 0
pi_w_bar 0
EIGENVALUES:
Modulus Real Imaginary
6.61e-20 6.61e-20 0
2.007e-17 -2.007e-17 0
3.596e-17 -3.596e-17 0
5.105e-17 5.105e-17 0
0.0879 -0.011 0.08721
0.0879 -0.011 -0.08721
0.1375 -0.1375 0
0.2524 -0.2524 0
0.2905 -0.2905 0
0.5427 0.5427 0
0.5693 -0.5693 0
0.6036 -0.06731 0.5999
0.6036 -0.06731 -0.5999
0.6057 -0.6057 0
0.6426 -0.07203 0.6385
0.6426 -0.07203 -0.6385
0.7834 0.6781 0.3923
0.7834 0.6781 -0.3923
0.864 0.864 0
0.8886 0.8886 0
0.9638 0.9638 0
0.9984 0.9979 0.02918
0.9984 0.9979 -0.02918
1 1 0
1 1 0
1 1 0
1.03 1.03 0.01891
1.03 1.03 -0.01891
1.124e+22 1.124e+22 0
Inf -Inf 0
There are 4 eigenvalue(s) larger than 1 in modulus
for 4 forward-looking variable(s)
The rank condition is verified.
MODEL SUMMARY
Number of variables: 28
Number of stochastic shocks: 10
Number of state variables: 26
Number of jumpers: 4
Number of static variables: 0
MATRIX OF COVARIANCE OF EXOGENOUS SHOCKS
Variables e_pi_p e_pi_w e_mu e_s e_g e_U e_x e_r e_pi_p_bar e_r_bar
e_pi_p 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
e_pi_w 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
e_mu 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
e_s 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
e_g 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
e_U 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
e_x 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
e_r 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
e_pi_p_bar 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
e_r_bar 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
POLICY AND TRANSITION FUNCTIONS
pi_p pi_w x_tilde pi_p_bar r r_bar mu_p g s U pi_w_bar
x_tilde(-1) 0 0 1.340068 0 0.527719 0 0 0 0 -0.433692 0
pi_p_bar(-1) 0.004455 0 0.051434 1.000000 -0.091220 0 -0.000011 0 -0.000014 -0.013203 1.000000
r(-1) 0 0 -0.093015 0 1.409234 0 0 0 0 0.023876 0
r_bar(-1) 0 0 0.056990 0 0.035760 1.000000 0 0 0 -0.014629 0
mu_p(-1) 0.002965 0 0 0 0.000340 0 0.999993 0 1.249991 0 0
g(-1) 0 0 0 0 0 0 -0.002500 1.000000 -0.003125 0 1.000000
s(-1) 0 -0.001406 0 0 0 0 -0.000004 0 -0.000004 0 0
U(-1) 0 -0.027309 0 0 0 0 -0.000068 0 -0.000085 0.290350 0
pi_w_bar(-1) 0 0 -0.005238 0 -0.001064 0 0 0 0 0.001345 0
EXPECTATION(-1)(...) 0.830111 0 0 0 0.095100 0 -0.002075 0 -0.002594 0 0
EXPECTATION(-1)(...) 0 0.787149 0 0 0 0 0.001968 0 0.002460 0 0
pi_p(-2) 0.003121 0 -0.108204 0 0.002327 0 -0.000008 0 -0.000010 0.027775 0
pi_p(-3) 0.003121 0 -0.149169 0 -0.024413 0 -0.000008 0 -0.000010 0.038290 0
pi_p(-4) 0.156069 0 0.251117 0 0.068870 0 -0.000390 0 -0.000488 -0.064459 0
pi_w(-2) 0 0.004016 0.184995 0 0.037564 0 0.000010 0 0.000013 -0.047486 0
pi_w(-3) 0 0.004016 -0.040391 0 -0.008202 0 0.000010 0 0.000013 0.010368 0
pi_w(-4) 0 0.200803 -0.042870 0 -0.008705 0 0.000502 0 0.000628 0.011004 0
U(-2) 0 0 0 0 0 0 0 0 0 0.137020 0
x_tilde(-2) 0 0 -0.148829 0 -0.518390 0 0 0 0 0.038203 0
x_tilde(-3) 0 0 -0.209489 0 -0.003615 0 0 0 0 0.053774 0
x_tilde(-4) 0 0 -0.009072 0 -0.001842 0 0 0 0 0.002329 0
r(-2) 0 0 -0.159634 0 -0.491283 0 0 0 0 0.040976 0
r(-3) 0 0 0.143998 0 0.053467 0 0 0 0 -0.036963 0
r(-4) 0 0 0.051661 0 -0.007178 0 0 0 0 -0.013261 0
pi_p(-1) 0.003121 0 -0.045178 0 -0.050665 0 -0.000008 0 -0.000010 0.011597 0
pi_w(-1) 0 0.004016 -0.096496 0 -0.019594 0 0.000010 0 0.000013 0.024770 0
e_pi_p 1.000000 0 0 0 0.114563 0 -0.002500 0 -0.003125 0 0
e_pi_w 0 1.000000 0 0 0 0 0.002500 0 0.003125 0 0
e_mu 0 0 0 0 0 0 1.000000 0 1.250000 0 0
e_s 0 0 0 0 0 0 0 0 1.000000 0 0
e_g 0 0 0 0 0 0 -0.002500 1.000000 -0.003125 0 1.000000
e_U 0 0 0 0 0 0 0 0 0 1.000000 0
e_x 0 0 1.000000 0 0.203054 0 0 0 0 -0.256690 0
e_r 0 0 0 0 1.000000 0 0 0 0 0 0
e_pi_p_bar 0 0 0 1.000000 0 0 0 0 0 0 1.000000
e_r_bar 0 0 0 0 0 1.000000 0 0 0 0 0
THEORETICAL MOMENTS
VARIABLE MEAN STD. DEV. VARIANCE
pi_p NaN NaN NaN
pi_w NaN NaN NaN
x_tilde 0.0000 0.0000 0.0000
pi_p_bar NaN NaN NaN
r NaN NaN NaN
r_bar NaN NaN NaN
mu_p 0.0000 0.0000 0.0000
g NaN NaN NaN
s 0.0000 0.0000 0.0000
U 0.0000 0.0000 0.0000
pi_w_bar NaN NaN NaN
VARIANCE DECOMPOSITION (in percent)
e_pi_p e_pi_w e_mu e_s e_g e_U e_x e_r e_pi_p_bar e_r_bar
x_tilde 3.01 0.11 67.89 0.00 0.56 0.00 4.66 15.12 6.15 2.50
mu_p 2.21 0.11 79.98 0.00 0.59 0.00 1.12 10.01 4.15 1.83
s 2.20 0.11 79.63 0.44 0.59 0.00 1.12 9.97 4.13 1.82
U 2.97 0.10 67.71 0.00 0.56 1.03 4.15 14.89 6.12 2.48
All endogenous are constant or non stationary, not displaying correlations and auto-correlations
Initial value of the log posterior (or likelihood): -14367.5389
-----------------
f at the beginning of new iteration, 14367.5389494184
Predicted improvement: 4016092.562939641
lambda = 1; f = 14379.9548579
lambda = 0.33333; f = 14372.1076016
lambda = 0.11111; f = 14368.0184782
lambda = 0.037037; f = 14367.5833973
lambda = 0.012346; f = 14367.5414771
lambda = 0.0041152; f = 3952.0910685
Norm of dx 28.341
----
Improvement on iteration 1 = 10415.447880911
-----------------
f at the beginning of new iteration, 3952.0910685071
Predicted improvement: 138257.007368678
lambda = 1; f = 817.6321078
lambda = 0.33333; f = 711.7981810
lambda = 0.11111; f = 1272.7366925
lambda = 0.037037; f = 1769.6341079
lambda = 0.012346; f = 2355.6478844
Norm of dx 5.3477
----
Improvement on iteration 2 = 3240.292887476
-----------------
f at the beginning of new iteration, 711.7981810308
Predicted improvement: 1967.982080679
lambda = 1; f = 712.0258953
lambda = 0.33333; f = 711.8077433
lambda = 0.11111; f = 398.7599433
lambda = 0.2148; f = 711.7991008
lambda = 0.14463; f = 378.0925309
Norm of dx 0.63002
----
Improvement on iteration 3 = 333.705650155
-----------------
f at the beginning of new iteration, 378.0925308754
Predicted improvement: 282.922235548
lambda = 1; f = 316.4481371
lambda = 0.33333; f = 1981.2660376
lambda = 0.11111; f = 311.6967155
lambda = 0.2148; f = 254.5868712
lambda = 0.41524; f = 438.7954071
lambda = 0.2796; f = 305.8602517
Norm of dx 0.24741
----
Improvement on iteration 4 = 123.505659678
-----------------
f at the beginning of new iteration, 254.5868711977
Predicted improvement: 2581.153460543
lambda = 1; f = 255.1250128
lambda = 0.33333; f = 449.3358267
lambda = 0.11111; f = 307.3968589
lambda = 0.037037; f = 472.7721171
lambda = 0.012346; f = 746.4584078
lambda = 0.0041152; f = 234.4439119
lambda = 0.0079555; f = 222.1115633
lambda = 0.015379; f = 432.9451416
lambda = 0.010356; f = 241.9299758
lambda = 0.008168; f = 222.0693996
lambda = 0.0094179; f = 226.3632141
Norm of dx 2.3752
----
Improvement on iteration 5 = 32.517471620
warning: possible inaccuracy in H matrix
-----------------
f at the beginning of new iteration, 222.0693995773
Predicted improvement: 520.147357514
lambda = 1; f = 796.4411412
lambda = 0.33333; f = 42590.4230452
lambda = 0.11111; f = 333.5406868
lambda = 0.037037; f = 31183.7124238
lambda = 0.012346; f = 214.8134985
Norm of dx 0.42273
----
Improvement on iteration 6 = 7.255901105
-----------------
f at the beginning of new iteration, 214.8134984720
Predicted improvement: 220.733463198
lambda = 1; f = -32.0839571
Norm of dx 0.24484
----
Improvement on iteration 7 = 246.897455556
-----------------
f at the beginning of new iteration, -32.0839570836
Predicted improvement: 137.189275950
lambda = 1; f = 2127621.2980580
lambda = 0.33333; f = -81.7728878
Norm of dx 0.25928
----
Improvement on iteration 8 = 49.688930696
-----------------
f at the beginning of new iteration, -81.7728877799
Predicted improvement: 357.366514839
lambda = 1; f = -81.5344614
lambda = 0.33333; f = -81.7721916
lambda = 0.11111; f = -123.9211199
Norm of dx 0.74814
----
Improvement on iteration 9 = 42.148232143
-----------------
f at the beginning of new iteration, -123.9211199232
Predicted improvement: 41.573413483
lambda = 1; f = -123.9208238
lambda = 0.33333; f = -146.0934429
lambda = 0.64439; f = -156.5321996
Norm of dx 0.15572
----
Improvement on iteration 10 = 32.611079679
-----------------
f at the beginning of new iteration, -156.5321996017
Predicted improvement: 8.587264578
lambda = 1; f = -156.5321835
lambda = 0.33333; f = -161.8579543
lambda = 0.64439; f = -166.0289817
lambda = 1.2457; f = -156.5320195
lambda = 0.8388; f = -168.1808886
lambda = 1.0635; f = -156.5321580
lambda = 0.92232; f = -156.5321985
lambda = 0.8468; f = -168.2611421
lambda = 0.89134; f = -168.6952910
lambda = 0.93822; f = -156.5321969
lambda = 0.9098; f = -156.5321993
lambda = 0.89317; f = -168.7126606
lambda = 0.90311; f = -156.5321995
Norm of dx 0.050388
Cliff. Perturbing search direction.
Predicted improvement: 33.544377869
lambda = 1; f = -156.5312502
lambda = 0.33333; f = -177.1973878
lambda = 0.64439; f = -156.5321407
lambda = 0.4339; f = -182.7279331
lambda = 0.55011; f = -156.5321972
lambda = 0.4771; f = -184.9959268
lambda = 0.51965; f = -187.1637838
lambda = 0.566; f = -156.5321930
lambda = 0.53772; f = -156.5321991
lambda = 0.52143; f = -187.2529611
lambda = 0.53114; f = -156.5321995
lambda = 0.52529; f = -187.4462227
Norm of dx 0.07969
Cliff again. Try traversing
Predicted improvement: 1382433.357464821
lambda = 1; f = 555047.0214850
lambda = 0.33333; f = 61077.2436681
lambda = 0.11111; f = 6496.7741013
lambda = 0.037037; f = 533.9645833
lambda = 0.012346; f = -94.6741074
lambda = 0.0041152; f = -153.2242841
lambda = 0.0013717; f = -156.4336620
lambda = 0.00045725; f = -156.5228427
lambda = 0.00015242; f = -156.5317294
lambda = 5.0805e-05; f = -266.1703704
lambda = 9.8216e-05; f = -156.5321300
lambda = 6.6133e-05; f = -156.5321994
lambda = 5.2163e-05; f = -268.1758499
lambda = 6.0144e-05; f = -279.1771655
lambda = 6.9347e-05; f = -156.5321981
lambda = 6.3669e-05; f = -283.6280117
lambda = 6.7018e-05; f = -156.5321992
lambda = 6.4988e-05; f = -156.5321996
lambda = 6.38e-05; f = -283.7884199
lambda = 6.451e-05; f = -156.5321996
Norm of dx 1662.8
----
Improvement on iteration 11 = 127.256220302
back and forth on step length never finished
-----------------
f at the beginning of new iteration, -283.7884199033
Predicted improvement: 8.587264578
lambda = 1; f = -283.7884040
lambda = 0.33333; f = -288.3916936
Norm of dx 0.050388
----
Improvement on iteration 12 = 4.603273739
-----------------
f at the beginning of new iteration, -288.3916936424
Predicted improvement: 14.193559398
lambda = 1; f = -288.3916595
lambda = 0.33333; f = -296.4271526
lambda = 0.64439; f = -301.4165606
lambda = 1.2457; f = -288.3915357
lambda = 0.8388; f = -288.3916916
lambda = 0.66161; f = -301.6247261
lambda = 0.76285; f = -302.7091371
Norm of dx 0.039464
----
Improvement on iteration 13 = 14.317443496
warning: possible inaccuracy in H matrix
-----------------
f at the beginning of new iteration, -302.7091371385
Predicted improvement: 9.098710584
lambda = 1; f = -302.7082341
lambda = 0.33333; f = -302.7090452
lambda = 0.11111; f = -302.7091295
lambda = 0.037037; f = -302.7091369
lambda = 0.012346; f = -302.9326595
lambda = 0.023866; f = -302.7091371
lambda = 0.01607; f = -302.9996572
lambda = 0.020374; f = -303.0768244
lambda = 0.025831; f = -302.7091371
lambda = 0.022403; f = -302.7091371
lambda = 0.020569; f = -303.0803020
lambda = 0.02165; f = -302.7091371
lambda = 0.020995; f = -303.0879251
lambda = 0.021386; f = -302.7091371
lambda = 0.02115; f = -302.7091371
Norm of dx 0.046267
Cliff. Perturbing search direction.
Predicted improvement: 17.476829023
lambda = 1; f = -302.7072057
lambda = 0.33333; f = -302.7089350
lambda = 0.11111; f = -302.7091186
lambda = 0.037037; f = -302.7091361
lambda = 0.012346; f = -303.1388428
lambda = 0.023866; f = -302.7091370
lambda = 0.01607; f = -302.7091371
lambda = 0.012676; f = -303.1502743
lambda = 0.014615; f = -302.7091371
lambda = 0.013418; f = -303.1760102
lambda = 0.014124; f = -303.2004475
lambda = 0.014867; f = -302.7091371
lambda = 0.014417; f = -303.2105752
lambda = 0.014685; f = -302.7091371
lambda = 0.014523; f = -302.7091371
Norm of dx 0.068455
Cliff again. Try traversing
Predicted improvement: 59472.426184416
lambda = 1; f = 98707.5283225
lambda = 0.33333; f = 10489.7065090
lambda = 0.11111; f = 828.2572269
lambda = 0.037037; f = -198.4040018
lambda = 0.012346; f = -296.8816198
lambda = 0.0041152; f = -302.6465630
lambda = 0.0013717; f = -302.7030762
lambda = 0.00045725; f = -302.7087202
lambda = 0.00015242; f = -302.7091355
lambda = 5.0805e-05; f = -307.9846674
lambda = 9.8216e-05; f = -311.5791739
lambda = 0.00018987; f = -302.7091240
lambda = 0.00012785; f = -313.1868301
Norm of dx 344.88
----
Improvement on iteration 14 = 10.477692915
-----------------
f at the beginning of new iteration, -313.1868300534
Predicted improvement: 31.765391815
lambda = 1; f = -313.1848546
lambda = 0.33333; f = -313.1866269
lambda = 0.11111; f = -313.1868125
lambda = 0.037037; f = -313.1868294
lambda = 0.012346; f = -313.9682311
lambda = 0.023866; f = -313.1868300
lambda = 0.01607; f = -314.2028396
lambda = 0.020374; f = -313.1868300
lambda = 0.01767; f = -314.3034683
lambda = 0.019246; f = -313.1868301
lambda = 0.018285; f = -314.3420715
lambda = 0.018856; f = -313.1868301
lambda = 0.018511; f = -314.3562912
lambda = 0.018717; f = -313.1868301
Norm of dx 0.1188
Cliff. Perturbing search direction.
Predicted improvement: 39.645444604
lambda = 1; f = -313.1830474
lambda = 0.33333; f = -313.1864325
lambda = 0.11111; f = -313.1867930
lambda = 0.037037; f = -313.1868279
lambda = 0.012346; f = -314.1611121
lambda = 0.023866; f = -313.1868296
lambda = 0.01607; f = -313.1868300
lambda = 0.012676; f = -314.1870170
lambda = 0.014615; f = -313.1868300
lambda = 0.013418; f = -314.2453341
lambda = 0.014124; f = -313.1868301
lambda = 0.013696; f = -313.1868301
lambda = 0.013446; f = -314.2474955
lambda = 0.013596; f = -313.1868301
Norm of dx 0.16383
Cliff again. Try traversing
Predicted improvement: 152917.738331425
lambda = 1; f = 6297.4357680
lambda = 0.33333; f = 411.3281416
lambda = 0.11111; f = -235.9723693
lambda = 0.037037; f = -305.6574991
lambda = 0.012346; f = -312.6545749
lambda = 0.0041152; f = -313.1834556
lambda = 0.0013717; f = -313.1866340
lambda = 0.00045725; f = -313.1868099
lambda = 0.00015242; f = -313.1868283
lambda = 5.0805e-05; f = -313.1868300
lambda = 1.6935e-05; f = -318.0911906
lambda = 3.2739e-05; f = -313.1868300
lambda = 2.2044e-05; f = -319.4336150
lambda = 2.7948e-05; f = -313.1868301
lambda = 2.4239e-05; f = -319.9892107
lambda = 2.6401e-05; f = -313.1868301
lambda = 2.5082e-05; f = -320.1990234
lambda = 2.5865e-05; f = -313.1868301
lambda = 2.5392e-05; f = -320.2758376
lambda = 2.5675e-05; f = -313.1868301
Norm of dx 553.02
----
Improvement on iteration 15 = 7.089007595
back and forth on step length never finished
-----------------
f at the beginning of new iteration, -320.2758376481
Predicted improvement: 31.765391815
lambda = 1; f = -320.2738620
lambda = 0.33333; f = -320.2756344
lambda = 0.11111; f = -320.2758201
lambda = 0.037037; f = -320.2758369
lambda = 0.012346; f = -320.9923961
Norm of dx 0.1188
----
Improvement on iteration 16 = 0.716558433
-----------------
f at the beginning of new iteration, -320.9923960811
Predicted improvement: 27.067243051
lambda = 1; f = -320.9905805
lambda = 0.33333; f = -320.9921996
lambda = 0.11111; f = -320.9923760
lambda = 0.037037; f = -320.9923944
lambda = 0.012346; f = -320.9923960
lambda = 0.0041152; f = -321.2149335
lambda = 0.0079555; f = -320.9923961
lambda = 0.0053568; f = -321.2819821
lambda = 0.0067914; f = -320.9923961
lambda = 0.0058901; f = -321.3107741
lambda = 0.0064154; f = -321.3391226
lambda = 0.0069876; f = -320.9923961
lambda = 0.0066385; f = -320.9923961
lambda = 0.0064374; f = -321.3403076
lambda = 0.0065573; f = -320.9923961
lambda = 0.0064851; f = -321.3428808
Norm of dx 0.10972
Cliff. Perturbing search direction.
Predicted improvement: 39.400676980
lambda = 1; f = -320.9894076
lambda = 0.33333; f = -320.9920708
lambda = 0.11111; f = -320.9923622
lambda = 0.037037; f = -320.9923930
lambda = 0.012346; f = -320.9923959
lambda = 0.0041152; f = -321.3162420
lambda = 0.0079555; f = -320.9923961
lambda = 0.0053568; f = -320.9923961
lambda = 0.0042252; f = -321.3248825
lambda = 0.0048717; f = -321.3756776
lambda = 0.0056171; f = -320.9923961
lambda = 0.0051572; f = -320.9923961
lambda = 0.0048995; f = -321.3778632
lambda = 0.0050525; f = -321.3898824
lambda = 0.0052103; f = -320.9923961
lambda = 0.0051151; f = -320.9923961
lambda = 0.0050588; f = -321.3903709
Norm of dx 0.15863
Cliff again. Try traversing
Predicted improvement: 86940.724186567
lambda = 1; f = 149437.9933819
lambda = 0.33333; f = 16073.3035228
lambda = 0.11111; f = 1419.9932632
lambda = 0.037037; f = -153.1561166
lambda = 0.012346; f = -309.6432936
lambda = 0.0041152; f = -320.9622110
lambda = 0.0013717; f = -164.0512694
lambda = 0.00045725; f = -351.5036232
Norm of dx 416.99
----
Improvement on iteration 17 = 30.511227160
-----------------
f at the beginning of new iteration, -351.5036232412
Predicted improvement: 35.329776669
lambda = 1; f = -351.5027079
lambda = 0.33333; f = -351.5035293
lambda = 0.11111; f = -351.5036152
lambda = 0.037037; f = -351.5036229
lambda = 0.012346; f = -352.3664277
lambda = 0.023866; f = -351.5036232
lambda = 0.01607; f = -352.6230253
lambda = 0.020374; f = -351.5036232
lambda = 0.01767; f = -352.7327418
lambda = 0.019246; f = -351.5036232
lambda = 0.018285; f = -352.7747767
lambda = 0.018856; f = -352.8138149
lambda = 0.019445; f = -351.5036232
lambda = 0.019089; f = -351.5036232
lambda = 0.018879; f = -352.8154007
Norm of dx 0.14764
Cliff. Perturbing search direction.
Predicted improvement: 54.771682328
lambda = 1; f = -351.5023824
lambda = 0.33333; f = -351.5034944
lambda = 0.11111; f = -351.5036117
lambda = 0.037037; f = -351.5036227
lambda = 0.012346; f = -352.8278136
lambda = 0.023866; f = -351.5036232
lambda = 0.01607; f = -353.2164774
lambda = 0.020374; f = -351.5036232
lambda = 0.01767; f = -351.5036232
lambda = 0.016224; f = -353.2323655
lambda = 0.017077; f = -351.5036232
lambda = 0.01656; f = -351.5036232
lambda = 0.016257; f = -353.2358161
lambda = 0.016438; f = -351.5036232
Norm of dx 0.22392
Cliff again. Try traversing
Predicted improvement: 111009.071089543
lambda = 1; f = 175264.1742428
lambda = 0.33333; f = 18929.5801148
lambda = 0.11111; f = 1714.5605120
lambda = 0.037037; f = -146.4283033
lambda = 0.012346; f = -335.9594346
lambda = 0.0041152; f = -351.2700133
lambda = 0.0013717; f = -157.7674737
lambda = 0.00045725; f = -353.0592334
lambda = 0.00015242; f = -366.3638263
Norm of dx 471.19
----
Improvement on iteration 18 = 14.860203075
-----------------
f at the beginning of new iteration, -366.3638263165
Predicted improvement: 18.062516653
lambda = 1; f = -366.3620045
lambda = 0.33333; f = -366.3636349
lambda = 0.11111; f = -366.3638085
lambda = 0.037037; f = -366.3638253
lambda = 0.012346; f = -366.8088926
lambda = 0.023866; f = -366.3638261
lambda = 0.01607; f = -366.3638263
lambda = 0.012676; f = -366.8207585
lambda = 0.014615; f = -366.3638263
lambda = 0.013418; f = -366.8474772
lambda = 0.014124; f = -366.3638263
lambda = 0.013696; f = -366.3638263
lambda = 0.013446; f = -366.8484676
lambda = 0.013596; f = -366.8538509
lambda = 0.013747; f = -366.3638263
Norm of dx 0.092285
Cliff. Perturbing search direction.
Predicted improvement: 23.326714097
lambda = 1; f = -366.3612448
lambda = 0.33333; f = -366.3635526
lambda = 0.11111; f = -366.3638001
lambda = 0.037037; f = -366.3638246
lambda = 0.012346; f = -366.3638263
lambda = 0.0041152; f = -366.5556456
lambda = 0.0079555; f = -366.7343485
lambda = 0.015379; f = -366.3638263
lambda = 0.010356; f = -366.8458889
lambda = 0.013129; f = -366.3638263
lambda = 0.011387; f = -366.8937729
lambda = 0.012402; f = -366.3638263
lambda = 0.011783; f = -366.3638263
lambda = 0.011426; f = -366.8955825
lambda = 0.011638; f = -366.3638263
lambda = 0.01151; f = -366.3638263
Norm of dx 0.11218
Cliff again. Try traversing
Predicted improvement: 106468.149555603
lambda = 1; f = 3715.0050621
lambda = 0.33333; f = 52.0206100
lambda = 0.11111; f = -330.2599960
lambda = 0.037037; f = -364.9057102
lambda = 0.012346; f = 6128.0103747
lambda = 0.0041152; f = 611.7627391
lambda = 0.0013717; f = 768167.3072312
lambda = 0.00045725; f = 1274821.1655296
lambda = 0.00015242; f = -369.7811270
lambda = 5.0805e-05; f = -375.2891838
lambda = 9.8216e-05; f = -378.3594861
Norm of dx 461.45
----
Improvement on iteration 19 = 11.995659794
-----------------
f at the beginning of new iteration, -378.3594861108
Predicted improvement: 13.075470877
lambda = 1; f = -378.3581666
lambda = 0.33333; f = -378.3593570
lambda = 0.11111; f = -378.3594769
lambda = 0.037037; f = -378.3594860
lambda = 0.012346; f = -378.6814021
lambda = 0.023866; f = -378.9801452
lambda = 0.046138; f = -378.3594857
lambda = 0.031067; f = -378.3594861
lambda = 0.024504; f = -378.9966331
lambda = 0.028254; f = -379.0934860
lambda = 0.032577; f = -378.3594861
lambda = 0.029909; f = -378.3594861
lambda = 0.028415; f = -379.0976505
lambda = 0.029302; f = -379.1205477
lambda = 0.030218; f = -378.3594861
lambda = 0.029665; f = -379.1299029
lambda = 0.029995; f = -378.3594861
Norm of dx 0.10468
Cliff. Perturbing search direction.
Predicted improvement: 15.196708111
lambda = 1; f = -378.3570951
lambda = 0.33333; f = -378.3592442
lambda = 0.11111; f = -378.3594664
lambda = 0.037037; f = -378.3594856
lambda = 0.012346; f = -378.7330378
lambda = 0.023866; f = -378.3594861
lambda = 0.01607; f = -378.8450862
lambda = 0.020374; f = -378.9741888
lambda = 0.025831; f = -378.3594861
lambda = 0.022403; f = -378.3594861
lambda = 0.020569; f = -378.9800080
lambda = 0.02165; f = -379.0123891
lambda = 0.022789; f = -378.3594861
lambda = 0.022099; f = -379.0258061
lambda = 0.02251; f = -378.3594861
lambda = 0.022263; f = -379.0307029
Norm of dx 0.11936
Cliff again. Try traversing
Predicted improvement: 17602.084384439
lambda = 1; f = 31480.4251436
lambda = 0.33333; f = 3068.0888398
lambda = 0.11111; f = -25.6831982
lambda = 0.037037; f = -348.3893905
lambda = 0.012346; f = -377.2760302
lambda = 0.0041152; f = -362.0027446
lambda = 0.0013717; f = 124.8769117
lambda = 0.00045725; f = -380.8291477
lambda = 0.00015242; f = -382.4894478
lambda = 0.00029465; f = -383.7057668
Norm of dx 187.63
----
Improvement on iteration 20 = 5.346280648
-----------------
f at the beginning of new iteration, -383.7057667591
Predicted improvement: 24.295082554
lambda = 1; f = -383.7042067
lambda = 0.33333; f = -383.7056172
lambda = 0.11111; f = -383.7057569
lambda = 0.037037; f = -383.7057667
lambda = 0.012346; f = -384.3023910
lambda = 0.023866; f = -384.8533170
lambda = 0.046138; f = -383.7057665
lambda = 0.031067; f = -385.1947967
lambda = 0.039387; f = -383.7057667
lambda = 0.03416; f = -385.3408304
lambda = 0.037207; f = -383.7057667
lambda = 0.035348; f = -383.7057668
lambda = 0.034277; f = -385.3463427
lambda = 0.034915; f = -383.7057668
lambda = 0.034531; f = -383.7057668
Norm of dx 0.19211
Cliff. Perturbing search direction.
Predicted improvement: 34.246394064
lambda = 1; f = -383.7040580
lambda = 0.33333; f = -383.7056018
lambda = 0.11111; f = -383.7057556
lambda = 0.037037; f = -383.7057667
lambda = 0.012346; f = -384.5440896
lambda = 0.023866; f = -385.3134073
lambda = 0.046138; f = -383.7057664
lambda = 0.031067; f = -385.7879122
lambda = 0.039387; f = -383.7057667
lambda = 0.03416; f = -383.7057668
lambda = 0.031363; f = -385.8073343
lambda = 0.033013; f = -383.7057668
lambda = 0.032013; f = -385.8498920
lambda = 0.032609; f = -385.8889126
lambda = 0.033216; f = -383.7057668
lambda = 0.032851; f = -383.7057668
Norm of dx 0.23194
Cliff again. Try traversing
Predicted improvement: 132610.379632158
lambda = 1; f = 220762.5815557
lambda = 0.33333; f = 23952.4705725
lambda = 0.11111; f = 2242.8604388
lambda = 0.037037; f = -116.8903456
lambda = 0.012346; f = -361.6518462
lambda = 0.0041152; f = -383.0387378
lambda = 0.0013717; f = -320.8339859
lambda = 0.00045725; f = 1126.6655818
lambda = 0.00015242; f = -396.6079721
Norm of dx 515
----
Improvement on iteration 21 = 12.902205390
-----------------
f at the beginning of new iteration, -396.6079721492
Predicted improvement: 33.886589090
lambda = 1; f = -396.6073481
lambda = 0.33333; f = -396.6079234
lambda = 0.11111; f = -396.6079712
lambda = 0.037037; f = -398.9886091
lambda = 0.071599; f = -400.9893223
lambda = 0.13841; f = -396.6079692
lambda = 0.0932; f = -396.6079719
lambda = 0.073512; f = -401.0942221
lambda = 0.084761; f = -396.6079721
lambda = 0.07782; f = -396.6079721
lambda = 0.073932; f = -401.1171666
lambda = 0.076241; f = -396.6079721
lambda = 0.074847; f = -396.6079721
lambda = 0.074023; f = -401.1221361
Norm of dx 0.18338
Cliff. Perturbing search direction.
Predicted improvement: 47.603397761
lambda = 1; f = -396.6069968
lambda = 0.33333; f = -396.6078900
lambda = 0.11111; f = -396.6079693
lambda = 0.037037; f = -399.9906457
lambda = 0.071599; f = -396.6079720
lambda = 0.048211; f = -400.9567129
lambda = 0.061123; f = -396.6079721
lambda = 0.053011; f = -401.3643299
lambda = 0.057739; f = -401.7614929
lambda = 0.062888; f = -396.6079721
lambda = 0.059746; f = -401.9288321
lambda = 0.061612; f = -396.6079721
lambda = 0.060486; f = -401.9902977
lambda = 0.061159; f = -396.6079721
Norm of dx 0.19122
Cliff again. Try traversing
Predicted improvement: 49471.186667053
lambda = 1; f = -352.0859696
lambda = 0.33333; f = -392.4880193
lambda = 0.11111; f = -396.3754599
lambda = 0.037037; f = -316.2973912
lambda = 0.012346; f = -391.0635956
lambda = 0.0041152; f = -396.5847599
lambda = 0.0013717; f = 290.1654938
lambda = 0.00045725; f = -284.6307862
lambda = 0.00015242; f = -400.3397547
lambda = 5.0805e-05; f = -400.0078569
Norm of dx 314.55
----
Improvement on iteration 22 = 5.382325564
back and forth on step length never finished
-----------------
f at the beginning of new iteration, -401.9902977135
Predicted improvement: 33.886589090
lambda = 1; f = -401.9895692
lambda = 0.33333; f = -401.9902168
lambda = 0.11111; f = -401.9902888
lambda = 0.037037; f = -401.9902967
lambda = 0.012346; f = -401.9902976
lambda = 0.0041152; f = -401.9902977
lambda = 0.0013717; f = -401.9902977
lambda = 0.00045725; f = -401.9902977
lambda = 0.00015242; f = -401.9988968
Norm of dx 0.18338
----
Improvement on iteration 23 = 0.008599045
-----------------
f at the beginning of new iteration, -401.9988967588
Predicted improvement: 23.773983614
lambda = 1; f = -401.9981442
lambda = 0.33333; f = -401.9988132
lambda = 0.11111; f = -401.9988875
lambda = 0.037037; f = -401.9988957
lambda = 0.012346; f = -401.9988966
lambda = 0.0041152; f = -401.9988967
lambda = 0.0013717; f = -401.9988968
lambda = 0.00045725; f = -401.9988968
lambda = 0.00015242; f = -401.9988968
lambda = 5.0805e-05; f = -402.0013131
lambda = 9.8216e-05; f = -401.9988968
lambda = 6.6133e-05; f = -402.0020421
lambda = 8.3845e-05; f = -402.0028844
lambda = 0.0001063; f = -401.9988968
lambda = 9.2193e-05; f = -401.9988968
lambda = 8.4644e-05; f = -402.0029224
lambda = 8.9096e-05; f = -401.9988968
lambda = 8.6398e-05; f = -402.0030058
lambda = 8.8007e-05; f = -401.9988968
lambda = 8.7038e-05; f = -402.0030362
Norm of dx 0.16037
Cliff. Perturbing search direction.
Predicted improvement: 41.206862133
lambda = 1; f = -401.9978971
lambda = 0.33333; f = -401.9987857
lambda = 0.11111; f = -401.9988844
lambda = 0.037037; f = -401.9988954
lambda = 0.012346; f = -401.9988966
lambda = 0.0041152; f = -401.9988967
lambda = 0.0013717; f = -401.9988968
lambda = 0.00045725; f = -401.9988968
lambda = 0.00015242; f = -401.9988968
lambda = 5.0805e-05; f = -402.0030848
lambda = 9.8216e-05; f = -401.9988968
lambda = 6.6133e-05; f = -402.0043482
lambda = 8.3845e-05; f = -401.9988968
lambda = 7.2718e-05; f = -402.0048909
lambda = 7.9203e-05; f = -401.9988968
lambda = 7.5246e-05; f = -402.0050992
lambda = 7.7596e-05; f = -401.9988968
lambda = 7.6177e-05; f = -401.9988968
lambda = 7.5338e-05; f = -402.0051069
Norm of dx 0.16687
Cliff again. Try traversing
Predicted improvement: 385696.893418075
lambda = 1; f = -22.2943827
lambda = 0.33333; f = -362.4391340
lambda = 0.11111; f = -398.3816762
lambda = 0.037037; f = -401.8059001
lambda = 0.012346; f = -346.8757554
lambda = 0.0041152; f = -398.6311532
lambda = 0.0013717; f = -10.9502383
lambda = 0.00045725; f = -393.2537030
lambda = 0.00015242; f = -447.1718850
Norm of dx 878.29
----
Improvement on iteration 24 = 45.172988207
-----------------
f at the beginning of new iteration, -447.1718849659
Predicted improvement: 10.856354821
lambda = 1; f = -447.1712678
lambda = 0.33333; f = -447.1718200
lambda = 0.11111; f = -447.1718789
lambda = 0.037037; f = -447.1718846
lambda = 0.012346; f = -447.4394525
lambda = 0.023866; f = -447.1718849
lambda = 0.01607; f = -447.1718850
lambda = 0.012676; f = -447.4465877
lambda = 0.014615; f = -447.1718850
lambda = 0.013418; f = -447.1718850
lambda = 0.012748; f = -447.4481537
lambda = 0.013146; f = -447.1718850
lambda = 0.012906; f = -447.4515667
lambda = 0.013049; f = -447.1718850
Norm of dx 0.075518
Cliff. Perturbing search direction.
Predicted improvement: 16.384388953
lambda = 1; f = -447.1709263
lambda = 0.33333; f = -447.1717829
lambda = 0.11111; f = -447.1718751
lambda = 0.037037; f = -447.1718843
lambda = 0.012346; f = -447.1718850
lambda = 0.0041152; f = -447.3066540
lambda = 0.0079555; f = -447.4322694
lambda = 0.015379; f = -447.1718849
lambda = 0.010356; f = -447.5107052
lambda = 0.013129; f = -447.1718850
lambda = 0.011387; f = -447.1718850
lambda = 0.010454; f = -447.5139315
lambda = 0.011004; f = -447.1718850
lambda = 0.010671; f = -447.1718850
lambda = 0.010476; f = -447.5146324
lambda = 0.010592; f = -447.1718850
Norm of dx 0.10392
Cliff again. Try traversing
Predicted improvement: 64736.156085619
lambda = 1; f = 74221.4579827
lambda = 0.33333; f = 7710.9924626
lambda = 0.11111; f = 413.9496335
lambda = 0.037037; f = -365.8299380
lambda = 0.012346; f = -442.1377909
lambda = 0.0041152; f = -447.1712846
lambda = 0.0013717; f = -447.1718218
lambda = 0.00045725; f = -447.1718791
lambda = 0.00015242; f = -447.1718846
lambda = 5.0805e-05; f = -452.8861803
lambda = 9.8216e-05; f = -447.1718849
lambda = 6.6133e-05; f = -447.1718850
lambda = 5.2163e-05; f = -453.0106175
lambda = 6.0144e-05; f = -447.1718850
lambda = 5.522e-05; f = -447.1718850
lambda = 5.2461e-05; f = -453.0376904
lambda = 5.4099e-05; f = -453.1850174
lambda = 5.5789e-05; f = -447.1718850
lambda = 5.4769e-05; f = -447.1718850
lambda = 5.4166e-05; f = -453.1909489
Norm of dx 359.82
----
Improvement on iteration 25 = 6.019063916
back and forth on step length never finished
-----------------
f at the beginning of new iteration, -453.1909488818
Predicted improvement: 10.856354821
lambda = 1; f = -453.1903154
lambda = 0.33333; f = -453.1908785
lambda = 0.11111; f = -453.1909411
lambda = 0.037037; f = -453.1909480
lambda = 0.012346; f = -453.1909488
lambda = 0.0041152; f = -453.1909489
lambda = 0.0013717; f = -453.1909489
lambda = 0.00045725; f = -453.1909489
lambda = 0.00015242; f = -453.1909489
lambda = 5.0805e-05; f = -453.1909489
lambda = 1.6935e-05; f = -453.1913252
Norm of dx 0.075518
----
Improvement on iteration 26 = 0.000376329
-----------------
f at the beginning of new iteration, -453.1913252109
Predicted improvement: 38.323576743
lambda = 1; f = -455.6568683
lambda = 0.33333; f = -456.4816987
lambda = 0.11111; f = -454.5609733
lambda = 0.037037; f = -453.6779177
lambda = 0.012346; f = -453.3568479
lambda = 0.0041152; f = -453.2468684
lambda = 0.0013717; f = -453.2098806
lambda = 0.00045725; f = -453.1975149
lambda = 0.00015242; f = -453.1933889
lambda = 5.0805e-05; f = -453.1920132
lambda = 1.6935e-05; f = -453.1915545
lambda = 5.645e-06; f = -453.1914017
lambda = 1.8817e-06; f = -453.1913507
lambda = 6.2723e-07; f = -453.1913337
lambda = 2.0908e-07; f = -453.1913280
lambda = 6.9692e-08; f = -453.1913262
lambda = 2.3231e-08; f = -453.1913255
lambda = 7.7435e-09; f = -453.1913253
lambda = 2.5812e-09; f = -453.1913252
Norm of dx 0.16592
----
Improvement on iteration 27 = 3.290373525
smallest step still improving too slow
-----------------
f at the beginning of new iteration, -456.4816987362
Predicted improvement: 12.935699230
lambda = 1; f = -476.9437416
lambda = 1.9332; f = -486.9000573
Norm of dx 0.073
----
Improvement on iteration 28 = 30.418358598
-----------------
f at the beginning of new iteration, -486.9000573344
Predicted improvement: 35.049024831
lambda = 1; f = -520.8478928
Norm of dx 0.28878
----
Improvement on iteration 29 = 33.947835511
-----------------
f at the beginning of new iteration, -520.8478928456
Predicted improvement: 12.294089648
lambda = 1; f = -540.7588411
lambda = 1.9332; f = -547.0314156
Norm of dx 0.15293
----
Improvement on iteration 30 = 26.183522747
-----------------
f at the beginning of new iteration, -547.0314155925
Predicted improvement: 12.842731818
lambda = 1; f = -562.5312316
Norm of dx 0.15526
----
Improvement on iteration 31 = 15.499816009
-----------------
f at the beginning of new iteration, -562.5312316012
Predicted improvement: 2.104971081
lambda = 1; f = -563.9456977
Norm of dx 0.058303
----
Improvement on iteration 32 = 1.414466096
-----------------
f at the beginning of new iteration, -563.9456976969
Predicted improvement: 1.884065616
lambda = 1; f = -564.5611116
lambda = 0.33333; f = -564.8129398
Norm of dx 0.09338
----
Improvement on iteration 33 = 0.867242055
-----------------
f at the beginning of new iteration, -564.8129397520
Predicted improvement: 0.619846732
lambda = 1; f = -565.7181161
lambda = 1.9332; f = -565.9446275
Norm of dx 0.051399
----
Improvement on iteration 34 = 1.131687738
-----------------
f at the beginning of new iteration, -565.9446274901
Predicted improvement: 0.319652733
lambda = 1; f = -566.5551117
lambda = 1.9332; f = -567.0734738
lambda = 3.7372; f = -565.9446259
lambda = 2.5164; f = -567.3724963
lambda = 3.1904; f = -565.9446275
lambda = 2.767; f = -567.4950895
lambda = 3.0137; f = -567.6123838
lambda = 3.2825; f = -565.9446274
lambda = 3.1185; f = -565.9446275
lambda = 3.0241; f = -567.6172129
lambda = 3.0804; f = -567.6434716
lambda = 3.1377; f = -565.9446275
lambda = 3.1032; f = -567.6540609
Norm of dx 0.022683
Cliff. Perturbing search direction.
Predicted improvement: 0.723764688
lambda = 1; f = -567.2037759
lambda = 1.9332; f = -565.9446273
lambda = 1.3017; f = -567.5054491
lambda = 1.6503; f = -567.8054612
lambda = 2.0923; f = -565.9446266
lambda = 1.8146; f = -567.9281267
lambda = 1.9765; f = -565.9446272
lambda = 1.8777; f = -565.9446275
lambda = 1.8209; f = -565.9446275
Norm of dx 0.037951
Cliff again. Try traversing
Predicted improvement: 38.382125640
lambda = 1; f = -554.6816588
lambda = 0.33333; f = -565.6110392
lambda = 0.11111; f = 325.3949289
lambda = 0.037037; f = -462.0555153
lambda = 0.012346; f = -562.9808494
lambda = 0.0041152; f = -566.2097390
lambda = 0.0079555; f = -565.4399110
lambda = 0.0053568; f = -566.0808274
Norm of dx 8.7615
----
Improvement on iteration 35 = 1.983499179
back and forth on step length never finished
-----------------
f at the beginning of new iteration, -567.9281266688
Predicted improvement: 0.319652733
lambda = 1; f = -567.9281226
lambda = 0.33333; f = -567.9281262
lambda = 0.11111; f = -567.9281266
lambda = 0.037037; f = -567.9281267
lambda = 0.012346; f = -567.9281267
lambda = 0.0041152; f = -567.9299858
Norm of dx 0.022683
----
Improvement on iteration 36 = 0.001859084
-----------------
f at the beginning of new iteration, -567.9299857525
Predicted improvement: 0.578851073
lambda = 1; f = -567.9299739
lambda = 0.33333; f = -567.9299845
lambda = 0.11111; f = -567.9299856
lambda = 0.037037; f = -567.9299857
lambda = 0.012346; f = -567.9299858
lambda = 0.0041152; f = -567.9299858
lambda = 0.0013717; f = -567.9315737
lambda = 0.0026518; f = -567.9330543
lambda = 0.0051265; f = -567.9299858
lambda = 0.0034519; f = -567.9339791
lambda = 0.0043764; f = -567.9299858
lambda = 0.0037956; f = -567.9299858
lambda = 0.0034848; f = -567.9340172
lambda = 0.0036681; f = -567.9299858
lambda = 0.003557; f = -567.9299858
lambda = 0.0034919; f = -567.9340255
lambda = 0.0035308; f = -567.9299858
Norm of dx 0.073708
Cliff. Perturbing search direction.
Predicted improvement: 0.944925245
lambda = 1; f = -567.9299668
lambda = 0.33333; f = -567.9299837
lambda = 0.11111; f = -567.9299855
lambda = 0.037037; f = -567.9299857
lambda = 0.012346; f = -567.9299858
lambda = 0.0041152; f = -567.9299858
lambda = 0.0013717; f = -567.9325781
lambda = 0.0026518; f = -567.9349958
lambda = 0.0051265; f = -567.9299858
lambda = 0.0034519; f = -567.9299858
lambda = 0.0027227; f = -567.9351295
lambda = 0.0031393; f = -567.9299858
lambda = 0.0028822; f = -567.9299858
lambda = 0.0027382; f = -567.9351589
lambda = 0.0028237; f = -567.9299858
lambda = 0.0027721; f = -567.9352229
lambda = 0.002803; f = -567.9299858
Norm of dx 0.11701
Cliff again. Try traversing
Predicted improvement: 132.607511485
lambda = 1; f = -567.9285197
lambda = 0.33333; f = -550.0353293
lambda = 0.11111; f = -567.5043327
lambda = 0.037037; f = -567.9299857
lambda = 0.012346; f = -567.9299857
lambda = 0.0041152; f = -567.9299858
lambda = 0.0013717; f = -568.2385571
lambda = 0.0026518; f = -567.9299858
lambda = 0.0017856; f = -567.9299858
lambda = 0.0014084; f = -568.2453402
lambda = 0.0016239; f = -567.9299858
lambda = 0.0014909; f = -567.9299858
lambda = 0.0014164; f = -568.2468191
lambda = 0.0014607; f = -568.2548891
lambda = 0.0015063; f = -567.9299858
lambda = 0.0014788; f = -567.9299858
lambda = 0.0014625; f = -567.9299858
Norm of dx 16.285
----
Improvement on iteration 37 = 0.324903315
back and forth on step length never finished
-----------------
f at the beginning of new iteration, -568.2548890676
Predicted improvement: 0.578851073
lambda = 1; f = -568.2548771
lambda = 0.33333; f = -568.2548877
lambda = 0.11111; f = -568.2548889
lambda = 0.037037; f = -568.2548891
lambda = 0.012346; f = -568.2548891
lambda = 0.0041152; f = -568.2548891
lambda = 0.0013717; f = -568.2548891
lambda = 0.00045725; f = -568.2548891
lambda = 0.00015242; f = -568.2548891
lambda = 5.0805e-05; f = -568.2548891
Norm of dx 0.073708
----
Improvement on iteration 38 = 0.000000000
improvement < crit termination
Final value of minus the log posterior (or likelihood):-568.254889
POSTERIOR KERNEL OPTIMIZATION PROBLEM!
(minus) the hessian matrix at the "mode" is not positive definite!
=> posterior variance of the estimated parameters are not positive.
You should try to change the initial values of the parameters using
the estimated_params_init block, or use another optimization routine.
[Warning: The results below are most likely wrong!]
[> In dynare_estimation_1 (line 324)
In dynare_estimation (line 105)
In nkpc_new.driver (line 575)
In dynare (line 293)]
[Warning: Matrix is singular, close to singular or badly scaled. Results may be inaccurate. RCOND = NaN.]
[> In dynare_estimation_1 (line 347)
In dynare_estimation (line 105)
In nkpc_new.driver (line 575)
In dynare (line 293)
]
RESULTS FROM MAXIMUM LIKELIHOOD ESTIMATION
parameters
Estimate s.d. t-stat
gamma_p 0.2839 NaN NaN
kappa_p 0.3086 NaN NaN
gamma_w 1.0000 NaN NaN
alpha_w -0.0092 NaN NaN
standard deviation of shocks
Estimate s.d. t-stat
e_pi_p 1.7371 NaN NaN
e_pi_w 0.9645 NaN NaN
e_mu 0.0062 NaN NaN
e_s 0.0253 NaN NaN
e_g 0.0839 NaN NaN
e_U 0.0432 NaN NaN
e_x 0.4011 NaN NaN
e_r 0.2757 NaN NaN
e_pi_p_bar 0.1280 NaN NaN
e_r_bar 0.1514 NaN NaN
correlation of shocks
Estimate s.d. t-stat
e_pi_p,e_pi_w -0.1130 NaN NaN
Total computing time : 0h00m07s
Note: warning(s) encountered in MATLAB/Octave code