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=66
TEMPERATURE T=1e-10
Nmax=66
Tmin=1e-10
Ninit=0
Lambda=2
strategy=kept
lastall=1
lastalloverride=0
z=1
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=	1.
>>  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=1
# keep=10000 keepenergy=10 Nmax=66 channels=1 dots=1 betabar=1

Number of invariant subspaces: 9
(-2 1) 1 states
0.0519779776193 
(-1 0) 2 states
0.0245921318578 0.0743154939419 
(-1 2) 2 states
0.0247536685875 0.0744480346552 
(0 -1) 1 states
0.0469296481804 
(0 1) 4 states
0 0.047076730208 0.0519779776193 0.0990546212151 
(0 3) 1 states
0.0472237256235 
(1 0) 2 states
0.0245921318578 0.0743154939419 
(1 2) 2 states
0.0247536685875 0.0744480346552 
(2 1) 1 states
0.0519779776193 
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>=-3.68459992152e-05

 ***** [NRG] Iteration 1/66 (scale 0.72134752044448169) ***** 
[1] xi(0)=0.755928946018 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=3.0340606947 nrkept=64 nrkeptmult=64 
Total energy=-1.09188092215  Egs=-1.4470918763
Kept: 64 out of 64, ratio=1
<SZd>=-5.21561119387e-05
Memory used: 11 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/66 (scale 0.51006972329839473) ***** 
[1] xi(1)=0.814613079963 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=5.49874995166 nrkept=256 nrkeptmult=256 
Total energy=-1.39987299315  Egs=-0.603823471442
Kept: 256 out of 256, ratio=1
<SZd>=-7.38523553166e-05
Memory used: 12 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/66 (scale 0.36067376022224085) ***** 
[1] xi(2)=0.892494029723 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Safeguard: keep additional 1 states
Emax=10.0070843598 nrkept=1014 nrkeptmult=1014 
Total energy=-1.9250089134  Egs=-1.45598593012
Kept: 1014 out of 1024, ratio=0.99
<SZd>=-0.000104622160554
Memory used: 16 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/66 (scale 0.25503486164919736) ***** 
[1] xi(3)=0.942104180991 zeta(4)=0
Stats: nr=49 min=1 max=394 sum3=221685426
Safeguard: keep additional 5 states
Emax=10.8970701777 nrkept=3046 nrkeptmult=3046 
Total energy=-2.12936322962  Egs=-0.801279930526
Kept: 3046 out of 4056, ratio=0.751
<SZd>=-0.000148311750674
Memory used: 46 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/66 (scale 0.18033688011112042) ***** 
[1] xi(4)=0.96993426226 zeta(5)=0
Stats: nr=64 min=1 max=892 sum3=4319913808
Safeguard: keep additional 9 states
Emax=10.0106390802 nrkept=4174 nrkeptmult=4174 
Total energy=-2.39724720701  Egs=-1.48546418915
Kept: 4174 out of 12184, ratio=0.343
<SZd>=-0.000210427087207
Memory used: 76 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 6/66 (scale 0.12751743082459868) ***** 
[1] xi(5)=0.984675515632 zeta(6)=0
Stats: nr=77 min=1 max=1324 sum3=10177111120
Safeguard: keep additional 8 states
Emax=10.0192634915 nrkept=5209 nrkeptmult=5209 
Total energy=-2.5088071214  Egs=-0.874860116526
Kept: 5209 out of 16696, ratio=0.312
<SZd>=-0.000299132521222
Memory used: 105 MB  Time elapsed: 16 s

 ***** [NRG] Iteration 7/66 (scale 0.090168440055560212) ***** 
[1] xi(6)=0.992263204714 zeta(7)=0
Stats: nr=80 min=1 max=1457 sum3=19510937620
Safeguard: keep additional 1 states
Emax=10.0005532337 nrkept=4202 nrkeptmult=4202 
Total energy=-2.64161896877  Egs=-1.47293052082
Kept: 4202 out of 20836, ratio=0.202
<SZd>=-0.000430277876909
Memory used: 93 MB  Time elapsed: 32 s

 ***** [NRG] Iteration 8/66 (scale 0.063758715412299341) ***** 
[1] xi(7)=0.996112749408 zeta(8)=0
Stats: nr=77 min=1 max=1350 sum3=10696177766
Safeguard: keep additional 3 states
Emax=10.0001515125 nrkept=4836 nrkeptmult=4836 
Total energy=-2.69987416287  Egs=-0.913682054556
Kept: 4836 out of 16808, ratio=0.288
<SZd>=-0.000607903108457
Memory used: 108 MB  Time elapsed: 43 s

 ***** [NRG] Iteration 9/66 (scale 0.045084220027780106) ***** 
[1] xi(8)=0.998051634085 zeta(9)=0
Stats: nr=80 min=1 max=1367 sum3=15978297264
Safeguard: keep additional 3 states
Emax=10.0127093817 nrkept=4096 nrkeptmult=4096 
Total energy=-2.76499638501  Egs=-1.44445710939
Kept: 4096 out of 19344, ratio=0.212
<SZd>=-0.0008707372809
Memory used: 105 MB  Time elapsed: 58 s

 ***** [NRG] Iteration 10/66 (scale 0.031879357706149671) ***** 
[1] xi(9)=0.99902462843 zeta(10)=0
Stats: nr=77 min=1 max=1312 sum3=9823752472
Safeguard: keep additional 2 states
Emax=10.0620759351 nrkept=4603 nrkeptmult=4603 
Total energy=-2.79508958235  Egs=-0.94397125636
Kept: 4603 out of 16384, ratio=0.281
<SZd>=-0.00124926183725
Memory used: 117 MB  Time elapsed: 68 s

 ***** [NRG] Iteration 11/66 (scale 0.022542110013890053) ***** 
[1] xi(10)=0.999512016628 zeta(11)=0
Stats: nr=76 min=1 max=1310 sum3=13968677740
Safeguard: keep additional 1 states
Emax=10.0042042285 nrkept=3834 nrkeptmult=3834 
Total energy=-2.82691125777  Egs=-1.41165469433
Kept: 3834 out of 18412, ratio=0.208
<SZd>=-0.00181553607099
Memory used: 107 MB  Time elapsed: 82 s

 ***** [NRG] Iteration 12/66 (scale 0.015939678853074835) ***** 
