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=67
TEMPERATURE T=1e-10
Nmax=67
Tmin=1e-10
Ninit=0
Lambda=2
strategy=kept
lastall=1
lastalloverride=0
z=0.5
discretization=Z
band=flat
betabar=1
T=1e-10
T_ratio=1
Tmin_ratio=-999
keep=10000
keepenergy=10
keepmin=0
keepblock=0
keepzrescale=0
log=
logall=0
logenumber=10
smooth=new
alpha=0.3
omega0=1e-99
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=0.1
broaden_min=1e-08
broaden_min_ratio=3
broaden_ratio=1.02
dumpenergies=0
saveall=0
dumpannotated=0
dumpscaled=1
dumpabs=0
dumpprecision=8
dumpgroups=1
grouptol=1e-06
dumpdiagonal=0
ops=A_d
specs=
specd=A_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=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.001	 Delta=	0	 t=	0.
>> 
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.5
>>  ops=	A_d
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=6
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||||> (z)

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

Number of invariant subspaces: 6
(-2 1) 1 states
0.0367538100397 
(-1 2) 2 states
0.0174462807917 0.0525956033849 
(0 1) 3 states
0 0.0367538100397 0.0700418841766 
(0 3) 1 states
0.0332880741369 
(1 2) 2 states
0.0174462807917 0.0525956033849 
(2 1) 1 states
0.0367538100397 
Number of states (multiplicity taken into account): 16

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

Computing the following spectra:

Before NRG iteration (N=-1)

 ***** [NRG] Iteration 1/67 (scale 1.0201394465967895) ***** 
[1] xi(0)=0.583569837945 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=2.34125880272 nrkept=35 nrkeptmult=64 
Total energy=-1.19170522967  Egs=-1.12110231009
Kept: 35 out of 35, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/67 (scale 0.72134752044448169) ***** 
[1] xi(1)=0.776754544846 zeta(2)=0
Stats: nr=15 min=1 max=20 sum3=30066
Emax=4.6306144545 nrkept=126 nrkeptmult=256 
Total energy=-1.66764112744  Egs=-0.659787251333
Kept: 126 out of 126, ratio=1
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/67 (scale 0.51006972329839473) ***** 
[1] xi(2)=0.992903283056 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=9.70884141923 nrkept=462 nrkeptmult=1024 
Total energy=-2.47359302813  Egs=-1.58008182778
Kept: 462 out of 462, ratio=1
Memory used: 5 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/67 (scale 0.36067376022224085) ***** 
[1] xi(3)=1.04266411315 zeta(4)=0
Stats: nr=28 min=1 max=210 sum3=37226316
Emax=10.2107510943 nrkept=1259 nrkeptmult=3041 
Total energy=-2.79981331329  Egs=-0.904474683605
Kept: 1259 out of 1716, ratio=0.734
Memory used: 13 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/67 (scale 0.25503486164919736) ***** 
[1] xi(4)=1.01328740971 zeta(5)=0
Stats: nr=36 min=1 max=561 sum3=559501320
Safeguard: keep additional 3 states
Emax=10.2083031786 nrkept=1883 nrkeptmult=4688 
Total energy=-3.18999951897  Egs=-1.52993282236
Kept: 1883 out of 4728, ratio=0.398
Memory used: 21 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 6/67 (scale 0.18033688011112042) ***** 
[1] xi(5)=0.999788616061 zeta(6)=0
Stats: nr=42 min=1 max=717 sum3=1710937002
Safeguard: keep additional 3 states
Emax=10.2375168905 nrkept=2124 nrkeptmult=5286 
Total energy=-3.34861375836  Egs=-0.879544102641
Kept: 2124 out of 7092, ratio=0.299
Memory used: 27 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 7/67 (scale 0.12751743082459868) ***** 
[1] xi(6)=0.998344002082 zeta(7)=0
Stats: nr=44 min=1 max=910 sum3=2520131436
Safeguard: keep additional 2 states
Emax=10.0192263144 nrkept=1921 nrkeptmult=4746 
Total energy=-3.5388071973  Egs=-1.49150933881
Kept: 1921 out of 8010, ratio=0.24
Memory used: 28 MB  Time elapsed: 13 s

 ***** [NRG] Iteration 8/67 (scale 0.090168440055560212) ***** 
[1] xi(7)=0.998902918682 zeta(8)=0
Stats: nr=42 min=1 max=745 sum3=1897784946
Safeguard: keep additional 2 states
Emax=10.0000486937 nrkept=1999 nrkeptmult=4973 
Total energy=-3.62044366844  Egs=-0.905377436836
Kept: 1999 out of 7230, ratio=0.276
Memory used: 28 MB  Time elapsed: 19 s

 ***** [NRG] Iteration 9/67 (scale 0.063758715412299341) ***** 
[1] xi(8)=0.999398551363 zeta(9)=0
Stats: nr=44 min=1 max=863 sum3=2129089633
Safeguard: keep additional 1 states
Emax=10.0019670946 nrkept=1664 nrkeptmult=4134 
Total energy=-3.71353337391  Egs=-1.46003106974
Kept: 1664 out of 7549, ratio=0.22
Memory used: 26 MB  Time elapsed: 24 s

 ***** [NRG] Iteration 10/67 (scale 0.045084220027780106) ***** 
[1] xi(9)=0.99968731683 zeta(10)=0
Stats: nr=42 min=1 max=652 sum3=1254088034
Emax=10.0473810799 nrkept=1873 nrkeptmult=4643 
Total energy=-3.7556615186  Egs=-0.934432150913
Kept: 1873 out of 6284, ratio=0.298
Memory used: 29 MB  Time elapsed: 28 s

 ***** [NRG] Iteration 11/67 (scale 0.031879357706149671) ***** 
[1] xi(10)=0.9998407918 zeta(11)=0
Stats: nr=42 min=1 max=815 sum3=1780810507
Safeguard: keep additional 7 states
Emax=10.0127857568 nrkept=1638 nrkeptmult=4060 
Total energy=-3.80107058612  Egs=-1.42440346325
Kept: 1638 out of 7057, ratio=0.232
Memory used: 28 MB  Time elapsed: 33 s

 ***** [NRG] Iteration 12/67 (scale 0.022542110013890053) ***** 
