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=66
TEMPERATURE T=1e-10
Nmax=66
Tmin=1e-10
Ninit=0
Lambda=2
strategy=kept
lastall=1
lastalloverride=0
z=0.75
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.75
>>  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.75
# keep=10000 keepenergy=10 Nmax=66 channels=1 dots=1 betabar=1

Number of invariant subspaces: 6
(-2 1) 1 states
0.0437078924027 
(-1 2) 2 states
0.0207472412478 0.0625470657632 
(0 1) 3 states
0 0.0437078924027 0.083294307011 
(0 3) 1 states
0.0395864146084 
(1 2) 2 states
0.0207472412478 0.0625470657632 
(2 1) 1 states
0.0437078924027 
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/66 (scale 0.85783160370214839) ***** 
[1] xi(0)=0.684778584642 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=2.7474472697 nrkept=35 nrkeptmult=64 
Total energy=-1.17592354873  Egs=-1.31482566553
Kept: 35 out of 35, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/66 (scale 0.6065785440939202) ***** 
[1] xi(1)=0.833755393818 zeta(2)=0
Stats: nr=15 min=1 max=20 sum3=30066
Emax=5.22193421287 nrkept=126 nrkeptmult=256 
Total energy=-1.5812566261  Egs=-0.668228511079
Kept: 126 out of 126, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/66 (scale 0.4289158018510742) ***** 
[1] xi(2)=0.947046165312 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=10.4050026003 nrkept=460 nrkeptmult=1022 
Total energy=-2.23141165172  Egs=-1.51581038241
Kept: 460 out of 462, ratio=0.996
Memory used: 5 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/66 (scale 0.3032892720469601) ***** 
[1] xi(3)=0.972572268687 zeta(4)=0
Stats: nr=28 min=1 max=209 sum3=36715794
Emax=10.730991904 nrkept=1259 nrkeptmult=3043 
Total energy=-2.48205848013  Egs=-0.826428270005
Kept: 1259 out of 1710, ratio=0.736
Memory used: 12 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/66 (scale 0.2144579009255371) ***** 
[1] xi(4)=0.980132362648 zeta(5)=0
Stats: nr=36 min=1 max=561 sum3=560127451
Safeguard: keep additional 7 states
Emax=10.429654582 nrkept=1730 nrkeptmult=4298 
Total energy=-2.80343808263  Egs=-1.49856732306
Kept: 1730 out of 4729, ratio=0.366
Memory used: 20 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 6/66 (scale 0.15164463602348005) ***** 
[1] xi(5)=0.988164304343 zeta(6)=0
Stats: nr=42 min=1 max=665 sum3=1351171682
Safeguard: keep additional 1 states
Emax=10.7114800551 nrkept=2122 nrkeptmult=5284 
Total energy=-2.93591339631  Egs=-0.873590501791
Kept: 2122 out of 6518, ratio=0.326
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 7/66 (scale 0.10722895046276855) ***** 
[1] xi(6)=0.993732382824 zeta(7)=0
Stats: nr=44 min=1 max=910 sum3=2517194209
Safeguard: keep additional 2 states
Emax=10.0274508866 nrkept=1817 nrkeptmult=4474 
Total energy=-3.09465533619  Egs=-1.48040188028
Kept: 1817 out of 8005, ratio=0.227
Memory used: 24 MB  Time elapsed: 11 s

 ***** [NRG] Iteration 8/66 (scale 0.075822318011740025) ***** 
[1] xi(7)=0.996801670793 zeta(8)=0
Stats: nr=42 min=1 max=705 sum3=1607566036
Safeguard: keep additional 1 states
Emax=10.0204790702 nrkept=1970 nrkeptmult=4888 
Total energy=-3.16363661765  Egs=-0.909775423186
Kept: 1970 out of 6832, ratio=0.288
Memory used: 26 MB  Time elapsed: 16 s

 ***** [NRG] Iteration 9/66 (scale 0.053614475231384275) ***** 
[1] xi(8)=0.998386814476 zeta(9)=0
Stats: nr=44 min=1 max=850 sum3=2036505503
Safeguard: keep additional 3 states
Emax=10.0056025545 nrkept=1654 nrkeptmult=4108 
Total energy=-3.24140446598  Egs=-1.45050096996
Kept: 1654 out of 7433, ratio=0.223
Memory used: 25 MB  Time elapsed: 21 s

 ***** [NRG] Iteration 10/66 (scale 0.037911159005870013) ***** 
[1] xi(9)=0.999190082691 zeta(10)=0
Stats: nr=42 min=1 max=646 sum3=1226297206
Safeguard: keep additional 2 states
Emax=10.0140218091 nrkept=1857 nrkeptmult=4609 
Total energy=-3.27708665009  Egs=-0.941205308618
Kept: 1857 out of 6244, ratio=0.297
Memory used: 27 MB  Time elapsed: 23 s

 ***** [NRG] Iteration 11/66 (scale 0.026807237615692137) ***** 
[1] xi(10)=0.999594228576 zeta(11)=0
Stats: nr=42 min=1 max=805 sum3=1725741631
Safeguard: keep additional 1 states
Emax=10.0000492513 nrkept=1568 nrkeptmult=3862 
Total energy=-3.31501685504  Egs=-1.41492404005
Kept: 1568 out of 6997, ratio=0.224
Memory used: 27 MB  Time elapsed: 28 s

 ***** [NRG] Iteration 12/66 (scale 0.018955579502935006) ***** 
