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=34
TEMPERATURE T=1e-09
Nmax=34
Tmin=1e-10
Ninit=0
Lambda=4
strategy=kept
lastall=0
lastalloverride=1
z=0.25
discretization=Z
band=asymode
betabar=1
T=1e-09
T_ratio=1
Tmin_ratio=-999
keep=8000
keepenergy=10
keepmin=500
keepblock=0
keepzrescale=0
log=
logall=0
logenumber=10
smooth=old
alpha=0.3
omega0=-1
omega0_ratio=1
brcut=0
loggauss_b=0.3
gauss_b=0.1
limitLL=0
eta=0
finite=0
bins=1000
discard_trim=1e-16
discard_immediately=1e-16
NN1=0
NN2even=1
NN2avg=0
NNtanh=0
goodE=2
savebins=1
broaden=0
broaden_max=4
broaden_min=1e-08
broaden_min_ratio=3
broaden_ratio=1.01
dumpenergies=0
saveall=0
dumpannotated=0
dumpscaled=1
dumpabs=0
dumpprecision=8
dumpgroups=1
grouptol=1e-06
dumpdiagonal=0
ops=A_d self_d n_d n_d_ud Himp Hhyb Hpot SdSk
specs=
specd=A_d-A_d self_d-A_d
spect=
specb=
reim=0
trsq=
specgt=
gtp=0.7
speci1t=
speci2t=
speckubo=
kubop=0.7
kubor=4
kubob=2
cond=
specchit=
chitp=1
chitp_ratio=-999
dm=1
T0opt=-1
dmnrg=0
cfs=0
fdm=1
fdmexpv=1
fdmexpvn=0
mats=0
nromegan=100
width_td=12
width_custom=20
prec_td=6
prec_custom=14
prec_xy=12
safeguard=1e-05
safeguardmax=200
fixeps=1e-10
diag=dsyev
diagratio=1
dsyevrlimit=100
restart=1
restartfactor=2
checkdiag=0
checkrho=0
calc0=1
tdht=0
spin=2
tri=cpp
preccpp=2000
checksumrules=0
showEs=0
globalB=0
globalBx=0
polarized=0
pol2x2=0
mpi=0
saveF=0
trunc=
noimag=1
dumpsubspaces=0
done=1
resume=0
forcestop=-1
stopafter=
removefiles=1
store=one

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QS
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	../model.m	 Variant: 		 Channels: 	1
>>  U=	0	 Gamma=	-999	 Delta=	0.024934636296339015	 t=	0.
>> 	 extraJkondo=	0.05
>>  Gamma^(1/2)=	0.05050305102594391
>>  Lambda=	4.	 BAND=	asymode
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d self_d n_d n_d_ud Himp Hhyb Hpot SdSk
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=8
Reading <||||> (e)
Reading <||Himp||> (s)
Reading <||Hhyb||> (s)
Reading <||Hpot||> (s)
Reading <||n_d||> (s)
Reading <||n_d_ud||> (s)
Reading <||A_d||> (d)
Reading <||self_d||> (d)
Reading <||||> (z)
Reading <||||> (T)
Tridiagonalisation, ch=0. Using GMP version 5.0.5
mMAX=80
Using precision of 2000 digits.
  xi(0)=            0.1105468817135617 -->            0.07224298735651534
zeta(0)=           0.02442593583648851
  xi(1)=            0.7234472345094897 -->             0.9455533906636359
zeta(1)=          -0.03089214812898047
  xi(2)=            0.5989449789337551 -->              1.565655181571762
zeta(2)=           0.01330425975981917
  xi(3)=            0.2679337207444153 -->              1.400769128904217
zeta(3)=           0.01201927351003285
  xi(4)=            0.3168906938433066 -->              3.313436621112528
zeta(4)=          -0.02853076788897458
  xi(5)=           0.08386979595491127 -->               1.75389974348474
zeta(5)=           0.03919705799554069
  xi(6)=           0.03033998727052085 -->              1.268950169372165
zeta(6)=          -0.01060564161230866
  xi(7)=           0.01092063679677809 -->             0.9134970156291409
zeta(7)=         -0.003263196940067326
  xi(8)=          0.004835265369878517 -->             0.8089272754610982
zeta(8)=         -0.003908710439885663
  xi(9)=          0.003587669339060117 -->              1.200415431914189
zeta(9)=          0.001824653764073181
  xi(10)=          0.001357410307157379 -->             0.9083648052013483
zeta(10)=         0.0005944939355874574
  xi(11)=         0.0008481319063189717 -->              1.135122033192516
zeta(11)=         0.0003196118558607157
  xi(12)=          0.000379846045821153 -->              1.016756031980918
zeta(12)=         0.0001913655618539865
  xi(13)=         0.0001940934653505897 -->              1.039082564288216
zeta(13)=         4.969797885117206e-05
  xi(14)=         9.533205531511519e-05 -->               1.02072345729691
zeta(14)=         2.223368999255758e-05
  xi(15)=         4.697306348067651e-05 -->              1.005884276749067
zeta(15)=         5.706003242836275e-06
  xi(16)=         2.335589526038261e-05 -->              1.000289360368479
zeta(16)=         1.457145079326178e-06
  xi(17)=         1.165825739413104e-05 -->             0.9986027683184079
zeta(17)=         3.787893320652774e-07
  xi(18)=         5.826201878451904e-06 -->             0.9981013676251282
