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 May 31 2013 at 17:57:13

Parallelization using MPI: Running on 1 processors.

Using Intel MKL library 10.3.9
Processor optimization: Intel(R) Core(TM) i7 Processor
max_threads=8 blas_max_threads=8 dynamic=1

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.5
omega0=1e-99
omega0_ratio=1
brcut=0
loggauss_b=0.3
gauss_b=0.1
limitLL=0
eta=0
finite=1
bins=100
discard_trim=1e-16
discard_immediately=1e-16
NN1=0
NN2even=1
NN2avg=0
NNtanh=0
goodE=2.1
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 SZd
specs=
specd=A_d-A_d-u A_d-A_d-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=1
cfs=1
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=1
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: QSZ
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 	MAGFIELD	 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 	 extraB=	0.0003
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.75
>>  ops=	A_d SZd
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=9
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||SZd||> (s)
Reading <||||> (z)

# $Id: nrg.cc,v 1.31 2012/10/26 10:18:22 rokzitko Exp rokzitko $ QSZ
# 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: 9
(-2 1) 1 states
0.0437080950522 
(-1 0) 2 states
0.0206794355226 0.0624916324536 
(-1 2) 2 states
0.0208152711796 0.0626030854643 
(0 -1) 1 states
0.039462972924 
(0 1) 4 states
0 0.0395866536738 0.0437080950522 0.0832946758941 
(0 3) 1 states
0.0397102615917 
(1 0) 2 states
0.0206794355226 0.0624916324536 
(1 2) 2 states
0.0208152711796 0.0626030854643 
(2 1) 1 states
0.0437080950522 
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-d FT
Spectrum: A_d-A_d-u FT

Before NRG iteration (N=-1)
<SZd>=-3.09725454829e-05

 ***** [NRG] Iteration 1/66 (scale 0.85783160370214839) ***** 
[1] xi(0)=0.684778584642 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=2.74762181274 nrkept=64 nrkeptmult=64 
Total energy=-1.17607327719  Egs=-1.3149999219
Kept: 64 out of 64, ratio=1
<SZd>=-4.38362170828e-05
Memory used: 11 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/66 (scale 0.6065785440939202) ***** 
[1] xi(1)=0.833755393818 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=5.22193534821 nrkept=256 nrkeptmult=256 
Total energy=-1.58125702224  Egs=-0.667982323129
Kept: 256 out of 256, ratio=1
<SZd>=-6.20604398371e-05
Memory used: 12 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/66 (scale 0.4289158018510742) ***** 
[1] xi(2)=0.947046165312 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Emax=10.4050027078 nrkept=1021 nrkeptmult=1021 
Total energy=-2.23156097179  Egs=-1.51615759258
Kept: 1021 out of 1024, ratio=0.997
<SZd>=-8.83195394462e-05
Memory used: 16 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/66 (scale 0.3032892720469601) ***** 
[1] xi(3)=0.972572268687 zeta(4)=0
Stats: nr=49 min=1 max=397 sum3=227060008
Emax=10.7309916699 nrkept=3041 nrkeptmult=3041 
Total energy=-2.48205908928  Egs=-0.825937942998
Kept: 3041 out of 4084, ratio=0.745
<SZd>=-0.000124556835405
Memory used: 47 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/66 (scale 0.2144579009255371) ***** 
[1] xi(4)=0.980132362648 zeta(5)=0
Stats: nr=64 min=1 max=890 sum3=4287484628
Safeguard: keep additional 7 states
Emax=10.4296513476 nrkept=4268 nrkeptmult=4268 
Total energy=-2.80358663989  Egs=-1.49925719327
Kept: 4268 out of 12164, ratio=0.351
<SZd>=-0.000179736619769
Memory used: 77 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 6/66 (scale 0.15164463602348005) ***** 
[1] xi(5)=0.988164304343 zeta(6)=0
Stats: nr=77 min=1 max=1352 sum3=10854929552
Safeguard: keep additional 3 states
Emax=10.7114802845 nrkept=5284 nrkeptmult=5284 
Total energy=-2.93591433664  Egs=-0.872617061957
Kept: 5284 out of 17072, ratio=0.31
<SZd>=-0.000250748904291
Memory used: 106 MB  Time elapsed: 16 s

 ***** [NRG] Iteration 7/66 (scale 0.10722895046276855) ***** 
[1] xi(6)=0.993732382824 zeta(7)=0
Stats: nr=80 min=1 max=1481 sum3=20456698096
Safeguard: keep additional 5 states
Emax=10.0288249044 nrkept=4474 nrkeptmult=4474 
Total energy=-3.09480252742  Egs=-1.48176579263
Kept: 4474 out of 21136, ratio=0.212
<SZd>=-0.000356500336755
Memory used: 102 MB  Time elapsed: 34 s

 ***** [NRG] Iteration 8/66 (scale 0.075822318011740025) ***** 
[1] xi(7)=0.996801670793 zeta(8)=0
Stats: nr=77 min=1 max=1434 sum3=12827421322
Safeguard: keep additional 1 states
Emax=10.0199069004 nrkept=4882 nrkeptmult=4882 
Total energy=-3.16363809853  Egs=-0.90785368893
Kept: 4882 out of 17896, ratio=0.273
<SZd>=-0.000509291678763
Memory used: 120 MB  Time elapsed: 46 s

 ***** [NRG] Iteration 9/66 (scale 0.053614475231384275) ***** 
[1] xi(8)=0.998386814476 zeta(9)=0
Stats: nr=80 min=1 max=1377 sum3=16368825400
Safeguard: keep additional 3 states
Emax=10.0055873511 nrkept=4096 nrkeptmult=4096 
Total energy=-3.2415493825  Egs=-1.45317628554
Kept: 4096 out of 19528, ratio=0.21
<SZd>=-0.000728048960907
Memory used: 108 MB  Time elapsed: 61 s

 ***** [NRG] Iteration 10/66 (scale 0.037911159005870013) ***** 
[1] xi(9)=0.999190082691 zeta(10)=0
Stats: nr=77 min=1 max=1312 sum3=9823752472
Safeguard: keep additional 2 states
Emax=10.0140824044 nrkept=4603 nrkeptmult=4603 
Total energy=-3.27708904555  Egs=-0.937445965475
Kept: 4603 out of 16384, ratio=0.281
<SZd>=-0.00104183923645
Memory used: 119 MB  Time elapsed: 71 s

 ***** [NRG] Iteration 11/66 (scale 0.026807237615692137) ***** 
