NRG Ljubljana 2.3.20 - (c) Rok Zitko <rok.zitko@ijs.si>
$Id: nrg.cc,v 1.31 2012/10/26 10:18:22 rokzitko Exp rokzitko $
Timestamp: Fri May 31 11:54:49 2013
Compiled on Jun  5 2013 at 18:59:39

No MPI: single node calculation.

workdir=.
Tmin=1e-10 -> Nmax=67
TEMPERATURE T=1e-10
Nmax=67
Tmin=1e-10
Ninit=0
Lambda=2
strategy=kept
lastall=1
lastalloverride=0
z=0.25
discretization=Z
band=flat
betabar=1
T=1e-10
T_ratio=1
Tmin_ratio=-999
keep=10000
keepenergy=10
keepmin=0
keepblock=0
keepzrescale=0
log=
logall=0
logenumber=10
smooth=new
alpha=0.3
omega0=1e-99
omega0_ratio=1
brcut=0
loggauss_b=0.3
gauss_b=0.1
limitLL=0
eta=0
finite=0
bins=100
discard_trim=1e-16
discard_immediately=1e-16
NN1=0
NN2even=1
NN2avg=0
NNtanh=0
goodE=2
savebins=0
broaden=1
broaden_max=0.1
broaden_min=1e-08
broaden_min_ratio=3
broaden_ratio=1.02
dumpenergies=0
saveall=0
dumpannotated=0
dumpscaled=1
dumpabs=0
dumpprecision=8
dumpgroups=1
grouptol=1e-06
dumpdiagonal=0
ops=A_d
specs=
specd=A_d-A_d
spect=
specb=
reim=0
trsq=
specgt=
gtp=0.7
speci1t=
speci2t=
speckubo=
kubop=0.7
kubor=4
kubob=2
cond=
specchit=
chitp=1
chitp_ratio=-999
dm=1
T0opt=-1
dmnrg=0
cfs=0
fdm=1
fdmexpv=0
fdmexpvn=0
mats=0
nromegan=100
width_td=12
width_custom=12
prec_td=6
prec_custom=6
prec_xy=6
safeguard=1e-05
safeguardmax=200
fixeps=1e-15
diag=dsyev
diagratio=1
dsyevrlimit=100
restart=1
restartfactor=2
checkdiag=0
checkrho=0
calc0=1
tdht=0
spin=2
tri=old
preccpp=2000
checksumrules=0
showEs=0
globalB=0
globalBx=0
polarized=0
pol2x2=0
mpi=0
saveF=0
trunc=
noimag=1
dumpsubspaces=0
done=0
resume=0
forcestop=-1
stopafter=
removefiles=1
store=one

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QS
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 		 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=6
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||||> (z)

# $Id: nrg.cc,v 1.31 2012/10/26 10:18:22 rokzitko Exp rokzitko $ QS
# disk=Z band=flat Lambda=2 z=0.25
# keep=10000 keepenergy=10 Nmax=67 channels=1 dots=1 betabar=1

Number of invariant subspaces: 6
(-2 1) 1 states
0.0309061471093 
(-1 2) 2 states
0.0146705149772 0.0442274543445 
(0 1) 3 states
0 0.0309061471093 0.0588979693217 
(0 3) 1 states
0.0279918222124 
(1 2) 2 states
0.0146705149772 0.0442274543445 
(2 1) 1 states
0.0309061471093 
Number of states (multiplicity taken into account): 16

Computing the following operators:
s=[]
p=[]
g=[]
d=[A_d ]
v=[]
t=[]

Computing the following spectra:

Before NRG iteration (N=-1)

 ***** [NRG] Iteration 1/67 (scale 1.2131570881878404) ***** 
[1] xi(0)=0.473132378817 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=1.89846285841 nrkept=35 nrkeptmult=64 
Total energy=-1.14906683667  Egs=-0.9075842757
Kept: 35 out of 35, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/67 (scale 0.85783160370214839) ***** 
[1] xi(1)=0.631518937152 zeta(2)=0
Stats: nr=15 min=1 max=20 sum3=30066
Emax=3.76220791894 nrkept=126 nrkeptmult=256 
Total energy=-1.61117042628  Egs=-0.538687998459
Kept: 126 out of 126, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/67 (scale 0.6065785440939202) ***** 
[1] xi(2)=0.974788923894 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=8.48694522066 nrkept=462 nrkeptmult=1024 
Total energy=-2.57149943788  Egs=-1.58318987861
Kept: 462 out of 462, ratio=1
Memory used: 5 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/67 (scale 0.4289158018510742) ***** 
[1] xi(3)=1.2027371305 zeta(4)=0
Stats: nr=28 min=1 max=210 sum3=37226316
Safeguard: keep additional 1 states
Emax=10.3163498334 nrkept=1468 nrkeptmult=3540 
Total energy=-3.05731539259  Egs=-1.1326604257
Kept: 1468 out of 1716, ratio=0.855
Memory used: 14 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/67 (scale 0.3032892720469601) ***** 
[1] xi(4)=1.12253968187 zeta(5)=0
Stats: nr=36 min=1 max=655 sum3=887844448
Emax=10.3386102527 nrkept=2232 nrkeptmult=5614 
Total energy=-3.55007274671  Egs=-1.62471079441
Kept: 2232 out of 5512, ratio=0.405
Memory used: 27 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 6/67 (scale 0.2144579009255371) ***** 
[1] xi(5)=1.03243378845 zeta(6)=0
Stats: nr=42 min=2 max=841 sum3=2779096102
Safeguard: keep additional 5 states
Emax=10.0287601055 nrkept=2242 nrkeptmult=5560 
Total energy=-3.74365639227  Egs=-0.902665020596
Kept: 2242 out of 8428, ratio=0.266
Memory used: 32 MB  Time elapsed: 10 s

 ***** [NRG] Iteration 7/67 (scale 0.15164463602348005) ***** 
[1] xi(6)=1.00829687067 zeta(7)=0
Stats: nr=44 min=1 max=963 sum3=2974659303
Safeguard: keep additional 3 states
Emax=10.0105258789 nrkept=2066 nrkeptmult=5134 
Total energy=-3.97203467541  Egs=-1.50600963631
Kept: 2066 out of 8445, ratio=0.245
Memory used: 31 MB  Time elapsed: 18 s

 ***** [NRG] Iteration 8/67 (scale 0.10722895046276855) ***** 
[1] xi(7)=1.00300214986 zeta(8)=0
Stats: nr=42 min=1 max=791 sum3=2303803142
Safeguard: keep additional 2 states
Emax=10.036799574 nrkept=2121 nrkeptmult=5273 
Total energy=-4.06890319315  Egs=-0.903380265514
Kept: 2121 out of 7778, ratio=0.273
Memory used: 32 MB  Time elapsed: 25 s

 ***** [NRG] Iteration 9/67 (scale 0.075822318011740025) ***** 
