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.55929062767e-10
Nmax=42
Tmin=1e-10
Ninit=0
Lambda=3
strategy=kept
lastall=0
lastalloverride=0
z=0.6
discretization=Z
band=flat
betabar=1
T=1.55929062767e-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.6
>>  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.6
# keep=5000 keepenergy=10 Nmax=42 channels=1 dots=1 betabar=1

Number of invariant subspaces: 6
(-2 1) 1 states
0.0255464132764 
(-1 2) 2 states
0.011933689913 0.0360936700531 
(0 1) 3 states
0 0.0255464132764 0.0480273599661 
(0 3) 1 states
0.0224809466897 
(1 2) 2 states
0.011933689913 0.0360936700531 
(2 1) 1 states
0.0255464132764 
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.94170047669666523) ***** 
[1] xi(0)=0.650627930866 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=2.60913427527 nrkept=35 nrkeptmult=64 
Total energy=-1.22601149539  Egs=-1.26297422383
Kept: 35 out of 35, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

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

 ***** [NRG] Iteration 3/42 (scale 0.31390015889888839) ***** 
[1] xi(2)=1.08675662916 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Safeguard: keep additional 7 states
Emax=10.3939832064 nrkept=399 nrkeptmult=904 
Total energy=-2.17510274429  Egs=-1.72913094967
Kept: 399 out of 462, ratio=0.864
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/42 (scale 0.1812303412389395) ***** 
[1] xi(3)=1.03462125972 zeta(4)=0
Stats: nr=28 min=1 max=176 sum3=22434060
Safeguard: keep additional 1 states
Emax=10.9599756687 nrkept=564 nrkeptmult=1320 
Total energy=-2.30611943129  Egs=-0.722929097312
Kept: 564 out of 1488, ratio=0.379
Memory used: 8 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/42 (scale 0.10463338629962946) ***** 
[1] xi(4)=1.00713692568 zeta(5)=0
Stats: nr=33 min=1 max=270 sum3=58136254
Safeguard: keep additional 3 states
Emax=10.4330551847 nrkept=431 nrkeptmult=976 
Total energy=-2.47473235992  Egs=-1.61146393698
Kept: 431 out of 2110, ratio=0.204
Memory used: 13 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 6/42 (scale 0.060410113746313165) ***** 
[1] xi(5)=1.00170828391 zeta(6)=0
Stats: nr=28 min=1 max=190 sum3=28423954
Safeguard: keep additional 7 states
Emax=10.2980600905 nrkept=420 nrkeptmult=976 
Total energy=-2.51893593405  Egs=-0.73172472933
Kept: 420 out of 1606, ratio=0.262
Memory used: 15 MB  Time elapsed: 0 s

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

 ***** [NRG] Iteration 8/42 (scale 0.020136704582104387) ***** 
[1] xi(7)=1.00016539487 zeta(8)=0
Stats: nr=28 min=1 max=173 sum3=21206352
Safeguard: keep additional 3 states
Emax=10.328834685 nrkept=416 nrkeptmult=972 
Total energy=-2.58862452428  Egs=-0.76776686935
Kept: 416 out of 1446, ratio=0.288
Memory used: 16 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 9/42 (scale 0.011625931811069938) ***** 
[1] xi(8)=1.00005458281 zeta(9)=0
Stats: nr=33 min=1 max=203 sum3=24107488
Safeguard: keep additional 1 states
Emax=10.1538557214 nrkept=357 nrkeptmult=802 
Total energy=-2.60605121589  Egs=-1.49895009616
Kept: 357 out of 1558, ratio=0.229
Memory used: 18 MB  Time elapsed: 0 s

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

 ***** [NRG] Iteration 11/42 (scale 0.0038753106036899816) ***** 
[1] xi(10)=1.00000603819 zeta(11)=0
Stats: nr=33 min=1 max=202 sum3=23580328
Safeguard: keep additional 2 states
Emax=10.0017186363 nrkept=330 nrkeptmult=752 
Total energy=-2.61703919436  Egs=-1.45172590911
Kept: 330 out of 1546, ratio=0.213
Memory used: 17 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 12/42 (scale 0.0022374116202338221) ***** 
[1] xi(11)=1.000002012 zeta(12)=0
Stats: nr=28 min=1 max=146 sum3=12815432
Safeguard: keep additional 1 states
Emax=10.1274402688 nrkept=334 nrkeptmult=772 
Total energy=-2.61888108011  Egs=-0.823221678872
Kept: 334 out of 1232, ratio=0.271
Memory used: 17 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 13/42 (scale 0.0012917702012299939) ***** 
[1] xi(12)=1.00000067058 zeta(13)=0
Stats: nr=31 min=1 max=169 sum3=13415530
Safeguard: keep additional 3 states
Emax=10.0071678104 nrkept=259 nrkeptmult=592 
Total energy=-2.62070169689  Egs=-1.40939679708
Kept: 259 out of 1246, ratio=0.208
Memory used: 17 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 14/42 (scale 0.00074580387341127404) ***** 
[1] xi(13)=1.00000022352 zeta(14)=0
Stats: nr=26 min=1 max=116 sum3=6559786
Safeguard: keep additional 2 states
Emax=10.0685823656 nrkept=235 nrkeptmult=563 
Total energy=-2.62133731841  Egs=-0.852263629625
Kept: 235 out of 964, ratio=0.244
Memory used: 17 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 15/42 (scale 0.0004305900670766646) ***** 
[1] xi(14)=1.00000007451 zeta(15)=0
Stats: nr=31 min=1 max=122 sum3=4917132
Emax=10.1081619623 nrkept=184 nrkeptmult=440 
Total energy=-2.6219227165  Egs=-1.35952531491
Kept: 184 out of 882, ratio=0.209
Memory used: 17 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 16/42 (scale 0.00024860129113709131) ***** 
[1] xi(15)=1.00000002483 zeta(16)=0
Stats: nr=26 min=1 max=85 sum3=2580766
Safeguard: keep additional 2 states
Emax=10.3619315297 nrkept=213 nrkeptmult=523 
Total energy=-2.6221450882  Egs=-0.894491302914
Kept: 213 out of 694, ratio=0.307
Memory used: 17 MB  Time elapsed: 1 s

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

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

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

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

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

 ***** [NRG] Iteration 22/42 (scale 9.2074552272996784e-06) ***** 