[1] xi(11)=0.999796913201 zeta(12)=0
Stats: nr=42 min=1 max=616 sum3=1053105550
Safeguard: keep additional 3 states
Emax=10.0098411107 nrkept=1668 nrkeptmult=4104 
Total energy=-3.33343022239  Egs=-0.971395643689
Kept: 1668 out of 5908, ratio=0.282
Memory used: 27 MB  Time elapsed: 31 s

 ***** [NRG] Iteration 13/66 (scale 0.013403618807846069) ***** 
[1] xi(12)=0.999898406582 zeta(13)=0
Stats: nr=42 min=1 max=732 sum3=1273918029
Emax=10.0091423014 nrkept=1547 nrkeptmult=3808 
Total energy=-3.35190041845  Egs=-1.37800069727
Kept: 1547 out of 6273, ratio=0.247
Memory used: 27 MB  Time elapsed: 35 s

 ***** [NRG] Iteration 14/66 (scale 0.0094777897514675032) ***** 
[1] xi(13)=0.999949190818 zeta(14)=0
Stats: nr=42 min=1 max=605 sum3=1003484280
Safeguard: keep additional 2 states
Emax=10.0045892942 nrkept=1547 nrkeptmult=3775 
Total energy=-3.36137895239  Egs=-1.00007851863
Kept: 1547 out of 5826, ratio=0.266
Memory used: 28 MB  Time elapsed: 38 s

 ***** [NRG] Iteration 15/66 (scale 0.0067018094039230343) ***** 
[1] xi(14)=0.999974592294 zeta(15)=0
Stats: nr=36 min=1 max=690 sum3=1048131907
Safeguard: keep additional 2 states
Emax=10.0361664552 nrkept=1473 nrkeptmult=3652 
Total energy=-3.370369783  Egs=-1.34155271613
Kept: 1473 out of 5821, ratio=0.253
Memory used: 27 MB  Time elapsed: 41 s

 ***** [NRG] Iteration 16/66 (scale 0.0047388948757337516) ***** 
[1] xi(15)=0.999987295369 zeta(16)=0
Stats: nr=40 min=1 max=578 sum3=885160508
Emax=10.0052154861 nrkept=1435 nrkeptmult=3467 
Total energy=-3.37524212767  Egs=-1.02816053096
Kept: 1435 out of 5546, ratio=0.259
Memory used: 28 MB  Time elapsed: 44 s

 ***** [NRG] Iteration 17/66 (scale 0.0033509047019615172) ***** 
[1] xi(16)=0.99999364749 zeta(17)=0
Stats: nr=36 min=1 max=647 sum3=851865453
Safeguard: keep additional 2 states
Emax=10.003231859 nrkept=1259 nrkeptmult=3110 
Total energy=-3.37961350333  Egs=-1.30453595506
Kept: 1259 out of 5397, ratio=0.233
Memory used: 28 MB  Time elapsed: 46 s

 ***** [NRG] Iteration 18/66 (scale 0.0023694474378668758) ***** 
[1] xi(17)=0.999996823696 zeta(18)=0
Stats: nr=40 min=1 max=501 sum3=569892192
Safeguard: keep additional 1 states
Emax=10.0007365047 nrkept=1234 nrkeptmult=3032 
Total energy=-3.38212244717  Egs=-1.05887296723
Kept: 1234 out of 4740, ratio=0.26
Memory used: 30 MB  Time elapsed: 48 s

 ***** [NRG] Iteration 19/66 (scale 0.0016754523509807586) ***** 
[1] xi(18)=0.999998411836 zeta(19)=0
Stats: nr=34 min=1 max=556 sum3=546872911
Safeguard: keep additional 1 states
Emax=10.0023682967 nrkept=1107 nrkeptmult=2758 
Total energy=-3.38424009245  Egs=-1.26392449908
Kept: 1107 out of 4645, ratio=0.238
Memory used: 29 MB  Time elapsed: 50 s

 ***** [NRG] Iteration 20/66 (scale 0.0011847237189334379) ***** 
[1] xi(19)=0.999999205915 zeta(20)=0
Stats: nr=40 min=1 max=444 sum3=396066248
Emax=10.0053941462 nrkept=1060 nrkeptmult=2616 
Total energy=-3.38553864112  Egs=-1.09607721341
Kept: 1060 out of 4178, ratio=0.254
Memory used: 30 MB  Time elapsed: 51 s

 ***** [NRG] Iteration 21/66 (scale 0.00083772617549037929) ***** 
[1] xi(20)=0.999999602957 zeta(21)=0
Stats: nr=34 min=1 max=484 sum3=359180251
Safeguard: keep additional 2 states
Emax=10.0020404209 nrkept=924 nrkeptmult=2332 
Total energy=-3.38655779869  Egs=-1.21657601508
Kept: 924 out of 3997, ratio=0.231
Memory used: 31 MB  Time elapsed: 53 s

 ***** [NRG] Iteration 22/66 (scale 0.00059236185946671895) ***** 
[1] xi(21)=0.999999801478 zeta(22)=0
Stats: nr=40 min=1 max=376 sum3=234913884
Safeguard: keep additional 1 states
Emax=10.0155455125 nrkept=866 nrkeptmult=2148 
Total energy=-3.38723476434  Egs=-1.14282450737
Kept: 866 out of 3492, ratio=0.248
Memory used: 31 MB  Time elapsed: 53 s

 ***** [NRG] Iteration 23/66 (scale 0.00041886308774518965) ***** 