[1] xi(8)=1.00131950393 zeta(9)=0
Stats: nr=44 min=1 max=911 sum3=2522727317
Safeguard: keep additional 1 states
Emax=10.0016424925 nrkept=1688 nrkeptmult=4186 
Total energy=-4.18036348403  Egs=-1.47001956414
Kept: 1688 out of 7997, ratio=0.211
Memory used: 29 MB  Time elapsed: 32 s

 ***** [NRG] Iteration 10/67 (scale 0.053614475231384275) ***** 
[1] xi(9)=1.00062206904 zeta(10)=0
Stats: nr=42 min=1 max=663 sum3=1318309146
Emax=10.0184681848 nrkept=1873 nrkeptmult=4643 
Total energy=-4.23013103141  Egs=-0.928248335288
Kept: 1873 out of 6372, ratio=0.294
Memory used: 31 MB  Time elapsed: 37 s

 ***** [NRG] Iteration 11/67 (scale 0.037911159005870013) ***** 
[1] xi(10)=1.00030223414 zeta(11)=0
Stats: nr=42 min=1 max=815 sum3=1780810507
Safeguard: keep additional 3 states
Emax=10.0197753408 nrkept=1654 nrkeptmult=4108 
Total energy=-4.28449472914  Egs=-1.43397614715
Kept: 1654 out of 7057, ratio=0.234
Memory used: 31 MB  Time elapsed: 42 s

 ***** [NRG] Iteration 12/67 (scale 0.026807237615692137) ***** 
[1] xi(11)=1.00014898132 zeta(12)=0
Stats: nr=42 min=1 max=646 sum3=1226297206
Safeguard: keep additional 1 states
Emax=10.000284816 nrkept=1794 nrkeptmult=4452 
Total energy=-4.31015161874  Egs=-0.957088155644
Kept: 1794 out of 6244, ratio=0.287
Memory used: 32 MB  Time elapsed: 46 s

 ***** [NRG] Iteration 13/67 (scale 0.018955579502935006) ***** 
[1] xi(12)=1.00007396421 zeta(13)=0
Stats: nr=42 min=1 max=780 sum3=1561236113
Emax=10.0000976608 nrkept=1555 nrkeptmult=3830 
Total energy=-4.33662645013  Egs=-1.39667749985
Kept: 1555 out of 6761, ratio=0.23
Memory used: 31 MB  Time elapsed: 51 s

 ***** [NRG] Iteration 14/67 (scale 0.013403618807846069) ***** 
[1] xi(13)=1.0000368514 zeta(14)=0
Stats: nr=42 min=1 max=609 sum3=1022420756
Emax=10.0001865563 nrkept=1644 nrkeptmult=4034 
Total energy=-4.34984357364  Egs=-0.986086197741
Kept: 1644 out of 5858, ratio=0.281
Memory used: 31 MB  Time elapsed: 54 s

 ***** [NRG] Iteration 15/67 (scale 0.0094777897514675032) ***** 
[1] xi(14)=1.00001839314 zeta(15)=0
Stats: nr=40 min=1 max=722 sum3=1223923710
Emax=10.0202490375 nrkept=1539 nrkeptmult=3784 
Total energy=-4.36273121425  Egs=-1.35977278997
Kept: 1539 out of 6180, ratio=0.249
Memory used: 30 MB  Time elapsed: 58 s

 ***** [NRG] Iteration 16/67 (scale 0.0067018094039230343) ***** 
[1] xi(15)=1.00000918844 zeta(16)=0
Stats: nr=42 min=1 max=602 sum3=986928322
Safeguard: keep additional 1 states
Emax=10.0030354167 nrkept=1494 nrkeptmult=3640 
Total energy=-4.36952735386  Egs=-1.01407533392
Kept: 1494 out of 5794, ratio=0.258
Memory used: 31 MB  Time elapsed: 62 s

 ***** [NRG] Iteration 17/67 (scale 0.0047388948757337516) ***** 
[1] xi(16)=1.00000459219 zeta(17)=0
Stats: nr=36 min=1 max=670 sum3=955230795
Safeguard: keep additional 3 states
Emax=10.0052483164 nrkept=1406 nrkeptmult=3476 
Total energy=-4.37579829384  Egs=-1.32329164125
Kept: 1406 out of 5625, ratio=0.25
Memory used: 30 MB  Time elapsed: 65 s

 ***** [NRG] Iteration 18/67 (scale 0.0033509047019615172) ***** 
[1] xi(17)=1.00000229559 zeta(18)=0
Stats: nr=40 min=1 max=551 sum3=765342560
Safeguard: keep additional 4 states
Emax=10.0073439161 nrkept=1334 nrkeptmult=3250 
Total energy=-4.37929311978  Egs=-1.04294996429
Kept: 1334 out of 5288, ratio=0.252
Memory used: 31 MB  Time elapsed: 68 s

 ***** [NRG] Iteration 19/67 (scale 0.0023694474378668758) ***** 
[1] xi(18)=1.00000114767 zeta(19)=0
Stats: nr=34 min=1 max=602 sum3=692347587
Safeguard: keep additional 2 states
Emax=10.0012259487 nrkept=1213 nrkeptmult=3004 
Total energy=-4.38233760343  Egs=-1.28489182941
Kept: 1213 out of 5019, ratio=0.242
Memory used: 31 MB  Time elapsed: 70 s

 ***** [NRG] Iteration 20/67 (scale 0.0016754523509807586) ***** 
[1] xi(19)=1.0000005738 zeta(20)=0
Stats: nr=40 min=1 max=482 sum3=507819832
Safeguard: keep additional 2 states
Emax=10.0056127819 nrkept=1120 nrkeptmult=2784 
Total energy=-4.3841411544  Egs=-1.07645613764
Kept: 1120 out of 4570, ratio=0.245
Memory used: 30 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 21/67 (scale 0.0011847237189334379) ***** 
[1] xi(20)=1.00000028689 zeta(21)=0
Stats: nr=34 min=1 max=511 sum3=422866429
Safeguard: keep additional 3 states
Emax=10.0081170502 nrkept=1003 nrkeptmult=2506 
Total energy=-4.3856116904  Egs=-1.24124804533
Kept: 1003 out of 4225, ratio=0.237
Memory used: 28 MB  Time elapsed: 73 s

 ***** [NRG] Iteration 22/67 (scale 0.00083772617549037929) ***** 