[1] xi(21)=1.00000000003 zeta(22)=0
Stats: nr=26 min=1 max=83 sum3=2328838
Safeguard: keep additional 2 states
Emax=10.0054955437 nrkept=200 nrkeptmult=480 
Total energy=-2.62252455864  Egs=-1.17039703945
Kept: 200 out of 664, ratio=0.301
Memory used: 18 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 23/42 (scale 5.3159267540328959e-06) ***** 
[1] xi(22)=1.00000000001 zeta(23)=0
Stats: nr=29 min=1 max=103 sum3=3124876
Safeguard: keep additional 1 states
Emax=10.0301186768 nrkept=171 nrkeptmult=402 
Total energy=-2.62252958284  Egs=-0.945122780404
Kept: 171 out of 754, ratio=0.227
Memory used: 18 MB  Time elapsed: 1 s

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

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

 ***** [NRG] Iteration 26/42 (scale 1.0230505808110754e-06) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=26 min=1 max=78 sum3=1930246
Safeguard: keep additional 2 states
Emax=10.0116233024 nrkept=159 nrkeptmult=365 
Total energy=-2.62253665441  Egs=-1.50952800423
Kept: 159 out of 622, ratio=0.256
Memory used: 18 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 27/42 (scale 5.9065852822587735e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=27 min=1 max=86 sum3=1665372
Safeguard: keep additional 1 states
Emax=10.1383844831 nrkept=157 nrkeptmult=358 
Total energy=-2.62253707117  Egs=-0.705585774671
Kept: 157 out of 594, ratio=0.264
Memory used: 18 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 28/42 (scale 3.4101686027035842e-07) ***** 
[1] xi(27)=1 zeta(28)=0
Stats: nr=25 min=1 max=73 sum3=1593760
Safeguard: keep additional 2 states
Emax=10.1597137707 nrkept=131 nrkeptmult=281 
Total energy=-2.62253762701  Egs=-1.6299507073
Kept: 131 out of 586, ratio=0.224
Memory used: 18 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 29/42 (scale 1.9688617607529248e-07) ***** 
[1] xi(28)=1 zeta(29)=0
Stats: nr=21 min=1 max=76 sum3=1088178
Safeguard: keep additional 2 states
Emax=10.1454006834 nrkept=148 nrkeptmult=324 
Total energy=-2.62253775648  Egs=-0.657602684267
Kept: 148 out of 486, ratio=0.305
Memory used: 18 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 30/42 (scale 1.1367228675678615e-07) ***** 
[1] xi(29)=1 zeta(30)=0
Stats: nr=25 min=1 max=71 sum3=1411012
Safeguard: keep additional 1 states
Emax=10.2134793718 nrkept=118 nrkeptmult=244 
Total energy=-2.62253794829  Egs=-1.68737474959
Kept: 118 out of 550, ratio=0.215
Memory used: 18 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 31/42 (scale 6.5628725358430826e-08) ***** 
[1] xi(30)=1 zeta(31)=0
Stats: nr=21 min=1 max=69 sum3=801675
Safeguard: keep additional 1 states
Emax=10.015068421 nrkept=120 nrkeptmult=268 
Total energy=-2.62253799021  Egs=-0.638761783911
Kept: 120 out of 435, ratio=0.276
Memory used: 18 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 32/42 (scale 3.7890762252262051e-08) ***** 
[1] xi(31)=1 zeta(32)=0
Stats: nr=25 min=1 max=58 sum3=785656
Safeguard: keep additional 1 states
Emax=10.1929636416 nrkept=108 nrkeptmult=220 
Total energy=-2.62253805498  Egs=-1.70937489063
Kept: 108 out of 448, ratio=0.241
Memory used: 18 MB  Time elapsed: 1 s

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

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

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

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

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

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

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

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

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

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

** Iteration completed.


Total energy:             -2.622538108021265
TKW=1.12748541366e-06
Grand canonical Z=              5.364016957209835
Zft (last shell)=               5.362124274382368

Memory usage report
===================
  after diag: 19064 kB
after recalc: 19084 kB
 after split: 19080 kB
  after trim: 19084 kB
 after trunc: 19068 kB
      diag_h: 19064 kB

Peak usage: 19084 kB

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