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=100
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
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.2	 Gamma=	-999	 Delta=	0.05899294080491952	 t=	0.
>> 
>>  Gamma^(1/2)=	0.05036083229510181
>>  Lambda=	4.	 BAND=	asymode
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d self_d n_d n_d_ud Himp Hhyb Hpot
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=6
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.1473951459775392 -->            0.09632352810147869
zeta(0)=           0.03923313656881788
  xi(1)=            0.5528372796341922 -->             0.7225643271658869
zeta(1)=           0.04145351061127775
  xi(2)=            0.7202045751962041 -->              1.882630399465026
zeta(2)=           -0.1511016234934212
  xi(3)=            0.3095918501556386 -->              1.618559638754981
zeta(3)=             0.195509932094321
  xi(4)=            0.1862408979018947 -->              1.947351006029083
zeta(4)=          -0.09864354100818817
  xi(5)=           0.07760953809630297 -->              1.622984143567646
zeta(5)=          -0.09432154359545168
  xi(6)=           0.02747473618755493 -->              1.149112912533363
zeta(6)=          0.005259024222632166
  xi(7)=           0.01319713069675421 -->              1.103922759331162
zeta(7)=         0.0005376054958964868
  xi(8)=          0.006264885822113212 -->              1.048099045551295
zeta(8)=        -0.0009187637532516218
  xi(9)=          0.003041828101908727 -->              1.017779803452588
zeta(9)=         -0.000229430482493302
  xi(10)=          0.001502730780634466 -->              1.005611748801046
zeta(10)=        -5.427394531333414e-05
  xi(11)=         0.0007485269308603607 -->               1.00181281393512
zeta(11)=        -1.349101337174052e-05
  xi(12)=          0.000373795543925092 -->              1.000560301192059
zeta(12)=        -3.544365190064809e-06
  xi(13)=         0.0001868289350476049 -->              1.000191781633889
zeta(13)=        -9.612227610605994e-07
  xi(14)=         9.340443968853119e-05 -->              1.000084413271239
zeta(14)=        -2.559538716603784e-07
  xi(15)=         4.670074139706885e-05 -->              1.000052753705543
zeta(15)=        -6.945679241567835e-08
  xi(16)=         2.335016282282465e-05 -->              1.000043850777248
zeta(16)=        -1.792582915136509e-08
  xi(17)=          1.16750461633892e-05 -->              1.000040831563288
zeta(17)=        -5.133014792920165e-09
  xi(18)=          5.83752409299739e-06 -->              1.000041004812125
zeta(18)=        -1.245381741137066e-09
  xi(19)=         2.918755013521892e-06 -->              1.000038595138065
zeta(19)=        -4.232113890869454e-10
  xi(20)=         1.459382938623613e-06 -->              1.000042317322549
zeta(20)=         -9.28722570582165e-11
  xi(21)=          7.29684866270723e-07 -->              1.000033267853308
zeta(21)=        -4.010085706158504e-11
  xi(22)=         3.648480057717875e-07 -->              1.000048542452175
zeta(22)=         -1.39276919509823e-11
  xi(23)=         1.824174536153641e-07 -->              1.000012639345436
zeta(23)=         7.690531750273912e-13
  xi(24)=         9.121408790641992e-08 -->              1.000071418441322
zeta(24)=        -8.281606464554863e-12
  xi(25)=         4.560120643731925e-08 -->             0.9999434133725656
zeta(25)=         6.841760032082915e-12
  xi(26)=         2.280345917165788e-08 -->              1.000068664067125
zeta(26)=         -6.95063461633226e-12
  xi(27)=         1.139617590324569e-08 -->             0.9995815393831131
zeta(27)=         8.281616416276373e-12
  xi(28)=         5.700379316255894e-09 -->             0.9999834997962618
zeta(28)=         -4.07567719778798e-12
  xi(29)=         2.848922945154022e-09 -->             0.9995390760122898
zeta(29)=          3.60010097459949e-12
  xi(30)=         1.424909669071249e-09 -->             0.9998535737493974
zeta(30)=         1.556659106443538e-12
  xi(31)=         7.124537624648508e-10 -->             0.9998520692136117
zeta(31)=         2.231023720558989e-13
  xi(32)=         3.562578540628786e-10 -->             0.9999390032723523
zeta(32)=         1.001854303300132e-13
  xi(33)=         1.781289130401012e-10 -->             0.9999389247310122
zeta(33)=        -1.915119214755541e-14
  xi(34)=         8.906639767522776e-11 -->               0.99996071834764