[1] xi(21)=1.00000014344 zeta(22)=0
Stats: nr=40 min=1 max=405 sum3=297823712
Safeguard: keep additional 4 states
Emax=10.0059396163 nrkept=927 nrkeptmult=2297 
Total energy=-4.38654836372  Egs=-1.11811395226
Kept: 927 out of 3782, ratio=0.245
Memory used: 32 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 23/67 (scale 0.00059236185946671895) ***** 
[1] xi(22)=1.00000007172 zeta(23)=0
Stats: nr=34 min=1 max=429 sum3=245882912
Safeguard: keep additional 2 states
Emax=10.0234973055 nrkept=883 nrkeptmult=2252 
Total energy=-4.38725308145  Egs=-1.18967437867
Kept: 883 out of 3500, ratio=0.252
Memory used: 32 MB  Time elapsed: 75 s

 ***** [NRG] Iteration 24/67 (scale 0.00041886308774518965) ***** 
[1] xi(23)=1.00000003586 zeta(24)=0
Stats: nr=40 min=1 max=363 sum3=205107938
Safeguard: keep additional 2 states
Emax=10.0031019786 nrkept=829 nrkeptmult=2051 
Total energy=-4.38774334409  Egs=-1.17046035829
Kept: 829 out of 3344, ratio=0.248
Memory used: 31 MB  Time elapsed: 76 s

 ***** [NRG] Iteration 25/67 (scale 0.00029618092973335947) ***** 
[1] xi(24)=1.00000001793 zeta(25)=0
Stats: nr=34 min=1 max=385 sum3=175545704
Safeguard: keep additional 2 states
Emax=10.0096480061 nrkept=827 nrkeptmult=2108 
Total energy=-4.38807772964  Egs=-1.12899081297
Kept: 827 out of 3122, ratio=0.265
Memory used: 32 MB  Time elapsed: 77 s

 ***** [NRG] Iteration 26/67 (scale 0.00020943154387259482) ***** 
[1] xi(25)=1.00000000897 zeta(26)=0
Stats: nr=40 min=1 max=337 sum3=166805448
Safeguard: keep additional 1 states
Emax=10.0049168764 nrkept=786 nrkeptmult=1962 
Total energy=-4.38833633899  Egs=-1.23481565453
Kept: 786 out of 3126, ratio=0.251
Memory used: 31 MB  Time elapsed: 78 s

 ***** [NRG] Iteration 27/67 (scale 0.00014809046486667974) ***** 
[1] xi(26)=1.00000000448 zeta(27)=0
Stats: nr=34 min=1 max=362 sum3=147227213
Safeguard: keep additional 2 states
Emax=10.0107454458 nrkept=805 nrkeptmult=2038 
Total energy=-4.38849339291  Egs=-1.06052693019
Kept: 805 out of 2963, ratio=0.272
Memory used: 31 MB  Time elapsed: 78 s

 ***** [NRG] Iteration 28/67 (scale 0.00010471577193629741) ***** 
[1] xi(27)=1.00000000224 zeta(28)=0
Stats: nr=38 min=1 max=330 sum3=155870198
Safeguard: keep additional 2 states
Emax=10.0240535905 nrkept=763 nrkeptmult=1925 
Total energy=-4.38863056571  Egs=-1.30995354964
Kept: 763 out of 3038, ratio=0.251
Memory used: 32 MB  Time elapsed: 79 s

 ***** [NRG] Iteration 29/67 (scale 7.4045232433339868e-05) ***** 
[1] xi(28)=1.00000000112 zeta(29)=0
Stats: nr=34 min=1 max=349 sum3=134658679
Safeguard: keep additional 3 states
Emax=10.0120592259 nrkept=776 nrkeptmult=1946 
Total energy=-4.388703841  Egs=-0.989601684964
Kept: 776 out of 2881, ratio=0.269
Memory used: 33 MB  Time elapsed: 79 s

 ***** [NRG] Iteration 30/67 (scale 5.2357885968148706e-05) ***** 
[1] xi(29)=1.00000000056 zeta(30)=0
Stats: nr=38 min=1 max=323 sum3=143658912
Safeguard: keep additional 1 states
Emax=10.0036920428 nrkept=723 nrkeptmult=1799 
Total energy=-4.38877660696  Egs=-1.38978046655
Kept: 723 out of 2928, ratio=0.247
Memory used: 33 MB  Time elapsed: 80 s

 ***** [NRG] Iteration 31/67 (scale 3.7022616216669934e-05) ***** 
[1] xi(30)=1.00000000028 zeta(31)=0
Stats: nr=33 min=1 max=334 sum3=116636744
Safeguard: keep additional 1 states
Emax=10.0076247786 nrkept=728 nrkeptmult=1800 
Total energy=-4.38881086325  Egs=-0.925280021212
Kept: 728 out of 2726, ratio=0.267
Memory used: 33 MB  Time elapsed: 81 s

 ***** [NRG] Iteration 32/67 (scale 2.6178942984074353e-05) ***** 
[1] xi(31)=1.00000000014 zeta(32)=0
Stats: nr=36 min=1 max=304 sum3=122567426
Safeguard: keep additional 1 states
Emax=10.0066276225 nrkept=658 nrkeptmult=1614 
Total energy=-4.38884916803  Egs=-1.46319053248
Kept: 658 out of 2744, ratio=0.24
Memory used: 33 MB  Time elapsed: 81 s

 ***** [NRG] Iteration 33/67 (scale 1.8511308108334967e-05) ***** 
[1] xi(32)=1.00000000007 zeta(33)=0
Stats: nr=33 min=1 max=311 sum3=91943565
Safeguard: keep additional 1 states
Emax=10.0034905278 nrkept=690 nrkeptmult=1694 
Total energy=-4.38886537884  Egs=-0.875724731773
Kept: 690 out of 2481, ratio=0.278
Memory used: 34 MB  Time elapsed: 82 s

 ***** [NRG] Iteration 34/67 (scale 1.3089471492037176e-05) ***** 
[1] xi(33)=1.00000000004 zeta(34)=0
Stats: nr=36 min=1 max=287 sum3=105360776
Safeguard: keep additional 2 states
Emax=10.0144525693 nrkept=607 nrkeptmult=1461 
Total energy=-4.38888527317  Egs=-1.51987282395
Kept: 607 out of 2600, ratio=0.233
Memory used: 33 MB  Time elapsed: 82 s

 ***** [NRG] Iteration 35/67 (scale 9.2556540541674836e-06) ***** 
[1] xi(34)=1.00000000002 zeta(35)=0
Stats: nr=33 min=1 max=292 sum3=74039870
Safeguard: keep additional 3 states
Emax=10.0116782702 nrkept=619 nrkeptmult=1500 
Total energy=-4.38889307567  Egs=-0.842997606196
Kept: 619 out of 2282, ratio=0.271
Memory used: 33 MB  Time elapsed: 83 s

 ***** [NRG] Iteration 36/67 (scale 6.5447357460185882e-06) ***** 