[1] xi(10)=0.999594228576 zeta(11)=0
Stats: nr=76 min=1 max=1310 sum3=13968677740
Safeguard: keep additional 1 states
Emax=10.0000011545 nrkept=3836 nrkeptmult=3836 
Total energy=-3.31515842375  Egs=-1.42011566974
Kept: 3836 out of 18412, ratio=0.208
<SZd>=-0.00151025754958
Memory used: 111 MB  Time elapsed: 84 s

 ***** [NRG] Iteration 12/66 (scale 0.018955579502935006) ***** 
[1] xi(11)=0.999796913201 zeta(12)=0
Stats: nr=77 min=1 max=1236 sum3=8171184212
Safeguard: keep additional 3 states
Emax=10.0022740135 nrkept=4084 nrkeptmult=4084 
Total energy=-3.33343419791  Egs=-0.964136926299
Kept: 4084 out of 15344, ratio=0.266
<SZd>=-0.0021792519847
Memory used: 113 MB  Time elapsed: 92 s

 ***** [NRG] Iteration 13/66 (scale 0.013403618807846069) ***** 
[1] xi(12)=0.999898406582 zeta(13)=0
Stats: nr=74 min=1 max=1179 sum3=10077672808
Safeguard: keep additional 2 states
Emax=10.0000460926 nrkept=3791 nrkeptmult=3791 
Total energy=-3.35203779515  Egs=-1.38795332126
Kept: 3791 out of 16336, ratio=0.232
<SZd>=-0.00318021461519
Memory used: 109 MB  Time elapsed: 103 s

 ***** [NRG] Iteration 14/66 (scale 0.0094777897514675032) ***** 
[1] xi(13)=0.999949190818 zeta(14)=0
Stats: nr=77 min=1 max=1219 sum3=7853778362
Safeguard: keep additional 2 states
Emax=10.0000467205 nrkept=3799 nrkeptmult=3799 
Total energy=-3.36138565415  Egs=-0.986291027572
Kept: 3799 out of 15164, ratio=0.251
<SZd>=-0.00469389872407
Memory used: 110 MB  Time elapsed: 111 s

 ***** [NRG] Iteration 15/66 (scale 0.0067018094039230343) ***** 
[1] xi(14)=0.999974592294 zeta(15)=0
Stats: nr=68 min=1 max=1110 sum3=8257321432
Emax=10.0349903173 nrkept=3641 nrkeptmult=3641 
Total energy=-3.37050264598  Egs=-1.36037766662
Kept: 3641 out of 15196, ratio=0.24
<SZd>=-0.00692470249648
Memory used: 108 MB  Time elapsed: 120 s

 ***** [NRG] Iteration 16/66 (scale 0.0047388948757337516) ***** 
[1] xi(15)=0.999987295369 zeta(16)=0
Stats: nr=73 min=1 max=1182 sum3=7101640322
Emax=10.003970589 nrkept=3469 nrkeptmult=3469 
Total energy=-3.37525339676  Egs=-1.00250182843
Kept: 3469 out of 14564, ratio=0.238
<SZd>=-0.0103117490022
Memory used: 106 MB  Time elapsed: 127 s

 ***** [NRG] Iteration 17/66 (scale 0.0033509047019615172) ***** 
[1] xi(16)=0.99999364749 zeta(17)=0
Stats: nr=64 min=1 max=1022 sum3=6469133410
Emax=10.0009156971 nrkept=3090 nrkeptmult=3090 
Total energy=-3.37974181122  Egs=-1.33946347445
Kept: 3090 out of 13876, ratio=0.223
<SZd>=-0.0153074751256
Memory used: 110 MB  Time elapsed: 134 s

 ***** [NRG] Iteration 18/66 (scale 0.0023694474378668758) ***** 
[1] xi(17)=0.999996823696 zeta(18)=0
Stats: nr=73 min=1 max=1019 sum3=4478809368
Emax=10.0017720061 nrkept=3036 nrkeptmult=3036 
Total energy=-3.38214093158  Egs=-1.01252314132
Kept: 3036 out of 12360, ratio=0.246
<SZd>=-0.022493646331
Memory used: 110 MB  Time elapsed: 140 s

 ***** [NRG] Iteration 19/66 (scale 0.0016754523509807586) ***** 
[1] xi(18)=0.999998411836 zeta(19)=0
Stats: nr=62 min=1 max=902 sum3=4370490330
Safeguard: keep additional 2 states
Emax=10.002610018 nrkept=2692 nrkeptmult=2692 
Total energy=-3.38436360587  Egs=-1.32661146084
Kept: 2692 out of 12144, ratio=0.222
<SZd>=-0.0324718972055
Memory used: 110 MB  Time elapsed: 145 s

 ***** [NRG] Iteration 20/66 (scale 0.0011847237189334379) ***** 
[1] xi(19)=0.999999205915 zeta(20)=0
Stats: nr=73 min=1 max=889 sum3=2991491638
Safeguard: keep additional 1 states
Emax=10.0009620502 nrkept=2586 nrkeptmult=2586 
Total energy=-3.38556750353  Egs=-1.01618431594
Kept: 2586 out of 10768, ratio=0.24
<SZd>=-0.0459323873429
Memory used: 108 MB  Time elapsed: 149 s

 ***** [NRG] Iteration 21/66 (scale 0.00083772617549037929) ***** 
[1] xi(20)=0.999999602957 zeta(21)=0
Stats: nr=62 min=1 max=782 sum3=2748041580
Safeguard: keep additional 1 states
Emax=10.003959981 nrkept=2228 nrkeptmult=2228 
Total energy=-3.38667594498  Egs=-1.32315484975
Kept: 2228 out of 10344, ratio=0.215
<SZd>=-0.0637573987131
Memory used: 106 MB  Time elapsed: 152 s

 ***** [NRG] Iteration 22/66 (scale 0.00059236185946671895) ***** 
[1] xi(21)=0.999999801478 zeta(22)=0
Stats: nr=71 min=1 max=743 sum3=1723641866
Safeguard: keep additional 1 states
Emax=10.0025791 nrkept=2103 nrkeptmult=2103 
Total energy=-3.38727654232  Egs=-1.01390278485
Kept: 2103 out of 8912, ratio=0.236
<SZd>=-0.0872266479405
Memory used: 104 MB  Time elapsed: 155 s

 ***** [NRG] Iteration 23/66 (scale 0.00041886308774518965) ***** 
[1] xi(22)=0.999999900739 zeta(23)=0
Stats: nr=62 min=1 max=654 sum3=1507450476
Safeguard: keep additional 1 states
Emax=10.0012705678 nrkept=1989 nrkeptmult=1989 
Total energy=-3.38783296781  Egs=-1.32841852794
Kept: 1989 out of 8412, ratio=0.236
<SZd>=-0.117486315936
Memory used: 104 MB  Time elapsed: 157 s

 ***** [NRG] Iteration 24/66 (scale 0.00029618092973335947) ***** 