zeta(34)=         6.265875226408689e-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: 6
(-2 1) 1 states
0.021328824436 
(-1 2) 2 states
0 0.0420780148386 
(0 1) 3 states
0.000124841655885 0.0522519818053 0.0752212906445 
(0 3) 1 states
0.0207491904026 
(1 2) 2 states
0.0275688473732 0.0787004422965 
(2 1) 1 states
0.0855200992671 
Number of states (multiplicity taken into account): 16

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.000828569261616
<Himp>=0.107621608959
<Hpot>=0.0482399310443
<n_d>=0.990520239267
<n_d_ud>=0.241199655221

 ***** [NRG] Iteration 1/34 (scale 1.5302091698951845) ***** 
[1] xi(0)=0.0963235281015 zeta(1)=0.0414535106113
Stats: nr=10 min=1 max=8 sum3=1055
Emax=0.502334330098 nrkept=35 nrkeptmult=64 
Total energy=-0.246208461644  Egs=-0.157952392628
Kept: 35 out of 35, ratio=1
<Hhyb>=-0.00165499758446
<Himp>=0.106279261991
<Hpot>=0.046521948326
<n_d>=0.98135996887
<n_d_ud>=0.23260974163
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/34 (scale 0.76510458494759226) ***** 
[1] xi(1)=0.722564327166 zeta(2)=-0.151101623493
Stats: nr=15 min=1 max=20 sum3=30066
Emax=3.36097198234 nrkept=126 nrkeptmult=256 
Total energy=-1.1644630165  Egs=-1.20016867357
Kept: 126 out of 126, ratio=1
<Hhyb>=-0.00329614426844
<Himp>=0.103685399937
<Hpot>=0.0432166371041
<n_d>=0.964010537275
<n_d_ud>=0.21608318552
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/34 (scale 0.38255229247379613) ***** 
[1] xi(2)=1.88263039947 zeta(3)=0.195509932094
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=10.0540362701 nrkept=433 nrkeptmult=970 
Total energy=-2.05329192041  Egs=-2.32341805655
Kept: 433 out of 462, ratio=0.937
<Hhyb>=-0.00649179938922
<Himp>=0.0988806145182
<Hpot>=0.0371491366631
<n_d>=0.93321791165
<n_d_ud>=0.185745683316
Memory used: 6 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/34 (scale 0.19127614623689806) ***** 
[1] xi(3)=1.61855963875 zeta(4)=-0.0986435410082
Stats: nr=28 min=1 max=195 sum3=29508117
Emax=10.0288136116 nrkept=507 nrkeptmult=1196 
Total energy=-2.39661989312  Egs=-1.79493355269
Kept: 507 out of 1611, ratio=0.315
<Hhyb>=-0.0123455437787
<Himp>=0.0908570996564
<Hpot>=0.0272125058995
<n_d>=0.886564580752
<n_d_ud>=0.136062529498
Memory used: 11 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/34 (scale 0.095638073118449032) ***** 
[1] xi(4)=1.94735100603 zeta(5)=-0.0943215435955
Stats: nr=33 min=1 max=243 sum3=42362374
Emax=10.0015020648 nrkept=315 nrkeptmult=710 
Total energy=-2.60082954935  Egs=-2.13523390398
Kept: 315 out of 1900, ratio=0.166
<Hhyb>=-0.021478498351
<Himp>=0.08049457485
<Hpot>=0.0149436471516
<n_d>=0.840076635043
<n_d_ud>=0.0747182357582
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 6/34 (scale 0.047819036559224516) ***** 
[1] xi(5)=1.62298414357 zeta(6)=0.00525902422263
Stats: nr=28 min=1 max=142 sum3=11075769
Emax=10.0027008327 nrkept=226 nrkeptmult=509 
Total energy=-2.68432305455  Egs=-1.74603068585
Kept: 226 out of 1173, ratio=0.193
<Hhyb>=-0.0323889939534
<Himp>=0.0728895330893
<Hpot>=0.00690552798942
<n_d>=0.829515589944
<n_d_ud>=0.0345276399471
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 7/34 (scale 0.023909518279612258) ***** 
[1] xi(6)=1.14911291253 zeta(7)=0.000537605495896
Stats: nr=28 min=1 max=123 sum3=4818734
Emax=10.0309581093 nrkept=170 nrkeptmult=388 
Total energy=-2.7063149334  Egs=-0.919795982248
Kept: 170 out of 842, ratio=0.202
<Hhyb>=-0.0421329441237
<Himp>=0.0710803212729
<Hpot>=0.00552960102215
<n_d>=0.84008169387
<n_d_ud>=0.0276480051107
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 8/34 (scale 0.011954759139806129) ***** 
[1] xi(7)=1.10392275933 zeta(8)=-0.000918763753252
Stats: nr=26 min=1 max=83 sum3=1974028
Emax=10.071389109 nrkept=106 nrkeptmult=228 
Total energy=-2.72512041983  Egs=-1.57305439645
Kept: 106 out of 634, ratio=0.167
<Hhyb>=-0.0504027776015
<Himp>=0.072165410374
<Hpot>=0.00619229471457
<n_d>=0.829781140333
<n_d_ud>=0.0309614735729
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 9/34 (scale 0.0059773795699030645) ***** 
[1] xi(8)=1.04809904555 zeta(9)=-0.000229430482493
Stats: nr=20 min=1 max=63 sum3=608589
Emax=10.0621705671 nrkept=101 nrkeptmult=225 
Total energy=-2.72923499876  Egs=-0.688358316602
Kept: 101 out of 393, ratio=0.257
<Hhyb>=-0.0560012729215
<Himp>=0.0735512388406
<Hpot>=0.00692497605023
<n_d>=0.81385346486
<n_d_ud>=0.0346248802512
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 10/34 (scale 0.0029886897849515323) ***** 
[1] xi(9)=1.01777980345 zeta(10)=-5.42739453133e-05
Stats: nr=24 min=1 max=52 sum3=475378
Emax=11.5482381535 nrkept=100 nrkeptmult=208 
Total energy=-2.73417772305  Egs=-1.65380974267
Kept: 100 out of 376, ratio=0.266
<Hhyb>=-0.0589591590818
<Himp>=0.0743931040334
<Hpot>=0.00739273200505
<n_d>=0.804730420064
<n_d_ud>=0.0369636600253
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 11/34 (scale 0.0014943448924757661) ***** 
[1] xi(10)=1.0056117488 zeta(11)=-1.34910133717e-05
Stats: nr=21 min=1 max=59 sum3=493749
Emax=11.8218116412 nrkept=100 nrkeptmult=219 
Total energy=-2.73512126714  Egs=-0.631409857786
Kept: 100 out of 369, ratio=0.271
<Hhyb>=-0.060003045812
<Himp>=0.0747004807357
<Hpot>=0.00756851398769
<n_d>=0.80152134528
<n_d_ud>=0.0378425699384
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 12/34 (scale 0.00074717244623788307) ***** 
[1] xi(11)=1.00181281394 zeta(12)=-3.54436519006e-06
Stats: nr=24 min=1 max=52 sum3=464795
Emax=13.0813325904 nrkept=100 nrkeptmult=206 
Total energy=-2.73633965849  Egs=-1.63066953119
Kept: 100 out of 371, ratio=0.27
<Hhyb>=-0.0603048221488
<Himp>=0.0747900363845
<Hpot>=0.00761980490482
<n_d>=0.800588219802
<n_d_ud>=0.0380990245241
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 13/34 (scale 0.00037358622311894153) ***** 
[1] xi(12)=1.00056030119 zeta(13)=-9.61222761061e-07
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.3064580664 nrkept=100 nrkeptmult=218 
Total energy=-2.73657170024  Egs=-0.621119663642
Kept: 100 out of 367, ratio=0.272
<Hhyb>=-0.0603816882873
<Himp>=0.0748128956454
<Hpot>=0.00763307490775
<n_d>=0.800354375723
<n_d_ud>=0.0381653745388
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 14/34 (scale 0.00018679311155947077) ***** 
[1] xi(13)=1.00019178163 zeta(14)=-2.5595387166e-07
Stats: nr=24 min=1 max=53 sum3=476423
Emax=13.410016781 nrkept=100 nrkeptmult=206 
Total energy=-2.73687545821  Egs=-1.62617332365
Kept: 100 out of 371, ratio=0.27
<Hhyb>=-0.0604015756619
<Himp>=0.0748188234637
<Hpot>=0.00763644700704
<n_d>=0.800292051844
<n_d_ud>=0.0381822350352
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 15/34 (scale 9.3396555779735383e-05) ***** 
[1] xi(14)=1.00008441327 zeta(15)=-6.94567924157e-08
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.5247401108 nrkept=100 nrkeptmult=216 
Total energy=-2.73693327702  Egs=-0.619067912176
Kept: 100 out of 367, ratio=0.272
<Hhyb>=-0.0604064054314
<Himp>=0.0748202567359
<Hpot>=0.00763729496491
<n_d>=0.800277778344
<n_d_ud>=0.0381864748246
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 16/34 (scale 4.6698277889867692e-05) ***** 
[1] xi(15)=1.00005275371 zeta(16)=-1.79258291514e-08
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4594218442 nrkept=100 nrkeptmult=206 
Total energy=-2.73700917591  Egs=-1.6253038577
Kept: 100 out of 371, ratio=0.27
<Hhyb>=-0.0604076538927
<Himp>=0.0748206302257
<Hpot>=0.00763750190589
<n_d>=0.800273716876
<n_d_ud>=0.0381875095295
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 17/34 (scale 2.3349138944933846e-05) ***** 
[1] xi(16)=1.00004385078 zeta(17)=-5.13301479292e-09
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.5933595346 nrkept=100 nrkeptmult=216 
Total energy=-2.73702361937  Egs=-0.618586226232
Kept: 100 out of 367, ratio=0.272
<Hhyb>=-0.0604079575825
<Himp>=0.0748207194945
<Hpot>=0.007637558534
<n_d>=0.800272920899
<n_d_ud>=0.03818779267
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 18/34 (scale 1.1674569472466923e-05) ***** 
[1] xi(17)=1.00004083156 zeta(18)=-1.24538174114e-09
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4712917096 nrkept=100 nrkeptmult=206 
Total energy=-2.73704259222  Egs=-1.62514394969
Kept: 100 out of 371, ratio=0.27
<Hhyb>=-0.0604080348017
<Himp>=0.0748207429244
<Hpot>=0.00763757004346
<n_d>=0.800272630205
<n_d_ud>=0.0381878502173
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 19/34 (scale 5.8372847362334614e-06) ***** 
[1] xi(18)=1.00004100481 zeta(19)=-4.23211389087e-10
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.6106847562 nrkept=100 nrkeptmult=216 
Total energy=-2.73704620233  Egs=-0.618456862438
Kept: 100 out of 367, ratio=0.272
<Hhyb>=-0.0604080543366
<Himp>=0.0748207484612
<Hpot>=0.00763757449419
<n_d>=0.800272603721
<n_d_ud>=0.0381878724709
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 20/34 (scale 2.9186423681167307e-06) ***** 
[1] xi(19)=1.00003859514 zeta(20)=-9.28722570582e-11
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4742510372 nrkept=100 nrkeptmult=206 
Total energy=-2.7370509455  Egs=-1.62512710842
Kept: 100 out of 371, ratio=0.27
<Hhyb>=-0.0604080589369
<Himp>=0.0748207499401
<Hpot>=0.00763757485257
<n_d>=0.800272576396
<n_d_ud>=0.0381878742629
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 21/34 (scale 1.4593211840583654e-06) ***** 
[1] xi(20)=1.00004231732 zeta(21)=-4.01008570616e-11
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.6149478231 nrkept=100 nrkeptmult=216 
Total energy=-2.73705184797  Egs=-0.618419084703
Kept: 100 out of 367, ratio=0.272
<Hhyb>=-0.0604080603001
<Himp>=0.0748207502768
<Hpot>=0.00763757535865
<n_d>=0.800272580526
<n_d_ud>=0.0381878767932
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 22/34 (scale 7.2966059202918268e-07) ***** 
[1] xi(21)=1.00003326785 zeta(22)=-1.3927691951e-11
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4750004527 nrkept=100 nrkeptmult=206 
Total energy=-2.73705303376  Egs=-1.625126492
Kept: 100 out of 371, ratio=0.27
<Hhyb>=-0.0604080605307
<Himp>=0.0748207503729
<Hpot>=0.0076375752902
<n_d>=0.800272576514
<n_d_ud>=0.038187876451
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 23/34 (scale 3.6483029601459134e-07) ***** 
[1] xi(22)=1.00004854245 zeta(23)=7.69053175027e-13
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.6159357835 nrkept=100 nrkeptmult=216 
Total energy=-2.73705325938  Egs=-0.618424484128
Kept: 100 out of 367, ratio=0.272
<Hhyb>=-0.0604080606511
<Himp>=0.0748207503917
<Hpot>=0.00763757537814
<n_d>=0.800272578201
<n_d_ud>=0.0381878768907
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 24/34 (scale 1.8241514800729567e-07) ***** 
[1] xi(23)=1.00001263935 zeta(24)=-8.28160646455e-12
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.475155117 nrkept=100 nrkeptmult=206 
Total energy=-2.73705355582  Egs=-1.6250760839
Kept: 100 out of 371, ratio=0.27
<Hhyb>=-0.0604080606509
<Himp>=0.0748207503986
<Hpot>=0.00763757535051
<n_d>=0.800272577358
<n_d_ud>=0.0381878767526
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 25/34 (scale 9.1207574003647835e-08) ***** 
[1] xi(24)=1.00007141844 zeta(25)=6.84176003208e-12
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.6161218977 nrkept=100 nrkeptmult=216 
Total energy=-2.73705361223  Egs=-0.618514416812
Kept: 100 out of 367, ratio=0.272
<Hhyb>=-0.0604080606668
<Himp>=0.0748207503992
<Hpot>=0.00763757536943
<n_d>=0.800272577803
<n_d_ud>=0.0381878768472
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 26/34 (scale 4.5603787001823918e-08) ***** 
[1] xi(25)=0.999943413373 zeta(26)=-6.95063461633e-12
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4750020434 nrkept=100 nrkeptmult=206 
Total energy=-2.73705368633  Egs=-1.62483182364
Kept: 100 out of 371, ratio=0.27
<Hhyb>=-0.0604080606628
<Himp>=0.0748207503999
<Hpot>=0.00763757536124
<n_d>=0.800272577587
<n_d_ud>=0.0381878768062
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 27/34 (scale 2.2801893500911959e-08) ***** 
[1] xi(26)=1.00006866407 zeta(27)=8.28161641628e-12
Stats: nr=21 min=1 max=57 sum3=462847
Safeguard: keep additional 1 states
Emax=12.6161210504 nrkept=101 nrkeptmult=219 
Total energy=-2.73705370044  Egs=-0.618816287186
Kept: 101 out of 367, ratio=0.275
<Hhyb>=-0.0604080606655
<Himp>=0.0748207503999
<Hpot>=0.00763757536521
<n_d>=0.800272577685
<n_d_ud>=0.0381878768261
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 28/34 (scale 1.1400946750455979e-08) ***** 
[1] xi(27)=0.999581539383 zeta(28)=-4.07567719779e-12
Stats: nr=24 min=1 max=54 sum3=516057
Emax=13.4737417373 nrkept=100 nrkeptmult=206 
Total energy=-2.73705371895  Egs=-1.62366137644
Kept: 100 out of 375, ratio=0.267
<Hhyb>=-0.0604080606639
<Himp>=0.0748207504
<Hpot>=0.00763757536251
<n_d>=0.800272577616
<n_d_ud>=0.0381878768125
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 29/34 (scale 5.7004733752279897e-09) ***** 
[1] xi(28)=0.999983499796 zeta(29)=3.6001009746e-12
Stats: nr=21 min=1 max=57 sum3=462847
Safeguard: keep additional 4 states
Emax=12.6152533817 nrkept=104 nrkeptmult=224 
Total energy=-2.73705372248  Egs=-0.619412040924
Kept: 104 out of 367, ratio=0.283
<Hhyb>=-0.0604080606642
<Himp>=0.0748207504
<Hpot>=0.00763757536297
<n_d>=0.800272577628
<n_d_ud>=0.0381878768149
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 30/34 (scale 2.8502366876139948e-09) ***** 
[1] xi(29)=0.999539076012 zeta(30)=1.55665910644e-12
Stats: nr=24 min=1 max=54 sum3=561463
Emax=13.4700925795 nrkept=100 nrkeptmult=206 
Total energy=-2.73705372711  Egs=-1.62248290496
Kept: 100 out of 385, ratio=0.26
<Hhyb>=-0.0604080606635
<Himp>=0.0748207504
<Hpot>=0.00763757536191
<n_d>=0.800272577601
<n_d_ud>=0.0381878768096
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 31/34 (scale 1.4251183438069974e-09) ***** 
[1] xi(30)=0.999853573749 zeta(31)=2.23102372056e-13
Stats: nr=21 min=1 max=57 sum3=462847
Safeguard: keep additional 4 states
Emax=12.6139419679 nrkept=104 nrkeptmult=224 
Total energy=-2.73705372799  Egs=-0.619609665657
Kept: 104 out of 367, ratio=0.283
<Hhyb>=-0.0604080606635
<Himp>=0.0748207504
<Hpot>=0.00763757536189
<n_d>=0.8002725776
<n_d_ud>=0.0381878768095
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 32/34 (scale 7.1255917190349871e-10) ***** 
[1] xi(31)=0.999852069214 zeta(32)=1.0018543033e-13
Stats: nr=24 min=1 max=54 sum3=561463
Emax=13.4665203455 nrkept=100 nrkeptmult=206 
Total energy=-2.73705372915  Egs=-1.62270602914
Kept: 100 out of 385, ratio=0.26
<Hhyb>=-0.0604080606634
<Himp>=0.0748207504
<Hpot>=0.00763757536167
<n_d>=0.800272577595
<n_d_ud>=0.0381878768084
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 33/34 (scale 3.5627958595174936e-10) ***** 
[1] xi(32)=0.999939003272 zeta(33)=-1.91511921476e-14
Stats: nr=21 min=1 max=57 sum3=462847
Safeguard: keep additional 4 states
Emax=12.6205627834 nrkept=104 nrkeptmult=224 
Total energy=-2.73705372937  Egs=-0.619563219322
Kept: 104 out of 367, ratio=0.283
<Hhyb>=-0.0604080606634
<Himp>=0.0748207504
<Hpot>=0.00763757536173
<n_d>=0.800272577596
<n_d_ud>=0.0381878768087
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 34/34 (scale 1.7813979297587468e-10) ***** 
[1] xi(33)=0.999938924731 zeta(34)=6.26587522641e-15
Stats: nr=24 min=1 max=54 sum3=561463
Emax=13.4691958153 nrkept=100 nrkeptmult=206 
Total energy=-2.73705372966  Egs=-1.622949931
Kept: 100 out of 385, ratio=0.26
<Hhyb>=-0.0604080606634
<Himp>=0.0748207504
<Hpot>=0.0076375753617
<n_d>=0.800272577596
<n_d_ud>=0.0381878768085
Memory used: 16 MB  Time elapsed: 1 s