zeta(18)=         1.023690165621956e-07
  xi(19)=         2.912744611189157e-06 -->             0.9979792807121554
zeta(19)=         2.654494031716698e-08
  xi(20)=         1.456269886949164e-06 -->             0.9979090983242525
zeta(20)=         7.365097918006369e-09
  xi(21)=         7.281803465825329e-07 -->             0.9979713232935696
zeta(21)=          1.72424644448412e-09
  xi(22)=         3.640466222485405e-07 -->             0.9978519498665223
zeta(22)=         5.252661340474692e-10
  xi(23)=         1.820741837979044e-07 -->             0.9981308339076226
zeta(23)=         7.485959323696829e-11
  xi(24)=         9.100355089213651e-08 -->             0.9977630902505622
zeta(24)=         3.904712346352647e-11
  xi(25)=         4.552586078924802e-08 -->             0.9982912337395843
zeta(25)=        -8.954542227472203e-12
  xi(26)=         2.282263155003591e-08 -->              1.000909488026646
zeta(26)=          2.44651171170514e-12
  xi(27)=         1.144669716871736e-08 -->              1.004012861323079
zeta(27)=        -2.049894656742349e-12
  xi(28)=         5.715567282279033e-09 -->              1.002647833970532
zeta(28)=         -1.06333982247079e-11
  xi(29)=         2.853268654161199e-09 -->              1.001063759567891
zeta(29)=         3.798313022109058e-12
  xi(30)=         1.427137177649404e-09 -->               1.00141660785659
zeta(30)=         2.096657167480259e-12
  xi(31)=         7.135666403588892e-10 -->              1.001413873394822
zeta(31)=         1.249719021397554e-13
  xi(32)=         3.568249620677095e-10 -->              1.001530753199073
zeta(32)=         1.340611871970969e-13
  xi(33)=         1.784123576832135e-10 -->              1.001530060761751
zeta(33)=        -6.939107657203172e-14
  xi(34)=         8.920881694138008e-11 -->              1.001559679071386
zeta(34)=         8.384760317852646e-15

# $Id: nrg.cc,v 1.31 2012/10/26 10:18:22 rokzitko Exp rokzitko $ QS
# disk=Z band=asymode Lambda=4 z=0.25
# keep=8000 keepenergy=10 Nmax=34 channels=1 dots=1 betabar=1

Number of invariant subspaces: 8
(-2 2) 1 states
0.0178591316151 
(-1 1) 2 states
0.00222301767681 0.0373707196932 
(-1 3) 2 states
0.0113273794983 0.0446039935961 
(0 2) 4 states
0 0.0261978645765 0.0339879879407 0.0675966598135 
(0 4) 1 states
0.0380722414794 
(1 1) 2 states
0.0182938663113 0.053557272925 
(1 3) 2 states
0.0274764948932 0.0607122800675 
(2 2) 1 states
0.0501165334814 
Number of states (multiplicity taken into account): 32

Computing the following operators:
s=[]
p=[]
g=[]
d=[A_d self_d ]
v=[]
t=[]

Computing the following spectra:

Before NRG iteration (N=-1)
<Hhyb>=0.0238253435111
<Himp>=0.0248330680231
<Hpot>=0
<n_d>=0.995926619019
<n_d_ud>=0.247964321549

 ***** [NRG] Iteration 1/34 (scale 1.5302091698951845) ***** 