[1] xi(11)=0.99991969299 zeta(12)=0
Stats: nr=42 min=1 max=640 sum3=1188522258
Safeguard: keep additional 2 states
Emax=10.0015985499 nrkept=1757 nrkeptmult=4357 
Total energy=-3.8228058567  Egs=-0.96420745744
Kept: 1757 out of 6180, ratio=0.284
Memory used: 29 MB  Time elapsed: 36 s

 ***** [NRG] Iteration 13/67 (scale 0.015939678853074835) ***** 
[1] xi(12)=0.999959672404 zeta(13)=0
Stats: nr=42 min=1 max=761 sum3=1454172497
Safeguard: keep additional 3 states
Emax=10.0035212517 nrkept=1554 nrkeptmult=3828 
Total energy=-3.84491889148  Egs=-1.38729487465
Kept: 1554 out of 6617, ratio=0.235
Memory used: 29 MB  Time elapsed: 41 s

 ***** [NRG] Iteration 14/67 (scale 0.011271055006945026) ***** 
[1] xi(13)=0.999979792879 zeta(14)=0
Stats: nr=42 min=1 max=608 sum3=1018483162
Emax=10.0015336857 nrkept=1597 nrkeptmult=3941 
Total energy=-3.85611220873  Egs=-0.993102884985
Kept: 1597 out of 5854, ratio=0.273
Memory used: 29 MB  Time elapsed: 44 s

 ***** [NRG] Iteration 15/67 (scale 0.0079698394265374176) ***** 
[1] xi(14)=0.999989885634 zeta(15)=0
Stats: nr=40 min=1 max=705 sum3=1134126634
Safeguard: keep additional 3 states
Emax=10.0011836559 nrkept=1516 nrkeptmult=3738 
Total energy=-3.86687668148  Egs=-1.35065114648
Kept: 1516 out of 6016, ratio=0.252
Memory used: 29 MB  Time elapsed: 47 s

 ***** [NRG] Iteration 16/67 (scale 0.0056355275034725132) ***** 
[1] xi(15)=0.999994940118 zeta(16)=0
Stats: nr=40 min=1 max=596 sum3=957480812
Safeguard: keep additional 2 states
Emax=10.0009225943 nrkept=1465 nrkeptmult=3555 
Total energy=-3.87263091092  Egs=-1.02106314521
Kept: 1465 out of 5708, ratio=0.257
Memory used: 28 MB  Time elapsed: 50 s

 ***** [NRG] Iteration 17/67 (scale 0.0039849197132687088) ***** 
[1] xi(16)=0.999997469385 zeta(17)=0
Stats: nr=36 min=1 max=658 sum3=902028083
Safeguard: keep additional 1 states
Emax=10.0002540846 nrkept=1327 nrkeptmult=3260 
Total energy=-3.87786708483  Egs=-1.31399734086
Kept: 1327 out of 5513, ratio=0.241
Memory used: 29 MB  Time elapsed: 53 s

 ***** [NRG] Iteration 18/67 (scale 0.0028177637517362566) ***** 
[1] xi(17)=0.999998734524 zeta(18)=0
Stats: nr=40 min=1 max=529 sum3=666590356
Safeguard: keep additional 1 states
Emax=10.0000835745 nrkept=1281 nrkeptmult=3147 
Total energy=-3.88082780332  Egs=-1.05073340162
Kept: 1281 out of 4990, ratio=0.257
Memory used: 29 MB  Time elapsed: 56 s

 ***** [NRG] Iteration 19/67 (scale 0.0019924598566343544) ***** 
[1] xi(18)=0.99999936722 zeta(19)=0
Stats: nr=34 min=1 max=579 sum3=616598005
Safeguard: keep additional 1 states
Emax=10.0009765181 nrkept=1154 nrkeptmult=2870 
Total energy=-3.8833673909  Egs=-1.27459911917
Kept: 1154 out of 4825, ratio=0.239
Memory used: 27 MB  Time elapsed: 57 s

 ***** [NRG] Iteration 20/67 (scale 0.0014088818758681283) ***** 
[1] xi(19)=0.999999683599 zeta(20)=0
Stats: nr=40 min=1 max=462 sum3=445382958
Safeguard: keep additional 1 states
Emax=10.0052397336 nrkept=1096 nrkeptmult=2722 
Total energy=-3.88489741659  Egs=-1.08598578724
Kept: 1096 out of 4356, ratio=0.252
Memory used: 29 MB  Time elapsed: 59 s

 ***** [NRG] Iteration 21/67 (scale 0.0009962299283171772) ***** 
[1] xi(20)=0.999999841797 zeta(21)=0
Stats: nr=34 min=1 max=500 sum3=395810408
Emax=10.0046911021 nrkept=963 nrkeptmult=2416 
Total energy=-3.88612196156  Egs=-1.22917906529
Kept: 963 out of 4136, ratio=0.233
Memory used: 26 MB  Time elapsed: 60 s

 ***** [NRG] Iteration 22/67 (scale 0.00070444093793406415) ***** 
[1] xi(21)=0.999999920898 zeta(22)=0
Stats: nr=40 min=1 max=391 sum3=266535976
Safeguard: keep additional 2 states
Emax=10.0086552413 nrkept=895 nrkeptmult=2217 
Total energy=-3.88691806198  Egs=-1.13011663131
Kept: 895 out of 3634, ratio=0.246
Memory used: 29 MB  Time elapsed: 61 s

 ***** [NRG] Iteration 23/67 (scale 0.0004981149641585886) ***** 
[1] xi(22)=0.999999960449 zeta(23)=0
Stats: nr=34 min=1 max=417 sum3=224530091
Safeguard: keep additional 2 states
Emax=10.0005993028 nrkept=868 nrkeptmult=2218 
Total energy=-3.88750352867  Egs=-1.17536458788
Kept: 868 out of 3377, ratio=0.257
Memory used: 30 MB  Time elapsed: 62 s

 ***** [NRG] Iteration 24/67 (scale 0.00035222046896703208) ***** 
[1] xi(23)=0.999999980224 zeta(24)=0
Stats: nr=40 min=1 max=354 sum3=191343178
Safeguard: keep additional 1 states
Emax=10.0076766313 nrkept=823 nrkeptmult=2039 
Total energy=-3.88792105272  Egs=-1.18540540022
Kept: 823 out of 3286, ratio=0.25
Memory used: 29 MB  Time elapsed: 63 s

 ***** [NRG] Iteration 25/67 (scale 0.0002490574820792943) ***** 
