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

Number of invariant subspaces: 6
(-2 1) 1 states
0.0318239574146 
(-1 2) 2 states
0.0148661667484 0.0449630015095 
(0 1) 3 states
0 0.0318239574146 0.0598291682578 
(0 3) 1 states
0.0280052108432 
(1 2) 2 states
0.0148661667484 0.0449630015095 
(2 1) 1 states
0.0318239574146 
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.7559421113738346) ***** 
[1] xi(0)=0.76200124248 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=3.05635805669 nrkept=35 nrkeptmult=64 
Total energy=-1.15271488125  Egs=-1.47636544875
Kept: 35 out of 35, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

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

 ***** [NRG] Iteration 3/42 (scale 0.25198070379127824) ***** 
[1] xi(2)=1.01778794749 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Safeguard: keep additional 1 states
Emax=10.3167138118 nrkept=353 nrkeptmult=794 
Total energy=-1.8679642602  Egs=-1.66067403885
Kept: 353 out of 462, ratio=0.764
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/42 (scale 0.14548112716448583) ***** 
[1] xi(3)=1.00124865449 zeta(4)=0
Stats: nr=28 min=1 max=156 sum3=15556360
Safeguard: keep additional 3 states
Emax=10.0689870954 nrkept=416 nrkeptmult=966 
Total energy=-1.96980713755  Egs=-0.700041849648
Kept: 416 out of 1312, ratio=0.317
Memory used: 7 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/42 (scale 0.083993567930426061) ***** 
[1] xi(4)=0.998675647559 zeta(5)=0
Stats: nr=33 min=1 max=203 sum3=23998447
Safeguard: keep additional 1 states
Emax=10.3291742474 nrkept=389 nrkeptmult=874 
Total energy=-2.10394778542  Egs=-1.59703476321
Kept: 389 out of 1555, ratio=0.25
Memory used: 10 MB  Time elapsed: 0 s

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

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

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

 ***** [NRG] Iteration 9/42 (scale 0.0093326186589362328) ***** 
[1] xi(8)=0.999975957408 zeta(9)=0
Stats: nr=33 min=1 max=202 sum3=23580328
Safeguard: keep additional 2 states
Emax=10.002440233 nrkept=355 nrkeptmult=800 
Total energy=-2.20930288491  Egs=-1.48863780012
Kept: 355 out of 1546, ratio=0.23
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 10/42 (scale 0.0053881898949809581) ***** 
[1] xi(9)=0.999991972975 zeta(10)=0
Stats: nr=28 min=1 max=161 sum3=16714434
Emax=10.4709729679 nrkept=413 nrkeptmult=963 
Total energy=-2.21363476941  Egs=-0.803959137242
Kept: 413 out of 1326, ratio=0.311
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 11/42 (scale 0.0031108728863120773) ***** 
[1] xi(10)=0.999997322906 zeta(11)=0
Stats: nr=33 min=1 max=202 sum3=23580328
Safeguard: keep additional 3 states
Emax=10.0605530849 nrkept=319 nrkeptmult=724 
Total energy=-2.21812457445  Egs=-1.4432621347
Kept: 319 out of 1546, ratio=0.206
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 12/42 (scale 0.001796063298326986) ***** 
[1] xi(11)=0.999999107478 zeta(12)=0
Stats: nr=28 min=1 max=142 sum3=11816520
Safeguard: keep additional 4 states
Emax=10.0029690871 nrkept=299 nrkeptmult=695 
Total energy=-2.21961221618  Egs=-0.828279122767
Kept: 299 out of 1188, ratio=0.252
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 13/42 (scale 0.0010369576287706924) ***** 
[1] xi(12)=0.999999702475 zeta(13)=0
Stats: nr=31 min=1 max=154 sum3=9949148
Safeguard: keep additional 2 states
Emax=10.0052215206 nrkept=238 nrkeptmult=552 
Total energy=-2.22106434435  Egs=-1.40037367953
Kept: 238 out of 1118, ratio=0.213
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 14/42 (scale 0.00059868776610899534) ***** 
[1] xi(13)=0.999999900823 zeta(14)=0
Stats: nr=26 min=1 max=109 sum3=5311408
Safeguard: keep additional 2 states
Emax=10.107607787 nrkept=219 nrkeptmult=531 
Total energy=-2.22157890419  Egs=-0.859479461348
Kept: 219 out of 892, ratio=0.246
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 15/42 (scale 0.00034565254292356413) ***** 
[1] xi(14)=0.999999966941 zeta(15)=0
Stats: nr=31 min=1 max=112 sum3=3879366
Emax=10.0217761966 nrkept=183 nrkeptmult=436 
Total energy=-2.22204478521  Egs=-1.34783044253
Kept: 183 out of 822, ratio=0.223
Memory used: 15 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 16/42 (scale 0.00019956258870299843) ***** 
[1] xi(15)=0.99999998898 zeta(16)=0
Stats: nr=26 min=1 max=85 sum3=2548242
Safeguard: keep additional 2 states
Emax=10.3933151013 nrkept=213 nrkeptmult=523 
Total energy=-2.22222538661  Egs=-0.904986258218
Kept: 213 out of 690, ratio=0.309
Memory used: 15 MB  Time elapsed: 1 s

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

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

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

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

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

 ***** [NRG] Iteration 22/42 (scale 7.3912069889999351e-06) ***** 