[1] xi(11)=0.999755933863 zeta(12)=0
Stats: nr=77 min=1 max=1236 sum3=8159976876
Safeguard: keep additional 3 states
Emax=10.0173842861 nrkept=4084 nrkeptmult=4084 
Total energy=-2.84237384334  Egs=-0.970068827073
Kept: 4084 out of 15336, ratio=0.266
<SZd>=-0.00262952398256
Memory used: 110 MB  Time elapsed: 91 s

 ***** [NRG] Iteration 13/66 (scale 0.011271055006945026) ***** 
[1] xi(12)=0.999877948312 zeta(13)=0
Stats: nr=74 min=1 max=1179 sum3=10077672808
Emax=10.0073376241 nrkept=3781 nrkeptmult=3781 
Total energy=-2.85793330157  Egs=-1.38047930914
Kept: 3781 out of 16336, ratio=0.231
<SZd>=-0.00384648209816
Memory used: 109 MB  Time elapsed: 102 s

 ***** [NRG] Iteration 14/66 (scale 0.0079698394265374176) ***** 
[1] xi(13)=0.9999389695 zeta(14)=0
Stats: nr=77 min=1 max=1216 sum3=7792061218
Safeguard: keep additional 1 states
Emax=10.000927538 nrkept=3728 nrkeptmult=3728 
Total energy=-2.86583060833  Egs=-0.990899105501
Kept: 3728 out of 15124, ratio=0.246
<SZd>=-0.00568986470024
Memory used: 110 MB  Time elapsed: 110 s

 ***** [NRG] Iteration 15/66 (scale 0.0056355275034725132) ***** 
[1] xi(14)=0.999969483586 zeta(15)=0
Stats: nr=64 min=1 max=1091 sum3=7900199396
Safeguard: keep additional 2 states
Emax=10.0010088696 nrkept=3578 nrkeptmult=3578 
Total energy=-2.87346366958  Egs=-1.35445373013
Kept: 3578 out of 14912, ratio=0.24
<SZd>=-0.00845559560354
Memory used: 110 MB  Time elapsed: 118 s

 ***** [NRG] Iteration 16/66 (scale 0.0039849197132687088) ***** 
[1] xi(15)=0.999984741502 zeta(16)=0
Stats: nr=73 min=1 max=1164 sum3=6769666658
Safeguard: keep additional 1 states
Emax=10.0025111732 nrkept=3380 nrkeptmult=3380 
Total energy=-2.8774708578  Egs=-1.00558819594
Kept: 3380 out of 14312, ratio=0.236
<SZd>=-0.0125786335441
Memory used: 108 MB  Time elapsed: 125 s

 ***** [NRG] Iteration 17/66 (scale 0.0028177637517362566) ***** 
[1] xi(16)=0.999992370678 zeta(17)=0
Stats: nr=64 min=1 max=996 sum3=5952281330
Emax=10.0050375005 nrkept=3053 nrkeptmult=3053 
Total energy=-2.88123377352  Egs=-1.33542626149
Kept: 3053 out of 13520, ratio=0.226
<SZd>=-0.0185858981023
Memory used: 107 MB  Time elapsed: 132 s

 ***** [NRG] Iteration 18/66 (scale 0.0019924598566343544) ***** 
[1] xi(17)=0.999996185321 zeta(18)=0
Stats: nr=73 min=1 max=1003 sum3=4285165010
Safeguard: keep additional 3 states
Emax=10.0000894037 nrkept=2920 nrkeptmult=2920 
Total energy=-2.88325422813  Egs=-1.01405034826
Kept: 2920 out of 12212, ratio=0.239
<SZd>=-0.027091671663
Memory used: 106 MB  Time elapsed: 137 s

 ***** [NRG] Iteration 19/66 (scale 0.0014088818758681283) ***** 
[1] xi(18)=0.999998092656 zeta(19)=0
Stats: nr=62 min=1 max=868 sum3=3877566166
Safeguard: keep additional 2 states
Emax=10.0034434081 nrkept=2541 nrkeptmult=2541 
Total energy=-2.88512074707  Egs=-1.32482287924
Kept: 2541 out of 11680, ratio=0.218
<SZd>=-0.0387153054535
Memory used: 101 MB  Time elapsed: 141 s

 ***** [NRG] Iteration 20/66 (scale 0.0009962299283171772) ***** 
[1] xi(19)=0.999999046327 zeta(20)=0
Stats: nr=73 min=1 max=851 sum3=2567408694
Safeguard: keep additional 2 states
Emax=10.0001586761 nrkept=2449 nrkeptmult=2449 
Total energy=-2.88613303855  Egs=-1.0161223306
Kept: 2449 out of 10164, ratio=0.241
<SZd>=-0.0542290282804
Memory used: 102 MB  Time elapsed: 145 s

 ***** [NRG] Iteration 21/66 (scale 0.00070444093793406415) ***** 
[1] xi(20)=0.999999523163 zeta(21)=0
Stats: nr=62 min=1 max=740 sum3=2311551112
Safeguard: keep additional 1 states
Emax=10.0017479225 nrkept=2143 nrkeptmult=2143 
Total energy=-2.88706554469  Egs=-1.3237534812
Kept: 2143 out of 9796, ratio=0.219
<SZd>=-0.0746863943143
Memory used: 101 MB  Time elapsed: 148 s

 ***** [NRG] Iteration 22/66 (scale 0.0004981149641585886) ***** 
[1] xi(21)=0.999999761581 zeta(22)=0
Stats: nr=71 min=1 max=713 sum3=1526010688
Safeguard: keep additional 1 states
Emax=10.0143505271 nrkept=2058 nrkeptmult=2058 
Total energy=-2.88756996979  Egs=-1.0126680284
Kept: 2058 out of 8572, ratio=0.24
<SZd>=-0.101430196599
Memory used: 101 MB  Time elapsed: 150 s

 ***** [NRG] Iteration 23/66 (scale 0.00035222046896703208) ***** 
