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.25
discretization=Z
band=flat
betabar=1
T=1e-10
T_ratio=1
Tmin_ratio=-999
keep=10000
keepenergy=10
keepmin=0
keepblock=0
keepzrescale=0
log=
logall=0
logenumber=10
smooth=new
alpha=0.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.25
>>  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.25
# keep=10000 keepenergy=10 Nmax=43 channels=1 dots=1 betabar=1

Number of invariant subspaces: 9
(-2 1) 1 states
0.0221315215209 
(-1 0) 2 states
0.0104939174264 0.0316613444408 
(-1 2) 2 states
0.0105168442488 0.0316801558595 
(0 -1) 1 states
0.0200237403463 
(0 1) 4 states
0 0.0200446115157 0.0221315215209 0.0421761289389 
(0 3) 1 states
0.0200654785874 
(1 0) 2 states
0.0104939174264 0.0316613444408 
(1 2) 2 states
0.0105168442488 0.0316801558595 
(2 1) 1 states
0.0221315215209 
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>=-5.2226274993e-06

 ***** [NRG] Iteration 1/43 (scale 1.3832644937419567) ***** 
[1] xi(0)=0.419770954375 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=1.68430550283 nrkept=64 nrkeptmult=64 
Total energy=-1.16244495302  Egs=-0.80564519025
Kept: 64 out of 64, ratio=1
<SZd>=-9.05243208127e-06
Memory used: 4 MB  Time elapsed: 0 s

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

 ***** [NRG] Iteration 3/43 (scale 0.46108816458065222) ***** 