[1] xi(24)=0.999999990112 zeta(25)=0
Stats: nr=34 min=1 max=383 sum3=172023467
Emax=10.0036971781 nrkept=816 nrkeptmult=2066 
Total energy=-3.88819812443  Egs=-1.11248096225
Kept: 816 out of 3101, ratio=0.263
Memory used: 30 MB  Time elapsed: 63 s

 ***** [NRG] Iteration 26/67 (scale 0.00017611023448351604) ***** 
[1] xi(25)=0.999999995056 zeta(26)=0
Stats: nr=38 min=1 max=333 sum3=161834170
Safeguard: keep additional 3 states
Emax=10.0237130106 nrkept=783 nrkeptmult=1957 
Total energy=-3.8884187413  Egs=-1.25272033038
Kept: 783 out of 3082, ratio=0.254
Memory used: 29 MB  Time elapsed: 64 s

 ***** [NRG] Iteration 27/67 (scale 0.00012452874103964715) ***** 
[1] xi(26)=0.999999997528 zeta(27)=0
Stats: nr=34 min=1 max=359 sum3=144739908
Emax=10.0104228215 nrkept=798 nrkeptmult=2024 
Total energy=-3.88854859037  Egs=-1.0427237059
Kept: 798 out of 2952, ratio=0.27
Memory used: 29 MB  Time elapsed: 65 s

 ***** [NRG] Iteration 28/67 (scale 8.8055117241758019e-05) ***** 
[1] xi(27)=0.999999998764 zeta(28)=0
Stats: nr=38 min=1 max=327 sum3=150905436
Safeguard: keep additional 1 states
Emax=10.0147851711 nrkept=758 nrkeptmult=1918 
Total energy=-3.88866568899  Egs=-1.32983321272
Kept: 758 out of 3012, ratio=0.252
Memory used: 29 MB  Time elapsed: 65 s

 ***** [NRG] Iteration 29/67 (scale 6.2264370519823575e-05) ***** 
[1] xi(28)=0.999999999382 zeta(29)=0
Stats: nr=34 min=1 max=346 sum3=131798132
Safeguard: keep additional 1 states
Emax=10.0123129494 nrkept=765 nrkeptmult=1910 
Total energy=-3.88872624242  Egs=-0.972521427668
Kept: 765 out of 2864, ratio=0.267
Memory used: 30 MB  Time elapsed: 66 s

 ***** [NRG] Iteration 30/67 (scale 4.402755862087901e-05) ***** 
[1] xi(29)=0.999999999691 zeta(30)=0
Stats: nr=38 min=1 max=316 sum3=137692612
Safeguard: keep additional 1 states
Emax=10.0127816539 nrkept=715 nrkeptmult=1779 
Total energy=-3.88878828659  Egs=-1.40921200186
Kept: 715 out of 2884, ratio=0.248
Memory used: 30 MB  Time elapsed: 66 s

 ***** [NRG] Iteration 31/67 (scale 3.1132185259911788e-05) ***** 
[1] xi(30)=0.999999999846 zeta(31)=0
Stats: nr=33 min=1 max=330 sum3=112532256
Safeguard: keep additional 2 states
Emax=10.0193565251 nrkept=722 nrkeptmult=1782 
Total energy=-3.88881665758  Egs=-0.911307522098
Kept: 722 out of 2694, ratio=0.268
Memory used: 31 MB  Time elapsed: 67 s

 ***** [NRG] Iteration 32/67 (scale 2.2013779310439505e-05) ***** 
[1] xi(31)=0.999999999923 zeta(32)=0
Stats: nr=36 min=1 max=301 sum3=119695024
Safeguard: keep additional 3 states
Emax=10.0171895907 nrkept=653 nrkeptmult=1597 
Total energy=-3.88884922005  Egs=-1.47918576948
Kept: 653 out of 2722, ratio=0.24
Memory used: 30 MB  Time elapsed: 67 s

 ***** [NRG] Iteration 33/67 (scale 1.5566092629955894e-05) ***** 
[1] xi(32)=0.999999999961 zeta(33)=0
Stats: nr=33 min=1 max=310 sum3=90015037
Safeguard: keep additional 4 states
Emax=10.0121834087 nrkept=679 nrkeptmult=1672 
Total energy=-3.88886270088  Egs=-0.866037886274
Kept: 679 out of 2461, ratio=0.276
Memory used: 30 MB  Time elapsed: 68 s

 ***** [NRG] Iteration 34/67 (scale 1.1006889655219752e-05) ***** 
[1] xi(33)=0.999999999981 zeta(34)=0
Stats: nr=36 min=1 max=282 sum3=98996162
Safeguard: keep additional 1 states
Emax=10.0086245136 nrkept=592 nrkeptmult=1416 
Total energy=-3.88887955151  Egs=-1.53091713637
Kept: 592 out of 2558, ratio=0.231
Memory used: 30 MB  Time elapsed: 68 s

 ***** [NRG] Iteration 35/67 (scale 7.7830463149779469e-06) ***** 
[1] xi(34)=0.99999999999 zeta(35)=0
Stats: nr=33 min=1 max=287 sum3=68884364
Safeguard: keep additional 1 states
Emax=10.0449491878 nrkept=603 nrkeptmult=1444 
Total energy=-3.88888606679  Egs=-0.837111993638
Kept: 603 out of 2222, ratio=0.271
Memory used: 31 MB  Time elapsed: 68 s

 ***** [NRG] Iteration 36/67 (scale 5.5034448276098762e-06) ***** 
[1] xi(35)=0.999999999995 zeta(36)=0
Stats: nr=36 min=1 max=256 sum3=72549386
Emax=10.0041334351 nrkept=535 nrkeptmult=1261 
Total energy=-3.88889467272  Egs=-1.56373386014
Kept: 535 out of 2264, ratio=0.236
Memory used: 31 MB  Time elapsed: 69 s

 ***** [NRG] Iteration 37/67 (scale 3.8915231574889735e-06) ***** 
