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=43
TEMPERATURE T=1e-10
Nmax=43
Tmin=1e-10
Ninit=0
Lambda=3
strategy=kept
lastall=1
lastalloverride=0
z=0.5
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=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 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=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: 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.0001
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	3.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.5
>>  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=3 z=0.5
# keep=10000 keepenergy=10 Nmax=43 channels=1 dots=1 betabar=1

Number of invariant subspaces: 9
(-2 1) 1 states
0.0291267203408 
(-1 0) 2 states
0.013810772019 0.0416686726337 
(-1 2) 2 states
0.0138409454141 0.041693429853 
(0 -1) 1 states
0.0263527243119 
(0 1) 4 states
0 0.0263801923155 0.0291267203408 0.0555069072635 
(0 3) 1 states
0.0264076549263 
(1 0) 2 states
0.013810772019 0.0416686726337 
(1 2) 2 states
0.0138409454141 0.041693429853 
(2 1) 1 states
0.0291267203408 
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)
<SZd>=-6.8755356617e-06

 ***** [NRG] Iteration 1/43 (scale 1.0510537250399226) ***** 
[1] xi(0)=0.586403356468 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=2.35236690813 nrkept=64 nrkeptmult=64 
Total energy=-1.23375695957  Egs=-1.12813680645
Kept: 64 out of 64, ratio=1
<SZd>=-1.19200052391e-05
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/43 (scale 0.60682615108455829) ***** 
[1] xi(1)=0.927938038869 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=5.61680651034 nrkept=256 nrkeptmult=256 
Total energy=-1.70171254342  Egs=-0.771152632452
Kept: 256 out of 256, ratio=1
<SZd>=-2.06783787318e-05
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/43 (scale 0.35035124167997428) ***** 
[1] xi(2)=1.1294419779 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Emax=10.0655735178 nrkept=945 nrkeptmult=945 
Total energy=-2.31793826276  Egs=-1.75887979271
Kept: 945 out of 1024, ratio=0.923
<SZd>=-3.67837922257e-05
Memory used: 7 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/43 (scale 0.20227538369485276) ***** 
[1] xi(3)=1.06460020271 zeta(4)=0
Stats: nr=49 min=1 max=357 sum3=170022450
Safeguard: keep additional 15 states
Emax=10.443625624 nrkept=1328 nrkeptmult=1328 
Total energy=-2.4731867701  Egs=-0.767510630826
Kept: 1328 out of 3780, ratio=0.351
<SZd>=-6.24816459204e-05
Memory used: 16 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/43 (scale 0.11678374722665809) ***** 
[1] xi(4)=1.01470147291 zeta(5)=0
Stats: nr=60 min=1 max=414 sum3=414425840
Safeguard: keep additional 5 states
Emax=10.0320024748 nrkept=998 nrkeptmult=998 
Total energy=-2.65841641869  Egs=-1.58609098438
Kept: 998 out of 5312, ratio=0.188
<SZd>=-0.000108950976629
Memory used: 20 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 6/43 (scale 0.067425127898284254) ***** 
[1] xi(5)=1.00360683793 zeta(6)=0
Stats: nr=49 min=1 max=378 sum3=202464302
Safeguard: keep additional 15 states
Emax=10.3050891487 nrkept=976 nrkeptmult=976 
Total energy=-2.70992699356  Egs=-0.763967032397
Kept: 976 out of 3992, ratio=0.244
<SZd>=-0.000190971043091
Memory used: 22 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 7/43 (scale 0.038927915742219356) ***** 
[1] xi(6)=1.00109067805 zeta(7)=0
Stats: nr=60 min=1 max=308 sum3=168775504
Safeguard: keep additional 1 states
Emax=10.6373111387 nrkept=866 nrkeptmult=866 
Total energy=-2.76859981596  Egs=-1.50721715463
Kept: 866 out of 3904, ratio=0.222
<SZd>=-0.000338377169056
Memory used: 23 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 8/43 (scale 0.022475042632761417) ***** 
[1] xi(7)=1.00035371682 zeta(8)=0
Stats: nr=49 min=1 max=332 sum3=134201072
Emax=10.3508170149 nrkept=961 nrkeptmult=961 
Total energy=-2.78695644832  Egs=-0.816756286225
Kept: 961 out of 3464, ratio=0.277
<SZd>=-0.000602689224295
Memory used: 24 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 9/43 (scale 0.01297597191407312) ***** 
[1] xi(8)=1.00011689308 zeta(9)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Safeguard: keep additional 1 states
Emax=10.1134176276 nrkept=802 nrkeptmult=802 
Total energy=-2.80544969838  Egs=-1.42519189994
Kept: 802 out of 3844, ratio=0.209
<SZd>=-0.00109570171926
Memory used: 26 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 10/43 (scale 0.0074916808775871389) ***** 
[1] xi(9)=1.00003885431 zeta(10)=0
Stats: nr=49 min=1 max=310 sum3=108330934
Emax=10.4857101621 nrkept=961 nrkeptmult=961 
Total energy=-2.81196857848  Egs=-0.870149196159
Kept: 961 out of 3208, ratio=0.3
<SZd>=-0.00203045441894
Memory used: 26 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 11/43 (scale 0.0043253239713577068) ***** 
[1] xi(10)=1.00001293929 zeta(11)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Emax=10.0919916224 nrkept=737 nrkeptmult=737 
Total energy=-2.81779533643  Egs=-1.34712636453
Kept: 737 out of 3844, ratio=0.192
<SZd>=-0.00381138314748
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 12/43 (scale 0.0024972269591957133) ***** 
[1] xi(11)=1.00000431175 zeta(12)=0
Stats: nr=49 min=1 max=280 sum3=81545942
Safeguard: keep additional 1 states
Emax=10.0893126107 nrkept=766 nrkeptmult=766 
Total energy=-2.82008311619  Egs=-0.916128086591
Kept: 766 out of 2948, ratio=0.26
<SZd>=-0.00707274342215
Memory used: 28 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 13/43 (scale 0.0014417746571192356) ***** 
[1] xi(12)=1.0000014371 zeta(13)=0
Stats: nr=56 min=1 max=250 sum3=87443998
Emax=10.0862504692 nrkept=577 nrkeptmult=577 
Total energy=-2.82192351942  Egs=-1.27648466043
Kept: 577 out of 3064, ratio=0.188
<SZd>=-0.012613925498
Memory used: 29 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 14/43 (scale 0.00083240898639857098) ***** 
[1] xi(13)=1.00000047902 zeta(14)=0
Stats: nr=47 min=1 max=228 sum3=42177400
Safeguard: keep additional 2 states
Emax=10.0233570292 nrkept=563 nrkeptmult=563 
Total energy=-2.82272334192  Egs=-0.960852783121
Kept: 563 out of 2308, ratio=0.244
<SZd>=-0.0214820243594
Memory used: 29 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 15/43 (scale 0.00048059155237307853) ***** 
[1] xi(14)=1.00000015967 zeta(15)=0
Stats: nr=56 min=1 max=188 sum3=35828120
Safeguard: keep additional 2 states
Emax=10.0060155992 nrkept=428 nrkeptmult=428 
Total energy=-2.82330626682  Egs=-1.21293206057
Kept: 428 out of 2252, ratio=0.19
<SZd>=-0.0354383943986
Memory used: 29 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 16/43 (scale 0.00027746966213285698) ***** 
[1] xi(15)=1.00000005322 zeta(16)=0
Stats: nr=47 min=1 max=171 sum3=17696450
Safeguard: keep additional 2 states
Emax=10.0103506131 nrkept=511 nrkeptmult=511 
Total energy=-2.82358447678  Egs=-1.00266805418
Kept: 511 out of 1712, ratio=0.298
<SZd>=-0.0570596340065
Memory used: 29 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 17/43 (scale 0.00016019718412435953) ***** 
[1] xi(16)=1.00000001774 zeta(17)=0
Stats: nr=54 min=1 max=169 sum3=25737772
Emax=10.0656156449 nrkept=404 nrkeptmult=404 
Total energy=-2.82377156941  Egs=-1.16788962724
Kept: 404 out of 2044, ratio=0.198
<SZd>=-0.0890460055273
Memory used: 29 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 18/43 (scale 9.2489887377619005e-05) ***** 
[1] xi(17)=1.00000000591 zeta(18)=0
Stats: nr=47 min=1 max=161 sum3=14891834
Safeguard: keep additional 1 states
Emax=10.0087667452 nrkept=415 nrkeptmult=415 
Total energy=-2.82386691582  Egs=-1.03088469171
Kept: 415 out of 1616, ratio=0.257
<SZd>=-0.132709522164
Memory used: 29 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 19/43 (scale 5.33990613747865e-05) ***** 
[1] xi(18)=1.00000000197 zeta(19)=0
Stats: nr=50 min=1 max=146 sum3=15268786
Safeguard: keep additional 1 states
Emax=10.002203057 nrkept=369 nrkeptmult=369 
Total energy=-2.82392843095  Egs=-1.15198894513
Kept: 369 out of 1660, ratio=0.222
<SZd>=-0.183355817136
Memory used: 29 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 20/43 (scale 3.0829962459206333e-05) ***** 
[1] xi(19)=1.00000000066 zeta(20)=0
Stats: nr=46 min=1 max=147 sum3=11253312
Safeguard: keep additional 3 states
Emax=10.0103009088 nrkept=345 nrkeptmult=345 
Total energy=-2.82396057816  Egs=-1.04272621685
Kept: 345 out of 1476, ratio=0.234
<SZd>=-0.22942349773
Memory used: 29 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 21/43 (scale 1.7799687124928835e-05) ***** 
[1] xi(20)=1.00000000022 zeta(21)=0
Stats: nr=44 min=1 max=132 sum3=9987912
Safeguard: keep additional 2 states
Emax=10.0181154563 nrkept=299 nrkeptmult=299 
Total energy=-2.82398111213  Egs=-1.15361423444
Kept: 299 out of 1380, ratio=0.217
<SZd>=-0.258086148397
Memory used: 29 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 22/43 (scale 1.0276654153068778e-05) ***** 
[1] xi(21)=1.00000000007 zeta(22)=0
Stats: nr=44 min=1 max=123 sum3=6487544
Safeguard: keep additional 2 states
Emax=10.0136193936 nrkept=271 nrkeptmult=271 
Total energy=-2.82399185632  Egs=-1.04549486394
Kept: 271 out of 1196, ratio=0.227
<SZd>=-0.27043662192
Memory used: 29 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 23/43 (scale 5.9332290416429446e-06) ***** 
[1] xi(22)=1.00000000002 zeta(23)=0
Stats: nr=42 min=1 max=114 sum3=5302348
Emax=10.0238366129 nrkept=246 nrkeptmult=246 
Total energy=-2.82399872197  Egs=-1.15715183328
Kept: 246 out of 1084, ratio=0.227
<SZd>=-0.274350511543
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 24/43 (scale 3.4255513843562592e-06) ***** 
[1] xi(23)=1.00000000001 zeta(24)=0
Stats: nr=42 min=1 max=102 sum3=3981612
Safeguard: keep additional 1 states
Emax=10.0024175957 nrkept=244 nrkeptmult=244 
Total energy=-2.82400230432  Egs=-1.04577274634
Kept: 244 out of 984, ratio=0.248
<SZd>=-0.27566540833
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 25/43 (scale 1.9777430138809818e-06) ***** 
[1] xi(24)=1 zeta(25)=0
Stats: nr=41 min=1 max=103 sum3=3938398
Safeguard: keep additional 1 states
Emax=10.0438195728 nrkept=232 nrkeptmult=232 
Total energy=-2.82400459632  Egs=-1.15889634788
Kept: 232 out of 976, ratio=0.238
<SZd>=-0.276034806639
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 26/43 (scale 1.1418504614520865e-06) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=40 min=1 max=93 sum3=3310666
Safeguard: keep additional 1 states
Emax=10.0173313404 nrkept=232 nrkeptmult=232 
Total energy=-2.82400579039  Egs=-1.04573118988
Kept: 232 out of 928, ratio=0.25
<SZd>=-0.27619159414
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 27/43 (scale 6.592476712936605e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0975842807 nrkept=231 nrkeptmult=231 
Total energy=-2.82400655482  Egs=-1.1595518798
Kept: 231 out of 928, ratio=0.249
<SZd>=-0.276222655134
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 28/43 (scale 3.8061682048402881e-07) ***** 
[1] xi(27)=1 zeta(28)=0
Stats: nr=40 min=1 max=93 sum3=3306120
Safeguard: keep additional 1 states
Emax=10.0740152889 nrkept=232 nrkeptmult=232 
Total energy=-2.82400695283  Egs=-1.04569862764
Kept: 232 out of 924, ratio=0.251
<SZd>=-0.276245279854
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 29/43 (scale 2.1974922376455349e-07) ***** 
[1] xi(28)=1 zeta(29)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Emax=10.0069987544 nrkept=229 nrkeptmult=229 
Total energy=-2.82400720769  Egs=-1.15977897828
Kept: 229 out of 928, ratio=0.247
<SZd>=-0.27624565734
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 30/43 (scale 1.2687227349467628e-07) ***** 
[1] xi(29)=1 zeta(30)=0
Stats: nr=40 min=1 max=91 sum3=3180808
Safeguard: keep additional 1 states
Emax=10.0868791085 nrkept=232 nrkeptmult=232 
Total energy=-2.82400734036  Egs=-1.04568553863
Kept: 232 out of 916, ratio=0.253
<SZd>=-0.276249941912
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 31/43 (scale 7.3249741254851168e-08) ***** 
[1] xi(30)=1 zeta(31)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0051609801 nrkept=228 nrkeptmult=228 
Total energy=-2.82400742532  Egs=-1.15985563447
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276248961084
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 32/43 (scale 4.2290757831558757e-08) ***** 
[1] xi(31)=1 zeta(32)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0905865944 nrkept=232 nrkeptmult=232 
Total energy=-2.82400746954  Egs=-1.04568093053
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276250037097
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 33/43 (scale 2.4416580418283724e-08) ***** 
[1] xi(32)=1 zeta(33)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0076441301 nrkept=228 nrkeptmult=228 
Total energy=-2.82400749786  Egs=-1.15988128312
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.2762495886
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 34/43 (scale 1.4096919277186252e-08) ***** 
[1] xi(33)=1 zeta(34)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0917926054 nrkept=232 nrkeptmult=232 
Total energy=-2.8240075126  Egs=-1.04567937669
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276249907474
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 35/43 (scale 8.1388601394279074e-09) ***** 
[1] xi(34)=1 zeta(35)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0084764394 nrkept=228 nrkeptmult=228 
Total energy=-2.82400752204  Egs=-1.15988983321
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276249744578
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 36/43 (scale 4.6989730923954177e-09) ***** 
[1] xi(35)=1 zeta(36)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0921919313 nrkept=232 nrkeptmult=232 
Total energy=-2.82400752695  Egs=-1.04567886643
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276249846451
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 37/43 (scale 2.7129533798093026e-09) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0087587617 nrkept=228 nrkeptmult=228 
Total energy=-2.8240075301  Egs=-1.15989267309
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276249790685
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 38/43 (scale 1.5663243641318058e-09) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 3 states
Emax=10.0923309525 nrkept=234 nrkeptmult=234 
Total energy=-2.82400753174  Egs=-1.04567870685
Kept: 234 out of 912, ratio=0.257
<SZd>=-0.276249824135
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 39/43 (scale 9.0431779326976753e-10) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=41 min=1 max=97 sum3=3398916
Safeguard: keep additional 1 states
Emax=10.0087662125 nrkept=228 nrkeptmult=228 
Total energy=-2.82400753279  Egs=-1.15989360842
Kept: 228 out of 936, ratio=0.244
<SZd>=-0.276249805401
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 40/43 (scale 5.2210812137726861e-10) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923703152 nrkept=238 nrkeptmult=238 
Total energy=-2.82400753333  Egs=-1.04567866443
Kept: 238 out of 912, ratio=0.261
<SZd>=-0.276249816485
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 41/43 (scale 3.0143926442325579e-10) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=41 min=1 max=99 sum3=3556156
Safeguard: keep additional 1 states
Emax=10.0087970965 nrkept=228 nrkeptmult=228 
Total energy=-2.82400753368  Egs=-1.15989390876
Kept: 228 out of 952, ratio=0.239
<SZd>=-0.276249810234
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 42/43 (scale 1.7403604045908955e-10) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923839835 nrkept=238 nrkeptmult=238 
Total energy=-2.82400753386  Egs=-1.04567866113
Kept: 238 out of 912, ratio=0.261
<SZd>=-0.27624981392
Memory used: 30 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 43/43 (scale 1.0047975480775194e-10) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=41 min=1 max=99 sum3=3556156
Safeguard: keep additional 1 states
Emax=10.008807421 nrkept=952 nrkeptmult=952 
Total energy=-2.82400753398  Egs=-1.15989399742
Kept: 952 out of 952, ratio=1
<SZd>=-0.276249811726
Memory used: 30 MB  Time elapsed: 8 s

** Iteration completed.


Total energy:             -2.824007533980854
TKW=0.000167991463601
Grand canonical Z=              5.779279112961445
Zft (last shell)=               5.779278469315136
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.0001
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	3.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.5
>>  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.0291267203408 
(-1 0) 2 states
0.013810772019 0.0416686726337 
(-1 2) 2 states
0.0138409454141 0.041693429853 
(0 -1) 1 states
0.0263527243119 
(0 1) 4 states
0 0.0263801923155 0.0291267203408 0.0555069072635 
(0 3) 1 states
0.0264076549263 
(1 0) 2 states
0.013810772019 0.0416686726337 
(1 2) 2 states
0.0138409454141 0.041693429853 
(2 1) 1 states
0.0291267203408 
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 FDM
Spectrum: A_d-A_d-u FDM

 ***** [DMNRG] Iteration 1/43 (scale 1.0510537250399226) ***** 
[1] xi(0)=0.586403356468 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=2.35236690813 nrkept=64 nrkeptmult=64 
Kept: 64 out of 64, ratio=1
Memory used: 36 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 2/43 (scale 0.60682615108455829) ***** 
[1] xi(1)=0.927938038869 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=5.61680651034 nrkept=256 nrkeptmult=256 
Kept: 256 out of 256, ratio=1
Memory used: 36 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 3/43 (scale 0.35035124167997428) ***** 
[1] xi(2)=1.1294419779 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Emax=10.0655735178 nrkept=945 nrkeptmult=945 
Kept: 945 out of 1024, ratio=0.923
Memory used: 34 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 4/43 (scale 0.20227538369485276) ***** 
[1] xi(3)=1.06460020271 zeta(4)=0
Stats: nr=49 min=1 max=357 sum3=170022450
Safeguard: keep additional 15 states
Emax=10.443625624 nrkept=1328 nrkeptmult=1328 
Kept: 1328 out of 3780, ratio=0.351
Memory used: 35 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 5/43 (scale 0.11678374722665809) ***** 
[1] xi(4)=1.01470147291 zeta(5)=0
Stats: nr=60 min=1 max=414 sum3=414425840
Safeguard: keep additional 5 states
Emax=10.0320024748 nrkept=998 nrkeptmult=998 
Kept: 998 out of 5312, ratio=0.188
Memory used: 36 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 6/43 (scale 0.067425127898284254) ***** 
[1] xi(5)=1.00360683793 zeta(6)=0
Stats: nr=49 min=1 max=378 sum3=202464302
Safeguard: keep additional 15 states
Emax=10.3050891487 nrkept=976 nrkeptmult=976 
Kept: 976 out of 3992, ratio=0.244
Memory used: 37 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 7/43 (scale 0.038927915742219356) ***** 
[1] xi(6)=1.00109067805 zeta(7)=0
Stats: nr=60 min=1 max=308 sum3=168775504
Safeguard: keep additional 1 states
Emax=10.6373111387 nrkept=866 nrkeptmult=866 
Kept: 866 out of 3904, ratio=0.222
Memory used: 36 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 8/43 (scale 0.022475042632761417) ***** 
[1] xi(7)=1.00035371682 zeta(8)=0
Stats: nr=49 min=1 max=332 sum3=134201072
Emax=10.3508170149 nrkept=961 nrkeptmult=961 
Kept: 961 out of 3464, ratio=0.277
Memory used: 37 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 9/43 (scale 0.01297597191407312) ***** 
[1] xi(8)=1.00011689308 zeta(9)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Safeguard: keep additional 1 states
Emax=10.1134176276 nrkept=802 nrkeptmult=802 
Kept: 802 out of 3844, ratio=0.209
Memory used: 37 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 10/43 (scale 0.0074916808775871389) ***** 
[1] xi(9)=1.00003885431 zeta(10)=0
Stats: nr=49 min=1 max=310 sum3=108330934
Emax=10.4857101621 nrkept=961 nrkeptmult=961 
Kept: 961 out of 3208, ratio=0.3
Memory used: 36 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 11/43 (scale 0.0043253239713577068) ***** 
[1] xi(10)=1.00001293929 zeta(11)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Emax=10.0919916224 nrkept=737 nrkeptmult=737 
Kept: 737 out of 3844, ratio=0.192
Memory used: 37 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 12/43 (scale 0.0024972269591957133) ***** 
[1] xi(11)=1.00000431175 zeta(12)=0
Stats: nr=49 min=1 max=280 sum3=81545942
Safeguard: keep additional 1 states
Emax=10.0893126107 nrkept=766 nrkeptmult=766 
Kept: 766 out of 2948, ratio=0.26
Memory used: 35 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 13/43 (scale 0.0014417746571192356) ***** 
[1] xi(12)=1.0000014371 zeta(13)=0
Stats: nr=56 min=1 max=250 sum3=87443998
Emax=10.0862504692 nrkept=577 nrkeptmult=577 
Kept: 577 out of 3064, ratio=0.188
Memory used: 37 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 14/43 (scale 0.00083240898639857098) ***** 
[1] xi(13)=1.00000047902 zeta(14)=0
Stats: nr=47 min=1 max=228 sum3=42177400
Safeguard: keep additional 2 states
Emax=10.0233570292 nrkept=563 nrkeptmult=563 
Kept: 563 out of 2308, ratio=0.244
Memory used: 35 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 15/43 (scale 0.00048059155237307853) ***** 
[1] xi(14)=1.00000015967 zeta(15)=0
Stats: nr=56 min=1 max=188 sum3=35828120
Safeguard: keep additional 2 states
Emax=10.0060155992 nrkept=428 nrkeptmult=428 
Kept: 428 out of 2252, ratio=0.19
Memory used: 35 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 16/43 (scale 0.00027746966213285698) ***** 
[1] xi(15)=1.00000005322 zeta(16)=0
Stats: nr=47 min=1 max=171 sum3=17696450
Safeguard: keep additional 2 states
Emax=10.0103506131 nrkept=511 nrkeptmult=511 
Kept: 511 out of 1712, ratio=0.298
Memory used: 35 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 17/43 (scale 0.00016019718412435953) ***** 
[1] xi(16)=1.00000001774 zeta(17)=0
Stats: nr=54 min=1 max=169 sum3=25737772
Emax=10.0656156449 nrkept=404 nrkeptmult=404 
Kept: 404 out of 2044, ratio=0.198
Memory used: 35 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 18/43 (scale 9.2489887377619005e-05) ***** 
[1] xi(17)=1.00000000591 zeta(18)=0
Stats: nr=47 min=1 max=161 sum3=14891834
Safeguard: keep additional 1 states
Emax=10.0087667452 nrkept=415 nrkeptmult=415 
Kept: 415 out of 1616, ratio=0.257
Memory used: 35 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 19/43 (scale 5.33990613747865e-05) ***** 
[1] xi(18)=1.00000000197 zeta(19)=0
Stats: nr=50 min=1 max=146 sum3=15268786
Safeguard: keep additional 1 states
Emax=10.002203057 nrkept=369 nrkeptmult=369 
Kept: 369 out of 1660, ratio=0.222
Memory used: 35 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 20/43 (scale 3.0829962459206333e-05) ***** 
[1] xi(19)=1.00000000066 zeta(20)=0
Stats: nr=46 min=1 max=147 sum3=11253312
Safeguard: keep additional 3 states
Emax=10.0103009088 nrkept=345 nrkeptmult=345 
Kept: 345 out of 1476, ratio=0.234
Memory used: 35 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 21/43 (scale 1.7799687124928835e-05) ***** 
[1] xi(20)=1.00000000022 zeta(21)=0
Stats: nr=44 min=1 max=132 sum3=9987912
Safeguard: keep additional 2 states
Emax=10.0181154563 nrkept=299 nrkeptmult=299 
Kept: 299 out of 1380, ratio=0.217
Memory used: 35 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 22/43 (scale 1.0276654153068778e-05) ***** 
[1] xi(21)=1.00000000007 zeta(22)=0
Stats: nr=44 min=1 max=123 sum3=6487544
Safeguard: keep additional 2 states
Emax=10.0136193936 nrkept=271 nrkeptmult=271 
Kept: 271 out of 1196, ratio=0.227
Memory used: 35 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 23/43 (scale 5.9332290416429446e-06) ***** 
[1] xi(22)=1.00000000002 zeta(23)=0
Stats: nr=42 min=1 max=114 sum3=5302348
Emax=10.0238366129 nrkept=246 nrkeptmult=246 
Kept: 246 out of 1084, ratio=0.227
Memory used: 35 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 24/43 (scale 3.4255513843562592e-06) ***** 
[1] xi(23)=1.00000000001 zeta(24)=0
Stats: nr=42 min=1 max=102 sum3=3981612
Safeguard: keep additional 1 states
Emax=10.0024175957 nrkept=244 nrkeptmult=244 
Kept: 244 out of 984, ratio=0.248
Memory used: 35 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 25/43 (scale 1.9777430138809818e-06) ***** 
[1] xi(24)=1 zeta(25)=0
Stats: nr=41 min=1 max=103 sum3=3938398
Safeguard: keep additional 1 states
Emax=10.0438195728 nrkept=232 nrkeptmult=232 
Kept: 232 out of 976, ratio=0.238
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 26/43 (scale 1.1418504614520865e-06) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=40 min=1 max=93 sum3=3310666
Safeguard: keep additional 1 states
Emax=10.0173313404 nrkept=232 nrkeptmult=232 
Kept: 232 out of 928, ratio=0.25
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 27/43 (scale 6.592476712936605e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0975842807 nrkept=231 nrkeptmult=231 
Kept: 231 out of 928, ratio=0.249
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 28/43 (scale 3.8061682048402881e-07) ***** 
[1] xi(27)=1 zeta(28)=0
Stats: nr=40 min=1 max=93 sum3=3306120
Safeguard: keep additional 1 states
Emax=10.0740152889 nrkept=232 nrkeptmult=232 
Kept: 232 out of 924, ratio=0.251
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 29/43 (scale 2.1974922376455349e-07) ***** 
[1] xi(28)=1 zeta(29)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Emax=10.0069987544 nrkept=229 nrkeptmult=229 
Kept: 229 out of 928, ratio=0.247
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 30/43 (scale 1.2687227349467628e-07) ***** 
[1] xi(29)=1 zeta(30)=0
Stats: nr=40 min=1 max=91 sum3=3180808
Safeguard: keep additional 1 states
Emax=10.0868791085 nrkept=232 nrkeptmult=232 
Kept: 232 out of 916, ratio=0.253
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 31/43 (scale 7.3249741254851168e-08) ***** 
[1] xi(30)=1 zeta(31)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0051609801 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 32/43 (scale 4.2290757831558757e-08) ***** 
[1] xi(31)=1 zeta(32)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0905865944 nrkept=232 nrkeptmult=232 
Kept: 232 out of 912, ratio=0.254
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 33/43 (scale 2.4416580418283724e-08) ***** 
[1] xi(32)=1 zeta(33)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0076441301 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 34/43 (scale 1.4096919277186252e-08) ***** 
[1] xi(33)=1 zeta(34)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0917926054 nrkept=232 nrkeptmult=232 
Kept: 232 out of 912, ratio=0.254
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 35/43 (scale 8.1388601394279074e-09) ***** 
[1] xi(34)=1 zeta(35)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0084764394 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 36/43 (scale 4.6989730923954177e-09) ***** 
[1] xi(35)=1 zeta(36)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0921919313 nrkept=232 nrkeptmult=232 
Kept: 232 out of 912, ratio=0.254
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 37/43 (scale 2.7129533798093026e-09) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0087587617 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 38/43 (scale 1.5663243641318058e-09) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 3 states
Emax=10.0923309525 nrkept=234 nrkeptmult=234 
Kept: 234 out of 912, ratio=0.257
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 39/43 (scale 9.0431779326976753e-10) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=41 min=1 max=97 sum3=3398916
Safeguard: keep additional 1 states
Emax=10.0087662125 nrkept=228 nrkeptmult=228 
Kept: 228 out of 936, ratio=0.244
Memory used: 35 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 40/43 (scale 5.2210812137726861e-10) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923703152 nrkept=238 nrkeptmult=238 
Kept: 238 out of 912, ratio=0.261
Memory used: 35 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 41/43 (scale 3.0143926442325579e-10) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=41 min=1 max=99 sum3=3556156
Safeguard: keep additional 1 states
Emax=10.0087970965 nrkept=228 nrkeptmult=228 
Kept: 228 out of 952, ratio=0.239
Memory used: 35 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 42/43 (scale 1.7403604045908955e-10) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923839835 nrkept=238 nrkeptmult=238 
Kept: 238 out of 912, ratio=0.261
Memory used: 35 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 43/43 (scale 1.0047975480775194e-10) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=41 min=1 max=99 sum3=3556156
Safeguard: keep additional 1 states
Emax=10.008807421 nrkept=952 nrkeptmult=952 
Kept: 952 out of 952, ratio=1
Memory used: 35 MB  Time elapsed: 16 s

** Iteration completed.

Spectrum: A_d-A_d-d FDM
[A_d-A_d-d] pos=0.223750183427 neg=0.776249816573 sum= 1
mu1=-0.00282282714984 mu2=0.000666637965285 mu3=-3.85005662435e-06 mu4=0.000243645610303

Spectrum: A_d-A_d-u FDM
[A_d-A_d-u] pos=0.776249816573 neg=0.223750183427 sum= 1
mu1=0.00282282714984 mu2=0.000666637965285 mu3=3.85005662434e-06 mu4=0.000243645610303


Memory usage report
===================
  after diag: 44052 kB
after recalc: 51316 kB
 after split: 49804 kB
  after trim: 38820 kB
 after trunc: 44184 kB
      diag_h: 32500 kB

Peak usage: 51316 kB

Timing report
=============
         All: 17 s
          DM: 0 s
         FDM: 0 s
     broaden: 0 s
        diag: 5 s
    recalc d: 1 s
    recalc f: 0 s
    recalc s: 0 s
      spec d: 4 s
       Other: 3 s