** Iteration completed.


Total energy:             -2.737053729657401
TKW=0.00742755509725
Grand canonical Z=               73.7148666608505
Zft (last shell)=                55.4642542448135
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.2	 Gamma=	-999	 Delta=	0.05899294080491952	 t=	0.
>> 
>>  Gamma^(1/2)=	0.05036083229510181
>>  Lambda=	4.	 BAND=	asymode
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d self_d n_d n_d_ud Himp Hhyb Hpot
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=6
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: 6
(-2 1) 1 states
0.021328824436 
(-1 2) 2 states
0 0.0420780148386 
(0 1) 3 states
0.000124841655885 0.0522519818053 0.0752212906445 
(0 3) 1 states
0.0207491904026 
(1 2) 2 states
0.0275688473732 0.0787004422965 
(2 1) 1 states
0.0855200992671 
Number of states (multiplicity taken into account): 16

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.0963235281015 zeta(1)=0.0414535106113
Stats: nr=10 min=1 max=8 sum3=1055
Emax=0.502334330098 nrkept=35 nrkeptmult=64 
<Hhyb>_fdm=-0.0604080606628
<Himp>_fdm=0.0748207504001
<Hpot>_fdm=0.00763757536072
<n_d>_fdm=0.800272577571
<n_d_ud>_fdm=0.0381878768036
Kept: 35 out of 35, ratio=1
Memory used: 18 MB  Time elapsed: 1 s

 ***** [DMNRG] Iteration 2/34 (scale 0.76510458494759226) ***** 