[1] xi(0)=0.0722429873565 zeta(1)=-0.030892148129
Stats: nr=13 min=1 max=9 sum3=2864
Emax=0.359128143583 nrkept=56 nrkeptmult=128 
Total energy=-0.257939342701  Egs=-0.116884015895
Kept: 56 out of 56, ratio=1
<Hhyb>=0.0227165485195
<Himp>=0.0247315519583
<Hpot>=0
<n_d>=0.991855331852
<n_d_ud>=0.245931676273
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/34 (scale 0.76510458494759226) ***** 
[1] xi(1)=0.945553390664 zeta(2)=0.0133042597598
Stats: nr=19 min=1 max=35 sum3=99630
Emax=4.09551767314 nrkept=210 nrkeptmult=512 
Total energy=-1.54240848213  Egs=-1.67881511195
Kept: 210 out of 210, ratio=1
<Hhyb>=0.0205028801504
<Himp>=0.0245288846727
<Hpot>=0
<n_d>=0.98372738953
<n_d_ud>=0.241879404359
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/34 (scale 0.38255229247379613) ***** 
[1] xi(2)=1.56565518157 zeta(3)=0.01201927351
Stats: nr=26 min=1 max=99 sum3=3786552
Emax=10.060076228 nrkept=741 nrkeptmult=1930 
Total energy=-2.06587338364  Egs=-1.36834862007
Kept: 741 out of 792, ratio=0.936
<Hhyb>=0.0161103552501
<Himp>=0.0241264331597
<Hpot>=0
<n_d>=0.967587129525
<n_d_ud>=0.233853182476
Memory used: 9 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/34 (scale 0.19127614623689806) ***** 
[1] xi(3)=1.4007691289 zeta(4)=-0.028530767889
Stats: nr=34 min=1 max=334 sum3=122312182
Emax=10.0296900901 nrkept=927 nrkeptmult=2490 
Total energy=-2.42564182124  Egs=-1.88088501717
Kept: 927 out of 2812, ratio=0.33
<Hhyb>=0.00756512048059
<Himp>=0.0233438600804
<Hpot>=0
<n_d>=0.936202148811
<n_d_ud>=0.218309014051
Memory used: 18 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/34 (scale 0.095638073118449032) ***** 
[1] xi(4)=3.31343662111 zeta(5)=0.0391970579955
Stats: nr=40 min=1 max=376 sum3=214773540
Emax=10.5642528426 nrkept=500 nrkeptmult=1295 
Total energy=-2.829124519  Egs=-4.2188501357
Kept: 500 out of 3528, ratio=0.142
<Hhyb>=-0.0079004238417
<Himp>=0.021938367618
<Hpot>=0
<n_d>=0.879835075886
<n_d_ud>=0.190436491893
Memory used: 21 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 6/34 (scale 0.047819036559224516) ***** 
[1] xi(5)=1.75389974348 zeta(6)=-0.0106056416123
Stats: nr=34 min=1 max=224 sum3=37036511
Emax=12.2628450576 nrkept=500 nrkeptmult=1288 
Total energy=-2.8955618913  Egs=-1.38934987152
Kept: 500 out of 1895, ratio=0.264
<Hhyb>=-0.0301285311405
<Himp>=0.0199915422414
<Hpot>=0
<n_d>=0.801757924352
<n_d_ud>=0.150621971085
Memory used: 21 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 7/34 (scale 0.023909518279612258) ***** 
[1] xi(6)=1.26895016937 zeta(7)=-0.00326319694007
Stats: nr=36 min=1 max=215 sum3=39814837
Emax=12.7226205001 nrkept=500 nrkeptmult=1282 
Total energy=-2.93130287976  Egs=-1.49484351975
Kept: 500 out of 1897, ratio=0.264
<Hhyb>=-0.0493592267583
<Himp>=0.0184976411502
<Hpot>=0
<n_d>=0.741845236097
<n_d_ud>=0.113324872871
Memory used: 21 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 8/34 (scale 0.011954759139806129) ***** 
[1] xi(7)=0.913497015629 zeta(8)=-0.00390871043989
Stats: nr=33 min=1 max=229 sum3=38648109
Emax=14.5602312891 nrkept=500 nrkeptmult=1271 
Total energy=-2.94047205892  Egs=-0.766989870469
Kept: 500 out of 1893, ratio=0.264
<Hhyb>=-0.0588498805901
<Himp>=0.0181423508066
<Hpot>=0
<n_d>=0.727596368
<n_d_ud>=0.0909722975412
Memory used: 21 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 9/34 (scale 0.0059773795699030645) ***** 
[1] xi(8)=0.808927275461 zeta(9)=0.00182465376407
Stats: nr=34 min=1 max=215 sum3=39307401
Emax=18.0453035361 nrkept=500 nrkeptmult=1244 
Total energy=-2.94730606996  Egs=-1.14331220861
Kept: 500 out of 1893, ratio=0.264
<Hhyb>=-0.0615492721598
<Himp>=0.018471880351
<Hpot>=0
<n_d>=0.740812102952
<n_d_ud>=0.080341517183
Memory used: 22 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 10/34 (scale 0.0029886897849515323) ***** 
[1] xi(9)=1.20041543191 zeta(10)=0.000594493935587
Stats: nr=32 min=1 max=235 sum3=41372763
Emax=21.2700767316 nrkept=500 nrkeptmult=1210 
Total energy=-2.95050585936  Egs=-1.07063282772
Kept: 500 out of 1887, ratio=0.265
<Hhyb>=-0.0621700782611
<Himp>=0.0190376167343
<Hpot>=0
<n_d>=0.763500879182
<n_d_ud>=0.0775224041173
Memory used: 22 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 11/34 (scale 0.0014943448924757661) ***** 
[1] xi(10)=0.908364805201 zeta(11)=0.000319611855861
Stats: nr=33 min=1 max=223 sum3=41933624
Emax=23.4351296394 nrkept=500 nrkeptmult=1230 
Total energy=-2.9523108312  Egs=-1.20786831522
Kept: 500 out of 1880, ratio=0.266
<Hhyb>=-0.0622845132778
<Himp>=0.0194521900205
<Hpot>=0
<n_d>=0.780127281156
<n_d_ud>=0.0781865880851
Memory used: 23 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 12/34 (scale 0.00074717244623788307) ***** 
[1] xi(11)=1.13512203319 zeta(12)=0.000191365561854
Stats: nr=32 min=1 max=236 sum3=40376669
Emax=25.2412035581 nrkept=500 nrkeptmult=1184 
Total energy=-2.95314144716  Egs=-1.11167905774
Kept: 500 out of 1883, ratio=0.266
<Hhyb>=-0.0622947105955
<Himp>=0.0196953158431
<Hpot>=0
<n_d>=0.789877807282
<n_d_ud>=0.078743309465
Memory used: 23 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 13/34 (scale 0.00037358622311894153) ***** 
[1] xi(12)=1.01675603198 zeta(13)=4.96979788512e-05
Stats: nr=33 min=1 max=230 sum3=42597639
Emax=27.1603339957 nrkept=500 nrkeptmult=1179 
Total energy=-2.95359318406  Egs=-1.20919045952
Kept: 500 out of 1875, ratio=0.267
<Hhyb>=-0.0622861480962
<Himp>=0.019808984227
<Hpot>=0
<n_d>=0.794436461458
<n_d_ud>=0.0789816216302
Memory used: 23 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 14/34 (scale 0.00018679311155947077) ***** 
[1] xi(13)=1.03908256429 zeta(14)=2.22336899926e-05
Stats: nr=32 min=1 max=246 sum3=43888731
Emax=27.3370054041 nrkept=500 nrkeptmult=1170 
Total energy=-2.95379605663  Egs=-1.08608165789
Kept: 500 out of 1875, ratio=0.267
<Hhyb>=-0.0622795863765
<Himp>=0.0198651244209
<Hpot>=0
<n_d>=0.796687955854
<n_d_ud>=0.0790548570096
Memory used: 24 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 15/34 (scale 9.3396555779735383e-05) ***** 
[1] xi(14)=1.0207234573 zeta(15)=5.70600324284e-06
Stats: nr=33 min=1 max=231 sum3=42628830
Emax=27.7614368862 nrkept=500 nrkeptmult=1181 
Total energy=-2.95389768268  Egs=-1.08811346205
Kept: 500 out of 1872, ratio=0.267
<Hhyb>=-0.0622767636612
<Himp>=0.019887141232
<Hpot>=0
<n_d>=0.797570936895
<n_d_ud>=0.0790646925105
Memory used: 25 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 16/34 (scale 4.6698277889867692e-05) ***** 
[1] xi(15)=1.00588427675 zeta(16)=1.45714507933e-06
Stats: nr=32 min=1 max=243 sum3=43287243
Emax=29.789170396 nrkept=500 nrkeptmult=1164 
Total energy=-2.95394808953  Egs=-1.07941562433
Kept: 500 out of 1875, ratio=0.267
<Hhyb>=-0.0622759580388
<Himp>=0.0198944669993
<Hpot>=0
<n_d>=0.797864735738
<n_d_ud>=0.0790616093672
Memory used: 25 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 17/34 (scale 2.3349138944933846e-05) ***** 
[1] xi(16)=1.00028936037 zeta(17)=3.78789332065e-07
Stats: nr=33 min=1 max=233 sum3=42510856
Emax=29.9960263743 nrkept=500 nrkeptmult=1169 
Total energy=-2.95397220099  Egs=-1.03264866476
Kept: 500 out of 1870, ratio=0.267
<Hhyb>=-0.0622757916012
<Himp>=0.0198963614921
<Hpot>=0
<n_d>=0.797940714098
<n_d_ud>=0.0790594283141
Memory used: 25 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 18/34 (scale 1.1674569472466923e-05) ***** 
[1] xi(17)=0.998602768318 zeta(18)=1.02369016562e-07
Stats: nr=33 min=1 max=239 sum3=42523332
Emax=31.9028978168 nrkept=500 nrkeptmult=1153 
Total energy=-2.95398487094  Egs=-1.08526045106
Kept: 500 out of 1872, ratio=0.267
<Hhyb>=-0.0622757272134
<Himp>=0.0198969073346
<Hpot>=0
<n_d>=0.797962605035
<n_d_ud>=0.0790591887142
Memory used: 25 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 19/34 (scale 5.8372847362334614e-06) ***** 
[1] xi(18)=0.998101367625 zeta(19)=2.65449403172e-08
Stats: nr=33 min=1 max=238 sum3=42308259
Emax=32.2985494507 nrkept=500 nrkeptmult=1149 
Total energy=-2.95399082397  Egs=-1.01982861829
Kept: 500 out of 1869, ratio=0.268
<Hhyb>=-0.062275724733
<Himp>=0.0198970096079
<Hpot>=0
<n_d>=0.797966706689
<n_d_ud>=0.0790589092808
Memory used: 25 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 20/34 (scale 2.9186423681167307e-06) ***** 
[1] xi(19)=0.997979280712 zeta(20)=7.36509791801e-09
Stats: nr=33 min=1 max=235 sum3=42007819
Emax=32.3188274108 nrkept=500 nrkeptmult=1145 
Total energy=-2.95399400072  Egs=-1.08843587158
Kept: 500 out of 1867, ratio=0.268
<Hhyb>=-0.0622757175437
<Himp>=0.0198970519997
<Hpot>=0
<n_d>=0.797968406808
<n_d_ud>=0.0790589430349
Memory used: 25 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 21/34 (scale 1.4593211840583654e-06) ***** 
[1] xi(20)=0.997909098324 zeta(21)=1.72424644448e-09
Stats: nr=33 min=1 max=237 sum3=41893189
Emax=32.7844257977 nrkept=500 nrkeptmult=1143 
Total energy=-2.95399548424  Egs=-1.01657673054
Kept: 500 out of 1867, ratio=0.268
<Hhyb>=-0.0622757189683
<Himp>=0.0198970543877
<Hpot>=0
<n_d>=0.797968502578
<n_d_ud>=0.0790589005647
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 22/34 (scale 7.2966059202918268e-07) ***** 
[1] xi(21)=0.997971323294 zeta(22)=5.25266134047e-10
Stats: nr=32 min=1 max=235 sum3=42477450
Emax=32.2998185656 nrkept=500 nrkeptmult=1147 
Total energy=-2.95399627922  Egs=-1.08951986885
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.0622757177453
<Himp>=0.0198970590046
<Hpot>=0
<n_d>=0.797968687737
<n_d_ud>=0.0790589148978
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 23/34 (scale 3.6483029601459134e-07) ***** 
[1] xi(22)=0.997851949867 zeta(23)=7.4859593237e-11
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=32.9013652009 nrkept=500 nrkeptmult=1141 
Total energy=-2.95399664976  Egs=-1.0156537891
Kept: 500 out of 1868, ratio=0.268
<Hhyb>=-0.062275718231
<Himp>=0.0198970581457
<Hpot>=0
<n_d>=0.797968653291
<n_d_ud>=0.0790589059743
Memory used: 26 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 24/34 (scale 1.8241514800729567e-07) ***** 
[1] xi(23)=0.998130833908 zeta(24)=3.90471234635e-11
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.296490827 nrkept=500 nrkeptmult=1147 
Total energy=-2.95399684861  Egs=-1.09010590117
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.0622757179685
<Himp>=0.0198970589072
<Hpot>=0
<n_d>=0.79796868383
<n_d_ud>=0.0790589098117
Memory used: 26 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 25/34 (scale 9.1207574003647835e-08) ***** 
[1] xi(24)=0.997763090251 zeta(25)=-8.95454222747e-12
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=32.9494976079 nrkept=500 nrkeptmult=1141 
Total energy=-2.9539969412  Egs=-1.01514279025
Kept: 500 out of 1868, ratio=0.268
<Hhyb>=-0.0622757181044
<Himp>=0.0198970585852
<Hpot>=0
<n_d>=0.797968670916
<n_d_ud>=0.079058907585
Memory used: 26 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 26/34 (scale 4.5603787001823918e-08) ***** 
[1] xi(25)=0.99829123374 zeta(26)=2.44651171171e-12
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.2976539594 nrkept=500 nrkeptmult=1147 
Total energy=-2.95399699093  Egs=-1.09064250052
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.0622757180469
<Himp>=0.019897058737
<Hpot>=0
<n_d>=0.797968677006
<n_d_ud>=0.0790589084734
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 27/34 (scale 2.2801893500911959e-08) ***** 
[1] xi(26)=1.00090948803 zeta(27)=-2.04989465674e-12
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=32.9646106193 nrkept=500 nrkeptmult=1141 
Total energy=-2.95399701418  Egs=-1.01959246692
Kept: 500 out of 1868, ratio=0.268
<Hhyb>=-0.0622757180832
<Himp>=0.0198970586464
<Hpot>=0
<n_d>=0.79796867337
<n_d_ud>=0.0790589078955
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 28/34 (scale 1.1400946750455979e-08) ***** 
[1] xi(27)=1.00401286132 zeta(28)=-1.06333982247e-11
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.3158008129 nrkept=500 nrkeptmult=1147 
Total energy=-2.95399702669  Egs=-1.09688275928
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.0622757180648
<Himp>=0.0198970586934
<Hpot>=0
<n_d>=0.797968675257
<n_d_ud>=0.0790589081847
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 29/34 (scale 5.7004733752279897e-09) ***** 
[1] xi(28)=1.00264783397 zeta(29)=3.79831302211e-12
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=32.9973772922 nrkept=500 nrkeptmult=1141 
Total energy=-2.9539970325  Egs=-1.0193249626
Kept: 500 out of 1868, ratio=0.268
<Hhyb>=-0.062275718072
<Himp>=0.0198970586752
<Hpot>=0
<n_d>=0.797968674529
<n_d_ud>=0.0790589080709
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 30/34 (scale 2.8502366876139948e-09) ***** 
[1] xi(29)=1.00106375957 zeta(30)=2.09665716748e-12
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.3506674786 nrkept=500 nrkeptmult=1147 
Total energy=-2.95399703562  Egs=-1.09378755138
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.0622757180713
<Himp>=0.019897058677
<Hpot>=0
<n_d>=0.797968674599
<n_d_ud>=0.0790589080815
Memory used: 27 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 31/34 (scale 1.4251183438069974e-09) ***** 
[1] xi(30)=1.00141660786 zeta(31)=1.2497190214e-13
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=33.0541851107 nrkept=500 nrkeptmult=1141 
Total energy=-2.95399703707  Egs=-1.01849119216
Kept: 500 out of 1868, ratio=0.268
<Hhyb>=-0.0622757180732
<Himp>=0.0198970586721
<Hpot>=0
<n_d>=0.797968674404
<n_d_ud>=0.0790589080513
Memory used: 27 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 32/34 (scale 7.1255917190349871e-10) ***** 
[1] xi(31)=1.00141387339 zeta(32)=1.34061187197e-13
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.4381978613 nrkept=500 nrkeptmult=1147 
Total energy=-2.95399703785  Egs=-1.09432605878
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.0622757180723
<Himp>=0.0198970586744
<Hpot>=0
<n_d>=0.797968674495
<n_d_ud>=0.0790589080654
Memory used: 28 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 33/34 (scale 3.5627958595174936e-10) ***** 
[1] xi(32)=1.0015307532 zeta(33)=-6.9391076572e-14
Stats: nr=33 min=1 max=238 sum3=42268052
Safeguard: keep additional 2 states
Emax=33.0979627004 nrkept=502 nrkeptmult=1147 
Total energy=-2.95399703821  Egs=-1.01859321122
Kept: 502 out of 1868, ratio=0.269
<Hhyb>=-0.0622757180727
<Himp>=0.0198970586735
<Hpot>=0
<n_d>=0.797968674459
<n_d_ud>=0.0790589080598
Memory used: 28 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 34/34 (scale 1.7813979297587468e-10) ***** 
[1] xi(33)=1.00153006076 zeta(34)=8.38476031785e-15
Stats: nr=32 min=1 max=235 sum3=42821078
Emax=32.3875797211 nrkept=500 nrkeptmult=1147 
Total energy=-2.95399703841  Egs=-1.09433495289
Kept: 500 out of 1874, ratio=0.267
<Hhyb>=-0.0622757180725
<Himp>=0.0198970586741
<Hpot>=0
<n_d>=0.797968674482
<n_d_ud>=0.0790589080633
Memory used: 28 MB  Time elapsed: 8 s