[1] xi(23)=0.99999995037 zeta(24)=0
Stats: nr=69 min=1 max=672 sum3=1257195180
Safeguard: keep additional 1 states
Emax=10.0145373054 nrkept=2013 nrkeptmult=2013 
Total energy=-3.38813156725  Egs=-1.00816565148
Kept: 2013 out of 7956, ratio=0.253
<SZd>=-0.155299043008
Memory used: 109 MB  Time elapsed: 159 s

 ***** [NRG] Iteration 25/66 (scale 0.00020943154387259482) ***** 
[1] xi(24)=0.999999975185 zeta(25)=0
Stats: nr=62 min=1 max=622 sum3=1320642150
Safeguard: keep additional 1 states
Emax=10.0181329769 nrkept=1766 nrkeptmult=1766 
Total energy=-3.38841177821  Egs=-1.33795966418
Kept: 1766 out of 8052, ratio=0.219
<SZd>=-0.200135587416
Memory used: 107 MB  Time elapsed: 161 s

 ***** [NRG] Iteration 26/66 (scale 0.00014809046486667974) ***** 
[1] xi(25)=0.999999987592 zeta(26)=0
Stats: nr=69 min=1 max=599 sum3=895388342
Safeguard: keep additional 1 states
Emax=10.0028225723 nrkept=1746 nrkeptmult=1746 
Total energy=-3.3885601964  Egs=-1.00221304607
Kept: 1746 out of 7064, ratio=0.247
<SZd>=-0.248203291997
Memory used: 110 MB  Time elapsed: 162 s

 ***** [NRG] Iteration 27/66 (scale 0.00010471577193629741) ***** 
[1] xi(26)=0.999999993796 zeta(27)=0
Stats: nr=61 min=1 max=566 sum3=902205654
Safeguard: keep additional 3 states
Emax=10.0015470658 nrkept=1536 nrkeptmult=1536 
Total energy=-3.38870121069  Egs=-1.34663845586
Kept: 1536 out of 6984, ratio=0.22
<SZd>=-0.293182758055
Memory used: 110 MB  Time elapsed: 164 s

 ***** [NRG] Iteration 28/66 (scale 7.4045232433339868e-05) ***** 
[1] xi(27)=0.999999996898 zeta(28)=0
Stats: nr=65 min=1 max=526 sum3=612952602
Safeguard: keep additional 1 states
Emax=10.0045042769 nrkept=1475 nrkeptmult=1475 
Total energy=-3.38877509103  Egs=-0.997773098668
Kept: 1475 out of 6144, ratio=0.24
<SZd>=-0.327632689926
Memory used: 113 MB  Time elapsed: 165 s

 ***** [NRG] Iteration 29/66 (scale 5.2357885968148706e-05) ***** 
[1] xi(28)=0.999999998449 zeta(29)=0
Stats: nr=59 min=1 max=500 sum3=567156254
Safeguard: keep additional 2 states
Emax=10.0016466226 nrkept=1302 nrkeptmult=1302 
Total energy=-3.38884589932  Egs=-1.3523901052
Kept: 1302 out of 5900, ratio=0.221
<SZd>=-0.348449098719
Memory used: 113 MB  Time elapsed: 166 s

 ***** [NRG] Iteration 30/66 (scale 3.7022616216669934e-05) ***** 
[1] xi(29)=0.999999999225 zeta(30)=0
Stats: nr=59 min=1 max=449 sum3=396478788
Safeguard: keep additional 1 states
Emax=10.0013765028 nrkept=1255 nrkeptmult=1255 
Total energy=-3.3888827372  Egs=-0.995010228603
Kept: 1255 out of 5208, ratio=0.241
<SZd>=-0.358160225136
Memory used: 115 MB  Time elapsed: 166 s

 ***** [NRG] Iteration 31/66 (scale 2.6178942984074353e-05) ***** 
[1] xi(30)=0.999999999612 zeta(31)=0
Stats: nr=57 min=1 max=449 sum3=372148378
Emax=10.004433867 nrkept=1122 nrkeptmult=1122 
Total energy=-3.38891822642  Egs=-1.35563981826
Kept: 1122 out of 5020, ratio=0.224
<SZd>=-0.361977546615
Memory used: 116 MB  Time elapsed: 167 s

 ***** [NRG] Iteration 32/66 (scale 1.8511308108334967e-05) ***** 
[1] xi(31)=0.999999999806 zeta(32)=0
Stats: nr=58 min=1 max=381 sum3=260161692
Safeguard: keep additional 3 states
Emax=10.0030030674 nrkept=1141 nrkeptmult=1141 
Total energy=-3.3889366167  Egs=-0.993462313621
Kept: 1141 out of 4488, ratio=0.254
<SZd>=-0.363436347712
Memory used: 114 MB  Time elapsed: 168 s

 ***** [NRG] Iteration 33/66 (scale 1.3089471492037176e-05) ***** 
[1] xi(32)=0.999999999903 zeta(33)=0
Stats: nr=57 min=1 max=418 sum3=287422162
Safeguard: keep additional 3 states
Emax=10.0020460748 nrkept=1083 nrkeptmult=1083 
Total energy=-3.38895438369  Egs=-1.35734978388
Kept: 1083 out of 4564, ratio=0.237
<SZd>=-0.364071638175
Memory used: 115 MB  Time elapsed: 168 s

 ***** [NRG] Iteration 34/66 (scale 9.2556540541674836e-06) ***** 
[1] xi(33)=0.999999999952 zeta(34)=0
Stats: nr=58 min=1 max=371 sum3=237920658
Safeguard: keep additional 1 states
Emax=10.018649788 nrkept=1116 nrkeptmult=1116 
Total energy=-3.38896357126  Egs=-0.992643228433
Kept: 1116 out of 4332, ratio=0.258
<SZd>=-0.364368496161
Memory used: 115 MB  Time elapsed: 169 s

 ***** [NRG] Iteration 35/66 (scale 6.5447357460185882e-06) ***** 
[1] xi(34)=0.999999999976 zeta(35)=0
Stats: nr=57 min=1 max=412 sum3=271853526
Safeguard: keep additional 1 states
Emax=10.0252201681 nrkept=1075 nrkeptmult=1075 
Total energy=-3.38897246047  Egs=-1.3582236732
Kept: 1075 out of 4464, ratio=0.241
<SZd>=-0.364513162597
Memory used: 116 MB  Time elapsed: 169 s

 ***** [NRG] Iteration 36/66 (scale 4.6278270270837418e-06) ***** 
[1] xi(35)=0.999999999988 zeta(36)=0
Stats: nr=58 min=1 max=371 sum3=234091648
Safeguard: keep additional 2 states
Emax=10.0535407837 nrkept=1117 nrkeptmult=1117 
Total energy=-3.3889770523  Egs=-0.992222128703
Kept: 1117 out of 4300, ratio=0.26
<SZd>=-0.36458457104
Memory used: 115 MB  Time elapsed: 170 s

 ***** [NRG] Iteration 37/66 (scale 3.2723678730092941e-06) ***** 
