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 May 31 2013 at 17:57:13

Parallelization using MPI: Running on 1 processors.

Using Intel MKL library 10.3.9
Processor optimization: Intel(R) Core(TM) i7 Processor
max_threads=8 blas_max_threads=8 dynamic=1

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.25
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.5
omega0=1e-99
omega0_ratio=1
brcut=0
loggauss_b=0.3
gauss_b=0.1
limitLL=0
eta=0
finite=1
bins=100
discard_trim=1e-16
discard_immediately=1e-16
NN1=0
NN2even=1
NN2avg=0
NNtanh=0
goodE=2.1
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 SZd
specs=
specd=A_d-A_d-u A_d-A_d-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=1
cfs=1
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=1
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: QSZ
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 	MAGFIELD	 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 	 extraB=	0.0003
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d SZd
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=9
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||SZd||> (s)
Reading <||||> (z)

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

Number of invariant subspaces: 9
(-2 1) 1 states
0.0309062904041 
(-1 0) 2 states
0.0146225690892 0.0441882570753 
(-1 2) 2 states
0.0147186194033 0.044267066255 
(0 -1) 1 states
0.0279045357604 
(0 1) 4 states
0 0.0279919912572 0.0309062904041 0.0588982301615 
(0 3) 1 states
0.0280793952542 
(1 0) 2 states
0.0146225690892 0.0441882570753 
(1 2) 2 states
0.0147186194033 0.044267066255 
(2 1) 1 states
0.0309062904041 
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-d FT
Spectrum: A_d-A_d-u FT

Before NRG iteration (N=-1)
<SZd>=-2.18884963691e-05

 ***** [NRG] Iteration 1/67 (scale 1.2131570881878404) ***** 
[1] xi(0)=0.473132378817 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=1.89858626849 nrkept=64 nrkeptmult=64 
Total energy=-1.14921655241  Egs=-0.907707483067
Kept: 64 out of 64, ratio=1
<SZd>=-3.0972594458e-05
Memory used: 11 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/67 (scale 0.85783160370214839) ***** 
[1] xi(1)=0.631518937152 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=3.76220868151 nrkept=256 nrkeptmult=256 
Total energy=-1.61117080017  Egs=-0.538513906195
Kept: 256 out of 256, ratio=1
<SZd>=-4.38361947131e-05
Memory used: 12 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/67 (scale 0.6065785440939202) ***** 
[1] xi(2)=0.974788923894 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Emax=8.48719165627 nrkept=1024 nrkeptmult=1024 
Total energy=-2.57164892023  Egs=-1.5834356975
Kept: 1024 out of 1024, ratio=1
<SZd>=-6.20603756185e-05
Memory used: 16 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/67 (scale 0.4289158018510742) ***** 
[1] xi(3)=1.2027371305 zeta(4)=0
Stats: nr=49 min=1 max=400 sum3=230553856
Safeguard: keep additional 3 states
Emax=10.316350082 nrkept=3540 nrkeptmult=3540 
Total energy=-3.05731589939  Egs=-1.13231309516
Kept: 3540 out of 4096, ratio=0.864
<SZd>=-8.79025854902e-05
Memory used: 55 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/67 (scale 0.3032892720469601) ***** 
[1] xi(4)=1.12253968187 zeta(5)=0
Stats: nr=64 min=1 max=1038 sum3=6797276016
Safeguard: keep additional 2 states
Emax=10.3387537306 nrkept=5615 nrkeptmult=5615 
Total energy=-3.55022171117  Egs=-1.62520028636
Kept: 5615 out of 14160, ratio=0.397
<SZd>=-0.000125699957759
Memory used: 106 MB  Time elapsed: 8 s

 ***** [NRG] Iteration 6/67 (scale 0.2144579009255371) ***** 
[1] xi(5)=1.03243378845 zeta(6)=0
Stats: nr=77 min=2 max=1725 sum3=19012513708
Safeguard: keep additional 5 states
Emax=10.0283233223 nrkept=5542 nrkeptmult=5542 
Total energy=-3.74365715335  Egs=-0.90197396015
Kept: 5542 out of 22460, ratio=0.247
<SZd>=-0.000178711579131
Memory used: 117 MB  Time elapsed: 27 s

 ***** [NRG] Iteration 7/67 (scale 0.15164463602348005) ***** 
[1] xi(6)=1.00829687067 zeta(7)=0
Stats: nr=80 min=1 max=1563 sum3=23906874634
Safeguard: keep additional 11 states
Emax=10.011517717 nrkept=5136 nrkeptmult=5136 
Total energy=-3.97218262747  Egs=-1.50698026722
Kept: 5136 out of 22168, ratio=0.232
<SZd>=-0.000250761994046
Memory used: 113 MB  Time elapsed: 47 s

 ***** [NRG] Iteration 8/67 (scale 0.10722895046276855) ***** 
[1] xi(7)=1.00300214986 zeta(8)=0
Stats: nr=77 min=1 max=1622 sum3=18823528974
Safeguard: keep additional 8 states
Emax=10.0368187077 nrkept=5273 nrkeptmult=5273 
Total energy=-4.06890437173  Egs=-0.902011479591
Kept: 5273 out of 20544, ratio=0.257
<SZd>=-0.000356707026497
Memory used: 121 MB  Time elapsed: 63 s

 ***** [NRG] Iteration 9/67 (scale 0.075822318011740025) ***** 
[1] xi(8)=1.00131950393 zeta(9)=0
Stats: nr=80 min=1 max=1481 sum3=20427425036
Safeguard: keep additional 1 states
Emax=10.0014819469 nrkept=4170 nrkeptmult=4170 
Total energy=-4.18050965839  Egs=-1.4719318742
Kept: 4170 out of 21092, ratio=0.198
<SZd>=-0.000512472170741
Memory used: 109 MB  Time elapsed: 79 s

 ***** [NRG] Iteration 10/67 (scale 0.053614475231384275) ***** 
[1] xi(9)=1.00062206904 zeta(10)=0
Stats: nr=77 min=1 max=1338 sum3=10427379846
Emax=10.018577483 nrkept=4641 nrkeptmult=4641 
Total energy=-4.23013290916  Egs=-0.925556961202
Kept: 4641 out of 16680, ratio=0.278
<SZd>=-0.000726798815512
Memory used: 121 MB  Time elapsed: 90 s

 ***** [NRG] Iteration 11/67 (scale 0.037911159005870013) ***** 
[1] xi(10)=1.00030223414 zeta(11)=0
Stats: nr=76 min=1 max=1324 sum3=14394708342
Safeguard: keep additional 2 states
Emax=10.0000431522 nrkept=4023 nrkeptmult=4023 
Total energy=-4.28463809571  Egs=-1.43770826254
Kept: 4023 out of 18564, ratio=0.217
<SZd>=-0.0010684782779
Memory used: 117 MB  Time elapsed: 103 s

 ***** [NRG] Iteration 12/67 (scale 0.026807237615692137) ***** 