[1] xi(22)=0.999999880791 zeta(23)=0
Stats: nr=62 min=1 max=639 sum3=1404724236
Safeguard: keep additional 1 states
Emax=10.0021578418 nrkept=1942 nrkeptmult=1942 
Total energy=-2.88803864074  Egs=-1.33061815194
Kept: 1942 out of 8232, ratio=0.236
<SZd>=-0.135460488292
Memory used: 99 MB  Time elapsed: 152 s

 ***** [NRG] Iteration 24/66 (scale 0.0002490574820792943) ***** 
[1] xi(23)=0.999999940395 zeta(24)=0
Stats: nr=69 min=1 max=656 sum3=1169945404
Safeguard: keep additional 3 states
Emax=10.0024201944 nrkept=1967 nrkeptmult=1967 
Total energy=-2.88828934434  Egs=-1.00660940139
Kept: 1967 out of 7768, ratio=0.253
<SZd>=-0.176917705471
Memory used: 103 MB  Time elapsed: 154 s

 ***** [NRG] Iteration 25/66 (scale 0.00017611023448351604) ***** 
[1] xi(24)=0.999999970198 zeta(25)=0
Stats: nr=61 min=1 max=612 sum3=1248954038
Emax=10.0035392215 nrkept=1713 nrkeptmult=1713 
Total energy=-2.88852539144  Egs=-1.34033716886
Kept: 1713 out of 7868, ratio=0.218
<SZd>=-0.224080733563
Memory used: 102 MB  Time elapsed: 156 s

 ***** [NRG] Iteration 26/66 (scale 0.00012452874103964715) ***** 
[1] xi(25)=0.999999985099 zeta(26)=0
Stats: nr=69 min=1 max=581 sum3=816220428
Safeguard: keep additional 4 states
Emax=10.0002797499 nrkept=1680 nrkeptmult=1680 
Total energy=-2.88865003838  Egs=-1.00094916459
Kept: 1680 out of 6852, ratio=0.245
<SZd>=-0.271524494967
Memory used: 102 MB  Time elapsed: 157 s

 ***** [NRG] Iteration 27/66 (scale 8.8055117241758019e-05) ***** 
[1] xi(26)=0.999999992549 zeta(27)=0
Stats: nr=61 min=1 max=548 sum3=798076332
Safeguard: keep additional 4 states
Emax=10.0018814628 nrkept=1477 nrkeptmult=1477 
Total energy=-2.88876876683  Egs=-1.34834251659
Kept: 1477 out of 6720, ratio=0.22
<SZd>=-0.3120507592
Memory used: 104 MB  Time elapsed: 158 s

 ***** [NRG] Iteration 28/66 (scale 6.2264370519823575e-05) ***** 
[1] xi(27)=0.999999996275 zeta(28)=0
Stats: nr=63 min=1 max=506 sum3=547457896
Safeguard: keep additional 3 states
Emax=10.0069144664 nrkept=1420 nrkeptmult=1420 
Total energy=-2.88883084147  Egs=-0.996952771157
Kept: 1420 out of 5908, ratio=0.24
<SZd>=-0.339652970868
Memory used: 104 MB  Time elapsed: 159 s

 ***** [NRG] Iteration 29/66 (scale 4.402755862087901e-05) ***** 
[1] xi(28)=0.999999998137 zeta(29)=0
Stats: nr=59 min=1 max=488 sum3=514851298
Safeguard: keep additional 3 states
Emax=10.0026862654 nrkept=1243 nrkeptmult=1243 
Total energy=-2.88889042758  Egs=-1.35338197266
Kept: 1243 out of 5680, ratio=0.219
<SZd>=-0.354347435795
Memory used: 105 MB  Time elapsed: 160 s

 ***** [NRG] Iteration 30/66 (scale 3.1132185259911788e-05) ***** 
[1] xi(29)=0.999999999069 zeta(30)=0
Stats: nr=58 min=1 max=428 sum3=348701728
Safeguard: keep additional 1 states
Emax=10.0005256635 nrkept=1214 nrkeptmult=1214 
Total energy=-2.88892138991  Egs=-0.994544148258
Kept: 1214 out of 4972, ratio=0.244
<SZd>=-0.36049349726
Memory used: 109 MB  Time elapsed: 161 s

 ***** [NRG] Iteration 31/66 (scale 2.2013779310439505e-05) ***** 
[1] xi(30)=0.999999999534 zeta(31)=0
Stats: nr=57 min=1 max=437 sum3=337641950
Safeguard: keep additional 1 states
Emax=10.0073109725 nrkept=1103 nrkeptmult=1103 
Total energy=-2.88895124416  Egs=-1.35616213781
Kept: 1103 out of 4856, ratio=0.227
<SZd>=-0.362844572835
Memory used: 109 MB  Time elapsed: 162 s

 ***** [NRG] Iteration 32/66 (scale 1.5566092629955894e-05) ***** 
[1] xi(31)=0.999999999767 zeta(32)=0
Stats: nr=58 min=1 max=375 sum3=247813640
Safeguard: keep additional 1 states
Emax=10.0026675186 nrkept=1132 nrkeptmult=1132 
Total energy=-2.8889667047  Egs=-0.993218642735
Kept: 1132 out of 4412, ratio=0.257
<SZd>=-0.363788600936
Memory used: 109 MB  Time elapsed: 162 s

 ***** [NRG] Iteration 33/66 (scale 1.1006889655219752e-05) ***** 
[1] xi(32)=0.999999999884 zeta(33)=0
Stats: nr=57 min=1 max=416 sum3=282675538
Safeguard: keep additional 1 states
Emax=10.0070238801 nrkept=1079 nrkeptmult=1079 
Total energy=-2.88898164777  Egs=-1.35761060639
Kept: 1079 out of 4528, ratio=0.238
<SZd>=-0.364220937071
Memory used: 110 MB  Time elapsed: 163 s

 ***** [NRG] Iteration 34/66 (scale 7.7830463149779469e-06) ***** 