[1] xi(2)=1.25250020132 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Safeguard: keep additional 1 states
Emax=10.5491930386 nrkept=1006 nrkeptmult=1006 
Total energy=-2.60872723536  Egs=-1.96710866831
Kept: 1006 out of 1024, ratio=0.982
<SZd>=-2.81298018408e-05
Memory used: 7 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/43 (scale 0.26620937594079003) ***** 
[1] xi(3)=1.24631181837 zeta(4)=0
Stats: nr=49 min=1 max=388 sum3=213989362
Safeguard: keep additional 3 states
Emax=10.3900401307 nrkept=1700 nrkeptmult=1700 
Total energy=-2.85582041165  Egs=-0.928191110551
Kept: 1700 out of 4024, ratio=0.422
<SZd>=-4.73500215468e-05
Memory used: 19 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/43 (scale 0.1536960548602174) ***** 
[1] xi(4)=1.06045468118 zeta(5)=0
Stats: nr=60 min=2 max=516 sum3=817346432
Safeguard: keep additional 3 states
Emax=10.10282057 nrkept=1108 nrkeptmult=1108 
Total energy=-3.10556547256  Egs=-1.62492824642
Kept: 1108 out of 6800, ratio=0.163
<SZd>=-8.62878541871e-05
Memory used: 22 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 6/43 (scale 0.088736458646930014) ***** 
[1] xi(5)=1.01240866565 zeta(6)=0
Stats: nr=53 min=1 max=420 sum3=276080440
Safeguard: keep additional 14 states
Emax=10.2982971075 nrkept=975 nrkeptmult=975 
Total energy=-3.17290764033  Egs=-0.758900781023
Kept: 975 out of 4432, ratio=0.22
<SZd>=-0.000144144748155
Memory used: 24 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 7/43 (scale 0.051232018286739137) ***** 
[1] xi(6)=1.00352829264 zeta(7)=0
Stats: nr=60 min=1 max=307 sum3=167640828
Safeguard: keep additional 1 states
Emax=10.6287998149 nrkept=866 nrkeptmult=866 
Total energy=-3.25118789774  Egs=-1.52795575944
Kept: 866 out of 3900, ratio=0.222
<SZd>=-0.000254576403959
Memory used: 24 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 8/43 (scale 0.029578819548976672) ***** 
[1] xi(7)=1.00113318518 zeta(8)=0
Stats: nr=49 min=1 max=332 sum3=134201072
Safeguard: keep additional 3 states
Emax=10.3342802339 nrkept=964 nrkeptmult=964 
Total energy=-3.27494729199  Egs=-0.803257013474
Kept: 964 out of 3464, ratio=0.278
<SZd>=-0.000449993521769
Memory used: 25 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 9/43 (scale 0.017077339428913048) ***** 
[1] xi(8)=1.00037363698 zeta(9)=0
Stats: nr=60 min=1 max=304 sum3=160889638
Safeguard: keep additional 3 states
Emax=10.0616915243 nrkept=804 nrkeptmult=804 
Total energy=-3.29963840816  Egs=-1.44584092097
Kept: 804 out of 3856, ratio=0.209
<SZd>=-0.000810694794586
Memory used: 26 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 10/43 (scale 0.0098596065163255569) ***** 
[1] xi(9)=1.00012410789 zeta(10)=0
Stats: nr=49 min=1 max=310 sum3=108932430
Emax=10.4427997471 nrkept=961 nrkeptmult=961 
Total energy=-3.30809403877  Egs=-0.857603251961
Kept: 961 out of 3216, ratio=0.299
<SZd>=-0.00148786106335
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 11/43 (scale 0.0056924464763043487) ***** 
[1] xi(10)=1.0000413212 zeta(11)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Safeguard: keep additional 2 states
Emax=10.0310857578 nrkept=751 nrkeptmult=751 
Total energy=-3.31586922948  Egs=-1.3658785802
Kept: 751 out of 3844, ratio=0.195
<SZd>=-0.0027838730001
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 12/43 (scale 0.0032865355054418526) ***** 
[1] xi(11)=1.00001376841 zeta(12)=0
Stats: nr=49 min=1 max=286 sum3=86414494
Safeguard: keep additional 1 states
Emax=10.0705916227 nrkept=774 nrkeptmult=774 
Total energy=-3.31884367362  Egs=-0.905039404544
Kept: 774 out of 3004, ratio=0.258
<SZd>=-0.00520900860123
Memory used: 28 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 13/43 (scale 0.0018974821587681161) ***** 
[1] xi(12)=1.00000458888 zeta(13)=0
Stats: nr=56 min=1 max=254 sum3=91449822
Emax=10.0403708897 nrkept=617 nrkeptmult=617 
Total energy=-3.32129824607  Egs=-1.29359448427
Kept: 617 out of 3096, ratio=0.199
<SZd>=-0.00950338009746
Memory used: 28 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 14/43 (scale 0.001095511835147284) ***** 
[1] xi(13)=1.00000152956 zeta(14)=0
Stats: nr=47 min=1 max=238 sum3=49756532
Safeguard: keep additional 1 states
Emax=10.0177108628 nrkept=589 nrkeptmult=589 
Total energy=-3.32233856983  Egs=-0.949623473478
Kept: 589 out of 2468, ratio=0.239
<SZd>=-0.0165427512835
Memory used: 28 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 15/43 (scale 0.00063249405292270544) ***** 
[1] xi(14)=1.00000050985 zeta(15)=0
Stats: nr=56 min=1 max=197 sum3=41125708
Safeguard: keep additional 2 states
Emax=10.0036182465 nrkept=459 nrkeptmult=459 
Total energy=-3.32311509854  Egs=-1.22772491647
Kept: 459 out of 2356, ratio=0.195
<SZd>=-0.0276778085784
Memory used: 28 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 16/43 (scale 0.00036517061171576138) ***** 
[1] xi(15)=1.00000016995 zeta(16)=0
Stats: nr=47 min=1 max=181 sum3=21435702
Safeguard: keep additional 2 states
Emax=10.2581407857 nrkept=523 nrkeptmult=523 
Total energy=-3.32347775264  Egs=-0.993108667193
Kept: 523 out of 1836, ratio=0.285
<SZd>=-0.0451026719294
Memory used: 28 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 17/43 (scale 0.00021083135097423514) ***** 
[1] xi(16)=1.00000005665 zeta(17)=0
Stats: nr=56 min=1 max=172 sum3=27536068
Safeguard: keep additional 1 states
Emax=10.0145318967 nrkept=406 nrkeptmult=406 
Total energy=-3.32372579347  Egs=-1.17648933663
Kept: 406 out of 2092, ratio=0.194
<SZd>=-0.0715914600797
Memory used: 29 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 18/43 (scale 0.00012172353723858711) ***** 
[1] xi(17)=1.00000001888 zeta(18)=0
Stats: nr=47 min=1 max=161 sum3=15051442
Emax=10.0176224282 nrkept=440 nrkeptmult=440 
Total energy=-3.32385062457  Egs=-1.02552969198
Kept: 440 out of 1624, ratio=0.271
<SZd>=-0.109547149792
Memory used: 28 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 19/43 (scale 7.0277116991411708e-05) ***** 
[1] xi(18)=1.00000000629 zeta(19)=0
Stats: nr=50 min=1 max=152 sum3=17846258
Safeguard: keep additional 1 states
Emax=10.0004151138 nrkept=376 nrkeptmult=376 
Total energy=-3.32393168837  Egs=-1.15348772594
Kept: 376 out of 1760, ratio=0.214
<SZd>=-0.15763185841
Memory used: 29 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 20/43 (scale 4.0574512412862373e-05) ***** 
[1] xi(19)=1.0000000021 zeta(20)=0
Stats: nr=47 min=1 max=150 sum3=11876536
Safeguard: keep additional 2 states
Emax=10.0055904109 nrkept=369 nrkeptmult=369 
Total energy=-3.32397392678  Egs=-1.04100860113
Kept: 369 out of 1504, ratio=0.245
<SZd>=-0.208009831983
Memory used: 29 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 21/43 (scale 2.3425705663803904e-05) ***** 
[1] xi(20)=1.0000000007 zeta(21)=0
Stats: nr=44 min=1 max=138 sum3=11927064
Safeguard: keep additional 1 states
Emax=10.0304094612 nrkept=320 nrkeptmult=320 
Total energy=-3.32400092711  Egs=-1.15259410086
Kept: 320 out of 1476, ratio=0.217
<SZd>=-0.246080148246
Memory used: 29 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 22/43 (scale 1.3524837470954124e-05) ***** 
[1] xi(21)=1.00000000023 zeta(22)=0
Stats: nr=46 min=1 max=130 sum3=7757846
Safeguard: keep additional 4 states
Emax=10.162253586 nrkept=295 nrkeptmult=295 
Total energy=-3.32401506396  Egs=-1.04525052074
Kept: 295 out of 1280, ratio=0.23
<SZd>=-0.265856467427
Memory used: 29 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 23/43 (scale 7.8085685546013012e-06) ***** 
[1] xi(22)=1.00000000008 zeta(23)=0
Stats: nr=44 min=1 max=121 sum3=6570736
Safeguard: keep additional 1 states
Emax=10.0442052242 nrkept=255 nrkeptmult=255 
Total energy=-3.32402409367  Egs=-1.15638480706
Kept: 255 out of 1180, ratio=0.216
<SZd>=-0.272869571726
Memory used: 29 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 24/43 (scale 4.5082791569847079e-06) ***** 
[1] xi(23)=1.00000000003 zeta(24)=0
Stats: nr=42 min=1 max=107 sum3=4407240
Safeguard: keep additional 3 states
Emax=10.0473306001 nrkept=250 nrkeptmult=250 
Total energy=-3.32402880848  Egs=-1.04581255783
Kept: 250 out of 1020, ratio=0.245
<SZd>=-0.275181673483
Memory used: 29 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 25/43 (scale 2.6028561848671004e-06) ***** 
[1] xi(24)=1.00000000001 zeta(25)=0
Stats: nr=41 min=1 max=105 sum3=4212304
Safeguard: keep additional 2 states
Emax=10.0135713067 nrkept=233 nrkeptmult=233 
Total energy=-3.32403182402  Egs=-1.15854877638
Kept: 233 out of 1000, ratio=0.233
<SZd>=-0.275845504234
Memory used: 29 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 26/43 (scale 1.5027597189949027e-06) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=42 min=1 max=91 sum3=3192902
Safeguard: keep additional 1 states
Emax=10.0073021329 nrkept=234 nrkeptmult=234 
Total energy=-3.32403339559  Egs=-1.04579069452
Kept: 234 out of 932, ratio=0.251
<SZd>=-0.276107990351
Memory used: 29 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 27/43 (scale 8.676187282890334e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=41 min=1 max=95 sum3=3355134
Safeguard: keep additional 1 states
Emax=10.0891281372 nrkept=231 nrkeptmult=231 
Total energy=-3.3240344015  Egs=-1.15939532776
Kept: 231 out of 936, ratio=0.247
<SZd>=-0.276167725174
Memory used: 29 MB  Time elapsed: 8 s

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

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

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

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

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

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

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

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

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

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

 ***** [NRG] Iteration 38/43 (scale 2.0613987914882068e-09) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0923112832 nrkept=232 nrkeptmult=232 
Total energy=-3.32403568721  Egs=-1.04572655467
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276212544128
Memory used: 30 MB  Time elapsed: 9 s

 ***** [NRG] Iteration 39/43 (scale 1.1901491471728853e-09) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0089838327 nrkept=228 nrkeptmult=228 
Total energy=-3.32403568859  Egs=-1.15984290943
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276212519483
Memory used: 30 MB  Time elapsed: 9 s

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

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

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

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

** Iteration completed.


Total energy:             -3.324035690158814
TKW=0.000170357849479
Grand canonical Z=              3.947590489392792
Zft (last shell)=               3.947593465479779
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.25
>>  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.0221315215209 
(-1 0) 2 states
0.0104939174264 0.0316613444408 
(-1 2) 2 states
0.0105168442488 0.0316801558595 
(0 -1) 1 states
0.0200237403463 
(0 1) 4 states
0 0.0200446115157 0.0221315215209 0.0421761289389 
(0 3) 1 states
0.0200654785874 
(1 0) 2 states
0.0104939174264 0.0316613444408 
(1 2) 2 states
0.0105168442488 0.0316801558595 
(2 1) 1 states
0.0221315215209 
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.3832644937419567) ***** 
[1] xi(0)=0.419770954375 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=1.68430550283 nrkept=64 nrkeptmult=64 
Kept: 64 out of 64, ratio=1
Memory used: 39 MB  Time elapsed: 10 s

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

 ***** [DMNRG] Iteration 3/43 (scale 0.46108816458065222) ***** 