[1] xi(36)=0.999999999998 zeta(37)=0
Stats: nr=33 min=1 max=262 sum3=51852835
Safeguard: keep additional 1 states
Emax=10.084200852 nrkept=569 nrkeptmult=1328 
Total energy=-3.88889786556  Egs=-0.820460290838
Kept: 569 out of 2005, ratio=0.284
Memory used: 31 MB  Time elapsed: 69 s

 ***** [NRG] Iteration 38/67 (scale 2.7517224138049381e-06) ***** 
[1] xi(37)=0.999999999999 zeta(38)=0
Stats: nr=34 min=1 max=245 sum3=62083694
Safeguard: keep additional 1 states
Emax=10.025184247 nrkept=505 nrkeptmult=1165 
Total energy=-3.88890222005  Egs=-1.58246120929
Kept: 505 out of 2132, ratio=0.237
Memory used: 31 MB  Time elapsed: 69 s

 ***** [NRG] Iteration 39/67 (scale 1.9457615787444867e-06) ***** 
[1] xi(38)=0.999999999999 zeta(39)=0
Stats: nr=33 min=1 max=249 sum3=44336516
Emax=10.4459103367 nrkept=567 nrkeptmult=1322 
Total energy=-3.88890379898  Egs=-0.811471726585
Kept: 567 out of 1886, ratio=0.301
Memory used: 32 MB  Time elapsed: 70 s

 ***** [NRG] Iteration 40/67 (scale 1.3758612069024691e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.1981441091 nrkept=499 nrkeptmult=1149 
Total energy=-3.88890599004  Egs=-1.59250011274
Kept: 499 out of 2124, ratio=0.235
Memory used: 31 MB  Time elapsed: 70 s

 ***** [NRG] Iteration 41/67 (scale 9.7288078937224336e-07) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.6685417454 nrkept=567 nrkeptmult=1322 
Total energy=-3.88890677495  Egs=-0.806793095503
Kept: 567 out of 1862, ratio=0.305
Memory used: 31 MB  Time elapsed: 70 s

 ***** [NRG] Iteration 42/67 (scale 6.8793060345123453e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.4002552911 nrkept=499 nrkeptmult=1149 
Total energy=-3.88890787406  Egs=-1.59770203897
Kept: 499 out of 2124, ratio=0.235
Memory used: 31 MB  Time elapsed: 70 s

 ***** [NRG] Iteration 43/67 (scale 4.8644039468612168e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.7537218724 nrkept=567 nrkeptmult=1322 
Total energy=-3.88890826536  Egs=-0.804404966245
Kept: 567 out of 1862, ratio=0.305
Memory used: 31 MB  Time elapsed: 71 s

 ***** [NRG] Iteration 44/67 (scale 3.4396530172561726e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5011034821 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890881582  Egs=-1.60035044726
Kept: 498 out of 2124, ratio=0.234
Memory used: 32 MB  Time elapsed: 71 s

 ***** [NRG] Iteration 45/67 (scale 2.4322019734306084e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.7976645703 nrkept=567 nrkeptmult=1322 
Total energy=-3.88890901118  Egs=-0.803198325593
Kept: 567 out of 1857, ratio=0.305
Memory used: 32 MB  Time elapsed: 71 s

 ***** [NRG] Iteration 46/67 (scale 1.7198265086280863e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5135166719 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890928664  Egs=-1.60168674875
Kept: 498 out of 2124, ratio=0.234
Memory used: 31 MB  Time elapsed: 71 s

 ***** [NRG] Iteration 47/67 (scale 1.2161009867153042e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8198704514 nrkept=567 nrkeptmult=1322 
Total energy=-3.88890938424  Egs=-0.8025918124
Kept: 567 out of 1857, ratio=0.305
Memory used: 31 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 48/67 (scale 8.5991325431404316e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5197766792 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890952203  Egs=-1.60235795384
Kept: 498 out of 2124, ratio=0.234
Memory used: 32 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 49/67 (scale 6.080504933576521e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8304354411 nrkept=567 nrkeptmult=1324 
Total energy=-3.88890957081  Egs=-0.802287751331
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 50/67 (scale 4.2995662715702158e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5229304884 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890963972  Egs=-1.60269432375
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 51/67 (scale 3.0402524667882605e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8340453147 nrkept=567 nrkeptmult=1324 
Total energy=-3.88890966411  Egs=-0.802135518889
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 52/67 (scale 2.1497831357851079e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5245060286 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890969857  Egs=-1.60286270102
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 73 s

 ***** [NRG] Iteration 53/67 (scale 1.5201262333941303e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8358467908 nrkept=567 nrkeptmult=1324 
Total energy=-3.88890971076  Egs=-0.802059352092
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 73 s

 ***** [NRG] Iteration 54/67 (scale 1.0748915678925539e-08) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5252946654 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890972799  Egs=-1.60294693779
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 73 s

 ***** [NRG] Iteration 55/67 (scale 7.6006311669706513e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8367466745 nrkept=567 nrkeptmult=1324 
Total energy=-3.88890973409  Egs=-0.802021256037
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 73 s

 ***** [NRG] Iteration 56/67 (scale 5.3744578394627697e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5256892009 nrkept=498 nrkeptmult=1142 
Total energy=-3.8889097427  Egs=-1.60298906821
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 57/67 (scale 3.8003155834853256e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8371964056 nrkept=567 nrkeptmult=1324 
Total energy=-3.88890974575  Egs=-0.802002204845
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 58/67 (scale 2.6872289197313849e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5258865229 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890975006  Egs=-1.60301013644
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 59/67 (scale 1.9001577917426628e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8374212188 nrkept=567 nrkeptmult=1324 
Total energy=-3.88890975158  Egs=-0.801992678457
Kept: 567 out of 1857, ratio=0.305
Memory used: 33 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 60/67 (scale 1.3436144598656924e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5259851976 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890975373  Egs=-1.6030206713
Kept: 498 out of 2124, ratio=0.234
Memory used: 33 MB  Time elapsed: 75 s

 ***** [NRG] Iteration 61/67 (scale 9.5007889587133141e-10) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 2 states
Emax=10.8375360849 nrkept=569 nrkeptmult=1330 
Total energy=-3.8889097545  Egs=-0.801987915066
Kept: 569 out of 1857, ratio=0.306
Memory used: 34 MB  Time elapsed: 75 s

 ***** [NRG] Iteration 62/67 (scale 6.7180722993284622e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=34 min=1 max=244 sum3=61697330
Emax=10.5260345383 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890975557  Egs=-1.60302593892
Kept: 498 out of 2132, ratio=0.234
Memory used: 34 MB  Time elapsed: 75 s

 ***** [NRG] Iteration 63/67 (scale 4.750394479356657e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8376027884 nrkept=575 nrkeptmult=1340 
Total energy=-3.88890975595  Egs=-0.801985533321
Kept: 575 out of 1857, ratio=0.31
Memory used: 34 MB  Time elapsed: 75 s

 ***** [NRG] Iteration 64/67 (scale 3.3590361496642311e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260447789 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890975649  Egs=-1.60302857278
Kept: 498 out of 2154, ratio=0.231
Memory used: 34 MB  Time elapsed: 76 s

 ***** [NRG] Iteration 65/67 (scale 2.3751972396783285e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8376184865 nrkept=575 nrkeptmult=1340 
Total energy=-3.88890975668  Egs=-0.801984342436
Kept: 575 out of 1857, ratio=0.31
Memory used: 34 MB  Time elapsed: 76 s

 ***** [NRG] Iteration 66/67 (scale 1.6795180748321155e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260571151 nrkept=498 nrkeptmult=1142 
Total energy=-3.88890975695  Egs=-1.60302988972
Kept: 498 out of 2154, ratio=0.231
Memory used: 34 MB  Time elapsed: 76 s

 ***** [NRG] Iteration 67/67 (scale 1.1875986198391643e-10) ***** 
[1] xi(66)=1 zeta(67)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 12 states
Emax=10.837639995 nrkept=1857 nrkeptmult=4568 
Total energy=-3.88890975705  Egs=-0.80198374699
Kept: 1857 out of 1857, ratio=1
Memory used: 37 MB  Time elapsed: 76 s

** Iteration completed.


Total energy:             -3.888909757046869
TKW=2.14348530017e-05
Grand canonical Z=              10.02886540076087
Zft (last shell)=               10.02884741207829
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: 	SIAM	 Variant: 		 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.5
>>  ops=	A_d
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=6
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||||> (z)
Number of invariant subspaces: 6
(-2 1) 1 states
0.0367538100397 
(-1 2) 2 states
0.0174462807917 0.0525956033849 
(0 1) 3 states
0 0.0367538100397 0.0700418841766 
(0 3) 1 states
0.0332880741369 
(1 2) 2 states
0.0174462807917 0.0525956033849 
(2 1) 1 states
0.0367538100397 
Number of states (multiplicity taken into account): 16

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

Computing the following spectra:
Spectrum: A_d-A_d FDM

 ***** [DMNRG] Iteration 1/67 (scale 1.0201394465967895) ***** 
[1] xi(0)=0.583569837945 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=2.34125880272 nrkept=35 nrkeptmult=64 
Kept: 35 out of 35, ratio=1
Memory used: 35 MB  Time elapsed: 83 s

 ***** [DMNRG] Iteration 2/67 (scale 0.72134752044448169) ***** 
[1] xi(1)=0.776754544846 zeta(2)=0
Stats: nr=15 min=1 max=20 sum3=30066
Emax=4.6306144545 nrkept=126 nrkeptmult=256 
Kept: 126 out of 126, ratio=1
Memory used: 35 MB  Time elapsed: 83 s

 ***** [DMNRG] Iteration 3/67 (scale 0.51006972329839473) ***** 
[1] xi(2)=0.992903283056 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=9.70884141923 nrkept=462 nrkeptmult=1024 
Kept: 462 out of 462, ratio=1
Memory used: 35 MB  Time elapsed: 83 s

 ***** [DMNRG] Iteration 4/67 (scale 0.36067376022224085) ***** 
[1] xi(3)=1.04266411315 zeta(4)=0
Stats: nr=28 min=1 max=210 sum3=37226316
Emax=10.2107510943 nrkept=1259 nrkeptmult=3041 
Kept: 1259 out of 1716, ratio=0.734
Memory used: 38 MB  Time elapsed: 84 s

 ***** [DMNRG] Iteration 5/67 (scale 0.25503486164919736) ***** 
[1] xi(4)=1.01328740971 zeta(5)=0
Stats: nr=36 min=1 max=561 sum3=559501320
Safeguard: keep additional 3 states
Emax=10.2083031786 nrkept=1883 nrkeptmult=4688 
Kept: 1883 out of 4728, ratio=0.398
Memory used: 46 MB  Time elapsed: 84 s

 ***** [DMNRG] Iteration 6/67 (scale 0.18033688011112042) ***** 
[1] xi(5)=0.999788616061 zeta(6)=0
Stats: nr=42 min=1 max=717 sum3=1710937002
Safeguard: keep additional 3 states
Emax=10.2375168905 nrkept=2124 nrkeptmult=5286 
Kept: 2124 out of 7092, ratio=0.299
Memory used: 48 MB  Time elapsed: 86 s

 ***** [DMNRG] Iteration 7/67 (scale 0.12751743082459868) ***** 
[1] xi(6)=0.998344002082 zeta(7)=0
Stats: nr=44 min=1 max=910 sum3=2520131436
Safeguard: keep additional 2 states
Emax=10.0192263144 nrkept=1921 nrkeptmult=4746 
Kept: 1921 out of 8010, ratio=0.24
Memory used: 49 MB  Time elapsed: 89 s

 ***** [DMNRG] Iteration 8/67 (scale 0.090168440055560212) ***** 
[1] xi(7)=0.998902918682 zeta(8)=0
Stats: nr=42 min=1 max=745 sum3=1897784946
Safeguard: keep additional 2 states
Emax=10.0000486937 nrkept=1999 nrkeptmult=4973 
Kept: 1999 out of 7230, ratio=0.276
Memory used: 49 MB  Time elapsed: 91 s

 ***** [DMNRG] Iteration 9/67 (scale 0.063758715412299341) ***** 
[1] xi(8)=0.999398551363 zeta(9)=0
Stats: nr=44 min=1 max=863 sum3=2129089633
Safeguard: keep additional 1 states
Emax=10.0019670946 nrkept=1664 nrkeptmult=4134 
Kept: 1664 out of 7549, ratio=0.22
Memory used: 45 MB  Time elapsed: 93 s

 ***** [DMNRG] Iteration 10/67 (scale 0.045084220027780106) ***** 
[1] xi(9)=0.99968731683 zeta(10)=0
Stats: nr=42 min=1 max=652 sum3=1254088034
Emax=10.0473810799 nrkept=1873 nrkeptmult=4643 
Kept: 1873 out of 6284, ratio=0.298
Memory used: 48 MB  Time elapsed: 95 s

 ***** [DMNRG] Iteration 11/67 (scale 0.031879357706149671) ***** 
[1] xi(10)=0.9998407918 zeta(11)=0
Stats: nr=42 min=1 max=815 sum3=1780810507
Safeguard: keep additional 7 states
Emax=10.0127857568 nrkept=1638 nrkeptmult=4060 
Kept: 1638 out of 7057, ratio=0.232
Memory used: 47 MB  Time elapsed: 97 s

 ***** [DMNRG] Iteration 12/67 (scale 0.022542110013890053) ***** 
[1] xi(11)=0.99991969299 zeta(12)=0
Stats: nr=42 min=1 max=640 sum3=1188522258
Safeguard: keep additional 2 states
Emax=10.0015985499 nrkept=1757 nrkeptmult=4357 
Kept: 1757 out of 6180, ratio=0.284
Memory used: 47 MB  Time elapsed: 98 s

 ***** [DMNRG] Iteration 13/67 (scale 0.015939678853074835) ***** 
[1] xi(12)=0.999959672404 zeta(13)=0
Stats: nr=42 min=1 max=761 sum3=1454172497
Safeguard: keep additional 3 states
Emax=10.0035212517 nrkept=1554 nrkeptmult=3828 
Kept: 1554 out of 6617, ratio=0.235
Memory used: 45 MB  Time elapsed: 100 s

 ***** [DMNRG] Iteration 14/67 (scale 0.011271055006945026) ***** 
[1] xi(13)=0.999979792879 zeta(14)=0
Stats: nr=42 min=1 max=608 sum3=1018483162
Emax=10.0015336857 nrkept=1597 nrkeptmult=3941 
Kept: 1597 out of 5854, ratio=0.273
Memory used: 45 MB  Time elapsed: 101 s

 ***** [DMNRG] Iteration 15/67 (scale 0.0079698394265374176) ***** 
[1] xi(14)=0.999989885634 zeta(15)=0
Stats: nr=40 min=1 max=705 sum3=1134126634
Safeguard: keep additional 3 states
Emax=10.0011836559 nrkept=1516 nrkeptmult=3738 
Kept: 1516 out of 6016, ratio=0.252
Memory used: 44 MB  Time elapsed: 103 s

 ***** [DMNRG] Iteration 16/67 (scale 0.0056355275034725132) ***** 
[1] xi(15)=0.999994940118 zeta(16)=0
Stats: nr=40 min=1 max=596 sum3=957480812
Safeguard: keep additional 2 states
Emax=10.0009225943 nrkept=1465 nrkeptmult=3555 
Kept: 1465 out of 5708, ratio=0.257
Memory used: 43 MB  Time elapsed: 104 s

 ***** [DMNRG] Iteration 17/67 (scale 0.0039849197132687088) ***** 
[1] xi(16)=0.999997469385 zeta(17)=0
Stats: nr=36 min=1 max=658 sum3=902028083
Safeguard: keep additional 1 states
Emax=10.0002540846 nrkept=1327 nrkeptmult=3260 
Kept: 1327 out of 5513, ratio=0.241
Memory used: 43 MB  Time elapsed: 105 s

 ***** [DMNRG] Iteration 18/67 (scale 0.0028177637517362566) ***** 
[1] xi(17)=0.999998734524 zeta(18)=0
Stats: nr=40 min=1 max=529 sum3=666590356
Safeguard: keep additional 1 states
Emax=10.0000835745 nrkept=1281 nrkeptmult=3147 
Kept: 1281 out of 4990, ratio=0.257
Memory used: 45 MB  Time elapsed: 106 s

 ***** [DMNRG] Iteration 19/67 (scale 0.0019924598566343544) ***** 
[1] xi(18)=0.99999936722 zeta(19)=0
Stats: nr=34 min=1 max=579 sum3=616598005
Safeguard: keep additional 1 states
Emax=10.0009765181 nrkept=1154 nrkeptmult=2870 
Kept: 1154 out of 4825, ratio=0.239
Memory used: 48 MB  Time elapsed: 107 s

 ***** [DMNRG] Iteration 20/67 (scale 0.0014088818758681283) ***** 
[1] xi(19)=0.999999683599 zeta(20)=0
Stats: nr=40 min=1 max=462 sum3=445382958
Safeguard: keep additional 1 states
Emax=10.0052397336 nrkept=1096 nrkeptmult=2722 
Kept: 1096 out of 4356, ratio=0.252
Memory used: 48 MB  Time elapsed: 108 s

 ***** [DMNRG] Iteration 21/67 (scale 0.0009962299283171772) ***** 
[1] xi(20)=0.999999841797 zeta(21)=0
Stats: nr=34 min=1 max=500 sum3=395810408
Emax=10.0046911021 nrkept=963 nrkeptmult=2416 
Kept: 963 out of 4136, ratio=0.233
Memory used: 46 MB  Time elapsed: 109 s

 ***** [DMNRG] Iteration 22/67 (scale 0.00070444093793406415) ***** 
[1] xi(21)=0.999999920898 zeta(22)=0
Stats: nr=40 min=1 max=391 sum3=266535976
Safeguard: keep additional 2 states
Emax=10.0086552413 nrkept=895 nrkeptmult=2217 
Kept: 895 out of 3634, ratio=0.246
Memory used: 48 MB  Time elapsed: 109 s

 ***** [DMNRG] Iteration 23/67 (scale 0.0004981149641585886) ***** 
[1] xi(22)=0.999999960449 zeta(23)=0
Stats: nr=34 min=1 max=417 sum3=224530091
Safeguard: keep additional 2 states
Emax=10.0005993028 nrkept=868 nrkeptmult=2218 
Kept: 868 out of 3377, ratio=0.257
Memory used: 46 MB  Time elapsed: 110 s

 ***** [DMNRG] Iteration 24/67 (scale 0.00035222046896703208) ***** 
[1] xi(23)=0.999999980224 zeta(24)=0
Stats: nr=40 min=1 max=354 sum3=191343178
Safeguard: keep additional 1 states
Emax=10.0076766313 nrkept=823 nrkeptmult=2039 
Kept: 823 out of 3286, ratio=0.25
Memory used: 46 MB  Time elapsed: 110 s

 ***** [DMNRG] Iteration 25/67 (scale 0.0002490574820792943) ***** 
[1] xi(24)=0.999999990112 zeta(25)=0
Stats: nr=34 min=1 max=383 sum3=172023467
Emax=10.0036971781 nrkept=816 nrkeptmult=2066 
Kept: 816 out of 3101, ratio=0.263
Memory used: 46 MB  Time elapsed: 111 s

 ***** [DMNRG] Iteration 26/67 (scale 0.00017611023448351604) ***** 
[1] xi(25)=0.999999995056 zeta(26)=0
Stats: nr=38 min=1 max=333 sum3=161834170
Safeguard: keep additional 3 states
Emax=10.0237130106 nrkept=783 nrkeptmult=1957 
Kept: 783 out of 3082, ratio=0.254
Memory used: 46 MB  Time elapsed: 111 s

 ***** [DMNRG] Iteration 27/67 (scale 0.00012452874103964715) ***** 
[1] xi(26)=0.999999997528 zeta(27)=0
Stats: nr=34 min=1 max=359 sum3=144739908
Emax=10.0104228215 nrkept=798 nrkeptmult=2024 
Kept: 798 out of 2952, ratio=0.27
Memory used: 46 MB  Time elapsed: 112 s

 ***** [DMNRG] Iteration 28/67 (scale 8.8055117241758019e-05) ***** 
[1] xi(27)=0.999999998764 zeta(28)=0
Stats: nr=38 min=1 max=327 sum3=150905436
Safeguard: keep additional 1 states
Emax=10.0147851711 nrkept=758 nrkeptmult=1918 
Kept: 758 out of 3012, ratio=0.252
Memory used: 46 MB  Time elapsed: 112 s

 ***** [DMNRG] Iteration 29/67 (scale 6.2264370519823575e-05) ***** 
[1] xi(28)=0.999999999382 zeta(29)=0
Stats: nr=34 min=1 max=346 sum3=131798132
Safeguard: keep additional 1 states
Emax=10.0123129494 nrkept=765 nrkeptmult=1910 
Kept: 765 out of 2864, ratio=0.267
Memory used: 46 MB  Time elapsed: 112 s

 ***** [DMNRG] Iteration 30/67 (scale 4.402755862087901e-05) ***** 
[1] xi(29)=0.999999999691 zeta(30)=0
Stats: nr=38 min=1 max=316 sum3=137692612
Safeguard: keep additional 1 states
Emax=10.0127816539 nrkept=715 nrkeptmult=1779 
Kept: 715 out of 2884, ratio=0.248
Memory used: 46 MB  Time elapsed: 113 s

 ***** [DMNRG] Iteration 31/67 (scale 3.1132185259911788e-05) ***** 
[1] xi(30)=0.999999999846 zeta(31)=0
Stats: nr=33 min=1 max=330 sum3=112532256
Safeguard: keep additional 2 states
Emax=10.0193565251 nrkept=722 nrkeptmult=1782 
Kept: 722 out of 2694, ratio=0.268
Memory used: 46 MB  Time elapsed: 113 s

 ***** [DMNRG] Iteration 32/67 (scale 2.2013779310439505e-05) ***** 
[1] xi(31)=0.999999999923 zeta(32)=0
Stats: nr=36 min=1 max=301 sum3=119695024
Safeguard: keep additional 3 states
Emax=10.0171895907 nrkept=653 nrkeptmult=1597 
Kept: 653 out of 2722, ratio=0.24
Memory used: 45 MB  Time elapsed: 113 s

 ***** [DMNRG] Iteration 33/67 (scale 1.5566092629955894e-05) ***** 
[1] xi(32)=0.999999999961 zeta(33)=0
Stats: nr=33 min=1 max=310 sum3=90015037
Safeguard: keep additional 4 states
Emax=10.0121834087 nrkept=679 nrkeptmult=1672 
Kept: 679 out of 2461, ratio=0.276
Memory used: 45 MB  Time elapsed: 114 s

 ***** [DMNRG] Iteration 34/67 (scale 1.1006889655219752e-05) ***** 
[1] xi(33)=0.999999999981 zeta(34)=0
Stats: nr=36 min=1 max=282 sum3=98996162
Safeguard: keep additional 1 states
Emax=10.0086245136 nrkept=592 nrkeptmult=1416 
Kept: 592 out of 2558, ratio=0.231
Memory used: 45 MB  Time elapsed: 114 s

 ***** [DMNRG] Iteration 35/67 (scale 7.7830463149779469e-06) ***** 
[1] xi(34)=0.99999999999 zeta(35)=0
Stats: nr=33 min=1 max=287 sum3=68884364
Safeguard: keep additional 1 states
Emax=10.0449491878 nrkept=603 nrkeptmult=1444 
Kept: 603 out of 2222, ratio=0.271
Memory used: 44 MB  Time elapsed: 114 s

 ***** [DMNRG] Iteration 36/67 (scale 5.5034448276098762e-06) ***** 
[1] xi(35)=0.999999999995 zeta(36)=0
Stats: nr=36 min=1 max=256 sum3=72549386
Emax=10.0041334351 nrkept=535 nrkeptmult=1261 
Kept: 535 out of 2264, ratio=0.236
Memory used: 44 MB  Time elapsed: 114 s

 ***** [DMNRG] Iteration 37/67 (scale 3.8915231574889735e-06) ***** 
[1] xi(36)=0.999999999998 zeta(37)=0
Stats: nr=33 min=1 max=262 sum3=51852835
Safeguard: keep additional 1 states
Emax=10.084200852 nrkept=569 nrkeptmult=1328 
Kept: 569 out of 2005, ratio=0.284
Memory used: 44 MB  Time elapsed: 115 s

 ***** [DMNRG] Iteration 38/67 (scale 2.7517224138049381e-06) ***** 
[1] xi(37)=0.999999999999 zeta(38)=0
Stats: nr=34 min=1 max=245 sum3=62083694
Safeguard: keep additional 1 states
Emax=10.025184247 nrkept=505 nrkeptmult=1165 
Kept: 505 out of 2132, ratio=0.237
Memory used: 44 MB  Time elapsed: 115 s

 ***** [DMNRG] Iteration 39/67 (scale 1.9457615787444867e-06) ***** 
[1] xi(38)=0.999999999999 zeta(39)=0
Stats: nr=33 min=1 max=249 sum3=44336516
Emax=10.4459103367 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1886, ratio=0.301
Memory used: 43 MB  Time elapsed: 115 s

 ***** [DMNRG] Iteration 40/67 (scale 1.3758612069024691e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.1981441091 nrkept=499 nrkeptmult=1149 
Kept: 499 out of 2124, ratio=0.235
Memory used: 44 MB  Time elapsed: 115 s

 ***** [DMNRG] Iteration 41/67 (scale 9.7288078937224336e-07) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.6685417454 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1862, ratio=0.305
Memory used: 43 MB  Time elapsed: 115 s

 ***** [DMNRG] Iteration 42/67 (scale 6.8793060345123453e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.4002552911 nrkept=499 nrkeptmult=1149 
Kept: 499 out of 2124, ratio=0.235
Memory used: 44 MB  Time elapsed: 116 s

 ***** [DMNRG] Iteration 43/67 (scale 4.8644039468612168e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.7537218724 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1862, ratio=0.305
Memory used: 43 MB  Time elapsed: 116 s

 ***** [DMNRG] Iteration 44/67 (scale 3.4396530172561726e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5011034821 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 116 s

 ***** [DMNRG] Iteration 45/67 (scale 2.4322019734306084e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.7976645703 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 116 s

 ***** [DMNRG] Iteration 46/67 (scale 1.7198265086280863e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5135166719 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 116 s

 ***** [DMNRG] Iteration 47/67 (scale 1.2161009867153042e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8198704514 nrkept=567 nrkeptmult=1322 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 116 s

 ***** [DMNRG] Iteration 48/67 (scale 8.5991325431404316e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5197766792 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 117 s

 ***** [DMNRG] Iteration 49/67 (scale 6.080504933576521e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8304354411 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 117 s

 ***** [DMNRG] Iteration 50/67 (scale 4.2995662715702158e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5229304884 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 117 s

 ***** [DMNRG] Iteration 51/67 (scale 3.0402524667882605e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8340453147 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 117 s

 ***** [DMNRG] Iteration 52/67 (scale 2.1497831357851079e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5245060286 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 117 s

 ***** [DMNRG] Iteration 53/67 (scale 1.5201262333941303e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8358467908 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 118 s

 ***** [DMNRG] Iteration 54/67 (scale 1.0748915678925539e-08) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5252946654 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 118 s

 ***** [DMNRG] Iteration 55/67 (scale 7.6006311669706513e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8367466745 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 118 s

 ***** [DMNRG] Iteration 56/67 (scale 5.3744578394627697e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5256892009 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 118 s

 ***** [DMNRG] Iteration 57/67 (scale 3.8003155834853256e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8371964056 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 43 MB  Time elapsed: 118 s

 ***** [DMNRG] Iteration 58/67 (scale 2.6872289197313849e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5258865229 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 119 s

 ***** [DMNRG] Iteration 59/67 (scale 1.9001577917426628e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8374212188 nrkept=567 nrkeptmult=1324 
Kept: 567 out of 1857, ratio=0.305
Memory used: 44 MB  Time elapsed: 119 s

 ***** [DMNRG] Iteration 60/67 (scale 1.3436144598656924e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5259851976 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2124, ratio=0.234
Memory used: 44 MB  Time elapsed: 119 s

 ***** [DMNRG] Iteration 61/67 (scale 9.5007889587133141e-10) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 2 states
Emax=10.8375360849 nrkept=569 nrkeptmult=1330 
Kept: 569 out of 1857, ratio=0.306
Memory used: 44 MB  Time elapsed: 119 s

 ***** [DMNRG] Iteration 62/67 (scale 6.7180722993284622e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=34 min=1 max=244 sum3=61697330
Emax=10.5260345383 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2132, ratio=0.234
Memory used: 45 MB  Time elapsed: 119 s

 ***** [DMNRG] Iteration 63/67 (scale 4.750394479356657e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8376027884 nrkept=575 nrkeptmult=1340 
Kept: 575 out of 1857, ratio=0.31
Memory used: 43 MB  Time elapsed: 120 s

 ***** [DMNRG] Iteration 64/67 (scale 3.3590361496642311e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260447789 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2154, ratio=0.231
Memory used: 44 MB  Time elapsed: 120 s

 ***** [DMNRG] Iteration 65/67 (scale 2.3751972396783285e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8376184865 nrkept=575 nrkeptmult=1340 
Kept: 575 out of 1857, ratio=0.31
Memory used: 44 MB  Time elapsed: 120 s

 ***** [DMNRG] Iteration 66/67 (scale 1.6795180748321155e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260571151 nrkept=498 nrkeptmult=1142 
Kept: 498 out of 2154, ratio=0.231
Memory used: 45 MB  Time elapsed: 120 s

 ***** [DMNRG] Iteration 67/67 (scale 1.1875986198391643e-10) ***** 
[1] xi(66)=1 zeta(67)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 12 states
Emax=10.837639995 nrkept=1857 nrkeptmult=4568 
Kept: 1857 out of 1857, ratio=1
Memory used: 47 MB  Time elapsed: 120 s

** Iteration completed.

Spectrum: A_d-A_d FDM
[A_d-A_d] pos=0.50000001681 neg=0.499999983189 sum= 0.999999999999
mu1=-2.55568136098e-15 mu2=0.000664290165804 mu3=1.54498809579e-18 mu4=0.000226976058064


Memory usage report
===================
  after diag: 90576 kB
after recalc: 102604 kB
 after split: 116012 kB
  after trim: 50920 kB
 after trunc: 90744 kB
      diag_h: 57388 kB

Peak usage: 116012 kB

Timing report
=============
         All: 121 s
          DM: 1 s
         FDM: 5 s
        diag: 65 s
    recalc d: 8 s
      spec d: 21 s
       Other: 19 s