[1] xi(22)=0.999999900739 zeta(23)=0
Stats: nr=34 min=1 max=402 sum3=201263845
Safeguard: keep additional 1 states
Emax=10.0026731864 nrkept=846 nrkeptmult=2158 
Total energy=-3.38772084336  Egs=-1.16047230505
Kept: 846 out of 3265, ratio=0.259
Memory used: 32 MB  Time elapsed: 54 s

 ***** [NRG] Iteration 24/66 (scale 0.00029618092973335947) ***** 
[1] xi(23)=0.99999995037 zeta(24)=0
Stats: nr=40 min=1 max=346 sum3=178056912
Emax=10.0177243709 nrkept=810 nrkeptmult=2008 
Total energy=-3.38807659184  Egs=-1.20111879018
Kept: 810 out of 3198, ratio=0.253
Memory used: 31 MB  Time elapsed: 55 s

 ***** [NRG] Iteration 25/66 (scale 0.00020943154387259482) ***** 
[1] xi(24)=0.999999975185 zeta(25)=0
Stats: nr=34 min=1 max=374 sum3=161496433
Emax=10.0276516003 nrkept=816 nrkeptmult=2066 
Total energy=-3.38830602758  Egs=-1.09551664446
Kept: 816 out of 3049, ratio=0.268
Memory used: 31 MB  Time elapsed: 55 s

 ***** [NRG] Iteration 26/66 (scale 0.00014809046486667974) ***** 
[1] xi(25)=0.999999987592 zeta(26)=0
Stats: nr=38 min=1 max=333 sum3=161834170
Emax=10.0100608586 nrkept=770 nrkeptmult=1936 
Total energy=-3.38849429027  Egs=-1.27126814307
Kept: 770 out of 3082, ratio=0.25
Memory used: 31 MB  Time elapsed: 56 s

 ***** [NRG] Iteration 27/66 (scale 0.00010471577193629741) ***** 
[1] xi(26)=0.999999993796 zeta(27)=0
Stats: nr=34 min=1 max=354 sum3=139173931
Safeguard: keep additional 2 states
Emax=10.000601265 nrkept=789 nrkeptmult=1982 
Total energy=-3.38860161043  Egs=-1.02487102765
Kept: 789 out of 2905, ratio=0.272
Memory used: 31 MB  Time elapsed: 57 s

 ***** [NRG] Iteration 28/66 (scale 7.4045232433339868e-05) ***** 
[1] xi(27)=0.999999996898 zeta(28)=0
Stats: nr=38 min=1 max=326 sum3=149194284
Emax=10.0159752848 nrkept=751 nrkeptmult=1897 
Total energy=-3.38870156286  Egs=-1.34988342693
Kept: 751 out of 2976, ratio=0.252
Memory used: 31 MB  Time elapsed: 57 s

 ***** [NRG] Iteration 29/66 (scale 5.2357885968148706e-05) ***** 
[1] xi(28)=0.999999998449 zeta(29)=0
Stats: nr=34 min=1 max=343 sum3=127606636
Safeguard: keep additional 2 states
Emax=10.0031835774 nrkept=753 nrkeptmult=1888 
Total energy=-3.38875161767  Egs=-0.956012727559
Kept: 753 out of 2836, ratio=0.266
Memory used: 31 MB  Time elapsed: 58 s

 ***** [NRG] Iteration 30/66 (scale 3.7022616216669934e-05) ***** 
[1] xi(29)=0.999999999225 zeta(30)=0
Stats: nr=38 min=1 max=314 sum3=132096960
Emax=10.0179729349 nrkept=693 nrkeptmult=1713 
Total energy=-3.38880448761  Egs=-1.4280445008
Kept: 693 out of 2844, ratio=0.244
Memory used: 32 MB  Time elapsed: 58 s

 ***** [NRG] Iteration 31/66 (scale 2.6178942984074353e-05) ***** 
[1] xi(30)=0.999999999612 zeta(31)=0
Stats: nr=33 min=1 max=325 sum3=106369948
Safeguard: keep additional 3 states
Emax=10.0067640763 nrkept=709 nrkeptmult=1744 
Total energy=-3.38882800589  Egs=-0.89836623853
Kept: 709 out of 2614, ratio=0.271
Memory used: 31 MB  Time elapsed: 59 s

 ***** [NRG] Iteration 32/66 (scale 1.8511308108334967e-05) ***** 
[1] xi(31)=0.999999999806 zeta(32)=0
Stats: nr=36 min=1 max=296 sum3=113592734
Emax=10.0011727924 nrkept=633 nrkeptmult=1541 
Total energy=-3.38885566173  Egs=-1.49399699642
Kept: 633 out of 2672, ratio=0.237
Memory used: 32 MB  Time elapsed: 59 s

 ***** [NRG] Iteration 33/66 (scale 1.3089471492037176e-05) ***** 
[1] xi(32)=0.999999999903 zeta(33)=0
Stats: nr=33 min=1 max=303 sum3=83225817
Safeguard: keep additional 1 states
Emax=10.018903115 nrkept=656 nrkeptmult=1588 
Total energy=-3.38886688455  Egs=-0.857392968351
Kept: 656 out of 2385, ratio=0.275
Memory used: 32 MB  Time elapsed: 60 s

 ***** [NRG] Iteration 34/66 (scale 9.2556540541674836e-06) ***** 