[1] xi(33)=0.999999999942 zeta(34)=0
Stats: nr=58 min=1 max=371 sum3=235967204
Safeguard: keep additional 1 states
Emax=10.030116036 nrkept=1116 nrkeptmult=1116 
Total energy=-2.88898937262  Egs=-0.992523672875
Kept: 1116 out of 4316, ratio=0.259
<SZd>=-0.364427551223
Memory used: 109 MB  Time elapsed: 163 s

 ***** [NRG] Iteration 35/66 (scale 5.5034448276098762e-06) ***** 
[1] xi(34)=0.999999999971 zeta(35)=0
Stats: nr=57 min=1 max=412 sum3=271853526
Safeguard: keep additional 1 states
Emax=10.0355827293 nrkept=1075 nrkeptmult=1075 
Total energy=-2.88899684822  Egs=-1.35834787778
Kept: 1075 out of 4464, ratio=0.241
<SZd>=-0.364529076165
Memory used: 110 MB  Time elapsed: 164 s

 ***** [NRG] Iteration 36/66 (scale 3.8915231574889735e-06) ***** 
[1] xi(35)=0.999999999985 zeta(36)=0
Stats: nr=58 min=1 max=371 sum3=234091648
Safeguard: keep additional 2 states
Emax=10.0584490488 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88900070926  Egs=-0.99216799527
Kept: 1117 out of 4300, ratio=0.26
<SZd>=-0.36457941443
Memory used: 110 MB  Time elapsed: 165 s

 ***** [NRG] Iteration 37/66 (scale 2.7517224138049381e-06) ***** 
[1] xi(36)=0.999999999993 zeta(37)=0
Stats: nr=57 min=1 max=413 sum3=273706324
Safeguard: keep additional 1 states
Emax=10.0017711281 nrkept=1067 nrkeptmult=1067 
Total energy=-2.88900444808  Egs=-1.35871952762
Kept: 1067 out of 4468, ratio=0.239
<SZd>=-0.364604362926
Memory used: 114 MB  Time elapsed: 165 s

 ***** [NRG] Iteration 38/66 (scale 1.9457615787444867e-06) ***** 
[1] xi(37)=0.999999999996 zeta(38)=0
Stats: nr=58 min=1 max=371 sum3=231996350
Emax=10.0694101108 nrkept=1115 nrkeptmult=1115 
Total energy=-2.88900637825  Egs=-0.991988101685
Kept: 1115 out of 4268, ratio=0.261
<SZd>=-0.364617019654
Memory used: 115 MB  Time elapsed: 166 s

 ***** [NRG] Iteration 39/66 (scale 1.3758612069024691e-06) ***** 
[1] xi(38)=0.999999999998 zeta(39)=0
Stats: nr=57 min=1 max=413 sum3=272139644
Emax=10.0056891737 nrkept=1060 nrkeptmult=1060 
Total energy=-2.88900824792  Egs=-1.358906072
Kept: 1060 out of 4460, ratio=0.238
<SZd>=-0.364623260937
Memory used: 116 MB  Time elapsed: 166 s

 ***** [NRG] Iteration 40/66 (scale 9.7288078937224336e-07) ***** 
[1] xi(39)=0.999999999999 zeta(40)=0
Stats: nr=58 min=1 max=371 sum3=231642436
Safeguard: keep additional 2 states
Emax=10.0731672522 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88900921292  Egs=-0.99189764121
Kept: 1117 out of 4240, ratio=0.263
<SZd>=-0.364626437449
Memory used: 115 MB  Time elapsed: 167 s

 ***** [NRG] Iteration 41/66 (scale 6.8793060345123453e-07) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Emax=10.014400013 nrkept=1060 nrkeptmult=1060 