[1] xi(21)=0.999999999985 zeta(22)=0
Stats: nr=26 min=1 max=83 sum3=2328838
Safeguard: keep additional 1 states
Emax=10.0256512845 nrkept=196 nrkeptmult=468 
Total energy=-2.22252876462  Egs=-1.20085859868
Kept: 196 out of 664, ratio=0.295
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 23/42 (scale 4.2673153447353554e-06) ***** 
[1] xi(22)=0.999999999995 zeta(23)=0
Stats: nr=29 min=1 max=101 sum3=2911002
Safeguard: keep additional 2 states
Emax=10.0867807809 nrkept=172 nrkeptmult=402 
Total energy=-2.22253267051  Egs=-0.915304915706
Kept: 172 out of 738, ratio=0.233
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 24/42 (scale 2.4637356629999781e-06) ***** 
[1] xi(23)=0.999999999998 zeta(24)=0
Stats: nr=26 min=1 max=81 sum3=2132230
Emax=10.0057575336 nrkept=169 nrkeptmult=387 
Total energy=-2.22253605279  Egs=-1.3728265918
Kept: 169 out of 646, ratio=0.262
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 25/42 (scale 1.4224384482451187e-06) ***** 
[1] xi(24)=0.999999999999 zeta(25)=0
Stats: nr=27 min=1 max=92 sum3=2041213
Emax=10.0075988698 nrkept=163 nrkeptmult=378 
Total energy=-2.22253716191  Egs=-0.779725810521
Kept: 163 out of 631, ratio=0.258
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 26/42 (scale 8.212452209999928e-07) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=26 min=1 max=76 sum3=1811914
Safeguard: keep additional 1 states
Emax=10.111918862 nrkept=153 nrkeptmult=341 
Total energy=-2.22253842592  Egs=-1.53914835829
Kept: 153 out of 610, ratio=0.251
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 27/42 (scale 4.7414614941503953e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=25 min=1 max=85 sum3=1573380
Safeguard: keep additional 2 states
Emax=10.0018258822 nrkept=154 nrkeptmult=350 
Total energy=-2.22253875428  Egs=-0.69252543131
Kept: 154 out of 570, ratio=0.27
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 28/42 (scale 2.7374840699999758e-07) ***** 
[1] xi(27)=1 zeta(28)=0
Stats: nr=25 min=1 max=71 sum3=1466038
Safeguard: keep additional 1 states
Emax=10.1101980409 nrkept=128 nrkeptmult=272 
Total energy=-2.2225392048  Egs=-1.64573509168
Kept: 128 out of 574, ratio=0.223
Memory used: 16 MB  Time elapsed: 1 s

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

 ***** [NRG] Iteration 30/42 (scale 9.1249468999999199e-08) ***** 
[1] xi(29)=1 zeta(30)=0
Stats: nr=25 min=1 max=69 sum3=1313500
Safeguard: keep additional 3 states
Emax=10.0256879862 nrkept=114 nrkeptmult=236 
Total energy=-2.22253946243  Egs=-1.69371699412
Kept: 114 out of 538, ratio=0.212
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 31/42 (scale 5.2682905490559947e-08) ***** 
[1] xi(30)=1 zeta(31)=0
Stats: nr=21 min=1 max=65 sum3=686951
Emax=10.0028321386 nrkept=115 nrkeptmult=254 
Total energy=-2.22253949598  Egs=-0.636807426267
Kept: 115 out of 419, ratio=0.274
Memory used: 16 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 32/42 (scale 3.04164896666664e-08) ***** 
[1] xi(31)=1 zeta(32)=0
Stats: nr=25 min=1 max=55 sum3=663626
Safeguard: keep additional 1 states
Emax=10.2683974409 nrkept=108 nrkeptmult=220 
Total energy=-2.22253954804  Egs=-1.7116332284
Kept: 108 out of 428, ratio=0.252
Memory used: 16 MB  Time elapsed: 1 s

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

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

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

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

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

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

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

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

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

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

** Iteration completed.


Total energy:             -2.222539590617925
TKW=1.14318196367e-06
Grand canonical Z=              5.363960727361752
Zft (last shell)=                5.36207100742527

Memory usage report
===================
  after diag: 16772 kB
after recalc: 16796 kB
 after split: 16792 kB
  after trim: 16796 kB
 after trunc: 16780 kB
      diag_h: 16772 kB

Peak usage: 16796 kB

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