** Iteration completed.


Total energy:             -2.953997038406395
TKW=0.000857469104549
Grand canonical Z=              75.22149714542778
Zft (last shell)=               73.90263638630256
Calculating density matrices.
Calculating full density matrices.

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QS
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	../model.m	 Variant: 		 Channels: 	1
>>  U=	0	 Gamma=	-999	 Delta=	0.024934636296339015	 t=	0.
>> 	 extraJkondo=	0.05
>>  Gamma^(1/2)=	0.05050305102594391
>>  Lambda=	4.	 BAND=	asymode
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d self_d n_d n_d_ud Himp Hhyb Hpot SdSk
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=8
Reading <||||> (e)
Reading <||Himp||> (s)
Reading <||Hhyb||> (s)
Reading <||Hpot||> (s)
Reading <||n_d||> (s)
Reading <||n_d_ud||> (s)
Reading <||A_d||> (d)
Reading <||self_d||> (d)
Reading <||||> (z)
Reading <||||> (T)
Number of invariant subspaces: 8
(-2 2) 1 states
0.0178591316151 
(-1 1) 2 states
0.00222301767681 0.0373707196932 
(-1 3) 2 states
0.0113273794983 0.0446039935961 
(0 2) 4 states
0 0.0261978645765 0.0339879879407 0.0675966598135 
(0 4) 1 states
0.0380722414794 
(1 1) 2 states
0.0182938663113 0.053557272925 
(1 3) 2 states
0.0274764948932 0.0607122800675 
(2 2) 1 states
0.0501165334814 
Number of states (multiplicity taken into account): 32