[1] xi(1)=0.722564327166 zeta(2)=-0.151101623493
Stats: nr=15 min=1 max=20 sum3=30066
Emax=3.36097198234 nrkept=126 nrkeptmult=256 
Kept: 126 out of 126, ratio=1
Memory used: 19 MB  Time elapsed: 1 s

 ***** [DMNRG] Iteration 3/34 (scale 0.38255229247379613) ***** 
[1] xi(2)=1.88263039947 zeta(3)=0.195509932094
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=10.0540362701 nrkept=433 nrkeptmult=970 
Kept: 433 out of 462, ratio=0.937
Memory used: 19 MB  Time elapsed: 1 s

 ***** [DMNRG] Iteration 4/34 (scale 0.19127614623689806) ***** 
[1] xi(3)=1.61855963875 zeta(4)=-0.0986435410082
Stats: nr=28 min=1 max=195 sum3=29508117
Emax=10.0288136116 nrkept=507 nrkeptmult=1196 
Kept: 507 out of 1611, ratio=0.315
Memory used: 24 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 5/34 (scale 0.095638073118449032) ***** 
[1] xi(4)=1.94735100603 zeta(5)=-0.0943215435955
Stats: nr=33 min=1 max=243 sum3=42362374
Emax=10.0015020648 nrkept=315 nrkeptmult=710 
Kept: 315 out of 1900, ratio=0.166
Memory used: 27 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 6/34 (scale 0.047819036559224516) ***** 
[1] xi(5)=1.62298414357 zeta(6)=0.00525902422263
Stats: nr=28 min=1 max=142 sum3=11075769
Emax=10.0027008327 nrkept=226 nrkeptmult=509 
Kept: 226 out of 1173, ratio=0.193
Memory used: 28 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 7/34 (scale 0.023909518279612258) ***** 
[1] xi(6)=1.14911291253 zeta(7)=0.000537605495896
Stats: nr=28 min=1 max=123 sum3=4818734
Emax=10.0309581093 nrkept=170 nrkeptmult=388 
Kept: 170 out of 842, ratio=0.202
Memory used: 29 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 8/34 (scale 0.011954759139806129) ***** 
[1] xi(7)=1.10392275933 zeta(8)=-0.000918763753252
Stats: nr=26 min=1 max=83 sum3=1974028
Emax=10.071389109 nrkept=106 nrkeptmult=228 
Kept: 106 out of 634, ratio=0.167
Memory used: 29 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 9/34 (scale 0.0059773795699030645) ***** 
[1] xi(8)=1.04809904555 zeta(9)=-0.000229430482493
Stats: nr=20 min=1 max=63 sum3=608589
Emax=10.0621705671 nrkept=101 nrkeptmult=225 
Kept: 101 out of 393, ratio=0.257
Memory used: 29 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 10/34 (scale 0.0029886897849515323) ***** 
[1] xi(9)=1.01777980345 zeta(10)=-5.42739453133e-05
Stats: nr=24 min=1 max=52 sum3=475378
Emax=11.5482381535 nrkept=100 nrkeptmult=208 
Kept: 100 out of 376, ratio=0.266
Memory used: 30 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 11/34 (scale 0.0014943448924757661) ***** 
[1] xi(10)=1.0056117488 zeta(11)=-1.34910133717e-05
Stats: nr=21 min=1 max=59 sum3=493749
Emax=11.8218116412 nrkept=100 nrkeptmult=219 
Kept: 100 out of 369, ratio=0.271
Memory used: 31 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 12/34 (scale 0.00074717244623788307) ***** 
[1] xi(11)=1.00181281394 zeta(12)=-3.54436519006e-06
Stats: nr=24 min=1 max=52 sum3=464795
Emax=13.0813325904 nrkept=100 nrkeptmult=206 
Kept: 100 out of 371, ratio=0.27
Memory used: 32 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 13/34 (scale 0.00037358622311894153) ***** 
[1] xi(12)=1.00056030119 zeta(13)=-9.61222761061e-07
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.3064580664 nrkept=100 nrkeptmult=218 
Kept: 100 out of 367, ratio=0.272
Memory used: 32 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 14/34 (scale 0.00018679311155947077) ***** 
[1] xi(13)=1.00019178163 zeta(14)=-2.5595387166e-07
Stats: nr=24 min=1 max=53 sum3=476423
Emax=13.410016781 nrkept=100 nrkeptmult=206 
Kept: 100 out of 371, ratio=0.27
Memory used: 33 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 15/34 (scale 9.3396555779735383e-05) ***** 
[1] xi(14)=1.00008441327 zeta(15)=-6.94567924157e-08
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.5247401108 nrkept=100 nrkeptmult=216 
Kept: 100 out of 367, ratio=0.272
Memory used: 34 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 16/34 (scale 4.6698277889867692e-05) ***** 
[1] xi(15)=1.00005275371 zeta(16)=-1.79258291514e-08
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4594218442 nrkept=100 nrkeptmult=206 
Kept: 100 out of 371, ratio=0.27
Memory used: 35 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 17/34 (scale 2.3349138944933846e-05) ***** 
[1] xi(16)=1.00004385078 zeta(17)=-5.13301479292e-09
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.5933595346 nrkept=100 nrkeptmult=216 
Kept: 100 out of 367, ratio=0.272
Memory used: 36 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 18/34 (scale 1.1674569472466923e-05) ***** 
[1] xi(17)=1.00004083156 zeta(18)=-1.24538174114e-09
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4712917096 nrkept=100 nrkeptmult=206 
Kept: 100 out of 371, ratio=0.27
Memory used: 36 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 19/34 (scale 5.8372847362334614e-06) ***** 
[1] xi(18)=1.00004100481 zeta(19)=-4.23211389087e-10
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.6106847562 nrkept=100 nrkeptmult=216 
Kept: 100 out of 367, ratio=0.272
Memory used: 37 MB  Time elapsed: 2 s

 ***** [DMNRG] Iteration 20/34 (scale 2.9186423681167307e-06) ***** 
