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=42
TEMPERATURE T=1e-10
Nmax=42
Tmin=1e-10
Ninit=0
Lambda=3
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=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.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=3 z=0.75
# keep=10000 keepenergy=10 Nmax=42 channels=1 dots=1 betabar=1

Number of invariant subspaces: 9
(-2 1) 1 states
0.038332919723 
(-1 0) 2 states
0.018175998153 0.0548390572074 
(-1 2) 2 states
0.0182157085742 0.0548716395404 
(0 -1) 1 states
0.0346821356374 
(0 1) 4 states
0 0.0347182855631 0.038332919723 0.0730511981888 
(0 3) 1 states
0.0347544283915 
(1 0) 2 states
0.018175998153 0.0548390572074 
(1 2) 2 states
0.0182157085742 0.0548716395404 
(2 1) 1 states
0.038332919723 
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>=-9.05242337574e-06

 ***** [NRG] Iteration 1/42 (scale 0.79862812782237003) ***** 
[1] xi(0)=0.737093468838 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=2.95738142489 nrkept=64 nrkeptmult=64 
Total energy=-1.17844895023  Egs=-1.4154577632
Kept: 64 out of 64, ratio=1
<SZd>=-1.56983514721e-05
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/42 (scale 0.46108816458065222) ***** 
[1] xi(1)=0.971787846328 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=6.55317114344 nrkept=256 nrkeptmult=256 
Total energy=-1.50829483403  Egs=-0.715364021764
Kept: 256 out of 256, ratio=1
<SZd>=-2.72448664142e-05
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/42 (scale 0.26620937594079003) ***** 
[1] xi(2)=1.03321268925 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Safeguard: keep additional 1 states
Emax=10.1724476763 nrkept=786 nrkeptmult=786 
Total energy=-1.94783641066  Egs=-1.65111230616
Kept: 786 out of 1024, ratio=0.768
<SZd>=-4.90980048713e-05
Memory used: 7 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/42 (scale 0.1536960548602174) ***** 
[1] xi(3)=1.00725946845 zeta(4)=0
Stats: nr=49 min=1 max=300 sum3=99197412
Emax=10.0038586857 nrkept=991 nrkeptmult=991 
Total energy=-2.05963992281  Egs=-0.727432543689
Kept: 991 out of 3144, ratio=0.315
<SZd>=-8.32789990176e-05
Memory used: 13 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/42 (scale 0.088736458646930014) ***** 
[1] xi(4)=1.00017398697 zeta(5)=0
Stats: nr=60 min=1 max=313 sum3=177097810
Safeguard: keep additional 1 states
Emax=10.1916004612 nrkept=866 nrkeptmult=866 
Total energy=-2.19823484608  Egs=-1.56187124643
Kept: 866 out of 3964, ratio=0.218
<SZd>=-0.000145922832772
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 6/42 (scale 0.051232018286739137) ***** 
[1] xi(5)=0.999821553706 zeta(6)=0
Stats: nr=49 min=1 max=332 sum3=134201072
Safeguard: keep additional 7 states
Emax=10.3128834409 nrkept=968 nrkeptmult=968 
Total energy=-2.2378787228  Egs=-0.773810559288
Kept: 968 out of 3464, ratio=0.279
<SZd>=-0.000253499695764
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 7/42 (scale 0.029578819548976672) ***** 
[1] xi(6)=0.999921934048 zeta(7)=0
Stats: nr=60 min=1 max=306 sum3=164094446
Safeguard: keep additional 1 states
Emax=10.0031991768 nrkept=842 nrkeptmult=842 
Total energy=-2.28184876492  Egs=-1.48653809678
Kept: 842 out of 3872, ratio=0.217
<SZd>=-0.000450125360123
Memory used: 18 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 8/42 (scale 0.017077339428913048) ***** 
[1] xi(7)=0.999972235408 zeta(8)=0
Stats: nr=49 min=1 max=324 sum3=125033756
Emax=10.3727745672 nrkept=961 nrkeptmult=961 
Total energy=-2.29603126098  Egs=-0.830486278135
Kept: 961 out of 3368, ratio=0.285
<SZd>=-0.000810719728246
Memory used: 21 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 9/42 (scale 0.0098596065163255569) ***** 
[1] xi(8)=0.999990560771 zeta(9)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Emax=10.0159086216 nrkept=785 nrkeptmult=785 
Total energy=-2.30988304929  Egs=-1.40490275038
Kept: 785 out of 3844, ratio=0.204
<SZd>=-0.00148896800299
Memory used: 22 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 10/42 (scale 0.0056924464763043487) ***** 
[1] xi(9)=0.999996833398 zeta(10)=0
Stats: nr=49 min=1 max=301 sum3=99841658
Safeguard: keep additional 1 states
Emax=10.0949797628 nrkept=938 nrkeptmult=938 
Total energy=-2.31490471794  Egs=-0.882163523986
Kept: 938 out of 3140, ratio=0.299
<SZd>=-0.00278072916331
Memory used: 23 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 11/42 (scale 0.0032865355054418526) ***** 
[1] xi(10)=0.999998942233 zeta(11)=0
Stats: nr=60 min=1 max=294 sum3=147273392
Emax=10.0238818668 nrkept=701 nrkeptmult=701 
Total energy=-2.31927212089  Egs=-1.3288774598
Kept: 701 out of 3752, ratio=0.187
<SZd>=-0.00521134602045
Memory used: 25 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 12/42 (scale 0.0018974821587681161) ***** 
[1] xi(11)=0.999999647163 zeta(12)=0
Stats: nr=49 min=1 max=271 sum3=72225008
Safeguard: keep additional 1 states
Emax=10.0008820373 nrkept=673 nrkeptmult=673 
Total energy=-2.32103147898  Egs=-0.927206657551
Kept: 673 out of 2804, ratio=0.24
<SZd>=-0.00950618929293
Memory used: 26 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 13/42 (scale 0.001095511835147284) ***** 
[1] xi(12)=0.99999988236 zeta(13)=0
Stats: nr=56 min=1 max=225 sum3=62504584
Emax=10.0207007127 nrkept=529 nrkeptmult=529 
Total energy=-2.32241153235  Egs=-1.25973387551
Kept: 529 out of 2692, ratio=0.197
<SZd>=-0.0165420947519
Memory used: 27 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 14/42 (scale 0.00063249405292270544) ***** 
[1] xi(13)=0.999999960784 zeta(14)=0
Stats: nr=47 min=1 max=209 sum3=32643730
Safeguard: keep additional 2 states
Emax=10.042447289 nrkept=531 nrkeptmult=531 
Total energy=-2.32302629119  Egs=-0.97195987873
Kept: 531 out of 2116, ratio=0.251
<SZd>=-0.027681803934
Memory used: 26 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 15/42 (scale 0.00036517061171576138) ***** 
[1] xi(14)=0.999999986928 zeta(15)=0
Stats: nr=56 min=1 max=176 sum3=29456196
Safeguard: keep additional 1 states
Emax=10.0129851248 nrkept=420 nrkeptmult=420 
Total energy=-2.32346423415  Egs=-1.19928313396
Kept: 420 out of 2124, ratio=0.198
<SZd>=-0.0451020550186
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 16/42 (scale 0.00021083135097423514) ***** 
[1] xi(15)=0.999999995642 zeta(16)=0
Stats: nr=47 min=1 max=168 sum3=16944024
Safeguard: keep additional 1 states
Emax=10.0002930586 nrkept=491 nrkeptmult=491 
Total energy=-2.32367745347  Egs=-1.01132646118
Kept: 491 out of 1680, ratio=0.292
<SZd>=-0.0716525808667
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 17/42 (scale 0.00012172353723858711) ***** 
[1] xi(16)=0.999999998547 zeta(17)=0
Stats: nr=54 min=1 max=164 sum3=23208272
Emax=10.0354985639 nrkept=397 nrkeptmult=397 
Total energy=-2.32381880651  Egs=-1.16126296773
Kept: 397 out of 1964, ratio=0.202
<SZd>=-0.109468522084
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 18/42 (scale 7.0277116991411708e-05) ***** 
[1] xi(17)=0.999999999516 zeta(18)=0
Stats: nr=47 min=1 max=156 sum3=13835500
Safeguard: keep additional 1 states
Emax=10.027215911 nrkept=406 nrkeptmult=406 
Total energy=-2.32389155472  Egs=-1.03516216689
Kept: 406 out of 1588, ratio=0.256
<SZd>=-0.157816494479
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 19/42 (scale 4.0574512412862373e-05) ***** 
[1] xi(18)=0.999999999839 zeta(19)=0
Stats: nr=50 min=1 max=144 sum3=14352334
Safeguard: keep additional 1 states
Emax=10.0074648036 nrkept=350 nrkeptmult=350 
Total energy=-2.32393827972  Egs=-1.15158504201
Kept: 350 out of 1624, ratio=0.216
<SZd>=-0.207898194222
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 20/42 (scale 2.3425705663803904e-05) ***** 
[1] xi(19)=0.999999999946 zeta(20)=0
Stats: nr=46 min=1 max=142 sum3=9936320
Safeguard: keep additional 4 states
Emax=10.0376039082 nrkept=329 nrkeptmult=329 
Total energy=-2.32396273377  Egs=-1.04389804528
Kept: 329 out of 1400, ratio=0.235
<SZd>=-0.246303441059
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 21/42 (scale 1.3524837470954124e-05) ***** 
[1] xi(20)=0.999999999982 zeta(21)=0
Stats: nr=44 min=1 max=128 sum3=8692940
Safeguard: keep additional 1 states
Emax=10.0244573243 nrkept=280 nrkeptmult=280 
Total energy=-2.32397835042  Egs=-1.15466448998
Kept: 280 out of 1316, ratio=0.213
<SZd>=-0.265796889759
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 22/42 (scale 7.8085685546013012e-06) ***** 
[1] xi(21)=0.999999999994 zeta(22)=0
Stats: nr=44 min=1 max=117 sum3=5525668
Emax=10.0132959916 nrkept=254 nrkeptmult=254 
Total energy=-2.32398651523  Egs=-1.04562117273
Kept: 254 out of 1120, ratio=0.227
<SZd>=-0.273028224731
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 23/42 (scale 4.5082791569847079e-06) ***** 
[1] xi(22)=0.999999999998 zeta(23)=0
Stats: nr=42 min=1 max=108 sum3=4542524
Safeguard: keep additional 3 states
Emax=10.0161503955 nrkept=243 nrkeptmult=243 
Total energy=-2.32399173482  Egs=-1.15777863884
Kept: 243 out of 1016, ratio=0.239
<SZd>=-0.275186633862
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 24/42 (scale 2.6028561848671004e-06) ***** 
[1] xi(23)=0.999999999999 zeta(24)=0
Stats: nr=42 min=1 max=99 sum3=3741714
Safeguard: keep additional 3 states
Emax=10.0068195492 nrkept=242 nrkeptmult=242 
Total energy=-2.32399445671  Egs=-1.04573326062
Kept: 242 out of 972, ratio=0.249
<SZd>=-0.275946742419
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 25/42 (scale 1.5027597189949027e-06) ***** 
[1] xi(24)=1 zeta(25)=0
Stats: nr=41 min=1 max=101 sum3=3790166
Emax=10.0447839993 nrkept=230 nrkeptmult=230 
Total energy=-2.32399619865  Egs=-1.15916398573
Kept: 230 out of 968, ratio=0.238
<SZd>=-0.276150571876
Memory used: 27 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 26/42 (scale 8.676187282890334e-07) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=40 min=1 max=91 sum3=3185354
Safeguard: keep additional 1 states
Emax=10.0383843568 nrkept=232 nrkeptmult=232 
Total energy=-2.32399710591  Egs=-1.0456865993
Kept: 232 out of 920, ratio=0.252
<SZd>=-0.276245175044
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 27/42 (scale 5.009199063316343e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.1040183156 nrkept=231 nrkeptmult=231 
Total energy=-2.32399768681  Egs=-1.15966922343
Kept: 231 out of 928, ratio=0.249
<SZd>=-0.276260486715
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 28/42 (scale 2.8920624276301115e-07) ***** 
[1] xi(27)=1 zeta(28)=0
Stats: nr=40 min=1 max=93 sum3=3306120
Safeguard: keep additional 1 states
Emax=10.0790872039 nrkept=232 nrkeptmult=232 
Total energy=-2.32399798922  Egs=-1.04565994466
Kept: 232 out of 924, ratio=0.251
<SZd>=-0.276274975695
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 29/42 (scale 1.6697330211054477e-07) ***** 
[1] xi(28)=1 zeta(29)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.000349477 nrkept=228 nrkeptmult=228 
Total energy=-2.32399818289  Egs=-1.15984261291
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276274335681
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 30/42 (scale 9.6402080921003717e-08) ***** 
[1] xi(29)=1 zeta(30)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0882271657 nrkept=232 nrkeptmult=232 
Total energy=-2.32399828369  Egs=-1.04564976204
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276277309577
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 31/42 (scale 5.5657767370181585e-08) ***** 
[1] xi(30)=1 zeta(31)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0059609497 nrkept=228 nrkeptmult=228 
Total energy=-2.32399834825  Egs=-1.15990095779
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.27627646376
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 32/42 (scale 3.213402697366791e-08) ***** 
[1] xi(31)=1 zeta(32)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0910171662 nrkept=232 nrkeptmult=232 
Total energy=-2.32399838185  Egs=-1.04564623048
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276277248207
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 33/42 (scale 1.8552589123393859e-08) ***** 
[1] xi(32)=1 zeta(33)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0078504861 nrkept=228 nrkeptmult=228 
Total energy=-2.32399840337  Egs=-1.15992045701
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276276896395
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 34/42 (scale 1.0711342324555969e-08) ***** 
[1] xi(33)=1 zeta(34)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0919312226 nrkept=232 nrkeptmult=232 
Total energy=-2.32399841457  Egs=-1.0456450476
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276277134977
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 35/42 (scale 6.1841963744646209e-09) ***** 
[1] xi(34)=1 zeta(35)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0084839201 nrkept=228 nrkeptmult=228 
Total energy=-2.32399842174  Egs=-1.1599269522
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276277010013
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 36/42 (scale 3.5704474415186564e-09) ***** 
[1] xi(35)=1 zeta(36)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0922343666 nrkept=232 nrkeptmult=232 
Total energy=-2.32399842547  Egs=-1.04564466234
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276277086995
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 37/42 (scale 2.0613987914882068e-09) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0087019684 nrkept=228 nrkeptmult=228 
Total energy=-2.32399842786  Egs=-1.15992910654
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.2762770445
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 38/42 (scale 1.1901491471728853e-09) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 3 states
Emax=10.0923380112 nrkept=234 nrkeptmult=234 
Total energy=-2.32399842911  Egs=-1.04564454459
Kept: 234 out of 912, ratio=0.257
<SZd>=-0.276277069864
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 39/42 (scale 6.8713293049606891e-10) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=41 min=1 max=97 sum3=3398916
Safeguard: keep additional 1 states
Emax=10.0086841986 nrkept=228 nrkeptmult=228 
Total energy=-2.32399842991  Egs=-1.15992981329
Kept: 228 out of 936, ratio=0.244
<SZd>=-0.276277055619
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 40/42 (scale 3.9671638239096179e-10) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923691959 nrkept=238 nrkeptmult=238 
Total energy=-2.32399843032  Egs=-1.04564451613
Kept: 238 out of 912, ratio=0.261
<SZd>=-0.276277064034
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 41/42 (scale 2.2904431016535628e-10) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=41 min=1 max=99 sum3=3556156
Safeguard: keep additional 1 states
Emax=10.0087076729 nrkept=228 nrkeptmult=228 
Total energy=-2.32399843059  Egs=-1.15993003742
Kept: 228 out of 952, ratio=0.239
<SZd>=-0.276277059284
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 42/42 (scale 1.3223879413032059e-10) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923803136 nrkept=912 nrkeptmult=912 
Total energy=-2.32399843073  Egs=-1.0456445175
Kept: 912 out of 912, ratio=1
<SZd>=-0.276277062262
Memory used: 28 MB  Time elapsed: 6 s

** Iteration completed.


Total energy:             -2.323998430725157
TKW=0.000171641554296
Grand canonical Z=              4.170721469916313
Zft (last shell)=                4.17071846568201
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.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.038332919723 
(-1 0) 2 states
0.018175998153 0.0548390572074 
(-1 2) 2 states
0.0182157085742 0.0548716395404 
(0 -1) 1 states
0.0346821356374 
(0 1) 4 states
0 0.0347182855631 0.038332919723 0.0730511981888 
(0 3) 1 states
0.0347544283915 
(1 0) 2 states
0.018175998153 0.0548390572074 
(1 2) 2 states
0.0182157085742 0.0548716395404 
(2 1) 1 states
0.038332919723 
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/42 (scale 0.79862812782237003) ***** 
[1] xi(0)=0.737093468838 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=2.95738142489 nrkept=64 nrkeptmult=64 
Kept: 64 out of 64, ratio=1
Memory used: 32 MB  Time elapsed: 7 s

 ***** [DMNRG] Iteration 2/42 (scale 0.46108816458065222) ***** 
[1] xi(1)=0.971787846328 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=6.55317114344 nrkept=256 nrkeptmult=256 
Kept: 256 out of 256, ratio=1
Memory used: 32 MB  Time elapsed: 7 s

 ***** [DMNRG] Iteration 3/42 (scale 0.26620937594079003) ***** 
[1] xi(2)=1.03321268925 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Safeguard: keep additional 1 states
Emax=10.1724476763 nrkept=786 nrkeptmult=786 
Kept: 786 out of 1024, ratio=0.768
Memory used: 33 MB  Time elapsed: 7 s

 ***** [DMNRG] Iteration 4/42 (scale 0.1536960548602174) ***** 
[1] xi(3)=1.00725946845 zeta(4)=0
Stats: nr=49 min=1 max=300 sum3=99197412
Emax=10.0038586857 nrkept=991 nrkeptmult=991 
Kept: 991 out of 3144, ratio=0.315
Memory used: 32 MB  Time elapsed: 8 s

 ***** [DMNRG] Iteration 5/42 (scale 0.088736458646930014) ***** 
[1] xi(4)=1.00017398697 zeta(5)=0
Stats: nr=60 min=1 max=313 sum3=177097810
Safeguard: keep additional 1 states
Emax=10.1916004612 nrkept=866 nrkeptmult=866 
Kept: 866 out of 3964, ratio=0.218
Memory used: 34 MB  Time elapsed: 8 s

 ***** [DMNRG] Iteration 6/42 (scale 0.051232018286739137) ***** 
[1] xi(5)=0.999821553706 zeta(6)=0
Stats: nr=49 min=1 max=332 sum3=134201072
Safeguard: keep additional 7 states
Emax=10.3128834409 nrkept=968 nrkeptmult=968 
Kept: 968 out of 3464, ratio=0.279
Memory used: 33 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 7/42 (scale 0.029578819548976672) ***** 
[1] xi(6)=0.999921934048 zeta(7)=0
Stats: nr=60 min=1 max=306 sum3=164094446
Safeguard: keep additional 1 states
Emax=10.0031991768 nrkept=842 nrkeptmult=842 
Kept: 842 out of 3872, ratio=0.217
Memory used: 35 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 8/42 (scale 0.017077339428913048) ***** 
[1] xi(7)=0.999972235408 zeta(8)=0
Stats: nr=49 min=1 max=324 sum3=125033756
Emax=10.3727745672 nrkept=961 nrkeptmult=961 
Kept: 961 out of 3368, ratio=0.285
Memory used: 35 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 9/42 (scale 0.0098596065163255569) ***** 
[1] xi(8)=0.999990560771 zeta(9)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Emax=10.0159086216 nrkept=785 nrkeptmult=785 
Kept: 785 out of 3844, ratio=0.204
Memory used: 36 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 10/42 (scale 0.0056924464763043487) ***** 
[1] xi(9)=0.999996833398 zeta(10)=0
Stats: nr=49 min=1 max=301 sum3=99841658
Safeguard: keep additional 1 states
Emax=10.0949797628 nrkept=938 nrkeptmult=938 
Kept: 938 out of 3140, ratio=0.299
Memory used: 36 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 11/42 (scale 0.0032865355054418526) ***** 
[1] xi(10)=0.999998942233 zeta(11)=0
Stats: nr=60 min=1 max=294 sum3=147273392
Emax=10.0238818668 nrkept=701 nrkeptmult=701 
Kept: 701 out of 3752, ratio=0.187
Memory used: 37 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 12/42 (scale 0.0018974821587681161) ***** 
[1] xi(11)=0.999999647163 zeta(12)=0
Stats: nr=49 min=1 max=271 sum3=72225008
Safeguard: keep additional 1 states
Emax=10.0008820373 nrkept=673 nrkeptmult=673 
Kept: 673 out of 2804, ratio=0.24
Memory used: 36 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 13/42 (scale 0.001095511835147284) ***** 
[1] xi(12)=0.99999988236 zeta(13)=0
Stats: nr=56 min=1 max=225 sum3=62504584
Emax=10.0207007127 nrkept=529 nrkeptmult=529 
Kept: 529 out of 2692, ratio=0.197
Memory used: 36 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 14/42 (scale 0.00063249405292270544) ***** 
[1] xi(13)=0.999999960784 zeta(14)=0
Stats: nr=47 min=1 max=209 sum3=32643730
Safeguard: keep additional 2 states
Emax=10.042447289 nrkept=531 nrkeptmult=531 
Kept: 531 out of 2116, ratio=0.251
Memory used: 34 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 15/42 (scale 0.00036517061171576138) ***** 
[1] xi(14)=0.999999986928 zeta(15)=0
Stats: nr=56 min=1 max=176 sum3=29456196
Safeguard: keep additional 1 states
Emax=10.0129851248 nrkept=420 nrkeptmult=420 
Kept: 420 out of 2124, ratio=0.198
Memory used: 34 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 16/42 (scale 0.00021083135097423514) ***** 
[1] xi(15)=0.999999995642 zeta(16)=0
Stats: nr=47 min=1 max=168 sum3=16944024
Safeguard: keep additional 1 states
Emax=10.0002930586 nrkept=491 nrkeptmult=491 
Kept: 491 out of 1680, ratio=0.292
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 17/42 (scale 0.00012172353723858711) ***** 
[1] xi(16)=0.999999998547 zeta(17)=0
Stats: nr=54 min=1 max=164 sum3=23208272
Emax=10.0354985639 nrkept=397 nrkeptmult=397 
Kept: 397 out of 1964, ratio=0.202
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 18/42 (scale 7.0277116991411708e-05) ***** 
[1] xi(17)=0.999999999516 zeta(18)=0
Stats: nr=47 min=1 max=156 sum3=13835500
Safeguard: keep additional 1 states
Emax=10.027215911 nrkept=406 nrkeptmult=406 
Kept: 406 out of 1588, ratio=0.256
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 19/42 (scale 4.0574512412862373e-05) ***** 
[1] xi(18)=0.999999999839 zeta(19)=0
Stats: nr=50 min=1 max=144 sum3=14352334
Safeguard: keep additional 1 states
Emax=10.0074648036 nrkept=350 nrkeptmult=350 
Kept: 350 out of 1624, ratio=0.216
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 20/42 (scale 2.3425705663803904e-05) ***** 
[1] xi(19)=0.999999999946 zeta(20)=0
Stats: nr=46 min=1 max=142 sum3=9936320
Safeguard: keep additional 4 states
Emax=10.0376039082 nrkept=329 nrkeptmult=329 
Kept: 329 out of 1400, ratio=0.235
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 21/42 (scale 1.3524837470954124e-05) ***** 
[1] xi(20)=0.999999999982 zeta(21)=0
Stats: nr=44 min=1 max=128 sum3=8692940
Safeguard: keep additional 1 states
Emax=10.0244573243 nrkept=280 nrkeptmult=280 
Kept: 280 out of 1316, ratio=0.213
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 22/42 (scale 7.8085685546013012e-06) ***** 
[1] xi(21)=0.999999999994 zeta(22)=0
Stats: nr=44 min=1 max=117 sum3=5525668
Emax=10.0132959916 nrkept=254 nrkeptmult=254 
Kept: 254 out of 1120, ratio=0.227
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 23/42 (scale 4.5082791569847079e-06) ***** 
[1] xi(22)=0.999999999998 zeta(23)=0
Stats: nr=42 min=1 max=108 sum3=4542524
Safeguard: keep additional 3 states
Emax=10.0161503955 nrkept=243 nrkeptmult=243 
Kept: 243 out of 1016, ratio=0.239
Memory used: 34 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 24/42 (scale 2.6028561848671004e-06) ***** 
[1] xi(23)=0.999999999999 zeta(24)=0
Stats: nr=42 min=1 max=99 sum3=3741714
Safeguard: keep additional 3 states
Emax=10.0068195492 nrkept=242 nrkeptmult=242 
Kept: 242 out of 972, ratio=0.249
Memory used: 34 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 25/42 (scale 1.5027597189949027e-06) ***** 
[1] xi(24)=1 zeta(25)=0
Stats: nr=41 min=1 max=101 sum3=3790166
Emax=10.0447839993 nrkept=230 nrkeptmult=230 
Kept: 230 out of 968, ratio=0.238
Memory used: 34 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 26/42 (scale 8.676187282890334e-07) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=40 min=1 max=91 sum3=3185354
Safeguard: keep additional 1 states
Emax=10.0383843568 nrkept=232 nrkeptmult=232 
Kept: 232 out of 920, ratio=0.252
Memory used: 34 MB  Time elapsed: 14 s

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

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

 ***** [DMNRG] Iteration 29/42 (scale 1.6697330211054477e-07) ***** 
[1] xi(28)=1 zeta(29)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.000349477 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 34 MB  Time elapsed: 14 s

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

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

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

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

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

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

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

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

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

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

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

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

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

** Iteration completed.

Spectrum: A_d-A_d-d FDM
[A_d-A_d-d] pos=0.223722947232 neg=0.776277052767 sum= 1
mu1=-0.00283315143311 mu2=0.000669572373817 mu3=-3.98214536918e-06 mu4=0.000225758041052

Spectrum: A_d-A_d-u FDM
[A_d-A_d-u] pos=0.776277052767 neg=0.223722947232 sum= 1
mu1=0.00283315143311 mu2=0.000669572373817 mu3=3.98214536918e-06 mu4=0.000225758041052


Memory usage report
===================
  after diag: 39972 kB
after recalc: 45472 kB
 after split: 42620 kB
  after trim: 38736 kB
 after trunc: 40004 kB
      diag_h: 29428 kB

Peak usage: 45472 kB

Timing report
=============
         All: 16 s
          DM: 0 s
         FDM: 0 s
     broaden: 0 s
        diag: 4 s
    recalc d: 1 s
    recalc f: 0 s
    recalc s: 0 s
      spec d: 4 s
       Other: 2 s