[1] xi(33)=0.999999999952 zeta(34)=0
Stats: nr=36 min=1 max=276 sum3=91791798
Safeguard: keep additional 2 states
Emax=10.052883151 nrkept=581 nrkeptmult=1383 
Total energy=-3.38888114523  Egs=-1.54075405073
Kept: 581 out of 2466, ratio=0.236
Memory used: 31 MB  Time elapsed: 60 s

 ***** [NRG] Iteration 35/66 (scale 6.5447357460185882e-06) ***** 
[1] xi(34)=0.999999999976 zeta(35)=0
Stats: nr=33 min=1 max=283 sum3=65536626
Safeguard: keep additional 1 states
Emax=10.0368229375 nrkept=589 nrkeptmult=1402 
Total energy=-3.38888659041  Egs=-0.831992849493
Kept: 589 out of 2178, ratio=0.27
Memory used: 32 MB  Time elapsed: 60 s

 ***** [NRG] Iteration 36/66 (scale 4.6278270270837418e-06) ***** 
[1] xi(35)=0.999999999988 zeta(36)=0
Stats: nr=36 min=1 max=251 sum3=67223648
Emax=10.0001206832 nrkept=523 nrkeptmult=1227 
Total energy=-3.38889385381  Egs=-1.56950685185
Kept: 523 out of 2210, ratio=0.237
Memory used: 32 MB  Time elapsed: 61 s

 ***** [NRG] Iteration 37/66 (scale 3.2723678730092941e-06) ***** 
[1] xi(36)=0.999999999994 zeta(37)=0
Stats: nr=33 min=1 max=257 sum3=48651854
Emax=10.101122394 nrkept=567 nrkeptmult=1322 
Total energy=-3.38889652947  Egs=-0.817651877635
Kept: 567 out of 1958, ratio=0.29
Memory used: 32 MB  Time elapsed: 61 s

 ***** [NRG] Iteration 38/66 (scale 2.3139135135418709e-06) ***** 
[1] xi(37)=0.999999999997 zeta(38)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.0464946451 nrkept=501 nrkeptmult=1157 
Total energy=-3.38890019842  Egs=-1.58560372751
Kept: 501 out of 2124, ratio=0.236
Memory used: 32 MB  Time elapsed: 61 s

 ***** [NRG] Iteration 39/66 (scale 1.6361839365046471e-06) ***** 
[1] xi(38)=0.999999999998 zeta(39)=0
Stats: nr=33 min=1 max=245 sum3=42665008
Emax=10.5590939931 nrkept=567 nrkeptmult=1322 
Total energy=-3.38890152373  Egs=-0.809997043932
Kept: 567 out of 1870, ratio=0.303
Memory used: 32 MB  Time elapsed: 61 s

 ***** [NRG] Iteration 40/66 (scale 1.1569567567709354e-06) ***** 
[1] xi(39)=0.999999999999 zeta(40)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.2591669732 nrkept=499 nrkeptmult=1149 
Total energy=-3.38890336808  Egs=-1.59414157977
Kept: 499 out of 2124, ratio=0.235
Memory used: 32 MB  Time elapsed: 62 s

 ***** [NRG] Iteration 41/66 (scale 8.1809196825232353e-07) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.695187534 nrkept=567 nrkeptmult=1322 