[1] xi(19)=1.00003859514 zeta(20)=-9.28722570582e-11
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4742510372 nrkept=100 nrkeptmult=206 
Kept: 100 out of 371, ratio=0.27
Memory used: 38 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 21/34 (scale 1.4593211840583654e-06) ***** 
[1] xi(20)=1.00004231732 zeta(21)=-4.01008570616e-11
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.6149478231 nrkept=100 nrkeptmult=216 
Kept: 100 out of 367, ratio=0.272
Memory used: 39 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 22/34 (scale 7.2966059202918268e-07) ***** 
[1] xi(21)=1.00003326785 zeta(22)=-1.3927691951e-11
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4750004527 nrkept=100 nrkeptmult=206 
Kept: 100 out of 371, ratio=0.27
Memory used: 40 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 23/34 (scale 3.6483029601459134e-07) ***** 
[1] xi(22)=1.00004854245 zeta(23)=7.69053175027e-13
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.6159357835 nrkept=100 nrkeptmult=216 
Kept: 100 out of 367, ratio=0.272
Memory used: 41 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 24/34 (scale 1.8241514800729567e-07) ***** 
[1] xi(23)=1.00001263935 zeta(24)=-8.28160646455e-12
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.475155117 nrkept=100 nrkeptmult=206 
Kept: 100 out of 371, ratio=0.27
Memory used: 41 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 25/34 (scale 9.1207574003647835e-08) ***** 
[1] xi(24)=1.00007141844 zeta(25)=6.84176003208e-12
Stats: nr=21 min=1 max=57 sum3=462847
Emax=12.6161218977 nrkept=100 nrkeptmult=216 
Kept: 100 out of 367, ratio=0.272
Memory used: 42 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 26/34 (scale 4.5603787001823918e-08) ***** 
[1] xi(25)=0.999943413373 zeta(26)=-6.95063461633e-12
Stats: nr=24 min=1 max=54 sum3=501233
Emax=13.4750020434 nrkept=100 nrkeptmult=206 
Kept: 100 out of 371, ratio=0.27
Memory used: 43 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 27/34 (scale 2.2801893500911959e-08) ***** 
[1] xi(26)=1.00006866407 zeta(27)=8.28161641628e-12
Stats: nr=21 min=1 max=57 sum3=462847
Safeguard: keep additional 1 states
Emax=12.6161210504 nrkept=101 nrkeptmult=219 
Kept: 101 out of 367, ratio=0.275
Memory used: 44 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 28/34 (scale 1.1400946750455979e-08) ***** 
[1] xi(27)=0.999581539383 zeta(28)=-4.07567719779e-12
Stats: nr=24 min=1 max=54 sum3=516057
Emax=13.4737417373 nrkept=100 nrkeptmult=206 
Kept: 100 out of 375, ratio=0.267
Memory used: 45 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 29/34 (scale 5.7004733752279897e-09) ***** 
[1] xi(28)=0.999983499796 zeta(29)=3.6001009746e-12
Stats: nr=21 min=1 max=57 sum3=462847
Safeguard: keep additional 4 states
Emax=12.6152533817 nrkept=104 nrkeptmult=224 
Kept: 104 out of 367, ratio=0.283
Memory used: 45 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 30/34 (scale 2.8502366876139948e-09) ***** 
[1] xi(29)=0.999539076012 zeta(30)=1.55665910644e-12
Stats: nr=24 min=1 max=54 sum3=561463
Emax=13.4700925795 nrkept=100 nrkeptmult=206 
Kept: 100 out of 385, ratio=0.26
Memory used: 46 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 31/34 (scale 1.4251183438069974e-09) ***** 
[1] xi(30)=0.999853573749 zeta(31)=2.23102372056e-13
Stats: nr=21 min=1 max=57 sum3=462847
Safeguard: keep additional 4 states
Emax=12.6139419679 nrkept=104 nrkeptmult=224 
Kept: 104 out of 367, ratio=0.283
Memory used: 47 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 32/34 (scale 7.1255917190349871e-10) ***** 
[1] xi(31)=0.999852069214 zeta(32)=1.0018543033e-13
Stats: nr=24 min=1 max=54 sum3=561463
Emax=13.4665203455 nrkept=100 nrkeptmult=206 
Kept: 100 out of 385, ratio=0.26
Memory used: 48 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 33/34 (scale 3.5627958595174936e-10) ***** 
[1] xi(32)=0.999939003272 zeta(33)=-1.91511921476e-14
Stats: nr=21 min=1 max=57 sum3=462847
Safeguard: keep additional 4 states
Emax=12.6205627834 nrkept=104 nrkeptmult=224 
Kept: 104 out of 367, ratio=0.283
Memory used: 49 MB  Time elapsed: 3 s

 ***** [DMNRG] Iteration 34/34 (scale 1.7813979297587468e-10) ***** 