[1] xi(35)=1.00000000001 zeta(36)=0
Stats: nr=36 min=1 max=260 sum3=76488430
Safeguard: keep additional 2 states
Emax=10.0253498898 nrkept=552 nrkeptmult=1306 
Total energy=-4.38890326636  Egs=-1.55708156998
Kept: 552 out of 2326, ratio=0.237
Memory used: 34 MB  Time elapsed: 83 s

 ***** [NRG] Iteration 37/67 (scale 4.6278270270837418e-06) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=33 min=1 max=269 sum3=57170652
Emax=10.0556843894 nrkept=570 nrkeptmult=1334 
Total energy=-4.38890707848  Egs=-0.823739620124
Kept: 570 out of 2070, ratio=0.275
Memory used: 34 MB  Time elapsed: 83 s

 ***** [NRG] Iteration 38/67 (scale 3.2723678730092941e-06) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=34 min=1 max=245 sum3=62091594
Emax=10.0036967653 nrkept=512 nrkeptmult=1190 
Total energy=-4.38891224485  Egs=-1.57878543072
Kept: 512 out of 2136, ratio=0.24
Memory used: 34 MB  Time elapsed: 84 s

 ***** [NRG] Iteration 39/67 (scale 2.3139135135418709e-06) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=33 min=1 max=251 sum3=45594660