Total energy=-2.88901014781  Egs=-1.35899951965
Kept: 1060 out of 4468, ratio=0.237
<SZd>=-0.3646279964
Memory used: 116 MB  Time elapsed: 167 s

 ***** [NRG] Iteration 42/66 (scale 4.8644039468612168e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=58 min=1 max=371 sum3=231642436
Safeguard: keep additional 2 states
Emax=10.0746925257 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88901063029  Egs=-0.991852282612
Kept: 1117 out of 4240, ratio=0.263
<SZd>=-0.364628809384
Memory used: 116 MB  Time elapsed: 168 s

 ***** [NRG] Iteration 43/66 (scale 3.4396530172561726e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0170631144 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901109776  Egs=-1.35904628678
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629205737
Memory used: 117 MB  Time elapsed: 168 s

 ***** [NRG] Iteration 44/66 (scale 2.4322019734306084e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0754553744 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88901133899  Egs=-0.991829571227
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.36462941262
Memory used: 116 MB  Time elapsed: 169 s

 ***** [NRG] Iteration 45/66 (scale 1.7198265086280863e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0181138061 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901157272  Egs=-1.35906968111
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629513434
Memory used: 118 MB  Time elapsed: 170 s

 ***** [NRG] Iteration 46/66 (scale 1.2161009867153042e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0758368642 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88901169334  Egs=-0.991818207506
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629568436
Memory used: 119 MB  Time elapsed: 170 s

 ***** [NRG] Iteration 47/66 (scale 8.5991325431404316e-08) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0186265535 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901181021  Egs=-1.35908138097
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629594522
Memory used: 120 MB  Time elapsed: 171 s

 ***** [NRG] Iteration 48/66 (scale 6.080504933576521e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0760276271 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88901187052  Egs=-0.991812523632
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629609919
Memory used: 120 MB  Time elapsed: 171 s

 ***** [NRG] Iteration 49/66 (scale 4.2995662715702158e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0188820554 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901192895  Egs=-1.35908723157
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629616879
Memory used: 121 MB  Time elapsed: 172 s

 ***** [NRG] Iteration 50/66 (scale 3.0402524667882605e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761230132 nrkept=1117 nrkeptmult=1117 
Total energy=-2.8890119591  Egs=-0.991809681184
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629621551
Memory used: 120 MB  Time elapsed: 172 s

 ***** [NRG] Iteration 51/66 (scale 2.1497831357851079e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190096923 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901198832  Egs=-1.35909015706
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.36462962351
Memory used: 121 MB  Time elapsed: 173 s

 ***** [NRG] Iteration 52/66 (scale 1.5201262333941303e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761707075 nrkept=1117 nrkeptmult=1117 
Total energy=-2.8890120034  Egs=-0.991808259824
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629625089
Memory used: 121 MB  Time elapsed: 173 s

 ***** [NRG] Iteration 53/66 (scale 1.0748915678925539e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.019073492 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901201801  Egs=-1.35909161985
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629625688
Memory used: 122 MB  Time elapsed: 174 s

 ***** [NRG] Iteration 54/66 (scale 7.6006311669706513e-09) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0761945549 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88901202555  Egs=-0.991807549102
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629626288
Memory used: 121 MB  Time elapsed: 175 s

 ***** [NRG] Iteration 55/66 (scale 5.3744578394627697e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191053882 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901203285  Egs=-1.35909235128
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629626492
Memory used: 122 MB  Time elapsed: 175 s

 ***** [NRG] Iteration 56/66 (scale 3.8003155834853256e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762064786 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88901203662  Egs=-0.991807193723
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629626745
Memory used: 122 MB  Time elapsed: 176 s

 ***** [NRG] Iteration 57/66 (scale 2.6872289197313849e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191213354 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901204027  Egs=-1.359092717
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629626823
Memory used: 123 MB  Time elapsed: 176 s

 ***** [NRG] Iteration 58/66 (scale 1.9001577917426628e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762124405 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88901204216  Egs=-0.991807016021
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629626938
Memory used: 122 MB  Time elapsed: 177 s

 ***** [NRG] Iteration 59/66 (scale 1.3436144598656924e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191293088 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901204398  Egs=-1.35909289988
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629626971
Memory used: 123 MB  Time elapsed: 177 s

 ***** [NRG] Iteration 60/66 (scale 9.5007889587133141e-10) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762154214 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88901204492  Egs=-0.991806927159
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629627025
Memory used: 122 MB  Time elapsed: 178 s

 ***** [NRG] Iteration 61/66 (scale 6.7180722993284622e-10) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191332955 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901204584  Egs=-1.35909299133
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629627041
Memory used: 123 MB  Time elapsed: 179 s

 ***** [NRG] Iteration 62/66 (scale 4.750394479356657e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762169119 nrkept=1117 nrkeptmult=1117 
Total energy=-2.88901204631  Egs=-0.991806882717
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629627067
Memory used: 123 MB  Time elapsed: 179 s

 ***** [NRG] Iteration 63/66 (scale 3.3590361496642311e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191352887 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901204677  Egs=-1.35909303706
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629627074
Memory used: 124 MB  Time elapsed: 180 s

 ***** [NRG] Iteration 64/66 (scale 2.3751972396783285e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762176571 nrkept=1117 nrkeptmult=1117 
Total energy=-2.889012047  Egs=-0.991806860486
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364629627087
Memory used: 123 MB  Time elapsed: 180 s

 ***** [NRG] Iteration 65/66 (scale 1.6795180748321155e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0191362853 nrkept=1061 nrkeptmult=1061 
Total energy=-2.88901204723  Egs=-1.35909305994
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364629627091
Memory used: 124 MB  Time elapsed: 181 s

 ***** [NRG] Iteration 66/66 (scale 1.1875986198391643e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0762180296 nrkept=4244 nrkeptmult=4244 
Total energy=-2.88901204735  Egs=-0.99180684936
Kept: 4244 out of 4244, ratio=1
<SZd>=-0.364629627099
Memory used: 155 MB  Time elapsed: 182 s

** Iteration completed.

Spectrum: A_d-A_d-d FT
[A_d-A_d-d] pos=0.103420951159 neg=0.907954492328 sum= 1.01137544349
mu1=-0.00427920873338 mu2=0.000662852847227 mu3=-6.64994679309e-06 mu4=0.000189918362647

Spectrum: A_d-A_d-u FT
[A_d-A_d-u] pos=0.907954492328 neg=0.103420951159 sum= 1.01137544349
mu1=0.00427920873338 mu2=0.000662852847227 mu3=6.64994679309e-06 mu4=0.000189918362647


Total energy:             -2.889012047346915
TKW=0
Grand canonical Z=              8.341331365068363
Zft (last shell)=               8.341305715665925
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=	1.
>>  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.0519779776193 
(-1 0) 2 states
0.0245921318578 0.0743154939419 
(-1 2) 2 states
0.0247536685875 0.0744480346552 
(0 -1) 1 states
0.0469296481804 
(0 1) 4 states
0 0.047076730208 0.0519779776193 0.0990546212151 
(0 3) 1 states
0.0472237256235 
(1 0) 2 states
0.0245921318578 0.0743154939419 
(1 2) 2 states
0.0247536685875 0.0744480346552 
(2 1) 1 states
0.0519779776193 
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/66 (scale 0.72134752044448169) ***** 
[1] xi(0)=0.755928946018 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=3.0340606947 nrkept=64 nrkeptmult=64 
Kept: 64 out of 64, ratio=1
Memory used: 191 MB  Time elapsed: 204 s

 ***** [DMNRG] Iteration 2/66 (scale 0.51006972329839473) ***** 
[1] xi(1)=0.814613079963 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=5.49874995166 nrkept=256 nrkeptmult=256 
Kept: 256 out of 256, ratio=1
Memory used: 191 MB  Time elapsed: 204 s

 ***** [DMNRG] Iteration 3/66 (scale 0.36067376022224085) ***** 
[1] xi(2)=0.892494029723 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Safeguard: keep additional 1 states
Emax=10.0070843598 nrkept=1014 nrkeptmult=1014 
Kept: 1014 out of 1024, ratio=0.99
Memory used: 192 MB  Time elapsed: 204 s

 ***** [DMNRG] Iteration 4/66 (scale 0.25503486164919736) ***** 
[1] xi(3)=0.942104180991 zeta(4)=0
Stats: nr=49 min=1 max=394 sum3=221685426
Safeguard: keep additional 5 states
Emax=10.8970701777 nrkept=3046 nrkeptmult=3046 
Kept: 3046 out of 4056, ratio=0.751
Memory used: 206 MB  Time elapsed: 204 s

 ***** [DMNRG] Iteration 5/66 (scale 0.18033688011112042) ***** 
[1] xi(4)=0.96993426226 zeta(5)=0
Stats: nr=64 min=1 max=892 sum3=4319913808
Safeguard: keep additional 9 states
Emax=10.0106390802 nrkept=4174 nrkeptmult=4174 
Kept: 4174 out of 12184, ratio=0.343
Memory used: 230 MB  Time elapsed: 206 s

 ***** [DMNRG] Iteration 6/66 (scale 0.12751743082459868) ***** 
[1] xi(5)=0.984675515632 zeta(6)=0
Stats: nr=77 min=1 max=1324 sum3=10177111120
Safeguard: keep additional 8 states
Emax=10.0192634915 nrkept=5209 nrkeptmult=5209 
Kept: 5209 out of 16696, ratio=0.312
Memory used: 264 MB  Time elapsed: 210 s

 ***** [DMNRG] Iteration 7/66 (scale 0.090168440055560212) ***** 
[1] xi(6)=0.992263204714 zeta(7)=0
Stats: nr=80 min=1 max=1457 sum3=19510937620
Safeguard: keep additional 1 states
Emax=10.0005532337 nrkept=4202 nrkeptmult=4202 
Kept: 4202 out of 20836, ratio=0.202
Memory used: 241 MB  Time elapsed: 216 s

 ***** [DMNRG] Iteration 8/66 (scale 0.063758715412299341) ***** 
[1] xi(7)=0.996112749408 zeta(8)=0
Stats: nr=77 min=1 max=1350 sum3=10696177766
Safeguard: keep additional 3 states
Emax=10.0001515125 nrkept=4836 nrkeptmult=4836 
Kept: 4836 out of 16808, ratio=0.288
Memory used: 260 MB  Time elapsed: 220 s

 ***** [DMNRG] Iteration 9/66 (scale 0.045084220027780106) ***** 
[1] xi(8)=0.998051634085 zeta(9)=0
Stats: nr=80 min=1 max=1367 sum3=15978297264
Safeguard: keep additional 3 states
Emax=10.0127093817 nrkept=4096 nrkeptmult=4096 
Kept: 4096 out of 19344, ratio=0.212
Memory used: 244 MB  Time elapsed: 225 s

 ***** [DMNRG] Iteration 10/66 (scale 0.031879357706149671) ***** 
[1] xi(9)=0.99902462843 zeta(10)=0
Stats: nr=77 min=1 max=1312 sum3=9823752472
Safeguard: keep additional 2 states
Emax=10.0620759351 nrkept=4603 nrkeptmult=4603 
Kept: 4603 out of 16384, ratio=0.281
Memory used: 256 MB  Time elapsed: 229 s

 ***** [DMNRG] Iteration 11/66 (scale 0.022542110013890053) ***** 
[1] xi(10)=0.999512016628 zeta(11)=0
Stats: nr=76 min=1 max=1310 sum3=13968677740
Safeguard: keep additional 1 states
Emax=10.0042042285 nrkept=3834 nrkeptmult=3834 
Kept: 3834 out of 18412, ratio=0.208
Memory used: 239 MB  Time elapsed: 233 s

 ***** [DMNRG] Iteration 12/66 (scale 0.015939678853074835) ***** 
[1] xi(11)=0.999755933863 zeta(12)=0
Stats: nr=77 min=1 max=1236 sum3=8159976876
Safeguard: keep additional 3 states
Emax=10.0173842861 nrkept=4084 nrkeptmult=4084 
Kept: 4084 out of 15336, ratio=0.266
Memory used: 245 MB  Time elapsed: 236 s

 ***** [DMNRG] Iteration 13/66 (scale 0.011271055006945026) ***** 
[1] xi(12)=0.999877948312 zeta(13)=0
Stats: nr=74 min=1 max=1179 sum3=10077672808
Emax=10.0073376241 nrkept=3781 nrkeptmult=3781 
Kept: 3781 out of 16336, ratio=0.231
Memory used: 239 MB  Time elapsed: 239 s

 ***** [DMNRG] Iteration 14/66 (scale 0.0079698394265374176) ***** 
[1] xi(13)=0.9999389695 zeta(14)=0
Stats: nr=77 min=1 max=1216 sum3=7792061218
Safeguard: keep additional 1 states
Emax=10.000927538 nrkept=3728 nrkeptmult=3728 
Kept: 3728 out of 15124, ratio=0.246
Memory used: 236 MB  Time elapsed: 242 s

 ***** [DMNRG] Iteration 15/66 (scale 0.0056355275034725132) ***** 
[1] xi(14)=0.999969483586 zeta(15)=0
Stats: nr=64 min=1 max=1091 sum3=7900199396
Safeguard: keep additional 2 states
Emax=10.0010088696 nrkept=3578 nrkeptmult=3578 
Kept: 3578 out of 14912, ratio=0.24
Memory used: 235 MB  Time elapsed: 245 s

 ***** [DMNRG] Iteration 16/66 (scale 0.0039849197132687088) ***** 
[1] xi(15)=0.999984741502 zeta(16)=0
Stats: nr=73 min=1 max=1164 sum3=6769666658
Safeguard: keep additional 1 states
Emax=10.0025111732 nrkept=3380 nrkeptmult=3380 
Kept: 3380 out of 14312, ratio=0.236
Memory used: 232 MB  Time elapsed: 248 s

 ***** [DMNRG] Iteration 17/66 (scale 0.0028177637517362566) ***** 
[1] xi(16)=0.999992370678 zeta(17)=0
Stats: nr=64 min=1 max=996 sum3=5952281330
Emax=10.0050375005 nrkept=3053 nrkeptmult=3053 
Kept: 3053 out of 13520, ratio=0.226
Memory used: 230 MB  Time elapsed: 250 s

 ***** [DMNRG] Iteration 18/66 (scale 0.0019924598566343544) ***** 
[1] xi(17)=0.999996185321 zeta(18)=0
Stats: nr=73 min=1 max=1003 sum3=4285165010
Safeguard: keep additional 3 states
Emax=10.0000894037 nrkept=2920 nrkeptmult=2920 
Kept: 2920 out of 12212, ratio=0.239
Memory used: 226 MB  Time elapsed: 252 s

 ***** [DMNRG] Iteration 19/66 (scale 0.0014088818758681283) ***** 
[1] xi(18)=0.999998092656 zeta(19)=0
Stats: nr=62 min=1 max=868 sum3=3877566166
Safeguard: keep additional 2 states
Emax=10.0034434081 nrkept=2541 nrkeptmult=2541 
Kept: 2541 out of 11680, ratio=0.218
Memory used: 220 MB  Time elapsed: 254 s

 ***** [DMNRG] Iteration 20/66 (scale 0.0009962299283171772) ***** 
[1] xi(19)=0.999999046327 zeta(20)=0
Stats: nr=73 min=1 max=851 sum3=2567408694
Safeguard: keep additional 2 states
Emax=10.0001586761 nrkept=2449 nrkeptmult=2449 
Kept: 2449 out of 10164, ratio=0.241
Memory used: 222 MB  Time elapsed: 255 s

 ***** [DMNRG] Iteration 21/66 (scale 0.00070444093793406415) ***** 
[1] xi(20)=0.999999523163 zeta(21)=0
Stats: nr=62 min=1 max=740 sum3=2311551112
Safeguard: keep additional 1 states
Emax=10.0017479225 nrkept=2143 nrkeptmult=2143 
Kept: 2143 out of 9796, ratio=0.219
Memory used: 219 MB  Time elapsed: 256 s

 ***** [DMNRG] Iteration 22/66 (scale 0.0004981149641585886) ***** 
[1] xi(21)=0.999999761581 zeta(22)=0
Stats: nr=71 min=1 max=713 sum3=1526010688
Safeguard: keep additional 1 states
Emax=10.0143505271 nrkept=2058 nrkeptmult=2058 
Kept: 2058 out of 8572, ratio=0.24
Memory used: 219 MB  Time elapsed: 257 s

 ***** [DMNRG] Iteration 23/66 (scale 0.00035222046896703208) ***** 
[1] xi(22)=0.999999880791 zeta(23)=0
Stats: nr=62 min=1 max=639 sum3=1404724236
Safeguard: keep additional 1 states
Emax=10.0021578418 nrkept=1942 nrkeptmult=1942 
Kept: 1942 out of 8232, ratio=0.236
Memory used: 217 MB  Time elapsed: 258 s

 ***** [DMNRG] Iteration 24/66 (scale 0.0002490574820792943) ***** 
[1] xi(23)=0.999999940395 zeta(24)=0
Stats: nr=69 min=1 max=656 sum3=1169945404
Safeguard: keep additional 3 states
Emax=10.0024201944 nrkept=1967 nrkeptmult=1967 
Kept: 1967 out of 7768, ratio=0.253
Memory used: 218 MB  Time elapsed: 259 s

 ***** [DMNRG] Iteration 25/66 (scale 0.00017611023448351604) ***** 
[1] xi(24)=0.999999970198 zeta(25)=0
Stats: nr=61 min=1 max=612 sum3=1248954038
Emax=10.0035392215 nrkept=1713 nrkeptmult=1713 
Kept: 1713 out of 7868, ratio=0.218
Memory used: 218 MB  Time elapsed: 260 s

 ***** [DMNRG] Iteration 26/66 (scale 0.00012452874103964715) ***** 
[1] xi(25)=0.999999985099 zeta(26)=0
Stats: nr=69 min=1 max=581 sum3=816220428
Safeguard: keep additional 4 states
Emax=10.0002797499 nrkept=1680 nrkeptmult=1680 
Kept: 1680 out of 6852, ratio=0.245
Memory used: 216 MB  Time elapsed: 260 s

 ***** [DMNRG] Iteration 27/66 (scale 8.8055117241758019e-05) ***** 
[1] xi(26)=0.999999992549 zeta(27)=0
Stats: nr=61 min=1 max=548 sum3=798076332
Safeguard: keep additional 4 states
Emax=10.0018814628 nrkept=1477 nrkeptmult=1477 
Kept: 1477 out of 6720, ratio=0.22
Memory used: 218 MB  Time elapsed: 261 s

 ***** [DMNRG] Iteration 28/66 (scale 6.2264370519823575e-05) ***** 
[1] xi(27)=0.999999996275 zeta(28)=0
Stats: nr=63 min=1 max=506 sum3=547457896
Safeguard: keep additional 3 states
Emax=10.0069144664 nrkept=1420 nrkeptmult=1420 
Kept: 1420 out of 5908, ratio=0.24
Memory used: 218 MB  Time elapsed: 261 s

 ***** [DMNRG] Iteration 29/66 (scale 4.402755862087901e-05) ***** 
[1] xi(28)=0.999999998137 zeta(29)=0
Stats: nr=59 min=1 max=488 sum3=514851298
Safeguard: keep additional 3 states
Emax=10.0026862654 nrkept=1243 nrkeptmult=1243 
Kept: 1243 out of 5680, ratio=0.219
Memory used: 223 MB  Time elapsed: 262 s

 ***** [DMNRG] Iteration 30/66 (scale 3.1132185259911788e-05) ***** 
[1] xi(29)=0.999999999069 zeta(30)=0
Stats: nr=58 min=1 max=428 sum3=348701728
Safeguard: keep additional 1 states
Emax=10.0005256635 nrkept=1214 nrkeptmult=1214 
Kept: 1214 out of 4972, ratio=0.244
Memory used: 223 MB  Time elapsed: 262 s

 ***** [DMNRG] Iteration 31/66 (scale 2.2013779310439505e-05) ***** 
[1] xi(30)=0.999999999534 zeta(31)=0
Stats: nr=57 min=1 max=437 sum3=337641950
Safeguard: keep additional 1 states
Emax=10.0073109725 nrkept=1103 nrkeptmult=1103 
Kept: 1103 out of 4856, ratio=0.227
Memory used: 223 MB  Time elapsed: 262 s

 ***** [DMNRG] Iteration 32/66 (scale 1.5566092629955894e-05) ***** 
[1] xi(31)=0.999999999767 zeta(32)=0
Stats: nr=58 min=1 max=375 sum3=247813640
Safeguard: keep additional 1 states
Emax=10.0026675186 nrkept=1132 nrkeptmult=1132 
Kept: 1132 out of 4412, ratio=0.257
Memory used: 223 MB  Time elapsed: 262 s

 ***** [DMNRG] Iteration 33/66 (scale 1.1006889655219752e-05) ***** 
[1] xi(32)=0.999999999884 zeta(33)=0
Stats: nr=57 min=1 max=416 sum3=282675538
Safeguard: keep additional 1 states
Emax=10.0070238801 nrkept=1079 nrkeptmult=1079 
Kept: 1079 out of 4528, ratio=0.238
Memory used: 225 MB  Time elapsed: 263 s

 ***** [DMNRG] Iteration 34/66 (scale 7.7830463149779469e-06) ***** 
[1] xi(33)=0.999999999942 zeta(34)=0
Stats: nr=58 min=1 max=371 sum3=235967204
Safeguard: keep additional 1 states
Emax=10.030116036 nrkept=1116 nrkeptmult=1116 
Kept: 1116 out of 4316, ratio=0.259
Memory used: 225 MB  Time elapsed: 263 s

 ***** [DMNRG] Iteration 35/66 (scale 5.5034448276098762e-06) ***** 
[1] xi(34)=0.999999999971 zeta(35)=0
Stats: nr=57 min=1 max=412 sum3=271853526
Safeguard: keep additional 1 states
Emax=10.0355827293 nrkept=1075 nrkeptmult=1075 
Kept: 1075 out of 4464, ratio=0.241
Memory used: 225 MB  Time elapsed: 263 s

 ***** [DMNRG] Iteration 36/66 (scale 3.8915231574889735e-06) ***** 
[1] xi(35)=0.999999999985 zeta(36)=0
Stats: nr=58 min=1 max=371 sum3=234091648
Safeguard: keep additional 2 states
Emax=10.0584490488 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4300, ratio=0.26
Memory used: 225 MB  Time elapsed: 263 s

 ***** [DMNRG] Iteration 37/66 (scale 2.7517224138049381e-06) ***** 
[1] xi(36)=0.999999999993 zeta(37)=0
Stats: nr=57 min=1 max=413 sum3=273706324
Safeguard: keep additional 1 states
Emax=10.0017711281 nrkept=1067 nrkeptmult=1067 
Kept: 1067 out of 4468, ratio=0.239
Memory used: 225 MB  Time elapsed: 264 s

 ***** [DMNRG] Iteration 38/66 (scale 1.9457615787444867e-06) ***** 
[1] xi(37)=0.999999999996 zeta(38)=0
Stats: nr=58 min=1 max=371 sum3=231996350
Emax=10.0694101108 nrkept=1115 nrkeptmult=1115 
Kept: 1115 out of 4268, ratio=0.261
Memory used: 227 MB  Time elapsed: 264 s

 ***** [DMNRG] Iteration 39/66 (scale 1.3758612069024691e-06) ***** 
[1] xi(38)=0.999999999998 zeta(39)=0
Stats: nr=57 min=1 max=413 sum3=272139644
Emax=10.0056891737 nrkept=1060 nrkeptmult=1060 
Kept: 1060 out of 4460, ratio=0.238
Memory used: 229 MB  Time elapsed: 264 s

 ***** [DMNRG] Iteration 40/66 (scale 9.7288078937224336e-07) ***** 
[1] xi(39)=0.999999999999 zeta(40)=0
Stats: nr=58 min=1 max=371 sum3=231642436
Safeguard: keep additional 2 states
Emax=10.0731672522 nrkept=1117 nrkeptmult=1117 
Kept: 1117 out of 4240, ratio=0.263
Memory used: 229 MB  Time elapsed: 264 s

 ***** [DMNRG] Iteration 41/66 (scale 6.8793060345123453e-07) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Emax=10.014400013 nrkept=1060 nrkeptmult=1060 
Kept: 1060 out of 4468, ratio=0.237
Memory used: 229 MB  Time elapsed: 265 s

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

** Iteration completed.

Spectrum: A_d-A_d-d DMNRG
[A_d-A_d-d] pos=0.135275812054 neg=0.864827687685 sum= 1.00010349974
mu1=-0.00379353215009 mu2=0.000662526497817 mu3=-5.0544834227e-06 mu4=0.000189844855361

Spectrum: A_d-A_d-d CFS
[A_d-A_d-d] pos=0.135370376023 neg=0.864629623976 sum= 0.999999999999
mu1=-0.00385029282228 mu2=0.000666150362682 mu3=-5.11294318925e-06 mu4=0.000191128693723

Spectrum: A_d-A_d-d FDM
[A_d-A_d-d] pos=0.135370376023 neg=0.864629623976 sum= 0.999999999998
mu1=-0.00385029282228 mu2=0.000666150362682 mu3=-5.11294318925e-06 mu4=0.000191128693723

Spectrum: A_d-A_d-u DMNRG
[A_d-A_d-u] pos=0.864827687685 neg=0.135275812054 sum= 1.00010349974
mu1=0.00379353215009 mu2=0.000662526497817 mu3=5.0544834227e-06 mu4=0.000189844855361

Spectrum: A_d-A_d-u CFS
[A_d-A_d-u] pos=0.864629623976 neg=0.135370376023 sum= 0.999999999999
mu1=0.00385029282228 mu2=0.000666150362682 mu3=5.11294318925e-06 mu4=0.000191128693723

Spectrum: A_d-A_d-u FDM
[A_d-A_d-u] pos=0.864629623976 neg=0.135370376023 sum= 0.999999999998
mu1=0.00385029282228 mu2=0.000666150362682 mu3=5.11294318925e-06 mu4=0.000191128693723

Master done. Terminating slave processes.
Master exiting.

Memory usage report (process 0)
===================
  after diag: 430120 kB
after recalc: 521064 kB
 after split: 489372 kB
  after trim: 271004 kB
 after trunc: 430508 kB
      diag_h: 342096 kB

Peak usage: 521064 kB

Timing report (process 0)
=============
         All: 272 s
          DM: 3 s
         FDM: 17 s
        diag: 159 s
    recalc d: 23 s
      spec d: 42 s
       Other: 25 s