[1] xi(33)=0.999938924731 zeta(34)=6.26587522641e-15
Stats: nr=24 min=1 max=54 sum3=561463
Emax=13.4691958153 nrkept=100 nrkeptmult=206 
Kept: 100 out of 385, ratio=0.26
Memory used: 50 MB  Time elapsed: 3 s

** Iteration completed.

Spectrum: A_d-A_d FDM
[A_d-A_d] pos=0.599863711154 neg=0.400136288846 sum= 1
mu1=0.0391963283322 mu2=0.0137734615837 mu3=0.00216065639962 mu4=0.000606706447566

Spectrum: self_d-A_d FDM
[self_d-A_d] pos=0.0477910356806 neg=-0.00877083711861 sum= 0.039020198562
mu1=0.0111729585981 mu2=0.00195133917083 mu3=0.000505988736441 mu4=0.000110486756641


FDM expectation values
<Hhyb>=-0.0604080606628
<Himp>=0.0748207504001
<Hpot>=0.00763757536072
<n_d>=0.800272577571
<n_d_ud>=0.0381878768036

Memory usage report
===================
  after diag: 50484 kB
after recalc: 51312 kB
 after split: 50484 kB
  after trim: 51312 kB
 after trunc: 50484 kB
      diag_h: 16524 kB

Peak usage: 51312 kB

Timing report
=============
         All: 3 s
          DM: 0 s
         FDM: 0 s
        diag: 0 s
    recalc d: 0 s
    recalc f: 0 s
    recalc s: 0 s
      spec d: 1 s
     tridiag: 0 s
       Other: 0 s