Emax=10.3327931792 nrkept=567 nrkeptmult=1322 
Total energy=-4.38891412654  Egs=-0.813208616745
Kept: 567 out of 1914, ratio=0.296
Memory used: 34 MB  Time elapsed: 84 s

 ***** [NRG] Iteration 40/67 (scale 1.6361839365046471e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.1296480344 nrkept=499 nrkeptmult=1149 
Total energy=-4.388916729  Egs=-1.59056469057
Kept: 499 out of 2124, ratio=0.235
Memory used: 34 MB  Time elapsed: 84 s

 ***** [NRG] Iteration 41/67 (scale 1.1569567567709354e-06) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.6374354844 nrkept=567 nrkeptmult=1322 
Total energy=-4.38891766346  Egs=-0.807687843727
Kept: 567 out of 1862, ratio=0.305
Memory used: 34 MB  Time elapsed: 84 s

 ***** [NRG] Iteration 42/67 (scale 8.1809196825232353e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.3597623579 nrkept=499 nrkeptmult=1149 
Total energy=-4.38891896971  Egs=-1.59670855676
Kept: 499 out of 2124, ratio=0.235
Memory used: 34 MB  Time elapsed: 85 s

 ***** [NRG] Iteration 43/67 (scale 5.7847837838546772e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.7372964209 nrkept=567 nrkeptmult=1322 
Total energy=-4.38891943531  Egs=-0.804859196512
Kept: 567 out of 1862, ratio=0.305
Memory used: 34 MB  Time elapsed: 85 s

 ***** [NRG] Iteration 44/67 (scale 4.0904598412616176e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.4899287193 nrkept=499 nrkeptmult=1149 
Total energy=-4.38892008972  Egs=-1.59984708178
Kept: 499 out of 2124, ratio=0.235
Memory used: 34 MB  Time elapsed: 85 s

 ***** [NRG] Iteration 45/67 (scale 2.8923918919273386e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.789385392 nrkept=567 nrkeptmult=1322 
Total energy=-4.3889203221  Egs=-0.803427192129
Kept: 567 out of 1862, ratio=0.305
Memory used: 34 MB  Time elapsed: 85 s

 ***** [NRG] Iteration 46/67 (scale 2.0452299206308088e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5111582642 nrkept=498 nrkeptmult=1142 
Total energy=-4.38892064963  Egs=-1.6014333867
Kept: 498 out of 2124, ratio=0.234
Memory used: 34 MB  Time elapsed: 86 s

 ***** [NRG] Iteration 47/67 (scale 1.4461959459636693e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8156857111 nrkept=567 nrkeptmult=1322 
Total energy=-4.38892076572  Egs=-0.802706688359
Kept: 567 out of 1857, ratio=0.305
Memory used: 34 MB  Time elapsed: 86 s

 ***** [NRG] Iteration 48/67 (scale 1.0226149603154044e-07) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5185900319 nrkept=498 nrkeptmult=1142 
Total energy=-4.38892092956  Egs=-1.60223085004
Kept: 498 out of 2124, ratio=0.234
Memory used: 34 MB  Time elapsed: 86 s

 ***** [NRG] Iteration 49/67 (scale 7.2309797298183465e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8286333834 nrkept=567 nrkeptmult=1322 
Total energy=-4.38892098758  Egs=-0.802345300596
Kept: 567 out of 1857, ratio=0.305
Memory used: 34 MB  Time elapsed: 86 s

 ***** [NRG] Iteration 50/67 (scale 5.113074801577022e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5223242957 nrkept=498 nrkeptmult=1142 
Total energy=-4.38892106952  Egs=-1.60263066579
Kept: 498 out of 2124, ratio=0.234
Memory used: 34 MB  Time elapsed: 87 s

 ***** [NRG] Iteration 51/67 (scale 3.6154898649091733e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8333633709 nrkept=567 nrkeptmult=1324 
Total energy=-4.38892109853  Egs=-0.802164321419
Kept: 567 out of 1857, ratio=0.305
Memory used: 35 MB  Time elapsed: 87 s

 ***** [NRG] Iteration 52/67 (scale 2.556537400788511e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5242078641 nrkept=498 nrkeptmult=1142 
Total energy=-4.3889211395  Egs=-1.60283084548
Kept: 498 out of 2124, ratio=0.234
Memory used: 35 MB  Time elapsed: 87 s

 ***** [NRG] Iteration 53/67 (scale 1.8077449324545866e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8355062475 nrkept=567 nrkeptmult=1324 
Total energy=-4.388921154  Egs=-0.802073760341
Kept: 567 out of 1857, ratio=0.305
Memory used: 35 MB  Time elapsed: 87 s

 ***** [NRG] Iteration 54/67 (scale 1.2782687003942555e-08) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5251454635 nrkept=498 nrkeptmult=1142 
Total energy=-4.38892117449  Egs=-1.60293100337
Kept: 498 out of 2124, ratio=0.234
Memory used: 36 MB  Time elapsed: 88 s

 ***** [NRG] Iteration 55/67 (scale 9.0387246622729332e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8365765199 nrkept=567 nrkeptmult=1324 
Total energy=-4.38892118174  Egs=-0.802028461909
Kept: 567 out of 1857, ratio=0.305
Memory used: 36 MB  Time elapsed: 88 s

 ***** [NRG] Iteration 56/67 (scale 6.3913435019712775e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.52561457 nrkept=498 nrkeptmult=1142 
Total energy=-4.38892119199  Egs=-1.60298109934
Kept: 498 out of 2124, ratio=0.234
Memory used: 36 MB  Time elapsed: 88 s

 ***** [NRG] Iteration 57/67 (scale 4.5193623311364666e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8371113572 nrkept=567 nrkeptmult=1324 
Total energy=-4.38892119561  Egs=-0.802005808218
Kept: 567 out of 1857, ratio=0.305
Memory used: 36 MB  Time elapsed: 88 s

 ***** [NRG] Iteration 58/67 (scale 3.1956717509856388e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5258492 nrkept=498 nrkeptmult=1142 
Total energy=-4.38892120073  Egs=-1.60300615159
Kept: 498 out of 2124, ratio=0.234
Memory used: 36 MB  Time elapsed: 89 s

 ***** [NRG] Iteration 59/67 (scale 2.2596811655682333e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8373787018 nrkept=567 nrkeptmult=1324 
Total energy=-4.38892120255  Egs=-0.801994480253
Kept: 567 out of 1857, ratio=0.305
Memory used: 36 MB  Time elapsed: 89 s

 ***** [NRG] Iteration 60/67 (scale 1.5978358754928194e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5259665342 nrkept=498 nrkeptmult=1142 
Total energy=-4.38892120511  Egs=-1.60301867878
Kept: 498 out of 2124, ratio=0.234
Memory used: 36 MB  Time elapsed: 89 s

 ***** [NRG] Iteration 61/67 (scale 1.1298405827841166e-09) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 2 states
Emax=10.8375174929 nrkept=569 nrkeptmult=1330 
Total energy=-4.38892120601  Egs=-0.801988815991
Kept: 569 out of 1857, ratio=0.306
Memory used: 36 MB  Time elapsed: 89 s

 ***** [NRG] Iteration 62/67 (scale 7.9891793774640969e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=34 min=1 max=244 sum3=61697330
Emax=10.5260252061 nrkept=498 nrkeptmult=1142 
Total energy=-4.3889212073  Egs=-1.60302494263
Kept: 498 out of 2132, ratio=0.234
Memory used: 36 MB  Time elapsed: 90 s

 ***** [NRG] Iteration 63/67 (scale 5.6492029139205832e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8375968195 nrkept=575 nrkeptmult=1340 
Total energy=-4.38892120775  Egs=-0.80198598379
Kept: 575 out of 1857, ratio=0.31
Memory used: 36 MB  Time elapsed: 90 s

 ***** [NRG] Iteration 64/67 (scale 3.9945896887320485e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260401124 nrkept=498 nrkeptmult=1142 
Total energy=-4.38892120839  Egs=-1.60302807463
Kept: 498 out of 2154, ratio=0.231
Memory used: 37 MB  Time elapsed: 90 s

 ***** [NRG] Iteration 65/67 (scale 2.8246014569602916e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8376155039 nrkept=575 nrkeptmult=1340 
Total energy=-4.38892120862  Egs=-0.801984567672
Kept: 575 out of 1857, ratio=0.31
Memory used: 36 MB  Time elapsed: 90 s

 ***** [NRG] Iteration 66/67 (scale 1.9972948443660242e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260547818 nrkept=498 nrkeptmult=1142 
Total energy=-4.38892120894  Egs=-1.60302964065
Kept: 498 out of 2154, ratio=0.231
Memory used: 37 MB  Time elapsed: 91 s

 ***** [NRG] Iteration 67/67 (scale 1.4123007284801458e-10) ***** 
[1] xi(66)=1 zeta(67)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 12 states
Emax=10.8376388455 nrkept=1857 nrkeptmult=4568 
Total energy=-4.38892120905  Egs=-0.801983859609
Kept: 1857 out of 1857, ratio=1
Memory used: 39 MB  Time elapsed: 91 s

** Iteration completed.


Total energy:             -4.388921209048739
TKW=2.16354528284e-05
Grand canonical Z=               7.81082756098191
Zft (last shell)=               7.810816704750114
Calculating density matrices.
Calculating full density matrices.

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QS
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 		 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=6
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||||> (z)
Number of invariant subspaces: 6
(-2 1) 1 states
0.0309061471093 
(-1 2) 2 states
0.0146705149772 0.0442274543445 
(0 1) 3 states
0 0.0309061471093 0.0588979693217 
(0 3) 1 states
0.0279918222124 
(1 2) 2 states
0.0146705149772 0.0442274543445 
(2 1) 1 states
0.0309061471093 
Number of states (multiplicity taken into account): 16

Computing the following operators:
s=[]
p=[]
g=[]
d=[A_d ]
v=[]
t=[]

Computing the following spectra:
Spectrum: A_d-A_d FDM

 ***** [DMNRG] Iteration 1/67 (scale 1.2131570881878404) ***** 
[1] xi(0)=0.473132378817 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=1.89846285841 nrkept=35 nrkeptmult=64 
Kept: 35 out of 35, ratio=1
Memory used: 38 MB  Time elapsed: 100 s

 ***** [DMNRG] Iteration 2/67 (scale 0.85783160370214839) ***** 
[1] xi(1)=0.631518937152 zeta(2)=0
Stats: nr=15 min=1 max=20 sum3=30066
Emax=3.76220791894 nrkept=126 nrkeptmult=256 
Kept: 126 out of 126, ratio=1
Memory used: 38 MB  Time elapsed: 100 s

 ***** [DMNRG] Iteration 3/67 (scale 0.6065785440939202) ***** 
[1] xi(2)=0.974788923894 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=8.48694522066 nrkept=462 nrkeptmult=1024 
Kept: 462 out of 462, ratio=1
Memory used: 38 MB  Time elapsed: 100 s

 ***** [DMNRG] Iteration 4/67 (scale 0.4289158018510742) ***** 
[1] xi(3)=1.2027371305 zeta(4)=0
Stats: nr=28 min=1 max=210 sum3=37226316
Safeguard: keep additional 1 states
Emax=10.3163498334 nrkept=1468 nrkeptmult=3540 
Kept: 1468 out of 1716, ratio=0.855
Memory used: 41 MB  Time elapsed: 100 s

 ***** [DMNRG] Iteration 5/67 (scale 0.3032892720469601) ***** 
[1] xi(4)=1.12253968187 zeta(5)=0
Stats: nr=36 min=1 max=655 sum3=887844448
Emax=10.3386102527 nrkept=2232 nrkeptmult=5614 
Kept: 2232 out of 5512, ratio=0.405
Memory used: 53 MB  Time elapsed: 101 s

 ***** [DMNRG] Iteration 6/67 (scale 0.2144579009255371) ***** 
[1] xi(5)=1.03243378845 zeta(6)=0
Stats: nr=42 min=2 max=841 sum3=2779096102
Safeguard: keep additional 5 states
Emax=10.0287601055 nrkept=2242 nrkeptmult=5560 
Kept: 2242 out of 8428, ratio=0.266
Memory used: 56 MB  Time elapsed: 104 s

 ***** [DMNRG] Iteration 7/67 (scale 0.15164463602348005) ***** 
[1] xi(6)=1.00829687067 zeta(7)=0
Stats: nr=44 min=1 max=963 sum3=2974659303
Safeguard: keep additional 3 states
Emax=10.0105258789 nrkept=2066 nrkeptmult=5134 
Kept: 2066 out of 8445, ratio=0.245
Memory used: 55 MB  Time elapsed: 107 s

 ***** [DMNRG] Iteration 8/67 (scale 0.10722895046276855) ***** 
[1] xi(7)=1.00300214986 zeta(8)=0
Stats: nr=42 min=1 max=791 sum3=2303803142
Safeguard: keep additional 2 states
Emax=10.036799574 nrkept=2121 nrkeptmult=5273 
Kept: 2121 out of 7778, ratio=0.273
Memory used: 53 MB  Time elapsed: 110 s

 ***** [DMNRG] Iteration 9/67 (scale 0.075822318011740025) ***** 
[1] xi(8)=1.00131950393 zeta(9)=0
Stats: nr=44 min=1 max=911 sum3=2522727317
Safeguard: keep additional 1 states
Emax=10.0016424925 nrkept=1688 nrkeptmult=4186 
Kept: 1688 out of 7997, ratio=0.211
Memory used: 49 MB  Time elapsed: 112 s

 ***** [DMNRG] Iteration 10/67 (scale 0.053614475231384275) ***** 
[1] xi(9)=1.00062206904 zeta(10)=0
Stats: nr=42 min=1 max=663 sum3=1318309146
Emax=10.0184681848 nrkept=1873 nrkeptmult=4643 
Kept: 1873 out of 6372, ratio=0.294
Memory used: 49 MB  Time elapsed: 114 s

 ***** [DMNRG] Iteration 11/67 (scale 0.037911159005870013) ***** 
[1] xi(10)=1.00030223414 zeta(11)=0
Stats: nr=42 min=1 max=815 sum3=1780810507
Safeguard: keep additional 3 states
Emax=10.0197753408 nrkept=1654 nrkeptmult=4108 
Kept: 1654 out of 7057, ratio=0.234
Memory used: 47 MB  Time elapsed: 116 s

 ***** [DMNRG] Iteration 12/67 (scale 0.026807237615692137) ***** 
[1] xi(11)=1.00014898132 zeta(12)=0
Stats: nr=42 min=1 max=646 sum3=1226297206
Safeguard: keep additional 1 states
Emax=10.000284816 nrkept=1794 nrkeptmult=4452 
Kept: 1794 out of 6244, ratio=0.287
Memory used: 48 MB  Time elapsed: 118 s

 ***** [DMNRG] Iteration 13/67 (scale 0.018955579502935006) ***** 
[1] xi(12)=1.00007396421 zeta(13)=0
Stats: nr=42 min=1 max=780 sum3=1561236113
Emax=10.0000976608 nrkept=1555 nrkeptmult=3830 
Kept: 1555 out of 6761, ratio=0.23
Memory used: 45 MB  Time elapsed: 120 s

 ***** [DMNRG] Iteration 14/67 (scale 0.013403618807846069) ***** 
[1] xi(13)=1.0000368514 zeta(14)=0
Stats: nr=42 min=1 max=609 sum3=1022420756
Emax=10.0001865563 nrkept=1644 nrkeptmult=4034 
Kept: 1644 out of 5858, ratio=0.281
Memory used: 46 MB  Time elapsed: 121 s

 ***** [DMNRG] Iteration 15/67 (scale 0.0094777897514675032) ***** 
[1] xi(14)=1.00001839314 zeta(15)=0
Stats: nr=40 min=1 max=722 sum3=1223923710
Emax=10.0202490375 nrkept=1539 nrkeptmult=3784 
Kept: 1539 out of 6180, ratio=0.249
Memory used: 45 MB  Time elapsed: 123 s

 ***** [DMNRG] Iteration 16/67 (scale 0.0067018094039230343) ***** 
[1] xi(15)=1.00000918844 zeta(16)=0
Stats: nr=42 min=1 max=602 sum3=986928322
Safeguard: keep additional 1 states
Emax=10.0030354167 nrkept=1494 nrkeptmult=3640 
Kept: 1494 out of 5794, ratio=0.258
Memory used: 45 MB  Time elapsed: 124 s

 ***** [DMNRG] Iteration 17/67 (scale 0.0047388948757337516) ***** 
[1] xi(16)=1.00000459219 zeta(17)=0
Stats: nr=36 min=1 max=670 sum3=955230795
Safeguard: keep additional 3 states
Emax=10.0052483164 nrkept=1406 nrkeptmult=3476 
Kept: 1406 out of 5625, ratio=0.25
Memory used: 44 MB  Time elapsed: 126 s

 ***** [DMNRG] Iteration 18/67 (scale 0.0033509047019615172) ***** 
[1] xi(17)=1.00000229559 zeta(18)=0
Stats: nr=40 min=1 max=551 sum3=765342560
Safeguard: keep additional 4 states
Emax=10.0073439161 nrkept=1334 nrkeptmult=3250 
Kept: 1334 out of 5288, ratio=0.252
Memory used: 45 MB  Time elapsed: 127 s

 ***** [DMNRG] Iteration 19/67 (scale 0.0023694474378668758) ***** 
[1] xi(18)=1.00000114767 zeta(19)=0
Stats: nr=34 min=1 max=602 sum3=692347587
Safeguard: keep additional 2 states
Emax=10.0012259487 nrkept=1213 nrkeptmult=3004 
Kept: 1213 out of 5019, ratio=0.242
Memory used: 45 MB  Time elapsed: 128 s

 ***** [DMNRG] Iteration 20/67 (scale 0.0016754523509807586) ***** 
[1] xi(19)=1.0000005738 zeta(20)=0
Stats: nr=40 min=1 max=482 sum3=507819832
Safeguard: keep additional 2 states
Emax=10.0056127819 nrkept=1120 nrkeptmult=2784 
Kept: 1120 out of 4570, ratio=0.245
Memory used: 48 MB  Time elapsed: 129 s

 ***** [DMNRG] Iteration 21/67 (scale 0.0011847237189334379) ***** 
[1] xi(20)=1.00000028689 zeta(21)=0
Stats: nr=34 min=1 max=511 sum3=422866429
Safeguard: keep additional 3 states
Emax=10.0081170502 nrkept=1003 nrkeptmult=2506 
Kept: 1003 out of 4225, ratio=0.237
Memory used: 46 MB  Time elapsed: 130 s

 ***** [DMNRG] Iteration 22/67 (scale 0.00083772617549037929) ***** 
[1] xi(21)=1.00000014344 zeta(22)=0
Stats: nr=40 min=1 max=405 sum3=297823712
Safeguard: keep additional 4 states
Emax=10.0059396163 nrkept=927 nrkeptmult=2297 
Kept: 927 out of 3782, ratio=0.245
Memory used: 48 MB  Time elapsed: 130 s

 ***** [DMNRG] Iteration 23/67 (scale 0.00059236185946671895) ***** 
[1] xi(22)=1.00000007172 zeta(23)=0
Stats: nr=34 min=1 max=429 sum3=245882912
Safeguard: keep additional 2 states
Emax=10.0234973055 nrkept=883 nrkeptmult=2252 
Kept: 883 out of 3500, ratio=0.252
Memory used: 47 MB  Time elapsed: 131 s

 ***** [DMNRG] Iteration 24/67 (scale 0.00041886308774518965) ***** 
[1] xi(23)=1.00000003586 zeta(24)=0
Stats: nr=40 min=1 max=363 sum3=205107938
Safeguard: keep additional 2 states
Emax=10.0031019786 nrkept=829 nrkeptmult=2051 
Kept: 829 out of 3344, ratio=0.248
Memory used: 45 MB  Time elapsed: 131 s

 ***** [DMNRG] Iteration 25/67 (scale 0.00029618092973335947) ***** 
[1] xi(24)=1.00000001793 zeta(25)=0
Stats: nr=34 min=1 max=385 sum3=175545704
Safeguard: keep additional 2 states
Emax=10.0096480061 nrkept=827 nrkeptmult=2108 
Kept: 827 out of 3122, ratio=0.265
Memory used: 46 MB  Time elapsed: 132 s

 ***** [DMNRG] Iteration 26/67 (scale 0.00020943154387259482) ***** 
[1] xi(25)=1.00000000897 zeta(26)=0
Stats: nr=40 min=1 max=337 sum3=166805448
Safeguard: keep additional 1 states
Emax=10.0049168764 nrkept=786 nrkeptmult=1962 
Kept: 786 out of 3126, ratio=0.251
Memory used: 45 MB  Time elapsed: 132 s

 ***** [DMNRG] Iteration 27/67 (scale 0.00014809046486667974) ***** 
[1] xi(26)=1.00000000448 zeta(27)=0
Stats: nr=34 min=1 max=362 sum3=147227213
Safeguard: keep additional 2 states
Emax=10.0107454458 nrkept=805 nrkeptmult=2038 
Kept: 805 out of 2963, ratio=0.272
Memory used: 45 MB  Time elapsed: 133 s

 ***** [DMNRG] Iteration 28/67 (scale 0.00010471577193629741) ***** 
[1] xi(27)=1.00000000224 zeta(28)=0
Stats: nr=38 min=1 max=330 sum3=155870198
Safeguard: keep additional 2 states
Emax=10.0240535905 nrkept=763 nrkeptmult=1925 
Kept: 763 out of 3038, ratio=0.251
Memory used: 45 MB  Time elapsed: 133 s

 ***** [DMNRG] Iteration 29/67 (scale 7.4045232433339868e-05) ***** 
[1] xi(28)=1.00000000112 zeta(29)=0
Stats: nr=34 min=1 max=349 sum3=134658679
Safeguard: keep additional 3 states
Emax=10.0120592259 nrkept=776 nrkeptmult=1946 
Kept: 776 out of 2881, ratio=0.269
Memory used: 45 MB  Time elapsed: 134 s

 ***** [DMNRG] Iteration 30/67 (scale 5.2357885968148706e-05) ***** 
[1] xi(29)=1.00000000056 zeta(30)=0
Stats: nr=38 min=1 max=323 sum3=143658912
Safeguard: keep additional 1 states
Emax=10.0036920428 nrkept=723 nrkeptmult=1799 
Kept: 723 out of 2928, ratio=0.247
Memory used: 45 MB  Time elapsed: 134 s

 ***** [DMNRG] Iteration 31/67 (scale 3.7022616216669934e-05) ***** 
[1] xi(30)=1.00000000028 zeta(31)=0
Stats: nr=33 min=1 max=334 sum3=116636744
Safeguard: keep additional 1 states
Emax=10.0076247786 nrkept=728 nrkeptmult=1800 
Kept: 728 out of 2726, ratio=0.267
Memory used: 45 MB  Time elapsed: 134 s

 ***** [DMNRG] Iteration 32/67 (scale 2.6178942984074353e-05) ***** 
[1] xi(31)=1.00000000014 zeta(32)=0
Stats: nr=36 min=1 max=304 sum3=122567426
Safeguard: keep additional 1 states
Emax=10.0066276225 nrkept=658 nrkeptmult=1614 
Kept: 658 out of 2744, ratio=0.24
Memory used: 45 MB  Time elapsed: 135 s

 ***** [DMNRG] Iteration 33/67 (scale 1.8511308108334967e-05) ***** 
[1] xi(32)=1.00000000007 zeta(33)=0
Stats: nr=33 min=1 max=311 sum3=91943565
Safeguard: keep additional 1 states
Emax=10.0034905278 nrkept=690 nrkeptmult=1694 
Kept: 690 out of 2481, ratio=0.278
Memory used: 44 MB  Time elapsed: 135 s

 ***** [DMNRG] Iteration 34/67 (scale 1.3089471492037176e-05) ***** 
[1] xi(33)=1.00000000004 zeta(34)=0
Stats: nr=36 min=1 max=287 sum3=105360776
Safeguard: keep additional 2 states
Emax=10.0144525693 nrkept=607 nrkeptmult=1461 
Kept: 607 out of 2600, ratio=0.233
Memory used: 44 MB  Time elapsed: 135 s

 ***** [DMNRG] Iteration 35/67 (scale 9.2556540541674836e-06) ***** 
[1] xi(34)=1.00000000002 zeta(35)=0
Stats: nr=33 min=1 max=292 sum3=74039870
Safeguard: keep additional 3 states
Emax=10.0116782702 nrkept=619 nrkeptmult=1500 
Kept: 619 out of 2282, ratio=0.271
Memory used: 44 MB  Time elapsed: 136 s

 ***** [DMNRG] Iteration 36/67 (scale 6.5447357460185882e-06) ***** 
[1] xi(35)=1.00000000001 zeta(36)=0
Stats: nr=36 min=1 max=260 sum3=76488430
Safeguard: keep additional 2 states
Emax=10.0253498898 nrkept=552 nrkeptmult=1306 
Kept: 552 out of 2326, ratio=0.237
Memory used: 44 MB  Time elapsed: 136 s

 ***** [DMNRG] Iteration 37/67 (scale 4.6278270270837418e-06) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=33 min=1 max=269 sum3=57170652
Emax=10.0556843894 nrkept=570 nrkeptmult=1334 
Kept: 570 out of 2070, ratio=0.275
Memory used: 43 MB  Time elapsed: 136 s

 ***** [DMNRG] Iteration 38/67 (scale 3.2723678730092941e-06) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=34 min=1 max=245 sum3=62091594
Emax=10.0036967653 nrkept=512 nrkeptmult=1190 
Kept: 512 out of 2136, ratio=0.24
Memory used: 43 MB  Time elapsed: 136 s

 ***** [DMNRG] Iteration 39/67 (scale 2.3139135135418709e-06) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=33 min=1 max=251 sum3=45594660
Emax=10.3327931792 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1914, ratio=0.296
Memory used: 43 MB  Time elapsed: 137 s

 ***** [DMNRG] Iteration 40/67 (scale 1.6361839365046471e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.1296480344 nrkept=499 nrkeptmult=1149 
Kept: 499 out of 2124, ratio=0.235
Memory used: 43 MB  Time elapsed: 137 s

 ***** [DMNRG] Iteration 41/67 (scale 1.1569567567709354e-06) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.6374354844 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1862, ratio=0.305
Memory used: 42 MB  Time elapsed: 137 s

 ***** [DMNRG] Iteration 42/67 (scale 8.1809196825232353e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.3597623579 nrkept=499 nrkeptmult=1149 
Kept: 499 out of 2124, ratio=0.235
Memory used: 43 MB  Time elapsed: 137 s

 ***** [DMNRG] Iteration 43/67 (scale 5.7847837838546772e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.7372964209 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1862, ratio=0.305
Memory used: 42 MB  Time elapsed: 137 s

 ***** [DMNRG] Iteration 44/67 (scale 4.0904598412616176e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.4899287193 nrkept=499 nrkeptmult=1149 
Kept: 499 out of 2124, ratio=0.235
Memory used: 43 MB  Time elapsed: 138 s

 ***** [DMNRG] Iteration 45/67 (scale 2.8923918919273386e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.789385392 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1862, ratio=0.305
Memory used: 42 MB  Time elapsed: 138 s

 ***** [DMNRG] Iteration 46/67 (scale 2.0452299206308088e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5111582642 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 43 MB  Time elapsed: 138 s

 ***** [DMNRG] Iteration 47/67 (scale 1.4461959459636693e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8156857111 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 138 s

 ***** [DMNRG] Iteration 48/67 (scale 1.0226149603154044e-07) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5185900319 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 43 MB  Time elapsed: 138 s

 ***** [DMNRG] Iteration 49/67 (scale 7.2309797298183465e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8286333834 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 139 s

 ***** [DMNRG] Iteration 50/67 (scale 5.113074801577022e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5223242957 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 139 s

 ***** [DMNRG] Iteration 51/67 (scale 3.6154898649091733e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8333633709 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 139 s

 ***** [DMNRG] Iteration 52/67 (scale 2.556537400788511e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5242078641 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 139 s

 ***** [DMNRG] Iteration 53/67 (scale 1.8077449324545866e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8355062475 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 140 s

 ***** [DMNRG] Iteration 54/67 (scale 1.2782687003942555e-08) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5251454635 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 140 s

 ***** [DMNRG] Iteration 55/67 (scale 9.0387246622729332e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8365765199 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 140 s

 ***** [DMNRG] Iteration 56/67 (scale 6.3913435019712775e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.52561457 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 140 s

 ***** [DMNRG] Iteration 57/67 (scale 4.5193623311364666e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8371113572 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 140 s

 ***** [DMNRG] Iteration 58/67 (scale 3.1956717509856388e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5258492 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 45 MB  Time elapsed: 141 s

 ***** [DMNRG] Iteration 59/67 (scale 2.2596811655682333e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8373787018 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 44 MB  Time elapsed: 141 s

 ***** [DMNRG] Iteration 60/67 (scale 1.5978358754928194e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5259665342 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 141 s

 ***** [DMNRG] Iteration 61/67 (scale 1.1298405827841166e-09) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 2 states
Emax=10.8375174929 nrkept=569 nrkeptmult=1330 
Kept: 569 out of 1857, ratio=0.306
Memory used: 44 MB  Time elapsed: 141 s

 ***** [DMNRG] Iteration 62/67 (scale 7.9891793774640969e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=34 min=1 max=244 sum3=61697330
Emax=10.5260252061 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2132, ratio=0.234
Memory used: 45 MB  Time elapsed: 142 s

 ***** [DMNRG] Iteration 63/67 (scale 5.6492029139205832e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8375968195 nrkept=575 nrkeptmult=1340 
Kept: 575 out of 1857, ratio=0.31
Memory used: 44 MB  Time elapsed: 142 s

 ***** [DMNRG] Iteration 64/67 (scale 3.9945896887320485e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260401124 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2154, ratio=0.231
Memory used: 45 MB  Time elapsed: 142 s

 ***** [DMNRG] Iteration 65/67 (scale 2.8246014569602916e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8376155039 nrkept=575 nrkeptmult=1340 
Kept: 575 out of 1857, ratio=0.31
Memory used: 44 MB  Time elapsed: 142 s

 ***** [DMNRG] Iteration 66/67 (scale 1.9972948443660242e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260547818 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2154, ratio=0.231
Memory used: 45 MB  Time elapsed: 142 s

 ***** [DMNRG] Iteration 67/67 (scale 1.4123007284801458e-10) ***** 
[1] xi(66)=1 zeta(67)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 12 states
Emax=10.8376388455 nrkept=1857 nrkeptmult=4568 
Kept: 1857 out of 1857, ratio=1
Memory used: 48 MB  Time elapsed: 143 s

** Iteration completed.

Spectrum: A_d-A_d FDM
[A_d-A_d] pos=0.500000025815 neg=0.499999974183 sum= 0.999999999998
mu1=3.35278679819e-15 mu2=0.000667030487483 mu3=2.40421831749e-17 mu4=0.000216439853132


Memory usage report
===================
  after diag: 101708 kB
after recalc: 116844 kB
 after split: 129772 kB
  after trim: 58200 kB
 after trunc: 101892 kB
      diag_h: 65192 kB

Peak usage: 129772 kB

Timing report
=============
         All: 143 s
          DM: 1 s
         FDM: 7 s
        diag: 79 s
    recalc d: 9 s
      spec d: 24 s
       Other: 21 s