[1] xi(36)=0.999999999994 zeta(37)=0
Stats: nr=57 min=1 max=413 sum3=273706324
Safeguard: keep additional 1 states
Emax=10.0005691486 nrkept=1069 nrkeptmult=1069 
Total energy=-3.38898149836  Egs=-1.35866494795
Kept: 1069 out of 4468, ratio=0.239
<SZd>=-0.364619964524
Memory used: 116 MB  Time elapsed: 171 s

 ***** [NRG] Iteration 38/66 (scale 2.3139135135418709e-06) ***** 
[1] xi(37)=0.999999999997 zeta(38)=0
Stats: nr=58 min=1 max=371 sum3=232214488
Emax=10.0679012837 nrkept=1115 nrkeptmult=1115 
Total energy=-3.38898379378  Egs=-0.992008674642
Kept: 1115 out of 4276, ratio=0.261
<SZd>=-0.364637798561
Memory used: 115 MB  Time elapsed: 171 s

 ***** [NRG] Iteration 39/66 (scale 1.6361839365046471e-06) ***** 
[1] xi(38)=0.999999999998 zeta(39)=0
Stats: nr=57 min=1 max=413 sum3=272139644
Emax=10.0023974172 nrkept=1060 nrkeptmult=1060 
Total energy=-3.38898601717  Egs=-1.35888661167
Kept: 1060 out of 4460, ratio=0.238
<SZd>=-0.364646627935
Memory used: 116 MB  Time elapsed: 172 s

 ***** [NRG] Iteration 40/66 (scale 1.1569567567709354e-06) ***** 
[1] xi(39)=0.999999999999 zeta(40)=0
Stats: nr=58 min=1 max=371 sum3=231642436
Safeguard: keep additional 2 states
Emax=10.0725678988 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38898716475  Egs=-0.991901221302
Kept: 1117 out of 4240, ratio=0.263
<SZd>=-0.364651104515
Memory used: 116 MB  Time elapsed: 172 s

 ***** [NRG] Iteration 41/66 (scale 8.1809196825232353e-07) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Emax=10.0128597108 nrkept=1060 nrkeptmult=1060 