Computing the following operators:
s=[]
p=[]
g=[]
d=[A_d self_d ]
v=[]
t=[]

Computing the following spectra:
Spectrum: A_d-A_d FDM
Spectrum: self_d-A_d FDM

 ***** [DMNRG] Iteration 1/34 (scale 1.5302091698951845) ***** 
[1] xi(0)=0.0722429873565 zeta(1)=-0.030892148129
Stats: nr=13 min=1 max=9 sum3=2864
Emax=0.359128143583 nrkept=56 nrkeptmult=128 
<Hhyb>_fdm=-0.0622757180666
<Himp>_fdm=0.0198970586891
<Hpot>_fdm=0
<n_d>_fdm=0.797968675084
<n_d_ud>_fdm=0.0790589081567
Kept: 56 out of 56, ratio=1
Memory used: 33 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 2/34 (scale 0.76510458494759226) ***** 
[1] xi(1)=0.945553390664 zeta(2)=0.0133042597598
Stats: nr=19 min=1 max=35 sum3=99630
Emax=4.09551767314 nrkept=210 nrkeptmult=512 
Kept: 210 out of 210, ratio=1
Memory used: 34 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 3/34 (scale 0.38255229247379613) ***** 
[1] xi(2)=1.56565518157 zeta(3)=0.01201927351
Stats: nr=26 min=1 max=99 sum3=3786552
Emax=10.060076228 nrkept=741 nrkeptmult=1930 
Kept: 741 out of 792, ratio=0.936
Memory used: 31 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 4/34 (scale 0.19127614623689806) ***** 
[1] xi(3)=1.4007691289 zeta(4)=-0.028530767889
Stats: nr=34 min=1 max=334 sum3=122312182
Emax=10.0296900901 nrkept=927 nrkeptmult=2490 
Kept: 927 out of 2812, ratio=0.33
Memory used: 37 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 5/34 (scale 0.095638073118449032) ***** 
[1] xi(4)=3.31343662111 zeta(5)=0.0391970579955
Stats: nr=40 min=1 max=376 sum3=214773540
Emax=10.5642528426 nrkept=500 nrkeptmult=1295 
Kept: 500 out of 3528, ratio=0.142
Memory used: 41 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 6/34 (scale 0.047819036559224516) ***** 
[1] xi(5)=1.75389974348 zeta(6)=-0.0106056416123
Stats: nr=34 min=1 max=224 sum3=37036511
Emax=12.2628450576 nrkept=500 nrkeptmult=1288 
Kept: 500 out of 1895, ratio=0.264
Memory used: 38 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 7/34 (scale 0.023909518279612258) ***** 
[1] xi(6)=1.26895016937 zeta(7)=-0.00326319694007
Stats: nr=36 min=1 max=215 sum3=39814837
Emax=12.7226205001 nrkept=500 nrkeptmult=1282 
Kept: 500 out of 1897, ratio=0.264
Memory used: 40 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 8/34 (scale 0.011954759139806129) ***** 
[1] xi(7)=0.913497015629 zeta(8)=-0.00390871043989
Stats: nr=33 min=1 max=229 sum3=38648109
Emax=14.5602312891 nrkept=500 nrkeptmult=1271 
Kept: 500 out of 1893, ratio=0.264
Memory used: 41 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 9/34 (scale 0.0059773795699030645) ***** 
[1] xi(8)=0.808927275461 zeta(9)=0.00182465376407
Stats: nr=34 min=1 max=215 sum3=39307401
Emax=18.0453035361 nrkept=500 nrkeptmult=1244 
Kept: 500 out of 1893, ratio=0.264
Memory used: 43 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 10/34 (scale 0.0029886897849515323) ***** 
[1] xi(9)=1.20041543191 zeta(10)=0.000594493935587
Stats: nr=32 min=1 max=235 sum3=41372763
Emax=21.2700767316 nrkept=500 nrkeptmult=1210 
Kept: 500 out of 1887, ratio=0.265
Memory used: 44 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 11/34 (scale 0.0014943448924757661) ***** 
[1] xi(10)=0.908364805201 zeta(11)=0.000319611855861
Stats: nr=33 min=1 max=223 sum3=41933624
Emax=23.4351296394 nrkept=500 nrkeptmult=1230 
Kept: 500 out of 1880, ratio=0.266
Memory used: 45 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 12/34 (scale 0.00074717244623788307) ***** 
[1] xi(11)=1.13512203319 zeta(12)=0.000191365561854
Stats: nr=32 min=1 max=236 sum3=40376669
Emax=25.2412035581 nrkept=500 nrkeptmult=1184 
Kept: 500 out of 1883, ratio=0.266
Memory used: 46 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 13/34 (scale 0.00037358622311894153) ***** 
[1] xi(12)=1.01675603198 zeta(13)=4.96979788512e-05
Stats: nr=33 min=1 max=230 sum3=42597639
Emax=27.1603339957 nrkept=500 nrkeptmult=1179 
Kept: 500 out of 1875, ratio=0.267
Memory used: 47 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 14/34 (scale 0.00018679311155947077) ***** 
[1] xi(13)=1.03908256429 zeta(14)=2.22336899926e-05
Stats: nr=32 min=1 max=246 sum3=43888731
Emax=27.3370054041 nrkept=500 nrkeptmult=1170 
Kept: 500 out of 1875, ratio=0.267
Memory used: 49 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 15/34 (scale 9.3396555779735383e-05) ***** 
[1] xi(14)=1.0207234573 zeta(15)=5.70600324284e-06
Stats: nr=33 min=1 max=231 sum3=42628830
Emax=27.7614368862 nrkept=500 nrkeptmult=1181 
Kept: 500 out of 1872, ratio=0.267
Memory used: 49 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 16/34 (scale 4.6698277889867692e-05) ***** 
[1] xi(15)=1.00588427675 zeta(16)=1.45714507933e-06
Stats: nr=32 min=1 max=243 sum3=43287243
Emax=29.789170396 nrkept=500 nrkeptmult=1164 
Kept: 500 out of 1875, ratio=0.267
Memory used: 51 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 17/34 (scale 2.3349138944933846e-05) ***** 
[1] xi(16)=1.00028936037 zeta(17)=3.78789332065e-07
Stats: nr=33 min=1 max=233 sum3=42510856
Emax=29.9960263743 nrkept=500 nrkeptmult=1169 
Kept: 500 out of 1870, ratio=0.267
Memory used: 51 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 18/34 (scale 1.1674569472466923e-05) ***** 
[1] xi(17)=0.998602768318 zeta(18)=1.02369016562e-07
Stats: nr=33 min=1 max=239 sum3=42523332
Emax=31.9028978168 nrkept=500 nrkeptmult=1153 
Kept: 500 out of 1872, ratio=0.267
Memory used: 53 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 19/34 (scale 5.8372847362334614e-06) ***** 
[1] xi(18)=0.998101367625 zeta(19)=2.65449403172e-08
Stats: nr=33 min=1 max=238 sum3=42308259
Emax=32.2985494507 nrkept=500 nrkeptmult=1149 
Kept: 500 out of 1869, ratio=0.268
Memory used: 54 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 20/34 (scale 2.9186423681167307e-06) ***** 
[1] xi(19)=0.997979280712 zeta(20)=7.36509791801e-09
Stats: nr=33 min=1 max=235 sum3=42007819
Emax=32.3188274108 nrkept=500 nrkeptmult=1145 
Kept: 500 out of 1867, ratio=0.268
Memory used: 55 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 21/34 (scale 1.4593211840583654e-06) ***** 
[1] xi(20)=0.997909098324 zeta(21)=1.72424644448e-09
Stats: nr=33 min=1 max=237 sum3=41893189
Emax=32.7844257977 nrkept=500 nrkeptmult=1143 
Kept: 500 out of 1867, ratio=0.268
Memory used: 57 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 22/34 (scale 7.2966059202918268e-07) ***** 
[1] xi(21)=0.997971323294 zeta(22)=5.25266134047e-10
Stats: nr=32 min=1 max=235 sum3=42477450
Emax=32.2998185656 nrkept=500 nrkeptmult=1147 
Kept: 500 out of 1866, ratio=0.268
Memory used: 57 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 23/34 (scale 3.6483029601459134e-07) ***** 
[1] xi(22)=0.997851949867 zeta(23)=7.4859593237e-11
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=32.9013652009 nrkept=500 nrkeptmult=1141 
Kept: 500 out of 1868, ratio=0.268
Memory used: 59 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 24/34 (scale 1.8241514800729567e-07) ***** 
[1] xi(23)=0.998130833908 zeta(24)=3.90471234635e-11
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.296490827 nrkept=500 nrkeptmult=1147 
Kept: 500 out of 1866, ratio=0.268
Memory used: 60 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 25/34 (scale 9.1207574003647835e-08) ***** 
[1] xi(24)=0.997763090251 zeta(25)=-8.95454222747e-12
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=32.9494976079 nrkept=500 nrkeptmult=1141 
Kept: 500 out of 1868, ratio=0.268
Memory used: 61 MB  Time elapsed: 15 s

 ***** [DMNRG] Iteration 26/34 (scale 4.5603787001823918e-08) ***** 