[1] xi(11)=1.00014898132 zeta(12)=0
Stats: nr=77 min=1 max=1288 sum3=9285140442
Safeguard: keep additional 4 states
Emax=10.0013605611 nrkept=4453 nrkeptmult=4453 
Total energy=-4.31015469625  Egs=-0.951854902395
Kept: 4453 out of 16092, ratio=0.277
<SZd>=-0.00150211753461
Memory used: 128 MB  Time elapsed: 113 s

 ***** [NRG] Iteration 13/67 (scale 0.018955579502935006) ***** 
[1] xi(12)=1.00007396421 zeta(13)=0
Stats: nr=76 min=1 max=1270 sum3=12681406924
Emax=10.0017093305 nrkept=3821 nrkeptmult=3821 
Total energy=-4.33676599163  Egs=-1.4038766459
Kept: 3821 out of 17812, ratio=0.215
<SZd>=-0.00217840522978
Memory used: 120 MB  Time elapsed: 124 s

 ***** [NRG] Iteration 14/67 (scale 0.013403618807846069) ***** 
[1] xi(13)=1.0000368514 zeta(14)=0
Stats: nr=77 min=1 max=1231 sum3=8065485290
Safeguard: keep additional 3 states
Emax=10.000017474 nrkept=3998 nrkeptmult=3998 
Total energy=-4.34984873004  Egs=-0.976060166882
Kept: 3998 out of 15284, ratio=0.262
<SZd>=-0.00320024310368
Memory used: 122 MB  Time elapsed: 133 s

 ***** [NRG] Iteration 15/67 (scale 0.0094777897514675032) ***** 
[1] xi(14)=1.00001839314 zeta(15)=0
Stats: nr=72 min=1 max=1155 sum3=9440440724
Emax=10.0200058669 nrkept=3781 nrkeptmult=3781 
Total energy=-4.36286633971  Egs=-1.37348580352
Kept: 3781 out of 15992, ratio=0.236
<SZd>=-0.00467118924566
Memory used: 121 MB  Time elapsed: 142 s

 ***** [NRG] Iteration 16/67 (scale 0.0067018094039230343) ***** 
[1] xi(15)=1.00000918844 zeta(16)=0
Stats: nr=77 min=1 max=1216 sum3=7792061218
Safeguard: keep additional 1 states
Emax=10.0010794918 nrkept=3636 nrkeptmult=3636 
Total energy=-4.36953605796  Egs=-0.995211569083
Kept: 3636 out of 15124, ratio=0.24
<SZd>=-0.00694203024837
Memory used: 117 MB  Time elapsed: 149 s

 ***** [NRG] Iteration 17/67 (scale 0.0047388948757337516) ***** 
[1] xi(16)=1.00000459219 zeta(17)=0
Stats: nr=64 min=1 max=1067 sum3=7356201924
Safeguard: keep additional 2 states
Emax=10.0028773308 nrkept=3404 nrkeptmult=3404 
Total energy=-4.37592888287  Egs=-1.34901175863
Kept: 3404 out of 14544, ratio=0.234
<SZd>=-0.0103200744649
Memory used: 116 MB  Time elapsed: 157 s

 ***** [NRG] Iteration 18/67 (scale 0.0033509047019615172) ***** 
[1] xi(17)=1.00000229559 zeta(18)=0
Stats: nr=73 min=1 max=1108 sum3=5829697196
Safeguard: keep additional 1 states
Emax=10.0040982903 nrkept=3245 nrkeptmult=3245 
Total energy=-4.37930761719  Egs=-1.00830510528
Kept: 3245 out of 13616, ratio=0.238
<SZd>=-0.0153115894459
Memory used: 117 MB  Time elapsed: 163 s

 ***** [NRG] Iteration 19/67 (scale 0.0023694474378668758) ***** 
[1] xi(18)=1.00000114767 zeta(19)=0
Stats: nr=62 min=1 max=961 sum3=5346491798
Safeguard: keep additional 3 states
Emax=10.005321246 nrkept=2931 nrkeptmult=2931 
Total energy=-4.38246356221  Egs=-1.33193290984
Kept: 2931 out of 12980, ratio=0.226
<SZd>=-0.0224900331387
Memory used: 118 MB  Time elapsed: 169 s

 ***** [NRG] Iteration 20/67 (scale 0.0016754523509807586) ***** 
[1] xi(19)=1.0000005738 zeta(20)=0
Stats: nr=73 min=1 max=967 sum3=3822280566
Safeguard: keep additional 1 states
Emax=10.0012636476 nrkept=2775 nrkeptmult=2775 
Total energy=-4.3841644307  Egs=-1.01516971845
Kept: 2775 out of 11724, ratio=0.237
<SZd>=-0.0324764646067
Memory used: 120 MB  Time elapsed: 173 s

 ***** [NRG] Iteration 21/67 (scale 0.0011847237189334379) ***** 
[1] xi(20)=1.00000028689 zeta(21)=0
Stats: nr=62 min=1 max=836 sum3=3396923856
Safeguard: keep additional 1 states
Emax=10.002424883 nrkept=2421 nrkeptmult=2421 
Total energy=-4.38573260056  Egs=-1.32365869667
Kept: 2421 out of 11100, ratio=0.218
<SZd>=-0.0459091603572
Memory used: 117 MB  Time elapsed: 177 s

 ***** [NRG] Iteration 22/67 (scale 0.00083772617549037929) ***** 
[1] xi(21)=1.00000014344 zeta(22)=0
Stats: nr=73 min=1 max=807 sum3=2211355914
Safeguard: keep additional 1 states
Emax=10.0001934725 nrkept=2254 nrkeptmult=2254 
Total energy=-4.38658347135  Egs=-1.01569082539
Kept: 2254 out of 9684, ratio=0.233
<SZd>=-0.06378240338
Memory used: 116 MB  Time elapsed: 180 s

 ***** [NRG] Iteration 23/67 (scale 0.00059236185946671895) ***** 
[1] xi(22)=1.00000007172 zeta(23)=0
Stats: nr=62 min=1 max=695 sum3=1844095000
Safeguard: keep additional 3 states
Emax=10.0056295614 nrkept=2102 nrkeptmult=2102 
Total energy=-4.38736828321  Egs=-1.32488587696
Kept: 2102 out of 9016, ratio=0.233
<SZd>=-0.0871816530402
Memory used: 119 MB  Time elapsed: 183 s

 ***** [NRG] Iteration 24/67 (scale 0.00041886308774518965) ***** 