Total energy=-3.38898827654  Egs=-1.35899769276
Kept: 1060 out of 4468, ratio=0.237
<SZd>=-0.364653309468
Memory used: 119 MB  Time elapsed: 173 s

 ***** [NRG] Iteration 42/66 (scale 5.7847837838546772e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=58 min=1 max=371 sum3=231642436
Safeguard: keep additional 2 states
Emax=10.0743813352 nrkept=1117 nrkeptmult=1117 
Total energy=-3.3889888503  Egs=-0.991847313135
Kept: 1117 out of 4240, ratio=0.263
<SZd>=-0.364654450898
Memory used: 119 MB  Time elapsed: 173 s

 ***** [NRG] Iteration 43/66 (scale 4.0904598412616176e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0165961783 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38898940622  Egs=-1.35905329467
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655009996
Memory used: 120 MB  Time elapsed: 174 s

 ***** [NRG] Iteration 44/66 (scale 2.8923918919273386e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.075288325 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38898969309  Egs=-0.991820313687
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655299687
Memory used: 119 MB  Time elapsed: 174 s

 ***** [NRG] Iteration 45/66 (scale 2.0452299206308088e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.017865681 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38898997105  Egs=-1.35908111086
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655441523
Memory used: 120 MB  Time elapsed: 175 s

 ***** [NRG] Iteration 46/66 (scale 1.4461959459636693e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0757419091 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38899011449  Egs=-0.991806802618
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655517856
Memory used: 120 MB  Time elapsed: 176 s

 ***** [NRG] Iteration 47/66 (scale 1.0226149603154044e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0184767017 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899025347  Egs=-1.35909502276
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655554364
Memory used: 121 MB  Time elapsed: 176 s

 ***** [NRG] Iteration 48/66 (scale 7.2309797298183465e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0759687261 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38899032519  Egs=-0.991800044239
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655575406
Memory used: 121 MB  Time elapsed: 177 s

 ***** [NRG] Iteration 49/66 (scale 5.113074801577022e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0187806375 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899039468  Egs=-1.35910197966
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655585057
Memory used: 122 MB  Time elapsed: 177 s

 ***** [NRG] Iteration 50/66 (scale 3.6154898649091733e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0760821412 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38899043054  Egs=-0.991796664332
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655591296
Memory used: 123 MB  Time elapsed: 178 s

 ***** [NRG] Iteration 51/66 (scale 2.556537400788511e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0189324201 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899046528  Egs=-1.35910545836
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655593969
Memory used: 124 MB  Time elapsed: 178 s

 ***** [NRG] Iteration 52/66 (scale 1.8077449324545866e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761388504 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38899048321  Egs=-0.991794974191
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655596018
Memory used: 124 MB  Time elapsed: 179 s

 ***** [NRG] Iteration 53/66 (scale 1.2782687003942555e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190082825 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899050059  Egs=-1.35910719779
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655596816
Memory used: 126 MB  Time elapsed: 180 s

 ***** [NRG] Iteration 54/66 (scale 9.0387246622729332e-09) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761672054 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38899050955  Egs=-0.991794129065
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655597573
Memory used: 125 MB  Time elapsed: 180 s

 ***** [NRG] Iteration 55/66 (scale 6.3913435019712775e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190462082 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899051824  Egs=-1.35910806752
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655597837
Memory used: 127 MB  Time elapsed: 181 s

 ***** [NRG] Iteration 56/66 (scale 4.5193623311364666e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.076181383 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38899052272  Egs=-0.991793706481
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655598148
Memory used: 128 MB  Time elapsed: 181 s

 ***** [NRG] Iteration 57/66 (scale 3.1956717509856388e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190651699 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899052706  Egs=-1.35910850241
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655598247
Memory used: 129 MB  Time elapsed: 182 s

 ***** [NRG] Iteration 58/66 (scale 2.2596811655682333e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761884718 nrkept=1117 nrkeptmult=1117 
Total energy=-3.3889905293  Egs=-0.991793495176
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655598386
Memory used: 128 MB  Time elapsed: 182 s

 ***** [NRG] Iteration 59/66 (scale 1.5978358754928194e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190746505 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899053148  Egs=-1.35910871986
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655598427
Memory used: 129 MB  Time elapsed: 183 s

 ***** [NRG] Iteration 60/66 (scale 1.1298405827841166e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761920162 nrkept=1117 nrkeptmult=1117 
Total energy=-3.3889905326  Egs=-0.991793389512
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655598492
Memory used: 129 MB  Time elapsed: 183 s

 ***** [NRG] Iteration 61/66 (scale 7.9891793774640969e-10) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190793906 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899053368  Egs=-1.3591088286
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.36465559851
Memory used: 130 MB  Time elapsed: 184 s

 ***** [NRG] Iteration 62/66 (scale 5.6492029139205832e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761937884 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38899053424  Egs=-0.991793336669
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655598542
Memory used: 132 MB  Time elapsed: 185 s

 ***** [NRG] Iteration 63/66 (scale 3.9945896887320485e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190817607 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899053478  Egs=-1.35910888299
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364655598551
Memory used: 132 MB  Time elapsed: 185 s

 ***** [NRG] Iteration 64/66 (scale 2.8246014569602916e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761946744 nrkept=1117 nrkeptmult=1117 
Total energy=-3.38899053506  Egs=-0.991793310238
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364655598566
Memory used: 132 MB  Time elapsed: 186 s

 ***** [NRG] Iteration 65/66 (scale 1.9972948443660242e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190829456 nrkept=1061 nrkeptmult=1061 
Total energy=-3.38899053534  Egs=-1.35910891019
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.36465559857
Memory used: 132 MB  Time elapsed: 186 s

 ***** [NRG] Iteration 66/66 (scale 1.4123007284801458e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761951174 nrkept=4244 nrkeptmult=4244 
Total energy=-3.38899053548  Egs=-0.991793297011
Kept: 4244 out of 4244, ratio=1
<SZd>=-0.364655598581
Memory used: 163 MB  Time elapsed: 187 s

** Iteration completed.

Spectrum: A_d-A_d-d FT
[A_d-A_d-d] pos=0.103451508816 neg=0.908152958751 sum= 1.01160446757
mu1=-0.00428034886173 mu2=0.000662892197833 mu3=-6.65106651056e-06 mu4=0.000220343164929

Spectrum: A_d-A_d-u FT
[A_d-A_d-u] pos=0.908152958752 neg=0.103451508816 sum= 1.01160446757
mu1=0.00428034886173 mu2=0.000662892197833 mu3=6.65106651054e-06 mu4=0.000220343164929


Total energy:               -3.3889905354765
TKW=0
Grand canonical Z=              6.274464397285286
Zft (last shell)=                6.27446353201182
Calculating density matrices.
Calculating full density matrices.

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QSZ
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 	MAGFIELD	 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 	 extraB=	0.0003
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.75
>>  ops=	A_d SZd
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=9
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||SZd||> (s)
Reading <||||> (z)
Number of invariant subspaces: 9
(-2 1) 1 states
0.0437080950522 
(-1 0) 2 states
0.0206794355226 0.0624916324536 
(-1 2) 2 states
0.0208152711796 0.0626030854643 
(0 -1) 1 states
0.039462972924 
(0 1) 4 states
0 0.0395866536738 0.0437080950522 0.0832946758941 
(0 3) 1 states
0.0397102615917 
(1 0) 2 states
0.0206794355226 0.0624916324536 
(1 2) 2 states
0.0208152711796 0.0626030854643 
(2 1) 1 states
0.0437080950522 
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-d DMNRG
Spectrum: A_d-A_d-d CFS
Spectrum: A_d-A_d-d FDM
Spectrum: A_d-A_d-u DMNRG
Spectrum: A_d-A_d-u CFS
Spectrum: A_d-A_d-u FDM

 ***** [DMNRG] Iteration 1/66 (scale 0.85783160370214839) ***** 
[1] xi(0)=0.684778584642 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=2.74762181274 nrkept=64 nrkeptmult=64 
Kept: 64 out of 64, ratio=1
Memory used: 190 MB  Time elapsed: 210 s

 ***** [DMNRG] Iteration 2/66 (scale 0.6065785440939202) ***** 
[1] xi(1)=0.833755393818 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=5.22193534821 nrkept=256 nrkeptmult=256 
Kept: 256 out of 256, ratio=1
Memory used: 190 MB  Time elapsed: 210 s

 ***** [DMNRG] Iteration 3/66 (scale 0.4289158018510742) ***** 
[1] xi(2)=0.947046165312 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Emax=10.4050027078 nrkept=1021 nrkeptmult=1021 
Kept: 1021 out of 1024, ratio=0.997
Memory used: 190 MB  Time elapsed: 210 s

 ***** [DMNRG] Iteration 4/66 (scale 0.3032892720469601) ***** 
[1] xi(3)=0.972572268687 zeta(4)=0
Stats: nr=49 min=1 max=397 sum3=227060008
Emax=10.7309916699 nrkept=3041 nrkeptmult=3041 
Kept: 3041 out of 4084, ratio=0.745
Memory used: 202 MB  Time elapsed: 211 s

 ***** [DMNRG] Iteration 5/66 (scale 0.2144579009255371) ***** 
[1] xi(4)=0.980132362648 zeta(5)=0
Stats: nr=64 min=1 max=890 sum3=4287484628
Safeguard: keep additional 7 states
Emax=10.4296513476 nrkept=4268 nrkeptmult=4268 
Kept: 4268 out of 12164, ratio=0.351
Memory used: 228 MB  Time elapsed: 213 s

 ***** [DMNRG] Iteration 6/66 (scale 0.15164463602348005) ***** 
[1] xi(5)=0.988164304343 zeta(6)=0
Stats: nr=77 min=1 max=1352 sum3=10854929552
Safeguard: keep additional 3 states
Emax=10.7114802845 nrkept=5284 nrkeptmult=5284 
Kept: 5284 out of 17072, ratio=0.31
Memory used: 257 MB  Time elapsed: 217 s

 ***** [DMNRG] Iteration 7/66 (scale 0.10722895046276855) ***** 
[1] xi(6)=0.993732382824 zeta(7)=0
Stats: nr=80 min=1 max=1481 sum3=20456698096
Safeguard: keep additional 5 states
Emax=10.0288249044 nrkept=4474 nrkeptmult=4474 
Kept: 4474 out of 21136, ratio=0.212
Memory used: 239 MB  Time elapsed: 223 s

 ***** [DMNRG] Iteration 8/66 (scale 0.075822318011740025) ***** 
[1] xi(7)=0.996801670793 zeta(8)=0
Stats: nr=77 min=1 max=1434 sum3=12827421322
Safeguard: keep additional 1 states
Emax=10.0199069004 nrkept=4882 nrkeptmult=4882 
Kept: 4882 out of 17896, ratio=0.273
Memory used: 249 MB  Time elapsed: 227 s

 ***** [DMNRG] Iteration 9/66 (scale 0.053614475231384275) ***** 
[1] xi(8)=0.998386814476 zeta(9)=0
Stats: nr=80 min=1 max=1377 sum3=16368825400
Safeguard: keep additional 3 states
Emax=10.0055873511 nrkept=4096 nrkeptmult=4096 
Kept: 4096 out of 19528, ratio=0.21
Memory used: 240 MB  Time elapsed: 232 s

 ***** [DMNRG] Iteration 10/66 (scale 0.037911159005870013) ***** 
[1] xi(9)=0.999190082691 zeta(10)=0
Stats: nr=77 min=1 max=1312 sum3=9823752472
Safeguard: keep additional 2 states
Emax=10.0140824044 nrkept=4603 nrkeptmult=4603 
Kept: 4603 out of 16384, ratio=0.281
Memory used: 250 MB  Time elapsed: 236 s

 ***** [DMNRG] Iteration 11/66 (scale 0.026807237615692137) ***** 
[1] xi(10)=0.999594228576 zeta(11)=0
Stats: nr=76 min=1 max=1310 sum3=13968677740
Safeguard: keep additional 1 states
Emax=10.0000011545 nrkept=3836 nrkeptmult=3836 
Kept: 3836 out of 18412, ratio=0.208
Memory used: 238 MB  Time elapsed: 240 s

 ***** [DMNRG] Iteration 12/66 (scale 0.018955579502935006) ***** 
[1] xi(11)=0.999796913201 zeta(12)=0
Stats: nr=77 min=1 max=1236 sum3=8171184212
Safeguard: keep additional 3 states
Emax=10.0022740135 nrkept=4084 nrkeptmult=4084 
Kept: 4084 out of 15344, ratio=0.266
Memory used: 239 MB  Time elapsed: 243 s

 ***** [DMNRG] Iteration 13/66 (scale 0.013403618807846069) ***** 
[1] xi(12)=0.999898406582 zeta(13)=0
Stats: nr=74 min=1 max=1179 sum3=10077672808
Safeguard: keep additional 2 states
Emax=10.0000460926 nrkept=3791 nrkeptmult=3791 
Kept: 3791 out of 16336, ratio=0.232
Memory used: 237 MB  Time elapsed: 247 s

 ***** [DMNRG] Iteration 14/66 (scale 0.0094777897514675032) ***** 
[1] xi(13)=0.999949190818 zeta(14)=0
Stats: nr=77 min=1 max=1219 sum3=7853778362
Safeguard: keep additional 2 states
Emax=10.0000467205 nrkept=3799 nrkeptmult=3799 
Kept: 3799 out of 15164, ratio=0.251
Memory used: 234 MB  Time elapsed: 250 s

 ***** [DMNRG] Iteration 15/66 (scale 0.0067018094039230343) ***** 
[1] xi(14)=0.999974592294 zeta(15)=0
Stats: nr=68 min=1 max=1110 sum3=8257321432
Emax=10.0349903173 nrkept=3641 nrkeptmult=3641 
Kept: 3641 out of 15196, ratio=0.24
Memory used: 233 MB  Time elapsed: 253 s

 ***** [DMNRG] Iteration 16/66 (scale 0.0047388948757337516) ***** 
[1] xi(15)=0.999987295369 zeta(16)=0
Stats: nr=73 min=1 max=1182 sum3=7101640322
Emax=10.003970589 nrkept=3469 nrkeptmult=3469 
Kept: 3469 out of 14564, ratio=0.238
Memory used: 227 MB  Time elapsed: 256 s

 ***** [DMNRG] Iteration 17/66 (scale 0.0033509047019615172) ***** 
[1] xi(16)=0.99999364749 zeta(17)=0
Stats: nr=64 min=1 max=1022 sum3=6469133410
Emax=10.0009156971 nrkept=3090 nrkeptmult=3090 
Kept: 3090 out of 13876, ratio=0.223
Memory used: 226 MB  Time elapsed: 258 s

 ***** [DMNRG] Iteration 18/66 (scale 0.0023694474378668758) ***** 
[1] xi(17)=0.999996823696 zeta(18)=0
Stats: nr=73 min=1 max=1019 sum3=4478809368
Emax=10.0017720061 nrkept=3036 nrkeptmult=3036 
Kept: 3036 out of 12360, ratio=0.246
Memory used: 223 MB  Time elapsed: 260 s

 ***** [DMNRG] Iteration 19/66 (scale 0.0016754523509807586) ***** 
[1] xi(18)=0.999998411836 zeta(19)=0
Stats: nr=62 min=1 max=902 sum3=4370490330
Safeguard: keep additional 2 states
Emax=10.002610018 nrkept=2692 nrkeptmult=2692 
Kept: 2692 out of 12144, ratio=0.222
Memory used: 219 MB  Time elapsed: 262 s

 ***** [DMNRG] Iteration 20/66 (scale 0.0011847237189334379) ***** 
[1] xi(19)=0.999999205915 zeta(20)=0
Stats: nr=73 min=1 max=889 sum3=2991491638
Safeguard: keep additional 1 states
Emax=10.0009620502 nrkept=2586 nrkeptmult=2586 
Kept: 2586 out of 10768, ratio=0.24
Memory used: 216 MB  Time elapsed: 263 s

 ***** [DMNRG] Iteration 21/66 (scale 0.00083772617549037929) ***** 
[1] xi(20)=0.999999602957 zeta(21)=0
Stats: nr=62 min=1 max=782 sum3=2748041580
Safeguard: keep additional 1 states
Emax=10.003959981 nrkept=2228 nrkeptmult=2228 
Kept: 2228 out of 10344, ratio=0.215
Memory used: 212 MB  Time elapsed: 264 s

 ***** [DMNRG] Iteration 22/66 (scale 0.00059236185946671895) ***** 
[1] xi(21)=0.999999801478 zeta(22)=0
Stats: nr=71 min=1 max=743 sum3=1723641866
Safeguard: keep additional 1 states
Emax=10.0025791 nrkept=2103 nrkeptmult=2103 
Kept: 2103 out of 8912, ratio=0.236
Memory used: 211 MB  Time elapsed: 265 s

 ***** [DMNRG] Iteration 23/66 (scale 0.00041886308774518965) ***** 
[1] xi(22)=0.999999900739 zeta(23)=0
Stats: nr=62 min=1 max=654 sum3=1507450476
Safeguard: keep additional 1 states
Emax=10.0012705678 nrkept=1989 nrkeptmult=1989 
Kept: 1989 out of 8412, ratio=0.236
Memory used: 208 MB  Time elapsed: 266 s

 ***** [DMNRG] Iteration 24/66 (scale 0.00029618092973335947) ***** 
[1] xi(23)=0.99999995037 zeta(24)=0
Stats: nr=69 min=1 max=672 sum3=1257195180
Safeguard: keep additional 1 states
Emax=10.0145373054 nrkept=2013 nrkeptmult=2013 
Kept: 2013 out of 7956, ratio=0.253
Memory used: 210 MB  Time elapsed: 267 s

 ***** [DMNRG] Iteration 25/66 (scale 0.00020943154387259482) ***** 
[1] xi(24)=0.999999975185 zeta(25)=0
Stats: nr=62 min=1 max=622 sum3=1320642150
Safeguard: keep additional 1 states
Emax=10.0181329769 nrkept=1766 nrkeptmult=1766 
Kept: 1766 out of 8052, ratio=0.219
Memory used: 210 MB  Time elapsed: 268 s

 ***** [DMNRG] Iteration 26/66 (scale 0.00014809046486667974) ***** 
[1] xi(25)=0.999999987592 zeta(26)=0
Stats: nr=69 min=1 max=599 sum3=895388342
Safeguard: keep additional 1 states
Emax=10.0028225723 nrkept=1746 nrkeptmult=1746 
Kept: 1746 out of 7064, ratio=0.247
Memory used: 210 MB  Time elapsed: 268 s

 ***** [DMNRG] Iteration 27/66 (scale 0.00010471577193629741) ***** 
[1] xi(26)=0.999999993796 zeta(27)=0
Stats: nr=61 min=1 max=566 sum3=902205654
Safeguard: keep additional 3 states
Emax=10.0015470658 nrkept=1536 nrkeptmult=1536 
Kept: 1536 out of 6984, ratio=0.22
Memory used: 212 MB  Time elapsed: 269 s

 ***** [DMNRG] Iteration 28/66 (scale 7.4045232433339868e-05) ***** 
[1] xi(27)=0.999999996898 zeta(28)=0
Stats: nr=65 min=1 max=526 sum3=612952602
Safeguard: keep additional 1 states
Emax=10.0045042769 nrkept=1475 nrkeptmult=1475 
Kept: 1475 out of 6144, ratio=0.24
Memory used: 212 MB  Time elapsed: 270 s

 ***** [DMNRG] Iteration 29/66 (scale 5.2357885968148706e-05) ***** 
[1] xi(28)=0.999999998449 zeta(29)=0
Stats: nr=59 min=1 max=500 sum3=567156254
Safeguard: keep additional 2 states
Emax=10.0016466226 nrkept=1302 nrkeptmult=1302 
Kept: 1302 out of 5900, ratio=0.221
Memory used: 216 MB  Time elapsed: 270 s

 ***** [DMNRG] Iteration 30/66 (scale 3.7022616216669934e-05) ***** 
[1] xi(29)=0.999999999225 zeta(30)=0
Stats: nr=59 min=1 max=449 sum3=396478788
Safeguard: keep additional 1 states
Emax=10.0013765028 nrkept=1255 nrkeptmult=1255 
Kept: 1255 out of 5208, ratio=0.241
Memory used: 215 MB  Time elapsed: 270 s

 ***** [DMNRG] Iteration 31/66 (scale 2.6178942984074353e-05) ***** 
[1] xi(30)=0.999999999612 zeta(31)=0
Stats: nr=57 min=1 max=449 sum3=372148378
Emax=10.004433867 nrkept=1122 nrkeptmult=1122 
Kept: 1122 out of 5020, ratio=0.224
Memory used: 221 MB  Time elapsed: 271 s

 ***** [DMNRG] Iteration 32/66 (scale 1.8511308108334967e-05) ***** 
[1] xi(31)=0.999999999806 zeta(32)=0
Stats: nr=58 min=1 max=381 sum3=260161692
Safeguard: keep additional 3 states
Emax=10.0030030674 nrkept=1141 nrkeptmult=1141 
Kept: 1141 out of 4488, ratio=0.254
Memory used: 222 MB  Time elapsed: 271 s

 ***** [DMNRG] Iteration 33/66 (scale 1.3089471492037176e-05) ***** 
[1] xi(32)=0.999999999903 zeta(33)=0
Stats: nr=57 min=1 max=418 sum3=287422162
Safeguard: keep additional 3 states
Emax=10.0020460748 nrkept=1083 nrkeptmult=1083 
Kept: 1083 out of 4564, ratio=0.237
Memory used: 223 MB  Time elapsed: 271 s

 ***** [DMNRG] Iteration 34/66 (scale 9.2556540541674836e-06) ***** 
[1] xi(33)=0.999999999952 zeta(34)=0
Stats: nr=58 min=1 max=371 sum3=237920658
Safeguard: keep additional 1 states
Emax=10.018649788 nrkept=1116 nrkeptmult=1116 
Kept: 1116 out of 4332, ratio=0.258
Memory used: 222 MB  Time elapsed: 271 s

 ***** [DMNRG] Iteration 35/66 (scale 6.5447357460185882e-06) ***** 
[1] xi(34)=0.999999999976 zeta(35)=0
Stats: nr=57 min=1 max=412 sum3=271853526
Safeguard: keep additional 1 states
Emax=10.0252201681 nrkept=1075 nrkeptmult=1075 
Kept: 1075 out of 4464, ratio=0.241
Memory used: 223 MB  Time elapsed: 272 s

 ***** [DMNRG] Iteration 36/66 (scale 4.6278270270837418e-06) ***** 
[1] xi(35)=0.999999999988 zeta(36)=0
Stats: nr=58 min=1 max=371 sum3=234091648
Safeguard: keep additional 2 states
Emax=10.0535407837 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4300, ratio=0.26
Memory used: 222 MB  Time elapsed: 272 s

 ***** [DMNRG] Iteration 37/66 (scale 3.2723678730092941e-06) ***** 
[1] xi(36)=0.999999999994 zeta(37)=0
Stats: nr=57 min=1 max=413 sum3=273706324
Safeguard: keep additional 1 states
Emax=10.0005691486 nrkept=1069 nrkeptmult=1069 
Kept: 1069 out of 4468, ratio=0.239
Memory used: 222 MB  Time elapsed: 272 s

 ***** [DMNRG] Iteration 38/66 (scale 2.3139135135418709e-06) ***** 
[1] xi(37)=0.999999999997 zeta(38)=0
Stats: nr=58 min=1 max=371 sum3=232214488
Emax=10.0679012837 nrkept=1115 nrkeptmult=1115 
Kept: 1115 out of 4276, ratio=0.261
Memory used: 222 MB  Time elapsed: 273 s

 ***** [DMNRG] Iteration 39/66 (scale 1.6361839365046471e-06) ***** 
[1] xi(38)=0.999999999998 zeta(39)=0
Stats: nr=57 min=1 max=413 sum3=272139644
Emax=10.0023974172 nrkept=1060 nrkeptmult=1060 
Kept: 1060 out of 4460, ratio=0.238
Memory used: 222 MB  Time elapsed: 273 s

 ***** [DMNRG] Iteration 40/66 (scale 1.1569567567709354e-06) ***** 
[1] xi(39)=0.999999999999 zeta(40)=0
Stats: nr=58 min=1 max=371 sum3=231642436
Safeguard: keep additional 2 states
Emax=10.0725678988 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4240, ratio=0.263
Memory used: 223 MB  Time elapsed: 273 s

 ***** [DMNRG] Iteration 41/66 (scale 8.1809196825232353e-07) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Emax=10.0128597108 nrkept=1060 nrkeptmult=1060 
Kept: 1060 out of 4468, ratio=0.237
Memory used: 223 MB  Time elapsed: 273 s

 ***** [DMNRG] Iteration 42/66 (scale 5.7847837838546772e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=58 min=1 max=371 sum3=231642436
Safeguard: keep additional 2 states
Emax=10.0743813352 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4240, ratio=0.263
Memory used: 223 MB  Time elapsed: 274 s

 ***** [DMNRG] Iteration 43/66 (scale 4.0904598412616176e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0165961783 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 223 MB  Time elapsed: 274 s

 ***** [DMNRG] Iteration 44/66 (scale 2.8923918919273386e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.075288325 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 223 MB  Time elapsed: 274 s

 ***** [DMNRG] Iteration 45/66 (scale 2.0452299206308088e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.017865681 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 223 MB  Time elapsed: 274 s

 ***** [DMNRG] Iteration 46/66 (scale 1.4461959459636693e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0757419091 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 223 MB  Time elapsed: 275 s

 ***** [DMNRG] Iteration 47/66 (scale 1.0226149603154044e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0184767017 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 225 MB  Time elapsed: 275 s

 ***** [DMNRG] Iteration 48/66 (scale 7.2309797298183465e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0759687261 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 225 MB  Time elapsed: 275 s

 ***** [DMNRG] Iteration 49/66 (scale 5.113074801577022e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0187806375 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 224 MB  Time elapsed: 275 s

 ***** [DMNRG] Iteration 50/66 (scale 3.6154898649091733e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0760821412 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 225 MB  Time elapsed: 276 s

 ***** [DMNRG] Iteration 51/66 (scale 2.556537400788511e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0189324201 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 224 MB  Time elapsed: 276 s

 ***** [DMNRG] Iteration 52/66 (scale 1.8077449324545866e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761388504 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 224 MB  Time elapsed: 276 s

 ***** [DMNRG] Iteration 53/66 (scale 1.2782687003942555e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190082825 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 225 MB  Time elapsed: 276 s

 ***** [DMNRG] Iteration 54/66 (scale 9.0387246622729332e-09) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761672054 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 226 MB  Time elapsed: 277 s

 ***** [DMNRG] Iteration 55/66 (scale 6.3913435019712775e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190462082 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 226 MB  Time elapsed: 277 s

 ***** [DMNRG] Iteration 56/66 (scale 4.5193623311364666e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.076181383 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 227 MB  Time elapsed: 277 s

 ***** [DMNRG] Iteration 57/66 (scale 3.1956717509856388e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190651699 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 226 MB  Time elapsed: 277 s

 ***** [DMNRG] Iteration 58/66 (scale 2.2596811655682333e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761884718 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 226 MB  Time elapsed: 278 s

 ***** [DMNRG] Iteration 59/66 (scale 1.5978358754928194e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190746505 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 227 MB  Time elapsed: 278 s

 ***** [DMNRG] Iteration 60/66 (scale 1.1298405827841166e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761920162 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 227 MB  Time elapsed: 278 s

 ***** [DMNRG] Iteration 61/66 (scale 7.9891793774640969e-10) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190793906 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 226 MB  Time elapsed: 279 s

 ***** [DMNRG] Iteration 62/66 (scale 5.6492029139205832e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761937884 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 226 MB  Time elapsed: 279 s

 ***** [DMNRG] Iteration 63/66 (scale 3.9945896887320485e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190817607 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 226 MB  Time elapsed: 279 s

 ***** [DMNRG] Iteration 64/66 (scale 2.8246014569602916e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761946744 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 228 MB  Time elapsed: 279 s

 ***** [DMNRG] Iteration 65/66 (scale 1.9972948443660242e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190829456 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 229 MB  Time elapsed: 280 s

 ***** [DMNRG] Iteration 66/66 (scale 1.4123007284801458e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761951174 nrkept=4244 nrkeptmult=4244 
Kept: 4244 out of 4244, ratio=1
Memory used: 247 MB  Time elapsed: 280 s

** Iteration completed.

Spectrum: A_d-A_d-d DMNRG
[A_d-A_d-d] pos=0.135159587535 neg=0.865013591765 sum= 1.0001731793
mu1=-0.00379526552288 mu2=0.000662543377272 mu3=-5.05338509112e-06 mu4=0.000220253407427

Spectrum: A_d-A_d-d CFS
[A_d-A_d-d] pos=0.135344405147 neg=0.864655594852 sum= 0.999999999998
mu1=-0.00382987307473 mu2=0.000665958031959 mu3=-5.12504589512e-06 mu4=0.000221622845747

Spectrum: A_d-A_d-d FDM
[A_d-A_d-d] pos=0.135344405147 neg=0.864655594852 sum= 0.999999999998
mu1=-0.00382987307473 mu2=0.000665958031959 mu3=-5.12504589512e-06 mu4=0.000221622845747

Spectrum: A_d-A_d-u DMNRG
[A_d-A_d-u] pos=0.865013591765 neg=0.135159587535 sum= 1.0001731793
mu1=0.00379526552288 mu2=0.000662543377272 mu3=5.05338509111e-06 mu4=0.000220253407427

Spectrum: A_d-A_d-u CFS
[A_d-A_d-u] pos=0.864655594852 neg=0.135344405147 sum= 0.999999999998
mu1=0.00382987307473 mu2=0.000665958031959 mu3=5.12504589511e-06 mu4=0.000221622845747

Spectrum: A_d-A_d-u FDM
[A_d-A_d-u] pos=0.864655594852 neg=0.135344405147 sum= 0.999999999998
mu1=0.00382987307473 mu2=0.000665958031959 mu3=5.12504589511e-06 mu4=0.000221622845747

Master done. Terminating slave processes.
Master exiting.

Memory usage report (process 0)
===================
  after diag: 427828 kB
after recalc: 525032 kB
 after split: 488668 kB
  after trim: 263368 kB
 after trunc: 428160 kB
      diag_h: 349120 kB

Peak usage: 525032 kB

Timing report (process 0)
=============
         All: 281 s
          DM: 4 s
         FDM: 18 s
        diag: 165 s
    recalc d: 24 s
      spec d: 43 s
       Other: 24 s