[1] xi(25)=0.99829123374 zeta(26)=2.44651171171e-12
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.2976539594 nrkept=500 nrkeptmult=1147 
Kept: 500 out of 1866, ratio=0.268
Memory used: 62 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 27/34 (scale 2.2801893500911959e-08) ***** 
[1] xi(26)=1.00090948803 zeta(27)=-2.04989465674e-12
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=32.9646106193 nrkept=500 nrkeptmult=1141 
Kept: 500 out of 1868, ratio=0.268
Memory used: 63 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 28/34 (scale 1.1400946750455979e-08) ***** 
[1] xi(27)=1.00401286132 zeta(28)=-1.06333982247e-11
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.3158008129 nrkept=500 nrkeptmult=1147 
Kept: 500 out of 1866, ratio=0.268
Memory used: 64 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 29/34 (scale 5.7004733752279897e-09) ***** 
[1] xi(28)=1.00264783397 zeta(29)=3.79831302211e-12
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=32.9973772922 nrkept=500 nrkeptmult=1141 
Kept: 500 out of 1868, ratio=0.268
Memory used: 65 MB  Time elapsed: 16 s

 ***** [DMNRG] Iteration 30/34 (scale 2.8502366876139948e-09) ***** 
[1] xi(29)=1.00106375957 zeta(30)=2.09665716748e-12
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.3506674786 nrkept=500 nrkeptmult=1147 
Kept: 500 out of 1866, ratio=0.268
Memory used: 66 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 31/34 (scale 1.4251183438069974e-09) ***** 
[1] xi(30)=1.00141660786 zeta(31)=1.2497190214e-13
Stats: nr=33 min=1 max=238 sum3=42268052
Emax=33.0541851107 nrkept=500 nrkeptmult=1141 
Kept: 500 out of 1868, ratio=0.268
Memory used: 67 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 32/34 (scale 7.1255917190349871e-10) ***** 
[1] xi(31)=1.00141387339 zeta(32)=1.34061187197e-13
Stats: nr=32 min=1 max=235 sum3=42624090
Emax=32.4381978613 nrkept=500 nrkeptmult=1147 
Kept: 500 out of 1866, ratio=0.268
Memory used: 68 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 33/34 (scale 3.5627958595174936e-10) ***** 
[1] xi(32)=1.0015307532 zeta(33)=-6.9391076572e-14
Stats: nr=33 min=1 max=238 sum3=42268052
Safeguard: keep additional 2 states
Emax=33.0979627004 nrkept=502 nrkeptmult=1147 
Kept: 502 out of 1868, ratio=0.269
Memory used: 69 MB  Time elapsed: 18 s

 ***** [DMNRG] Iteration 34/34 (scale 1.7813979297587468e-10) ***** 