Total energy=-3.38890402749  Egs=-0.806036905086
Kept: 567 out of 1862, ratio=0.305
Memory used: 32 MB  Time elapsed: 62 s

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

 ***** [NRG] Iteration 43/66 (scale 4.0904598412616176e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.7676439049 nrkept=567 nrkeptmult=1322 
Total energy=-3.38890528109  Egs=-0.804021988569
Kept: 567 out of 1862, ratio=0.305
Memory used: 33 MB  Time elapsed: 62 s

 ***** [NRG] Iteration 44/66 (scale 2.8923918919273386e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5050375077 nrkept=498 nrkeptmult=1142 
Total energy=-3.3889057441  Egs=-1.60077471498
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 63 s

 ***** [NRG] Iteration 45/66 (scale 2.0452299206308088e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8047350699 nrkept=567 nrkeptmult=1322 
Total energy=-3.38890590833  Egs=-0.803005593677
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 63 s

 ***** [NRG] Iteration 46/66 (scale 1.4461959459636693e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5155044897 nrkept=498 nrkeptmult=1142 
Total energy=-3.38890614  Egs=-1.60190007302
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 63 s

 ***** [NRG] Iteration 47/66 (scale 1.0226149603154044e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8233626321 nrkept=567 nrkeptmult=1322 
Total energy=-3.38890622206  Egs=-0.802495132785
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 63 s

 ***** [NRG] Iteration 48/66 (scale 7.2309797298183465e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5207756558 nrkept=498 nrkeptmult=1142 
Total energy=-3.38890633794  Egs=-1.60246491533
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 64 s

 ***** [NRG] Iteration 49/66 (scale 5.113074801577022e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8315849979 nrkept=567 nrkeptmult=1324 
Total energy=-3.38890637896  Egs=-0.802239332744
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 64 s

 ***** [NRG] Iteration 50/66 (scale 3.6154898649091733e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5234315011 nrkept=498 nrkeptmult=1142 
Total energy=-3.3889064369  Egs=-1.60274787955
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 64 s

 ***** [NRG] Iteration 51/66 (scale 2.556537400788511e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8346187093 nrkept=567 nrkeptmult=1324 
Total energy=-3.38890645741  Egs=-0.802111289854
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 64 s

 ***** [NRG] Iteration 52/66 (scale 1.8077449324545866e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5247568732 nrkept=498 nrkeptmult=1142 
Total energy=-3.38890648639  Egs=-1.60288949771
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 64 s

 ***** [NRG] Iteration 53/66 (scale 1.2782687003942555e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8361331529 nrkept=567 nrkeptmult=1324 
Total energy=-3.38890649664  Egs=-0.802047232633
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 65 s

 ***** [NRG] Iteration 54/66 (scale 9.0387246622729332e-09) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5254201724 nrkept=498 nrkeptmult=1142 
Total energy=-3.38890651113  Egs=-1.60296034084
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 65 s

 ***** [NRG] Iteration 55/66 (scale 6.3913435019712775e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.836889773 nrkept=567 nrkeptmult=1324 
Total energy=-3.38890651625  Egs=-0.802015195072
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 65 s

 ***** [NRG] Iteration 56/66 (scale 4.5193623311364666e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5257519756 nrkept=498 nrkeptmult=1142 
Total energy=-3.3889065235  Egs=-1.60299577091
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 65 s

 ***** [NRG] Iteration 57/66 (scale 3.1956717509856388e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8372679344 nrkept=567 nrkeptmult=1324 
Total energy=-3.38890652606  Egs=-0.801999174053
Kept: 567 out of 1857, ratio=0.305
Memory used: 34 MB  Time elapsed: 66 s

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

 ***** [NRG] Iteration 59/66 (scale 1.5978358754928194e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8374569781 nrkept=567 nrkeptmult=1324 
Total energy=-3.38890653097  Egs=-0.801991162984
Kept: 567 out of 1857, ratio=0.305
Memory used: 34 MB  Time elapsed: 66 s

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

 ***** [NRG] Iteration 61/66 (scale 7.9891793774640969e-10) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 2 states
Emax=10.837551717 nrkept=569 nrkeptmult=1330 
Total energy=-3.38890653342  Egs=-0.80198715731
Kept: 569 out of 1857, ratio=0.306
Memory used: 34 MB  Time elapsed: 67 s

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

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

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

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

 ***** [NRG] Iteration 66/66 (scale 1.4123007284801458e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260590775 nrkept=2154 nrkeptmult=5360 
Total energy=-3.38890653548  Egs=-1.60303009922
Kept: 2154 out of 2154, ratio=1
Memory used: 40 MB  Time elapsed: 68 s

** Iteration completed.


Total energy:             -3.388906535482716
TKW=2.16102831736e-05
Grand canonical Z=              4.891300339574657
Zft (last shell)=               4.891310464705425
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.75
>>  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.0437078924027 
(-1 2) 2 states
0.0207472412478 0.0625470657632 
(0 1) 3 states
0 0.0437078924027 0.083294307011 
(0 3) 1 states
0.0395864146084 
(1 2) 2 states
0.0207472412478 0.0625470657632 
(2 1) 1 states
0.0437078924027 
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/66 (scale 0.85783160370214839) ***** 
[1] xi(0)=0.684778584642 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=2.7474472697 nrkept=35 nrkeptmult=64 
Kept: 35 out of 35, ratio=1
Memory used: 36 MB  Time elapsed: 75 s

 ***** [DMNRG] Iteration 2/66 (scale 0.6065785440939202) ***** 
[1] xi(1)=0.833755393818 zeta(2)=0
Stats: nr=15 min=1 max=20 sum3=30066
Emax=5.22193421287 nrkept=126 nrkeptmult=256 
Kept: 126 out of 126, ratio=1
Memory used: 36 MB  Time elapsed: 75 s

 ***** [DMNRG] Iteration 3/66 (scale 0.4289158018510742) ***** 
[1] xi(2)=0.947046165312 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=10.4050026003 nrkept=460 nrkeptmult=1022 
Kept: 460 out of 462, ratio=0.996
Memory used: 36 MB  Time elapsed: 75 s

 ***** [DMNRG] Iteration 4/66 (scale 0.3032892720469601) ***** 
[1] xi(3)=0.972572268687 zeta(4)=0
Stats: nr=28 min=1 max=209 sum3=36715794
Emax=10.730991904 nrkept=1259 nrkeptmult=3043 
Kept: 1259 out of 1710, ratio=0.736
Memory used: 40 MB  Time elapsed: 75 s

 ***** [DMNRG] Iteration 5/66 (scale 0.2144579009255371) ***** 
[1] xi(4)=0.980132362648 zeta(5)=0
Stats: nr=36 min=1 max=561 sum3=560127451
Safeguard: keep additional 7 states
Emax=10.429654582 nrkept=1730 nrkeptmult=4298 
Kept: 1730 out of 4729, ratio=0.366
Memory used: 45 MB  Time elapsed: 76 s

 ***** [DMNRG] Iteration 6/66 (scale 0.15164463602348005) ***** 
[1] xi(5)=0.988164304343 zeta(6)=0
Stats: nr=42 min=1 max=665 sum3=1351171682
Safeguard: keep additional 1 states
Emax=10.7114800551 nrkept=2122 nrkeptmult=5284 
Kept: 2122 out of 6518, ratio=0.326
Memory used: 50 MB  Time elapsed: 77 s

 ***** [DMNRG] Iteration 7/66 (scale 0.10722895046276855) ***** 
[1] xi(6)=0.993732382824 zeta(7)=0
Stats: nr=44 min=1 max=910 sum3=2517194209
Safeguard: keep additional 2 states
Emax=10.0274508866 nrkept=1817 nrkeptmult=4474 
Kept: 1817 out of 8005, ratio=0.227
Memory used: 48 MB  Time elapsed: 80 s

 ***** [DMNRG] Iteration 8/66 (scale 0.075822318011740025) ***** 
[1] xi(7)=0.996801670793 zeta(8)=0
Stats: nr=42 min=1 max=705 sum3=1607566036
Safeguard: keep additional 1 states
Emax=10.0204790702 nrkept=1970 nrkeptmult=4888 
Kept: 1970 out of 6832, ratio=0.288
Memory used: 49 MB  Time elapsed: 82 s

 ***** [DMNRG] Iteration 9/66 (scale 0.053614475231384275) ***** 
[1] xi(8)=0.998386814476 zeta(9)=0
Stats: nr=44 min=1 max=850 sum3=2036505503
Safeguard: keep additional 3 states
Emax=10.0056025545 nrkept=1654 nrkeptmult=4108 
Kept: 1654 out of 7433, ratio=0.223
Memory used: 48 MB  Time elapsed: 84 s

 ***** [DMNRG] Iteration 10/66 (scale 0.037911159005870013) ***** 
[1] xi(9)=0.999190082691 zeta(10)=0
Stats: nr=42 min=1 max=646 sum3=1226297206
Safeguard: keep additional 2 states
Emax=10.0140218091 nrkept=1857 nrkeptmult=4609 
Kept: 1857 out of 6244, ratio=0.297
Memory used: 49 MB  Time elapsed: 85 s

 ***** [DMNRG] Iteration 11/66 (scale 0.026807237615692137) ***** 
[1] xi(10)=0.999594228576 zeta(11)=0
Stats: nr=42 min=1 max=805 sum3=1725741631
Safeguard: keep additional 1 states
Emax=10.0000492513 nrkept=1568 nrkeptmult=3862 
Kept: 1568 out of 6997, ratio=0.224
Memory used: 46 MB  Time elapsed: 87 s

 ***** [DMNRG] Iteration 12/66 (scale 0.018955579502935006) ***** 
[1] xi(11)=0.999796913201 zeta(12)=0
Stats: nr=42 min=1 max=616 sum3=1053105550
Safeguard: keep additional 3 states
Emax=10.0098411107 nrkept=1668 nrkeptmult=4104 
Kept: 1668 out of 5908, ratio=0.282
Memory used: 46 MB  Time elapsed: 89 s

 ***** [DMNRG] Iteration 13/66 (scale 0.013403618807846069) ***** 
[1] xi(12)=0.999898406582 zeta(13)=0
Stats: nr=42 min=1 max=732 sum3=1273918029
Emax=10.0091423014 nrkept=1547 nrkeptmult=3808 
Kept: 1547 out of 6273, ratio=0.247
Memory used: 46 MB  Time elapsed: 90 s

 ***** [DMNRG] Iteration 14/66 (scale 0.0094777897514675032) ***** 
[1] xi(13)=0.999949190818 zeta(14)=0
Stats: nr=42 min=1 max=605 sum3=1003484280
Safeguard: keep additional 2 states
Emax=10.0045892942 nrkept=1547 nrkeptmult=3775 
Kept: 1547 out of 5826, ratio=0.266
Memory used: 46 MB  Time elapsed: 92 s

 ***** [DMNRG] Iteration 15/66 (scale 0.0067018094039230343) ***** 
[1] xi(14)=0.999974592294 zeta(15)=0
Stats: nr=36 min=1 max=690 sum3=1048131907
Safeguard: keep additional 2 states
Emax=10.0361664552 nrkept=1473 nrkeptmult=3652 
Kept: 1473 out of 5821, ratio=0.253
Memory used: 44 MB  Time elapsed: 93 s

 ***** [DMNRG] Iteration 16/66 (scale 0.0047388948757337516) ***** 
[1] xi(15)=0.999987295369 zeta(16)=0
Stats: nr=40 min=1 max=578 sum3=885160508
Emax=10.0052154861 nrkept=1435 nrkeptmult=3467 
Kept: 1435 out of 5546, ratio=0.259
Memory used: 45 MB  Time elapsed: 94 s

 ***** [DMNRG] Iteration 17/66 (scale 0.0033509047019615172) ***** 
[1] xi(16)=0.99999364749 zeta(17)=0
Stats: nr=36 min=1 max=647 sum3=851865453
Safeguard: keep additional 2 states
Emax=10.003231859 nrkept=1259 nrkeptmult=3110 
Kept: 1259 out of 5397, ratio=0.233
Memory used: 43 MB  Time elapsed: 95 s

 ***** [DMNRG] Iteration 18/66 (scale 0.0023694474378668758) ***** 
[1] xi(17)=0.999996823696 zeta(18)=0
Stats: nr=40 min=1 max=501 sum3=569892192
Safeguard: keep additional 1 states
Emax=10.0007365047 nrkept=1234 nrkeptmult=3032 
Kept: 1234 out of 4740, ratio=0.26
Memory used: 45 MB  Time elapsed: 96 s

 ***** [DMNRG] Iteration 19/66 (scale 0.0016754523509807586) ***** 
[1] xi(18)=0.999998411836 zeta(19)=0
Stats: nr=34 min=1 max=556 sum3=546872911
Safeguard: keep additional 1 states
Emax=10.0023682967 nrkept=1107 nrkeptmult=2758 
Kept: 1107 out of 4645, ratio=0.238
Memory used: 47 MB  Time elapsed: 97 s

 ***** [DMNRG] Iteration 20/66 (scale 0.0011847237189334379) ***** 
[1] xi(19)=0.999999205915 zeta(20)=0
Stats: nr=40 min=1 max=444 sum3=396066248
Emax=10.0053941462 nrkept=1060 nrkeptmult=2616 
Kept: 1060 out of 4178, ratio=0.254
Memory used: 47 MB  Time elapsed: 98 s

 ***** [DMNRG] Iteration 21/66 (scale 0.00083772617549037929) ***** 
[1] xi(20)=0.999999602957 zeta(21)=0
Stats: nr=34 min=1 max=484 sum3=359180251
Safeguard: keep additional 2 states
Emax=10.0020404209 nrkept=924 nrkeptmult=2332 
Kept: 924 out of 3997, ratio=0.231
Memory used: 47 MB  Time elapsed: 99 s

 ***** [DMNRG] Iteration 22/66 (scale 0.00059236185946671895) ***** 
[1] xi(21)=0.999999801478 zeta(22)=0
Stats: nr=40 min=1 max=376 sum3=234913884
Safeguard: keep additional 1 states
Emax=10.0155455125 nrkept=866 nrkeptmult=2148 
Kept: 866 out of 3492, ratio=0.248
Memory used: 45 MB  Time elapsed: 99 s

 ***** [DMNRG] Iteration 23/66 (scale 0.00041886308774518965) ***** 
[1] xi(22)=0.999999900739 zeta(23)=0
Stats: nr=34 min=1 max=402 sum3=201263845
Safeguard: keep additional 1 states
Emax=10.0026731864 nrkept=846 nrkeptmult=2158 
Kept: 846 out of 3265, ratio=0.259
Memory used: 44 MB  Time elapsed: 100 s

 ***** [DMNRG] Iteration 24/66 (scale 0.00029618092973335947) ***** 
[1] xi(23)=0.99999995037 zeta(24)=0
Stats: nr=40 min=1 max=346 sum3=178056912
Emax=10.0177243709 nrkept=810 nrkeptmult=2008 
Kept: 810 out of 3198, ratio=0.253
Memory used: 44 MB  Time elapsed: 100 s

 ***** [DMNRG] Iteration 25/66 (scale 0.00020943154387259482) ***** 
[1] xi(24)=0.999999975185 zeta(25)=0
Stats: nr=34 min=1 max=374 sum3=161496433
Emax=10.0276516003 nrkept=816 nrkeptmult=2066 
Kept: 816 out of 3049, ratio=0.268
Memory used: 43 MB  Time elapsed: 100 s

 ***** [DMNRG] Iteration 26/66 (scale 0.00014809046486667974) ***** 
[1] xi(25)=0.999999987592 zeta(26)=0
Stats: nr=38 min=1 max=333 sum3=161834170
Emax=10.0100608586 nrkept=770 nrkeptmult=1936 
Kept: 770 out of 3082, ratio=0.25
Memory used: 44 MB  Time elapsed: 101 s

 ***** [DMNRG] Iteration 27/66 (scale 0.00010471577193629741) ***** 
[1] xi(26)=0.999999993796 zeta(27)=0
Stats: nr=34 min=1 max=354 sum3=139173931
Safeguard: keep additional 2 states
Emax=10.000601265 nrkept=789 nrkeptmult=1982 
Kept: 789 out of 2905, ratio=0.272
Memory used: 43 MB  Time elapsed: 101 s

 ***** [DMNRG] Iteration 28/66 (scale 7.4045232433339868e-05) ***** 
[1] xi(27)=0.999999996898 zeta(28)=0
Stats: nr=38 min=1 max=326 sum3=149194284
Emax=10.0159752848 nrkept=751 nrkeptmult=1897 
Kept: 751 out of 2976, ratio=0.252
Memory used: 43 MB  Time elapsed: 102 s

 ***** [DMNRG] Iteration 29/66 (scale 5.2357885968148706e-05) ***** 
[1] xi(28)=0.999999998449 zeta(29)=0
Stats: nr=34 min=1 max=343 sum3=127606636
Safeguard: keep additional 2 states
Emax=10.0031835774 nrkept=753 nrkeptmult=1888 
Kept: 753 out of 2836, ratio=0.266
Memory used: 44 MB  Time elapsed: 102 s

 ***** [DMNRG] Iteration 30/66 (scale 3.7022616216669934e-05) ***** 
[1] xi(29)=0.999999999225 zeta(30)=0
Stats: nr=38 min=1 max=314 sum3=132096960
Emax=10.0179729349 nrkept=693 nrkeptmult=1713 
Kept: 693 out of 2844, ratio=0.244
Memory used: 43 MB  Time elapsed: 102 s

 ***** [DMNRG] Iteration 31/66 (scale 2.6178942984074353e-05) ***** 
[1] xi(30)=0.999999999612 zeta(31)=0
Stats: nr=33 min=1 max=325 sum3=106369948
Safeguard: keep additional 3 states
Emax=10.0067640763 nrkept=709 nrkeptmult=1744 
Kept: 709 out of 2614, ratio=0.271
Memory used: 43 MB  Time elapsed: 103 s

 ***** [DMNRG] Iteration 32/66 (scale 1.8511308108334967e-05) ***** 
[1] xi(31)=0.999999999806 zeta(32)=0
Stats: nr=36 min=1 max=296 sum3=113592734
Emax=10.0011727924 nrkept=633 nrkeptmult=1541 
Kept: 633 out of 2672, ratio=0.237
Memory used: 43 MB  Time elapsed: 103 s

 ***** [DMNRG] Iteration 33/66 (scale 1.3089471492037176e-05) ***** 
[1] xi(32)=0.999999999903 zeta(33)=0
Stats: nr=33 min=1 max=303 sum3=83225817
Safeguard: keep additional 1 states
Emax=10.018903115 nrkept=656 nrkeptmult=1588 
Kept: 656 out of 2385, ratio=0.275
Memory used: 42 MB  Time elapsed: 103 s

 ***** [DMNRG] Iteration 34/66 (scale 9.2556540541674836e-06) ***** 
[1] xi(33)=0.999999999952 zeta(34)=0
Stats: nr=36 min=1 max=276 sum3=91791798
Safeguard: keep additional 2 states
Emax=10.052883151 nrkept=581 nrkeptmult=1383 
Kept: 581 out of 2466, ratio=0.236
Memory used: 42 MB  Time elapsed: 103 s

 ***** [DMNRG] Iteration 35/66 (scale 6.5447357460185882e-06) ***** 
[1] xi(34)=0.999999999976 zeta(35)=0
Stats: nr=33 min=1 max=283 sum3=65536626
Safeguard: keep additional 1 states
Emax=10.0368229375 nrkept=589 nrkeptmult=1402 
Kept: 589 out of 2178, ratio=0.27
Memory used: 42 MB  Time elapsed: 104 s

 ***** [DMNRG] Iteration 36/66 (scale 4.6278270270837418e-06) ***** 
[1] xi(35)=0.999999999988 zeta(36)=0
Stats: nr=36 min=1 max=251 sum3=67223648
Emax=10.0001206832 nrkept=523 nrkeptmult=1227 
Kept: 523 out of 2210, ratio=0.237
Memory used: 42 MB  Time elapsed: 104 s

 ***** [DMNRG] Iteration 37/66 (scale 3.2723678730092941e-06) ***** 
[1] xi(36)=0.999999999994 zeta(37)=0
Stats: nr=33 min=1 max=257 sum3=48651854
Emax=10.101122394 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1958, ratio=0.29
Memory used: 41 MB  Time elapsed: 104 s

 ***** [DMNRG] Iteration 38/66 (scale 2.3139135135418709e-06) ***** 
[1] xi(37)=0.999999999997 zeta(38)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.0464946451 nrkept=501 nrkeptmult=1157 
Kept: 501 out of 2124, ratio=0.236
Memory used: 41 MB  Time elapsed: 104 s

 ***** [DMNRG] Iteration 39/66 (scale 1.6361839365046471e-06) ***** 
[1] xi(38)=0.999999999998 zeta(39)=0
Stats: nr=33 min=1 max=245 sum3=42665008
Emax=10.5590939931 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1870, ratio=0.303
Memory used: 41 MB  Time elapsed: 104 s

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

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

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

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

 ***** [DMNRG] Iteration 44/66 (scale 2.8923918919273386e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5050375077 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 42 MB  Time elapsed: 105 s

 ***** [DMNRG] Iteration 45/66 (scale 2.0452299206308088e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8047350699 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1857, ratio=0.305
Memory used: 41 MB  Time elapsed: 105 s

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

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

 ***** [DMNRG] Iteration 48/66 (scale 7.2309797298183465e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5207756558 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 42 MB  Time elapsed: 106 s

 ***** [DMNRG] Iteration 49/66 (scale 5.113074801577022e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8315849979 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 41 MB  Time elapsed: 106 s

 ***** [DMNRG] Iteration 50/66 (scale 3.6154898649091733e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5234315011 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 42 MB  Time elapsed: 106 s

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

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

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

 ***** [DMNRG] Iteration 54/66 (scale 9.0387246622729332e-09) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5254201724 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 43 MB  Time elapsed: 107 s

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

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

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

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

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

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

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

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

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

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

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

 ***** [DMNRG] Iteration 66/66 (scale 1.4123007284801458e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260590775 nrkept=2154 nrkeptmult=5360 
Kept: 2154 out of 2154, ratio=1
Memory used: 49 MB  Time elapsed: 111 s

** Iteration completed.

Spectrum: A_d-A_d FDM
[A_d-A_d] pos=0.500000000001 neg=0.500000000001 sum= 1
mu1=9.35015953551e-16 mu2=0.000664752065275 mu3=3.78928659284e-17 mu4=0.000221607098398


Memory usage report
===================
  after diag: 91288 kB
after recalc: 102264 kB
 after split: 114652 kB
  after trim: 51208 kB
 after trunc: 91320 kB
      diag_h: 54992 kB

Peak usage: 114652 kB

Timing report
=============
         All: 111 s
          DM: 1 s
         FDM: 5 s
        diag: 58 s
    recalc d: 7 s
      spec d: 20 s
       Other: 18 s
