NRG Ljubljana 2.3.20 - (c) Rok Zitko <rok.zitko@ijs.si>
$Id: nrg.cc,v 1.31 2012/10/26 10:18:22 rokzitko Exp rokzitko $
Timestamp: Fri May 31 11:54:49 2013
Compiled on Jun  5 2013 at 18:59:39

No MPI: single node calculation.

workdir=.
Tmin=1e-10 -> Nmax=42
TEMPERATURE T=1.3223879413e-10
Nmax=42
Tmin=1e-10
Ninit=0
Lambda=3
strategy=kept
lastall=0
lastalloverride=0
z=0.75
discretization=Z
band=flat
betabar=1
T=1.3223879413e-10
T_ratio=1
Tmin_ratio=-999
keep=5000
keepenergy=10
keepmin=0
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=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=2.5
broaden_min=-99
broaden_min_ratio=3
broaden_ratio=1.05
dumpenergies=0
saveall=0
dumpannotated=0
dumpscaled=1
dumpabs=0
dumpprecision=8
dumpgroups=1
grouptol=1e-06
dumpdiagonal=0
ops=
specs=
specd=
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=0
T0opt=-1
dmnrg=0
cfs=0
fdm=0
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: QS
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 		 Channels: 	1
>>  U=	0.01	 Gamma=	0.0006	 Delta=	0	 t=	0.
>> 
>>  Gamma^(1/2)=	0.019544100476116797
>>  Lambda=	3.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.75
>>  ops=	
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=6
Reading <||||> (e)
Reading <||||> (z)

# $Id: nrg.cc,v 1.31 2012/10/26 10:18:22 rokzitko Exp rokzitko $ QS
# disk=Z band=flat Lambda=3 z=0.75
# keep=5000 keepenergy=10 Nmax=42 channels=1 dots=1 betabar=1

Number of invariant subspaces: 6
(-2 1) 1 states
0.0301229930704 
(-1 2) 2 states
0.0140715823652 0.0425597660672 
(0 1) 3 states
0 0.0301229930704 0.0566313484323 
(0 3) 1 states
0.0265083553619 
(1 2) 2 states
0.0140715823652 0.0425597660672 
(2 1) 1 states
0.0301229930704 
Number of states (multiplicity taken into account): 16

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

Computing the following spectra:

Before NRG iteration (N=-1)

 ***** [NRG] Iteration 1/42 (scale 0.79862812782237003) ***** 