[1] xi(23)=1.00000003586 zeta(24)=0
Stats: nr=71 min=1 max=695 sum3=1431129056
Safeguard: keep additional 1 states
Emax=10.0056758785 nrkept=2035 nrkeptmult=2035 
Total energy=-4.38779186111  Egs=-1.01125620868
Kept: 2035 out of 8408, ratio=0.242
<SZd>=-0.11746425168
Memory used: 119 MB  Time elapsed: 185 s

 ***** [NRG] Iteration 25/67 (scale 0.00029618092973335947) ***** 
[1] xi(24)=1.00000001793 zeta(25)=0
Stats: nr=62 min=1 max=630 sum3=1364926552
Safeguard: keep additional 1 states
Emax=10.0036573789 nrkept=1893 nrkeptmult=1893 
Total energy=-4.38818667341  Egs=-1.33301053968
Kept: 1893 out of 8140, ratio=0.233
<SZd>=-0.155360262611
Memory used: 118 MB  Time elapsed: 187 s

 ***** [NRG] Iteration 26/67 (scale 0.00020943154387259482) ***** 
[1] xi(25)=1.00000000897 zeta(26)=0
Stats: nr=69 min=1 max=641 sum3=1086480156
Safeguard: keep additional 2 states
Emax=10.0133600457 nrkept=1900 nrkeptmult=1900 
Total energy=-4.38839716726  Egs=-1.00507233932
Kept: 1900 out of 7572, ratio=0.251
<SZd>=-0.200016972709
Memory used: 119 MB  Time elapsed: 188 s

 ***** [NRG] Iteration 27/67 (scale 0.00014809046486667974) ***** 
[1] xi(26)=1.00000000448 zeta(27)=0
Stats: nr=61 min=1 max=599 sum3=1148042122
Emax=10.0006479317 nrkept=1660 nrkeptmult=1660 
Total energy=-4.38859599461  Egs=-1.34260738169
Kept: 1660 out of 7600, ratio=0.218
<SZd>=-0.248212643382
Memory used: 119 MB  Time elapsed: 190 s

 ***** [NRG] Iteration 28/67 (scale 0.00010471577193629741) ***** 
[1] xi(27)=1.00000000224 zeta(28)=0
Stats: nr=67 min=1 max=561 sum3=742975444
Safeguard: keep additional 3 states
Emax=10.0014214886 nrkept=1629 nrkeptmult=1629 
Total energy=-4.38870068726  Egs=-0.999779178318
Kept: 1629 out of 6640, ratio=0.245
<SZd>=-0.293089985991
Memory used: 120 MB  Time elapsed: 191 s

 ***** [NRG] Iteration 29/67 (scale 7.4045232433339868e-05) ***** 
[1] xi(28)=1.00000000112 zeta(29)=0
Stats: nr=61 min=1 max=533 sum3=734075190
Emax=10.01068027 nrkept=1421 nrkeptmult=1421 
Total energy=-4.38880063884  Egs=-1.3498719313
Kept: 1421 out of 6516, ratio=0.218
<SZd>=-0.327647757117
Memory used: 119 MB  Time elapsed: 192 s

 ***** [NRG] Iteration 30/67 (scale 5.2357885968148706e-05) ***** 
[1] xi(29)=1.00000000056 zeta(30)=0
Stats: nr=59 min=1 max=487 sum3=494943710
Emax=10.0158508548 nrkept=1349 nrkeptmult=1349 
Total energy=-4.38885279864  Egs=-0.996216610899
Kept: 1349 out of 5684, ratio=0.237
<SZd>=-0.348402531819
Memory used: 119 MB  Time elapsed: 193 s

 ***** [NRG] Iteration 31/67 (scale 3.7022616216669934e-05) ***** 
[1] xi(30)=1.00000000028 zeta(31)=0
Stats: nr=57 min=1 max=474 sum3=452028422
Safeguard: keep additional 1 states
Emax=10.0077809736 nrkept=1190 nrkeptmult=1190 
Total energy=-4.3889029366  Egs=-1.35425229831
Kept: 1190 out of 5396, ratio=0.221
<SZd>=-0.358167305283
Memory used: 120 MB  Time elapsed: 194 s

 ***** [NRG] Iteration 32/67 (scale 2.6178942984074353e-05) ***** 
[1] xi(31)=1.00000000014 zeta(32)=0
Stats: nr=58 min=1 max=407 sum3=310550402
Safeguard: keep additional 1 states
Emax=10.0061179257 nrkept=1180 nrkeptmult=1180 
Total energy=-4.38892896183  Egs=-0.994128547913
Kept: 1180 out of 4760, ratio=0.248
<SZd>=-0.361948355024
Memory used: 122 MB  Time elapsed: 195 s

 ***** [NRG] Iteration 33/67 (scale 1.8511308108334967e-05) ***** 
[1] xi(32)=1.00000000007 zeta(33)=0
Stats: nr=57 min=1 max=432 sum3=318384184
Emax=10.0022900413 nrkept=1094 nrkeptmult=1094 
Total energy=-4.38895407469  Egs=-1.35662259834
Kept: 1094 out of 4720, ratio=0.232
<SZd>=-0.363419000719
Memory used: 122 MB  Time elapsed: 195 s

 ***** [NRG] Iteration 34/67 (scale 1.3089471492037176e-05) ***** 
[1] xi(33)=1.00000000004 zeta(34)=0
Stats: nr=58 min=1 max=373 sum3=241218758
Safeguard: keep additional 2 states
Emax=10.0033430926 nrkept=1125 nrkeptmult=1125 
Total energy=-4.3889670725  Egs=-0.992996886873
Kept: 1125 out of 4376, ratio=0.257
<SZd>=-0.364050260551
Memory used: 123 MB  Time elapsed: 196 s

 ***** [NRG] Iteration 35/67 (scale 9.2556540541674836e-06) ***** 
[1] xi(34)=1.00000000002 zeta(35)=0
Stats: nr=57 min=1 max=415 sum3=278950578
Safeguard: keep additional 1 states
Emax=10.000505498 nrkept=1075 nrkeptmult=1075 
Total energy=-4.38897964027  Egs=-1.35784779152
Kept: 1075 out of 4500, ratio=0.239
<SZd>=-0.364348238061
Memory used: 124 MB  Time elapsed: 197 s

 ***** [NRG] Iteration 36/67 (scale 6.5447357460185882e-06) ***** 
[1] xi(35)=1.00000000001 zeta(36)=0
Stats: nr=58 min=1 max=371 sum3=234091648
Safeguard: keep additional 1 states
Emax=10.0395198524 nrkept=1116 nrkeptmult=1116 
Total energy=-4.38898613531  Egs=-0.992408123421
Kept: 1116 out of 4300, ratio=0.26
<SZd>=-0.364492837418
Memory used: 124 MB  Time elapsed: 197 s

 ***** [NRG] Iteration 37/67 (scale 4.6278270270837418e-06) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=57 min=1 max=412 sum3=271853526