[1] xi(33)=1.00153006076 zeta(34)=8.38476031785e-15
Stats: nr=32 min=1 max=235 sum3=42821078
Emax=32.3875797211 nrkept=500 nrkeptmult=1147 
Kept: 500 out of 1874, ratio=0.267
Memory used: 70 MB  Time elapsed: 18 s

** Iteration completed.

Spectrum: A_d-A_d FDM
[A_d-A_d] pos=0.601015662912 neg=0.398984337088 sum= 1
mu1=0.0251731708805 mu2=0.00418964596767 mu3=0.000277912529574 mu4=6.2624831011e-05

Spectrum: self_d-A_d FDM
[self_d-A_d] pos=0.0248603221007 neg=7.43141955261e-05 sum= 0.0249346362963
mu1=0.00159347366818 mu2=0.000145039530508 mu3=1.38677446832e-05 mu4=2.15405611506e-06


FDM expectation values
<Hhyb>=-0.0622757180666
<Himp>=0.0198970586891
<Hpot>=0
<n_d>=0.797968675084
<n_d_ud>=0.0790589081567

Memory usage report
===================
  after diag: 71940 kB
after recalc: 74512 kB
 after split: 71940 kB
  after trim: 72508 kB
 after trunc: 71940 kB
      diag_h: 29660 kB

Peak usage: 74512 kB

Timing report
=============
         All: 18 s
          DM: 0 s
         FDM: 0 s
        diag: 5 s
    recalc d: 2 s
    recalc s: 0 s
      spec d: 6 s
       Other: 3 s