[1] xi(0)=0.737093468838 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=2.9562454861 nrkept=35 nrkeptmult=64 
Total energy=-1.17797039898  Egs=-1.42907855666
Kept: 35 out of 35, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/42 (scale 0.46108816458065222) ***** 
[1] xi(1)=0.971787846328 zeta(2)=0
Stats: nr=15 min=1 max=20 sum3=30066
Emax=6.52249718023 nrkept=126 nrkeptmult=256 
Total energy=-1.50122312666  Egs=-0.701064899325
Kept: 126 out of 126, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/42 (scale 0.26620937594079003) ***** 
[1] xi(2)=1.03321268925 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Safeguard: keep additional 3 states
Emax=10.1694928073 nrkept=355 nrkeptmult=800 
Total energy=-1.94709286135  Egs=-1.67488366297
Kept: 355 out of 462, ratio=0.768
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/42 (scale 0.1536960548602174) ***** 
[1] xi(3)=1.00725946845 zeta(4)=0
Stats: nr=28 min=1 max=157 sum3=15891426
Emax=10.009984317 nrkept=433 nrkeptmult=1009 
Total energy=-2.0552323377  Egs=-0.703593052185
Kept: 433 out of 1320, ratio=0.328
Memory used: 7 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/42 (scale 0.088736458646930014) ***** 
[1] xi(4)=1.00017398697 zeta(5)=0
Stats: nr=33 min=1 max=213 sum3=27438589
Safeguard: keep additional 1 states
Emax=10.1851436223 nrkept=389 nrkeptmult=874 
Total energy=-2.19724245547  Egs=-1.60035818346
Kept: 389 out of 1621, ratio=0.24
Memory used: 10 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 6/42 (scale 0.051232018286739137) ***** 
[1] xi(5)=0.999821553706 zeta(6)=0
Stats: nr=28 min=1 max=173 sum3=21206352
Safeguard: keep additional 7 states
Emax=10.3011233798 nrkept=420 nrkeptmult=976 
Total energy=-2.23492538393  Egs=-0.73553472464
Kept: 420 out of 1446, ratio=0.29
Memory used: 11 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 7/42 (scale 0.029578819548976672) ***** 
[1] xi(6)=0.999921934048 zeta(7)=0
Stats: nr=33 min=1 max=205 sum3=24752824
Safeguard: keep additional 1 states
Emax=10.6288070329 nrkept=389 nrkeptmult=874 
Total energy=-2.28065283331  Egs=-1.54595247807
Kept: 389 out of 1570, ratio=0.248
Memory used: 12 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 8/42 (scale 0.017077339428913048) ***** 
[1] xi(7)=0.999972235408 zeta(8)=0
Stats: nr=28 min=1 max=173 sum3=21206352
Emax=10.3372914404 nrkept=413 nrkeptmult=963 
Total energy=-2.29385126403  Egs=-0.772862235162
Kept: 413 out of 1446, ratio=0.286
Memory used: 13 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 9/42 (scale 0.0098596065163255569) ***** 
[1] xi(8)=0.999990560771 zeta(9)=0
Stats: nr=33 min=1 max=202 sum3=23580328
Safeguard: keep additional 1 states
Emax=10.1922848206 nrkept=357 nrkeptmult=802 
Total energy=-2.30855364675  Egs=-1.49117337444
Kept: 357 out of 1546, ratio=0.231
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 10/42 (scale 0.0056924464763043487) ***** 
[1] xi(9)=0.999996833398 zeta(10)=0
Stats: nr=28 min=1 max=162 sum3=17013396
Emax=10.4666373137 nrkept=413 nrkeptmult=963 
Total energy=-2.31312301359  Egs=-0.802707036411
Kept: 413 out of 1332, ratio=0.31
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 11/42 (scale 0.0032865355054418526) ***** 
[1] xi(10)=0.999998942233 zeta(11)=0
Stats: nr=33 min=1 max=202 sum3=23580328
Emax=10.0428071846 nrkept=322 nrkeptmult=730 
Total energy=-2.31787325015  Egs=-1.44536292188
Kept: 322 out of 1546, ratio=0.208
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 12/42 (scale 0.0018974821587681161) ***** 
[1] xi(11)=0.999999647163 zeta(12)=0
Stats: nr=28 min=1 max=144 sum3=12231774
Safeguard: keep additional 1 states
Emax=10.1442797246 nrkept=320 nrkeptmult=740 
Total energy=-2.31944244689  Egs=-0.826988927604
Kept: 320 out of 1200, ratio=0.267
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 13/42 (scale 0.001095511835147284) ***** 
[1] xi(12)=0.99999988236 zeta(13)=0
Stats: nr=31 min=1 max=161 sum3=11616834
Safeguard: keep additional 1 states
Emax=10.0709709947 nrkept=245 nrkeptmult=570 
Total energy=-2.32097908342  Egs=-1.40266538702
Kept: 245 out of 1194, ratio=0.205
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 14/42 (scale 0.00063249405292270544) ***** 
[1] xi(13)=0.999999960784 zeta(14)=0
Stats: nr=26 min=1 max=111 sum3=5696496
Safeguard: keep additional 2 states
Emax=10.0043240756 nrkept=219 nrkeptmult=531 
Total energy=-2.32152152447  Egs=-0.857622373791
Kept: 219 out of 918, ratio=0.239
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 15/42 (scale 0.00036517061171576138) ***** 
[1] xi(14)=0.999999986928 zeta(15)=0
Stats: nr=31 min=1 max=112 sum3=3879366
Safeguard: keep additional 1 states
Emax=10.1239244114 nrkept=185 nrkeptmult=438 
Total energy=-2.32201480376  Egs=-1.35081868697
Kept: 185 out of 822, ratio=0.225
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 16/42 (scale 0.00021083135097423514) ***** 
[1] xi(15)=0.999999995642 zeta(16)=0
Stats: nr=26 min=1 max=86 sum3=2628288
Safeguard: keep additional 2 states
Emax=10.3860159004 nrkept=213 nrkeptmult=523 
Total energy=-2.32220503455  Egs=-0.902288946219
Kept: 213 out of 696, ratio=0.306
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 17/42 (scale 0.00012172353723858711) ***** 
[1] xi(16)=0.999999998547 zeta(17)=0
Stats: nr=31 min=1 max=108 sum3=3544234
Safeguard: keep additional 2 states
Emax=10.0379175393 nrkept=178 nrkeptmult=424 
Total energy=-2.32236107098  Egs=-1.28189200056
Kept: 178 out of 802, ratio=0.222
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 18/42 (scale 7.0277116991411708e-05) ***** 
[1] xi(17)=0.999999999516 zeta(18)=0
Stats: nr=26 min=1 max=82 sum3=2302678
Safeguard: keep additional 2 states
Emax=10.4987901163 nrkept=213 nrkeptmult=523 
Total energy=-2.32242901848  Egs=-0.966850999709
Kept: 213 out of 670, ratio=0.318
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 19/42 (scale 4.0574512412862373e-05) ***** 
[1] xi(18)=0.999999999839 zeta(19)=0
Stats: nr=31 min=1 max=108 sum3=3544234
Safeguard: keep additional 1 states
Emax=10.0984657895 nrkept=177 nrkeptmult=414 
Total energy=-2.32247727794  Egs=-1.1894031765
Kept: 177 out of 802, ratio=0.221
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 20/42 (scale 2.3425705663803904e-05) ***** 
[1] xi(19)=0.999999999946 zeta(20)=0
Stats: nr=26 min=1 max=83 sum3=2328838
Emax=10.0164930687 nrkept=207 nrkeptmult=513 
Total energy=-2.32250211495  Egs=-1.06024592788
Kept: 207 out of 664, ratio=0.312
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 21/42 (scale 1.3524837470954124e-05) ***** 
[1] xi(20)=0.999999999982 zeta(21)=0
Stats: nr=31 min=1 max=104 sum3=3210253
Safeguard: keep additional 1 states
Emax=10.0620295648 nrkept=177 nrkeptmult=414 
Total energy=-2.32251655884  Egs=-1.06795341928
Kept: 177 out of 781, ratio=0.227
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 22/42 (scale 7.8085685546013012e-06) ***** 
[1] xi(21)=0.999999999994 zeta(22)=0
Stats: nr=26 min=1 max=83 sum3=2328838
Emax=10.0152253907 nrkept=195 nrkeptmult=467 
Total energy=-2.32252587511  Egs=-1.19308279184
Kept: 195 out of 664, ratio=0.294
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 23/42 (scale 4.5082791569847079e-06) ***** 
[1] xi(22)=0.999999999998 zeta(23)=0
Stats: nr=29 min=1 max=100 sum3=2849246
Safeguard: keep additional 1 states
Emax=10.0772147585 nrkept=171 nrkeptmult=402 
Total energy=-2.3225300351  Egs=-0.92274366943
Kept: 171 out of 734, ratio=0.233
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 24/42 (scale 2.6028561848671004e-06) ***** 
[1] xi(23)=0.999999999999 zeta(24)=0
Stats: nr=26 min=1 max=80 sum3=2095848
Safeguard: keep additional 1 states
Emax=10.1145340544 nrkept=171 nrkeptmult=395 
Total energy=-2.32253358493  Egs=-1.36382213974
Kept: 171 out of 642, ratio=0.266
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 25/42 (scale 1.5027597189949027e-06) ***** 
[1] xi(24)=1 zeta(25)=0
Stats: nr=27 min=1 max=92 sum3=2081511
Safeguard: keep additional 2 states
Emax=10.0879517067 nrkept=166 nrkeptmult=384 
Total energy=-2.32253476541  Egs=-0.785543455452
Kept: 166 out of 639, ratio=0.26
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 26/42 (scale 8.676187282890334e-07) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=26 min=1 max=78 sum3=1930246
Safeguard: keep additional 1 states
Emax=10.079898437 nrkept=153 nrkeptmult=341 
Total energy=-2.32253609459  Egs=-1.53197819929
Kept: 153 out of 622, ratio=0.246
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 27/42 (scale 5.009199063316343e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=25 min=1 max=85 sum3=1573380
Emax=10.093732146 nrkept=155 nrkeptmult=352 
Total energy=-2.32253644303  Egs=-0.695603386132
Kept: 155 out of 570, ratio=0.272
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 28/42 (scale 2.8920624276301115e-07) ***** 
[1] xi(27)=1 zeta(28)=0
Stats: nr=25 min=1 max=72 sum3=1520960
Safeguard: keep additional 1 states
Emax=10.0807999316 nrkept=128 nrkeptmult=272 
Total energy=-2.32253691791  Egs=-1.64202382757
Kept: 128 out of 578, ratio=0.221
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 29/42 (scale 1.6697330211054477e-07) ***** 
[1] xi(28)=1 zeta(29)=0
Stats: nr=21 min=1 max=75 sum3=1027062
Safeguard: keep additional 3 states
Emax=10.0202124527 nrkept=145 nrkeptmult=318 
Total energy=-2.32253702702  Egs=-0.653458337819
Kept: 145 out of 474, ratio=0.306
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 30/42 (scale 9.6402080921003717e-08) ***** 
[1] xi(29)=1 zeta(30)=0
Stats: nr=25 min=1 max=69 sum3=1313500
Safeguard: keep additional 1 states
Emax=10.0435068638 nrkept=116 nrkeptmult=240 
Total energy=-2.32253719016  Egs=-1.69224698101
Kept: 116 out of 538, ratio=0.216
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 31/42 (scale 5.5657767370181585e-08) ***** 
[1] xi(30)=1 zeta(31)=0
Stats: nr=21 min=1 max=67 sum3=742771
Safeguard: keep additional 1 states
Emax=10.0170050723 nrkept=118 nrkeptmult=262 
Total energy=-2.32253722563  Egs=-0.637258395039
Kept: 118 out of 427, ratio=0.276
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 32/42 (scale 3.213402697366791e-08) ***** 
[1] xi(31)=1 zeta(32)=0
Stats: nr=25 min=1 max=57 sum3=739814
Safeguard: keep additional 1 states
Emax=10.2506778579 nrkept=108 nrkeptmult=220 
Total energy=-2.32253728061  Egs=-1.71111253538
Kept: 108 out of 440, ratio=0.245
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 33/42 (scale 1.8552589123393859e-08) ***** 
[1] xi(32)=1 zeta(33)=0
Stats: nr=21 min=1 max=63 sum3=612299
Emax=10.1630613674 nrkept=112 nrkeptmult=244 
Total energy=-2.32253729233  Egs=-0.631563097403
Kept: 112 out of 395, ratio=0.284
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 34/42 (scale 1.0711342324555969e-08) ***** 
[1] xi(33)=1 zeta(34)=0
Stats: nr=25 min=1 max=54 sum3=614330
Safeguard: keep additional 1 states
Emax=10.4789542459 nrkept=108 nrkeptmult=220 
Total energy=-2.32253731073  Egs=-1.71766835809
Kept: 108 out of 416, ratio=0.26
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 35/42 (scale 6.1841963744646209e-09) ***** 
[1] xi(34)=1 zeta(35)=0
Stats: nr=21 min=1 max=63 sum3=612299
Emax=10.2437698997 nrkept=112 nrkeptmult=244 
Total energy=-2.32253731462  Egs=-0.629628956716
Kept: 112 out of 395, ratio=0.284
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 36/42 (scale 3.5704474415186564e-09) ***** 
[1] xi(35)=1 zeta(36)=0
Stats: nr=25 min=1 max=54 sum3=614330
Safeguard: keep additional 1 states
Emax=10.5590855971 nrkept=108 nrkeptmult=220 
Total energy=-2.32253732076  Egs=-1.71988460013
Kept: 108 out of 416, ratio=0.26
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 37/42 (scale 2.0613987914882068e-09) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=21 min=1 max=63 sum3=612299
Emax=10.2701837067 nrkept=112 nrkeptmult=244 
Total energy=-2.32253732206  Egs=-0.62898015354
Kept: 112 out of 395, ratio=0.284
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 38/42 (scale 1.1901491471728853e-09) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=25 min=1 max=54 sum3=614330
Safeguard: keep additional 1 states
Emax=10.5855978927 nrkept=108 nrkeptmult=220 
Total energy=-2.32253732411  Egs=-1.72062683528
Kept: 108 out of 416, ratio=0.26
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 39/42 (scale 6.8713293049606891e-10) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=21 min=1 max=63 sum3=612299
Emax=10.2789737272 nrkept=112 nrkeptmult=244 
Total energy=-2.32253732454  Egs=-0.628763426999
Kept: 112 out of 395, ratio=0.284
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 40/42 (scale 3.9671638239096179e-10) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=25 min=1 max=54 sum3=614330
Safeguard: keep additional 1 states
Emax=10.5911210617 nrkept=108 nrkeptmult=224 
Total energy=-2.32253732522  Egs=-1.72087463105
Kept: 108 out of 416, ratio=0.26
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 41/42 (scale 2.2904431016535628e-10) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=21 min=1 max=61 sum3=575735
Emax=10.2818991272 nrkept=112 nrkeptmult=244 
Total energy=-2.32253732536  Egs=-0.628691133278
Kept: 112 out of 395, ratio=0.284
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 42/42 (scale 1.3223879413032059e-10) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=25 min=1 max=54 sum3=614330
Safeguard: keep additional 1 states
Emax=10.5921984509 nrkept=108 nrkeptmult=224 
Total energy=-2.32253732559  Egs=-1.72095726143
Kept: 108 out of 416, ratio=0.26
Memory used: 15 MB  Time elapsed: 1 s

** Iteration completed.


Total energy:             -2.322537325592553
TKW=1.14488068857e-06
Grand canonical Z=              5.363972596059725
Zft (last shell)=                5.36208325489966

Memory usage report
===================
  after diag: 16240 kB
after recalc: 16260 kB
 after split: 16256 kB
  after trim: 16260 kB
 after trunc: 16244 kB
      diag_h: 16240 kB

Peak usage: 16260 kB

Timing report
=============
         All: 1 s
        diag: 1 s
    recalc f: 0 s
       Other: 0 s