Safeguard: keep additional 3 states
Emax=10.000534209 nrkept=1073 nrkeptmult=1073 
Total energy=-4.38899242208  Egs=-1.35846939527
Kept: 1073 out of 4464, ratio=0.24
<SZd>=-0.364564208639
Memory used: 124 MB  Time elapsed: 198 s

 ***** [NRG] Iteration 38/67 (scale 3.2723678730092941e-06) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=58 min=1 max=371 sum3=232464242
Safeguard: keep additional 2 states
Emax=10.0626513977 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38899566862  Egs=-0.992107940214
Kept: 1117 out of 4292, ratio=0.26
<SZd>=-0.364599770394
Memory used: 124 MB  Time elapsed: 198 s

 ***** [NRG] Iteration 39/67 (scale 2.3139135135418709e-06) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=57 min=1 max=413 sum3=273706324
Safeguard: keep additional 1 states
Emax=10.0034837992 nrkept=1065 nrkeptmult=1065 
Total energy=-4.38899881272  Egs=-1.3587822985
Kept: 1065 out of 4468, ratio=0.238
<SZd>=-0.3646173543
Memory used: 125 MB  Time elapsed: 199 s

 ***** [NRG] Iteration 40/67 (scale 1.6361839365046471e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=58 min=1 max=371 sum3=231782964
Emax=10.070661664 nrkept=1115 nrkeptmult=1115 
Total energy=-4.38900043575  Egs=-0.991956395649
Kept: 1115 out of 4260, ratio=0.262
<SZd>=-0.364626349441
Memory used: 124 MB  Time elapsed: 199 s

 ***** [NRG] Iteration 41/67 (scale 1.1569567567709354e-06) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=57 min=1 max=413 sum3=272139644
Emax=10.0085129322 nrkept=1060 nrkeptmult=1060 
Total energy=-4.38900200798  Egs=-1.35893925491
Kept: 1060 out of 4460, ratio=0.238
<SZd>=-0.364630762817
Memory used: 125 MB  Time elapsed: 200 s

 ***** [NRG] Iteration 42/67 (scale 8.1809196825232353e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=58 min=1 max=371 sum3=231642436
Safeguard: keep additional 2 states
Emax=10.0736476623 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900281943  Egs=-0.991880260207
Kept: 1117 out of 4240, ratio=0.263
<SZd>=-0.36463301959
Memory used: 125 MB  Time elapsed: 201 s

 ***** [NRG] Iteration 43/67 (scale 5.7847837838546772e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0155357331 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900360559  Egs=-1.35901785667
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364634132954
Memory used: 126 MB  Time elapsed: 201 s

 ***** [NRG] Iteration 44/67 (scale 4.0904598412616176e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0749303156 nrkept=1117 nrkeptmult=1117 
Total energy=-4.3890040113  Egs=-0.991842101731
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364634703113
Memory used: 125 MB  Time elapsed: 202 s

 ***** [NRG] Iteration 45/67 (scale 2.8923918919273386e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0173921736 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900440439  Egs=-1.35905718812
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364634984378
Memory used: 128 MB  Time elapsed: 202 s

 ***** [NRG] Iteration 46/67 (scale 2.0452299206308088e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0755718035 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900460724  Egs=-0.991822999804
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635132479
Memory used: 127 MB  Time elapsed: 203 s

 ***** [NRG] Iteration 47/67 (scale 1.4461959459636693e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0182657594 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900480379  Egs=-1.35907686145
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364635204273
Memory used: 128 MB  Time elapsed: 203 s

 ***** [NRG] Iteration 48/67 (scale 1.0226149603154044e-07) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0758925952 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900490522  Egs=-0.991813443161
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635244067
Memory used: 127 MB  Time elapsed: 204 s

 ***** [NRG] Iteration 49/67 (scale 7.2309797298183465e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0186963227 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900500349  Egs=-1.35908670003
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364635262758
Memory used: 128 MB  Time elapsed: 204 s

 ***** [NRG] Iteration 50/67 (scale 5.113074801577022e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.076053004 nrkept=1117 nrkeptmult=1117 
Total energy=-4.3890050542  Egs=-0.991808663413
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635274091
Memory used: 128 MB  Time elapsed: 205 s

 ***** [NRG] Iteration 51/67 (scale 3.6154898649091733e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.018911094 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900510334  Egs=-1.3590916198
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364635279133
Memory used: 128 MB  Time elapsed: 206 s

 ***** [NRG] Iteration 52/67 (scale 2.556537400788511e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761332117 nrkept=1117 nrkeptmult=1117 
Total energy=-4.3890051287  Egs=-0.991806273174
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635282658
Memory used: 128 MB  Time elapsed: 206 s

 ***** [NRG] Iteration 53/67 (scale 1.8077449324545866e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190184086 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900515327  Egs=-1.35909407981
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364635284102
Memory used: 129 MB  Time elapsed: 207 s

 ***** [NRG] Iteration 54/67 (scale 1.2782687003942555e-08) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761733163 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900516595  Egs=-0.991805077956
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635285329
Memory used: 129 MB  Time elapsed: 207 s

 ***** [NRG] Iteration 55/67 (scale 9.0387246622729332e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190720535 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900517823  Egs=-1.35909530986
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364635285782
Memory used: 130 MB  Time elapsed: 208 s

 ***** [NRG] Iteration 56/67 (scale 6.3913435019712775e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761933689 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900518457  Egs=-0.991804480314
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635286261
Memory used: 129 MB  Time elapsed: 208 s

 ***** [NRG] Iteration 57/67 (scale 4.5193623311364666e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190988735 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900519071  Egs=-1.35909592491
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364635286421
Memory used: 130 MB  Time elapsed: 209 s

 ***** [NRG] Iteration 58/67 (scale 3.1956717509856388e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762033952 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900519388  Egs=-0.991804181477
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635286627
Memory used: 130 MB  Time elapsed: 209 s

 ***** [NRG] Iteration 59/67 (scale 2.2596811655682333e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191122829 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900519695  Egs=-1.35909623244
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.36463528669
Memory used: 131 MB  Time elapsed: 210 s

 ***** [NRG] Iteration 60/67 (scale 1.5978358754928194e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762084083 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900519854  Egs=-0.991804032047
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635286784
Memory used: 130 MB  Time elapsed: 210 s

 ***** [NRG] Iteration 61/67 (scale 1.1298405827841166e-09) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191189874 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900520007  Egs=-1.35909638622
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364635286811
Memory used: 131 MB  Time elapsed: 211 s

 ***** [NRG] Iteration 62/67 (scale 7.9891793774640969e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762109149 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900520086  Egs=-0.991803957321
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635286857
Memory used: 131 MB  Time elapsed: 211 s

 ***** [NRG] Iteration 63/67 (scale 5.6492029139205832e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191223396 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900520163  Egs=-1.35909646313
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364635286869
Memory used: 132 MB  Time elapsed: 212 s

 ***** [NRG] Iteration 64/67 (scale 3.9945896887320485e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762121681 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900520203  Egs=-0.991803919947
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635286891
Memory used: 131 MB  Time elapsed: 213 s

 ***** [NRG] Iteration 65/67 (scale 2.8246014569602916e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191240157 nrkept=1061 nrkeptmult=1061 
Total energy=-4.38900520241  Egs=-1.35909650159
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364635286897
Memory used: 132 MB  Time elapsed: 213 s

 ***** [NRG] Iteration 66/67 (scale 1.9972948443660242e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762127947 nrkept=1117 nrkeptmult=1117 
Total energy=-4.38900520261  Egs=-0.99180390125
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364635286908
Memory used: 132 MB  Time elapsed: 214 s

 ***** [NRG] Iteration 67/67 (scale 1.4123007284801458e-10) ***** 
[1] xi(66)=1 zeta(67)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191248537 nrkept=4468 nrkeptmult=4468 
Total energy=-4.3890052028  Egs=-1.35909652083
Kept: 4468 out of 4468, ratio=1
<SZd>=-0.364635286919
Memory used: 167 MB  Time elapsed: 214 s

** Iteration completed.

Spectrum: A_d-A_d-d FT
[A_d-A_d-d] pos=0.103397107979 neg=0.907834099675 sum= 1.01123120765
mu1=-0.00427819853991 mu2=0.000663186560933 mu3=-6.65271886567e-06 mu4=0.000210628538824

Spectrum: A_d-A_d-u FT
[A_d-A_d-u] pos=0.907834099675 neg=0.103397107979 sum= 1.01123120765
mu1=0.00427819853991 mu2=0.000663186560933 mu3=6.65271886568e-06 mu4=0.000210628538824


Total energy:              -4.38900520280194
TKW=0
Grand canonical Z=              5.063116181487962
Zft (last shell)=               5.063114728277493
Calculating density matrices.
Calculating full density matrices.

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QSZ
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 	MAGFIELD	 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 	 extraB=	0.0003
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d SZd
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=9
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||SZd||> (s)
Reading <||||> (z)
Number of invariant subspaces: 9
(-2 1) 1 states
0.0309062904041 
(-1 0) 2 states
0.0146225690892 0.0441882570753 
(-1 2) 2 states
0.0147186194033 0.044267066255 
(0 -1) 1 states
0.0279045357604 
(0 1) 4 states
0 0.0279919912572 0.0309062904041 0.0588982301615 
(0 3) 1 states
0.0280793952542 
(1 0) 2 states
0.0146225690892 0.0441882570753 
(1 2) 2 states
0.0147186194033 0.044267066255 
(2 1) 1 states
0.0309062904041 
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-d DMNRG
Spectrum: A_d-A_d-d CFS
Spectrum: A_d-A_d-d FDM
Spectrum: A_d-A_d-u DMNRG
Spectrum: A_d-A_d-u CFS
Spectrum: A_d-A_d-u FDM

 ***** [DMNRG] Iteration 1/67 (scale 1.2131570881878404) ***** 
[1] xi(0)=0.473132378817 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=1.89858626849 nrkept=64 nrkeptmult=64 
Kept: 64 out of 64, ratio=1
Memory used: 188 MB  Time elapsed: 241 s

 ***** [DMNRG] Iteration 2/67 (scale 0.85783160370214839) ***** 
[1] xi(1)=0.631518937152 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=3.76220868151 nrkept=256 nrkeptmult=256 
Kept: 256 out of 256, ratio=1
Memory used: 188 MB  Time elapsed: 241 s

 ***** [DMNRG] Iteration 3/67 (scale 0.6065785440939202) ***** 
[1] xi(2)=0.974788923894 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Emax=8.48719165627 nrkept=1024 nrkeptmult=1024 
Kept: 1024 out of 1024, ratio=1
Memory used: 188 MB  Time elapsed: 241 s

 ***** [DMNRG] Iteration 4/67 (scale 0.4289158018510742) ***** 
[1] xi(3)=1.2027371305 zeta(4)=0
Stats: nr=49 min=1 max=400 sum3=230553856
Safeguard: keep additional 3 states
Emax=10.316350082 nrkept=3540 nrkeptmult=3540 
Kept: 3540 out of 4096, ratio=0.864
Memory used: 210 MB  Time elapsed: 241 s

 ***** [DMNRG] Iteration 5/67 (scale 0.3032892720469601) ***** 
[1] xi(4)=1.12253968187 zeta(5)=0
Stats: nr=64 min=1 max=1038 sum3=6797276016
Safeguard: keep additional 2 states
Emax=10.3387537306 nrkept=5615 nrkeptmult=5615 
Kept: 5615 out of 14160, ratio=0.397
Memory used: 258 MB  Time elapsed: 244 s

 ***** [DMNRG] Iteration 6/67 (scale 0.2144579009255371) ***** 
[1] xi(5)=1.03243378845 zeta(6)=0
Stats: nr=77 min=2 max=1725 sum3=19012513708
Safeguard: keep additional 5 states
Emax=10.0283233223 nrkept=5542 nrkeptmult=5542 
Kept: 5542 out of 22460, ratio=0.247
Memory used: 262 MB  Time elapsed: 252 s

 ***** [DMNRG] Iteration 7/67 (scale 0.15164463602348005) ***** 
[1] xi(6)=1.00829687067 zeta(7)=0
Stats: nr=80 min=1 max=1563 sum3=23906874634
Safeguard: keep additional 11 states
Emax=10.011517717 nrkept=5136 nrkeptmult=5136 
Kept: 5136 out of 22168, ratio=0.232
Memory used: 259 MB  Time elapsed: 258 s

 ***** [DMNRG] Iteration 8/67 (scale 0.10722895046276855) ***** 
[1] xi(7)=1.00300214986 zeta(8)=0
Stats: nr=77 min=1 max=1622 sum3=18823528974
Safeguard: keep additional 8 states
Emax=10.0368187077 nrkept=5273 nrkeptmult=5273 
Kept: 5273 out of 20544, ratio=0.257
Memory used: 258 MB  Time elapsed: 264 s

 ***** [DMNRG] Iteration 9/67 (scale 0.075822318011740025) ***** 
[1] xi(8)=1.00131950393 zeta(9)=0
Stats: nr=80 min=1 max=1481 sum3=20427425036
Safeguard: keep additional 1 states
Emax=10.0014819469 nrkept=4170 nrkeptmult=4170 
Kept: 4170 out of 21092, ratio=0.198
Memory used: 238 MB  Time elapsed: 269 s

 ***** [DMNRG] Iteration 10/67 (scale 0.053614475231384275) ***** 
[1] xi(9)=1.00062206904 zeta(10)=0
Stats: nr=77 min=1 max=1338 sum3=10427379846
Emax=10.018577483 nrkept=4641 nrkeptmult=4641 
Kept: 4641 out of 16680, ratio=0.278
Memory used: 248 MB  Time elapsed: 273 s

 ***** [DMNRG] Iteration 11/67 (scale 0.037911159005870013) ***** 
[1] xi(10)=1.00030223414 zeta(11)=0
Stats: nr=76 min=1 max=1324 sum3=14394708342
Safeguard: keep additional 2 states
Emax=10.0000431522 nrkept=4023 nrkeptmult=4023 
Kept: 4023 out of 18564, ratio=0.217
Memory used: 235 MB  Time elapsed: 277 s

 ***** [DMNRG] Iteration 12/67 (scale 0.026807237615692137) ***** 
[1] xi(11)=1.00014898132 zeta(12)=0
Stats: nr=77 min=1 max=1288 sum3=9285140442
Safeguard: keep additional 4 states
Emax=10.0013605611 nrkept=4453 nrkeptmult=4453 
Kept: 4453 out of 16092, ratio=0.277
Memory used: 242 MB  Time elapsed: 280 s

 ***** [DMNRG] Iteration 13/67 (scale 0.018955579502935006) ***** 
[1] xi(12)=1.00007396421 zeta(13)=0
Stats: nr=76 min=1 max=1270 sum3=12681406924
Emax=10.0017093305 nrkept=3821 nrkeptmult=3821 
Kept: 3821 out of 17812, ratio=0.215
Memory used: 229 MB  Time elapsed: 284 s

 ***** [DMNRG] Iteration 14/67 (scale 0.013403618807846069) ***** 
[1] xi(13)=1.0000368514 zeta(14)=0
Stats: nr=77 min=1 max=1231 sum3=8065485290
Safeguard: keep additional 3 states
Emax=10.000017474 nrkept=3998 nrkeptmult=3998 
Kept: 3998 out of 15284, ratio=0.262
Memory used: 233 MB  Time elapsed: 287 s

 ***** [DMNRG] Iteration 15/67 (scale 0.0094777897514675032) ***** 
[1] xi(14)=1.00001839314 zeta(15)=0
Stats: nr=72 min=1 max=1155 sum3=9440440724
Emax=10.0200058669 nrkept=3781 nrkeptmult=3781 
Kept: 3781 out of 15992, ratio=0.236
Memory used: 234 MB  Time elapsed: 290 s

 ***** [DMNRG] Iteration 16/67 (scale 0.0067018094039230343) ***** 
[1] xi(15)=1.00000918844 zeta(16)=0
Stats: nr=77 min=1 max=1216 sum3=7792061218
Safeguard: keep additional 1 states
Emax=10.0010794918 nrkept=3636 nrkeptmult=3636 
Kept: 3636 out of 15124, ratio=0.24
Memory used: 229 MB  Time elapsed: 293 s

 ***** [DMNRG] Iteration 17/67 (scale 0.0047388948757337516) ***** 
[1] xi(16)=1.00000459219 zeta(17)=0
Stats: nr=64 min=1 max=1067 sum3=7356201924
Safeguard: keep additional 2 states
Emax=10.0028773308 nrkept=3404 nrkeptmult=3404 
Kept: 3404 out of 14544, ratio=0.234
Memory used: 227 MB  Time elapsed: 295 s

 ***** [DMNRG] Iteration 18/67 (scale 0.0033509047019615172) ***** 
[1] xi(17)=1.00000229559 zeta(18)=0
Stats: nr=73 min=1 max=1108 sum3=5829697196
Safeguard: keep additional 1 states
Emax=10.0040982903 nrkept=3245 nrkeptmult=3245 
Kept: 3245 out of 13616, ratio=0.238
Memory used: 227 MB  Time elapsed: 298 s

 ***** [DMNRG] Iteration 19/67 (scale 0.0023694474378668758) ***** 
[1] xi(18)=1.00000114767 zeta(19)=0
Stats: nr=62 min=1 max=961 sum3=5346491798
Safeguard: keep additional 3 states
Emax=10.005321246 nrkept=2931 nrkeptmult=2931 
Kept: 2931 out of 12980, ratio=0.226
Memory used: 222 MB  Time elapsed: 300 s

 ***** [DMNRG] Iteration 20/67 (scale 0.0016754523509807586) ***** 
[1] xi(19)=1.0000005738 zeta(20)=0
Stats: nr=73 min=1 max=967 sum3=3822280566
Safeguard: keep additional 1 states
Emax=10.0012636476 nrkept=2775 nrkeptmult=2775 
Kept: 2775 out of 11724, ratio=0.237
Memory used: 218 MB  Time elapsed: 301 s

 ***** [DMNRG] Iteration 21/67 (scale 0.0011847237189334379) ***** 
[1] xi(20)=1.00000028689 zeta(21)=0
Stats: nr=62 min=1 max=836 sum3=3396923856
Safeguard: keep additional 1 states
Emax=10.002424883 nrkept=2421 nrkeptmult=2421 
Kept: 2421 out of 11100, ratio=0.218
Memory used: 213 MB  Time elapsed: 303 s

 ***** [DMNRG] Iteration 22/67 (scale 0.00083772617549037929) ***** 
[1] xi(21)=1.00000014344 zeta(22)=0
Stats: nr=73 min=1 max=807 sum3=2211355914
Safeguard: keep additional 1 states
Emax=10.0001934725 nrkept=2254 nrkeptmult=2254 
Kept: 2254 out of 9684, ratio=0.233
Memory used: 213 MB  Time elapsed: 304 s

 ***** [DMNRG] Iteration 23/67 (scale 0.00059236185946671895) ***** 
[1] xi(22)=1.00000007172 zeta(23)=0
Stats: nr=62 min=1 max=695 sum3=1844095000
Safeguard: keep additional 3 states
Emax=10.0056295614 nrkept=2102 nrkeptmult=2102 
Kept: 2102 out of 9016, ratio=0.233
Memory used: 211 MB  Time elapsed: 305 s

 ***** [DMNRG] Iteration 24/67 (scale 0.00041886308774518965) ***** 
[1] xi(23)=1.00000003586 zeta(24)=0
Stats: nr=71 min=1 max=695 sum3=1431129056
Safeguard: keep additional 1 states
Emax=10.0056758785 nrkept=2035 nrkeptmult=2035 
Kept: 2035 out of 8408, ratio=0.242
Memory used: 211 MB  Time elapsed: 306 s

 ***** [DMNRG] Iteration 25/67 (scale 0.00029618092973335947) ***** 
[1] xi(24)=1.00000001793 zeta(25)=0
Stats: nr=62 min=1 max=630 sum3=1364926552
Safeguard: keep additional 1 states
Emax=10.0036573789 nrkept=1893 nrkeptmult=1893 
Kept: 1893 out of 8140, ratio=0.233
Memory used: 209 MB  Time elapsed: 307 s

 ***** [DMNRG] Iteration 26/67 (scale 0.00020943154387259482) ***** 
[1] xi(25)=1.00000000897 zeta(26)=0
Stats: nr=69 min=1 max=641 sum3=1086480156
Safeguard: keep additional 2 states
Emax=10.0133600457 nrkept=1900 nrkeptmult=1900 
Kept: 1900 out of 7572, ratio=0.251
Memory used: 212 MB  Time elapsed: 308 s

 ***** [DMNRG] Iteration 27/67 (scale 0.00014809046486667974) ***** 
[1] xi(26)=1.00000000448 zeta(27)=0
Stats: nr=61 min=1 max=599 sum3=1148042122
Emax=10.0006479317 nrkept=1660 nrkeptmult=1660 
Kept: 1660 out of 7600, ratio=0.218
Memory used: 210 MB  Time elapsed: 308 s

 ***** [DMNRG] Iteration 28/67 (scale 0.00010471577193629741) ***** 
[1] xi(27)=1.00000000224 zeta(28)=0
Stats: nr=67 min=1 max=561 sum3=742975444
Safeguard: keep additional 3 states
Emax=10.0014214886 nrkept=1629 nrkeptmult=1629 
Kept: 1629 out of 6640, ratio=0.245
Memory used: 211 MB  Time elapsed: 309 s

 ***** [DMNRG] Iteration 29/67 (scale 7.4045232433339868e-05) ***** 
[1] xi(28)=1.00000000112 zeta(29)=0
Stats: nr=61 min=1 max=533 sum3=734075190
Emax=10.01068027 nrkept=1421 nrkeptmult=1421 
Kept: 1421 out of 6516, ratio=0.218
Memory used: 211 MB  Time elapsed: 309 s

 ***** [DMNRG] Iteration 30/67 (scale 5.2357885968148706e-05) ***** 
[1] xi(29)=1.00000000056 zeta(30)=0
Stats: nr=59 min=1 max=487 sum3=494943710
Emax=10.0158508548 nrkept=1349 nrkeptmult=1349 
Kept: 1349 out of 5684, ratio=0.237
Memory used: 214 MB  Time elapsed: 310 s

 ***** [DMNRG] Iteration 31/67 (scale 3.7022616216669934e-05) ***** 
[1] xi(30)=1.00000000028 zeta(31)=0
Stats: nr=57 min=1 max=474 sum3=452028422
Safeguard: keep additional 1 states
Emax=10.0077809736 nrkept=1190 nrkeptmult=1190 
Kept: 1190 out of 5396, ratio=0.221
Memory used: 216 MB  Time elapsed: 310 s

 ***** [DMNRG] Iteration 32/67 (scale 2.6178942984074353e-05) ***** 
[1] xi(31)=1.00000000014 zeta(32)=0
Stats: nr=58 min=1 max=407 sum3=310550402
Safeguard: keep additional 1 states
Emax=10.0061179257 nrkept=1180 nrkeptmult=1180 
Kept: 1180 out of 4760, ratio=0.248
Memory used: 216 MB  Time elapsed: 310 s

 ***** [DMNRG] Iteration 33/67 (scale 1.8511308108334967e-05) ***** 
[1] xi(32)=1.00000000007 zeta(33)=0
Stats: nr=57 min=1 max=432 sum3=318384184
Emax=10.0022900413 nrkept=1094 nrkeptmult=1094 
Kept: 1094 out of 4720, ratio=0.232
Memory used: 218 MB  Time elapsed: 311 s

 ***** [DMNRG] Iteration 34/67 (scale 1.3089471492037176e-05) ***** 
[1] xi(33)=1.00000000004 zeta(34)=0
Stats: nr=58 min=1 max=373 sum3=241218758
Safeguard: keep additional 2 states
Emax=10.0033430926 nrkept=1125 nrkeptmult=1125 
Kept: 1125 out of 4376, ratio=0.257
Memory used: 217 MB  Time elapsed: 311 s

 ***** [DMNRG] Iteration 35/67 (scale 9.2556540541674836e-06) ***** 
[1] xi(34)=1.00000000002 zeta(35)=0
Stats: nr=57 min=1 max=415 sum3=278950578
Safeguard: keep additional 1 states
Emax=10.000505498 nrkept=1075 nrkeptmult=1075 
Kept: 1075 out of 4500, ratio=0.239
Memory used: 216 MB  Time elapsed: 311 s

 ***** [DMNRG] Iteration 36/67 (scale 6.5447357460185882e-06) ***** 
[1] xi(35)=1.00000000001 zeta(36)=0
Stats: nr=58 min=1 max=371 sum3=234091648
Safeguard: keep additional 1 states
Emax=10.0395198524 nrkept=1116 nrkeptmult=1116 
Kept: 1116 out of 4300, ratio=0.26
Memory used: 216 MB  Time elapsed: 311 s

 ***** [DMNRG] Iteration 37/67 (scale 4.6278270270837418e-06) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=57 min=1 max=412 sum3=271853526
Safeguard: keep additional 3 states
Emax=10.000534209 nrkept=1073 nrkeptmult=1073 
Kept: 1073 out of 4464, ratio=0.24
Memory used: 216 MB  Time elapsed: 312 s

 ***** [DMNRG] Iteration 38/67 (scale 3.2723678730092941e-06) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=58 min=1 max=371 sum3=232464242
Safeguard: keep additional 2 states
Emax=10.0626513977 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4292, ratio=0.26
Memory used: 217 MB  Time elapsed: 312 s

 ***** [DMNRG] Iteration 39/67 (scale 2.3139135135418709e-06) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=57 min=1 max=413 sum3=273706324
Safeguard: keep additional 1 states
Emax=10.0034837992 nrkept=1065 nrkeptmult=1065 
Kept: 1065 out of 4468, ratio=0.238
Memory used: 222 MB  Time elapsed: 312 s

 ***** [DMNRG] Iteration 40/67 (scale 1.6361839365046471e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=58 min=1 max=371 sum3=231782964
Emax=10.070661664 nrkept=1115 nrkeptmult=1115 
Kept: 1115 out of 4260, ratio=0.262
Memory used: 223 MB  Time elapsed: 312 s

 ***** [DMNRG] Iteration 41/67 (scale 1.1569567567709354e-06) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=57 min=1 max=413 sum3=272139644
Emax=10.0085129322 nrkept=1060 nrkeptmult=1060 
Kept: 1060 out of 4460, ratio=0.238
Memory used: 222 MB  Time elapsed: 313 s

 ***** [DMNRG] Iteration 42/67 (scale 8.1809196825232353e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=58 min=1 max=371 sum3=231642436
Safeguard: keep additional 2 states
Emax=10.0736476623 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4240, ratio=0.263
Memory used: 223 MB  Time elapsed: 313 s

 ***** [DMNRG] Iteration 43/67 (scale 5.7847837838546772e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0155357331 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 223 MB  Time elapsed: 313 s

 ***** [DMNRG] Iteration 44/67 (scale 4.0904598412616176e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0749303156 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 223 MB  Time elapsed: 313 s

 ***** [DMNRG] Iteration 45/67 (scale 2.8923918919273386e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0173921736 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 223 MB  Time elapsed: 314 s

 ***** [DMNRG] Iteration 46/67 (scale 2.0452299206308088e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0755718035 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 223 MB  Time elapsed: 314 s

 ***** [DMNRG] Iteration 47/67 (scale 1.4461959459636693e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0182657594 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 223 MB  Time elapsed: 314 s

 ***** [DMNRG] Iteration 48/67 (scale 1.0226149603154044e-07) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0758925952 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 224 MB  Time elapsed: 314 s

 ***** [DMNRG] Iteration 49/67 (scale 7.2309797298183465e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0186963227 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 225 MB  Time elapsed: 315 s

 ***** [DMNRG] Iteration 50/67 (scale 5.113074801577022e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.076053004 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 225 MB  Time elapsed: 315 s

 ***** [DMNRG] Iteration 51/67 (scale 3.6154898649091733e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.018911094 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 226 MB  Time elapsed: 315 s

 ***** [DMNRG] Iteration 52/67 (scale 2.556537400788511e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761332117 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 229 MB  Time elapsed: 315 s

 ***** [DMNRG] Iteration 53/67 (scale 1.8077449324545866e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190184086 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 229 MB  Time elapsed: 316 s

 ***** [DMNRG] Iteration 54/67 (scale 1.2782687003942555e-08) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761733163 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 229 MB  Time elapsed: 316 s

 ***** [DMNRG] Iteration 55/67 (scale 9.0387246622729332e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190720535 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 229 MB  Time elapsed: 316 s

 ***** [DMNRG] Iteration 56/67 (scale 6.3913435019712775e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761933689 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 231 MB  Time elapsed: 316 s

 ***** [DMNRG] Iteration 57/67 (scale 4.5193623311364666e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190988735 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 233 MB  Time elapsed: 317 s

 ***** [DMNRG] Iteration 58/67 (scale 3.1956717509856388e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762033952 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 234 MB  Time elapsed: 317 s

 ***** [DMNRG] Iteration 59/67 (scale 2.2596811655682333e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191122829 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 234 MB  Time elapsed: 317 s

 ***** [DMNRG] Iteration 60/67 (scale 1.5978358754928194e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762084083 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 234 MB  Time elapsed: 317 s

 ***** [DMNRG] Iteration 61/67 (scale 1.1298405827841166e-09) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191189874 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 234 MB  Time elapsed: 318 s

 ***** [DMNRG] Iteration 62/67 (scale 7.9891793774640969e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762109149 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 234 MB  Time elapsed: 318 s

 ***** [DMNRG] Iteration 63/67 (scale 5.6492029139205832e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191223396 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 234 MB  Time elapsed: 318 s

 ***** [DMNRG] Iteration 64/67 (scale 3.9945896887320485e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762121681 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 234 MB  Time elapsed: 319 s

 ***** [DMNRG] Iteration 65/67 (scale 2.8246014569602916e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191240157 nrkept=1061 nrkeptmult=1061 
Kept: 1061 out of 4468, ratio=0.237
Memory used: 234 MB  Time elapsed: 319 s

 ***** [DMNRG] Iteration 66/67 (scale 1.9972948443660242e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762127947 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4244, ratio=0.263
Memory used: 234 MB  Time elapsed: 319 s

 ***** [DMNRG] Iteration 67/67 (scale 1.4123007284801458e-10) ***** 
[1] xi(66)=1 zeta(67)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191248537 nrkept=4468 nrkeptmult=4468 
Kept: 4468 out of 4468, ratio=1
Memory used: 257 MB  Time elapsed: 320 s

** Iteration completed.

Spectrum: A_d-A_d-d DMNRG
[A_d-A_d-d] pos=0.135070888675 neg=0.864582736942 sum= 0.999653625617
mu1=-0.00379354085288 mu2=0.000662888625065 mu3=-5.05951939501e-06 mu4=0.000210577392707

Spectrum: A_d-A_d-d CFS
[A_d-A_d-d] pos=0.13536471167 neg=0.864635288328 sum= 0.999999999998
mu1=-0.0038204175935 mu2=0.000668203765323 mu3=-5.12413325071e-06 mu4=0.000216423444471

Spectrum: A_d-A_d-d FDM
[A_d-A_d-d] pos=0.13536471167 neg=0.864635288328 sum= 0.999999999998
mu1=-0.0038204175935 mu2=0.000668203765323 mu3=-5.12413325071e-06 mu4=0.000216423444471

Spectrum: A_d-A_d-u DMNRG
[A_d-A_d-u] pos=0.864582736942 neg=0.135070888675 sum= 0.999653625617
mu1=0.00379354085288 mu2=0.000662888625065 mu3=5.05951939503e-06 mu4=0.000210577392707

Spectrum: A_d-A_d-u CFS
[A_d-A_d-u] pos=0.864635288328 neg=0.13536471167 sum= 0.999999999998
mu1=0.0038204175935 mu2=0.000668203765323 mu3=5.12413325061e-06 mu4=0.000216423444471

Spectrum: A_d-A_d-u FDM
[A_d-A_d-u] pos=0.864635288328 neg=0.13536471167 sum= 0.999999999998
mu1=0.0038204175935 mu2=0.000668203765323 mu3=5.12413325061e-06 mu4=0.000216423444471

Master done. Terminating slave processes.
Master exiting.

Memory usage report (process 0)
===================
  after diag: 470988 kB
after recalc: 597288 kB
 after split: 521240 kB
  after trim: 269072 kB
 after trunc: 471432 kB
      diag_h: 423900 kB

Peak usage: 597288 kB

Timing report (process 0)
=============
         All: 321 s
          DM: 4 s
         FDM: 20 s
        diag: 189 s
    recalc d: 27 s
      spec d: 49 s
       Other: 29 s