[1] xi(2)=1.25250020132 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Safeguard: keep additional 1 states
Emax=10.5491930386 nrkept=1006 nrkeptmult=1006 
Kept: 1006 out of 1024, ratio=0.982
Memory used: 36 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 4/43 (scale 0.26620937594079003) ***** 
[1] xi(3)=1.24631181837 zeta(4)=0
Stats: nr=49 min=1 max=388 sum3=213989362
Safeguard: keep additional 3 states
Emax=10.3900401307 nrkept=1700 nrkeptmult=1700 
Kept: 1700 out of 4024, ratio=0.422
Memory used: 39 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 5/43 (scale 0.1536960548602174) ***** 
[1] xi(4)=1.06045468118 zeta(5)=0
Stats: nr=60 min=2 max=516 sum3=817346432
Safeguard: keep additional 3 states
Emax=10.10282057 nrkept=1108 nrkeptmult=1108 
Kept: 1108 out of 6800, ratio=0.163
Memory used: 40 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 6/43 (scale 0.088736458646930014) ***** 
[1] xi(5)=1.01240866565 zeta(6)=0
Stats: nr=53 min=1 max=420 sum3=276080440
Safeguard: keep additional 14 states
Emax=10.2982971075 nrkept=975 nrkeptmult=975 
Kept: 975 out of 4432, ratio=0.22
Memory used: 41 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 7/43 (scale 0.051232018286739137) ***** 
[1] xi(6)=1.00352829264 zeta(7)=0
Stats: nr=60 min=1 max=307 sum3=167640828
Safeguard: keep additional 1 states
Emax=10.6287998149 nrkept=866 nrkeptmult=866 
Kept: 866 out of 3900, ratio=0.222
Memory used: 40 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 8/43 (scale 0.029578819548976672) ***** 
[1] xi(7)=1.00113318518 zeta(8)=0
Stats: nr=49 min=1 max=332 sum3=134201072
Safeguard: keep additional 3 states
Emax=10.3342802339 nrkept=964 nrkeptmult=964 
Kept: 964 out of 3464, ratio=0.278
Memory used: 40 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 9/43 (scale 0.017077339428913048) ***** 
[1] xi(8)=1.00037363698 zeta(9)=0
Stats: nr=60 min=1 max=304 sum3=160889638
Safeguard: keep additional 3 states
Emax=10.0616915243 nrkept=804 nrkeptmult=804 
Kept: 804 out of 3856, ratio=0.209
Memory used: 40 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 10/43 (scale 0.0098596065163255569) ***** 
[1] xi(9)=1.00012410789 zeta(10)=0
Stats: nr=49 min=1 max=310 sum3=108932430
Emax=10.4427997471 nrkept=961 nrkeptmult=961 
Kept: 961 out of 3216, ratio=0.299
Memory used: 39 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 11/43 (scale 0.0056924464763043487) ***** 
[1] xi(10)=1.0000413212 zeta(11)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Safeguard: keep additional 2 states
Emax=10.0310857578 nrkept=751 nrkeptmult=751 
Kept: 751 out of 3844, ratio=0.195
Memory used: 40 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 12/43 (scale 0.0032865355054418526) ***** 
[1] xi(11)=1.00001376841 zeta(12)=0
Stats: nr=49 min=1 max=286 sum3=86414494
Safeguard: keep additional 1 states
Emax=10.0705916227 nrkept=774 nrkeptmult=774 
Kept: 774 out of 3004, ratio=0.258
Memory used: 39 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 13/43 (scale 0.0018974821587681161) ***** 
[1] xi(12)=1.00000458888 zeta(13)=0
Stats: nr=56 min=1 max=254 sum3=91449822
Emax=10.0403708897 nrkept=617 nrkeptmult=617 
Kept: 617 out of 3096, ratio=0.199
Memory used: 40 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 14/43 (scale 0.001095511835147284) ***** 
[1] xi(13)=1.00000152956 zeta(14)=0
Stats: nr=47 min=1 max=238 sum3=49756532
Safeguard: keep additional 1 states
Emax=10.0177108628 nrkept=589 nrkeptmult=589 
Kept: 589 out of 2468, ratio=0.239
Memory used: 39 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 15/43 (scale 0.00063249405292270544) ***** 
[1] xi(14)=1.00000050985 zeta(15)=0
Stats: nr=56 min=1 max=197 sum3=41125708
Safeguard: keep additional 2 states
Emax=10.0036182465 nrkept=459 nrkeptmult=459 
Kept: 459 out of 2356, ratio=0.195
Memory used: 39 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 16/43 (scale 0.00036517061171576138) ***** 
[1] xi(15)=1.00000016995 zeta(16)=0
Stats: nr=47 min=1 max=181 sum3=21435702
Safeguard: keep additional 2 states
Emax=10.2581407857 nrkept=523 nrkeptmult=523 
Kept: 523 out of 1836, ratio=0.285
Memory used: 38 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 17/43 (scale 0.00021083135097423514) ***** 
[1] xi(16)=1.00000005665 zeta(17)=0
Stats: nr=56 min=1 max=172 sum3=27536068
Safeguard: keep additional 1 states
Emax=10.0145318967 nrkept=406 nrkeptmult=406 
Kept: 406 out of 2092, ratio=0.194
Memory used: 37 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 18/43 (scale 0.00012172353723858711) ***** 
[1] xi(17)=1.00000001888 zeta(18)=0
Stats: nr=47 min=1 max=161 sum3=15051442
Emax=10.0176224282 nrkept=440 nrkeptmult=440 
Kept: 440 out of 1624, ratio=0.271
Memory used: 37 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 19/43 (scale 7.0277116991411708e-05) ***** 
[1] xi(18)=1.00000000629 zeta(19)=0
Stats: nr=50 min=1 max=152 sum3=17846258
Safeguard: keep additional 1 states
Emax=10.0004151138 nrkept=376 nrkeptmult=376 
Kept: 376 out of 1760, ratio=0.214
Memory used: 37 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 20/43 (scale 4.0574512412862373e-05) ***** 
[1] xi(19)=1.0000000021 zeta(20)=0
Stats: nr=47 min=1 max=150 sum3=11876536
Safeguard: keep additional 2 states
Emax=10.0055904109 nrkept=369 nrkeptmult=369 
Kept: 369 out of 1504, ratio=0.245
Memory used: 37 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 21/43 (scale 2.3425705663803904e-05) ***** 
[1] xi(20)=1.0000000007 zeta(21)=0
Stats: nr=44 min=1 max=138 sum3=11927064
Safeguard: keep additional 1 states
Emax=10.0304094612 nrkept=320 nrkeptmult=320 
Kept: 320 out of 1476, ratio=0.217
Memory used: 37 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 22/43 (scale 1.3524837470954124e-05) ***** 
[1] xi(21)=1.00000000023 zeta(22)=0
Stats: nr=46 min=1 max=130 sum3=7757846
Safeguard: keep additional 4 states
Emax=10.162253586 nrkept=295 nrkeptmult=295 
Kept: 295 out of 1280, ratio=0.23
Memory used: 37 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 23/43 (scale 7.8085685546013012e-06) ***** 
[1] xi(22)=1.00000000008 zeta(23)=0
Stats: nr=44 min=1 max=121 sum3=6570736
Safeguard: keep additional 1 states
Emax=10.0442052242 nrkept=255 nrkeptmult=255 
Kept: 255 out of 1180, ratio=0.216
Memory used: 37 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 24/43 (scale 4.5082791569847079e-06) ***** 
[1] xi(23)=1.00000000003 zeta(24)=0
Stats: nr=42 min=1 max=107 sum3=4407240
Safeguard: keep additional 3 states
Emax=10.0473306001 nrkept=250 nrkeptmult=250 
Kept: 250 out of 1020, ratio=0.245
Memory used: 37 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 25/43 (scale 2.6028561848671004e-06) ***** 
[1] xi(24)=1.00000000001 zeta(25)=0
Stats: nr=41 min=1 max=105 sum3=4212304
Safeguard: keep additional 2 states
Emax=10.0135713067 nrkept=233 nrkeptmult=233 
Kept: 233 out of 1000, ratio=0.233
Memory used: 37 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 26/43 (scale 1.5027597189949027e-06) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=42 min=1 max=91 sum3=3192902
Safeguard: keep additional 1 states
Emax=10.0073021329 nrkept=234 nrkeptmult=234 
Kept: 234 out of 932, ratio=0.251
Memory used: 37 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 27/43 (scale 8.676187282890334e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=41 min=1 max=95 sum3=3355134
Safeguard: keep additional 1 states
Emax=10.0891281372 nrkept=231 nrkeptmult=231 
Kept: 231 out of 936, ratio=0.247
Memory used: 37 MB  Time elapsed: 17 s

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

** Iteration completed.

Spectrum: A_d-A_d-d FDM
[A_d-A_d-d] pos=0.223787464121 neg=0.776212535879 sum= 1
mu1=-0.00282263819061 mu2=0.000667635832865 mu3=-3.87564242737e-06 mu4=0.000216804712389

Spectrum: A_d-A_d-u FDM
[A_d-A_d-u] pos=0.776212535879 neg=0.223787464121 sum= 1
mu1=0.00282263819061 mu2=0.000667635832865 mu3=3.87564242736e-06 mu4=0.000216804712389


Memory usage report
===================
  after diag: 54908 kB
after recalc: 66168 kB
 after split: 64624 kB
  after trim: 42468 kB
 after trunc: 55084 kB
      diag_h: 34968 kB

Peak usage: 66168 kB

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