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.5
discretization=Z
band=flat
betabar=1
T=1e-10
T_ratio=1
Tmin_ratio=-999
keep=10000
keepenergy=10
keepmin=0
keepblock=0
keepzrescale=0
log=
logall=0
logenumber=10
smooth=new
alpha=0.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.5
>>  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.5
# keep=10000 keepenergy=10 Nmax=67 channels=1 dots=1 betabar=1

Number of invariant subspaces: 9
(-2 1) 1 states
0.036753980447 
(-1 0) 2 states
0.0173892632005 0.0525489897136 
(-1 2) 2 states
0.0175034869175 0.0526427101507 
(0 -1) 1 states
0.0331842724671 
(0 1) 4 states
0 0.0332882751662 0.036753980447 0.0700421943691 
(0 3) 1 states
0.0333922166212 
(1 0) 2 states
0.0173892632005 0.0525489897136 
(1 2) 2 states
0.0175034869175 0.0526427101507 
(2 1) 1 states
0.036753980447 
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.60367320798e-05

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

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

 ***** [NRG] Iteration 3/67 (scale 0.51006972329839473) ***** 
[1] xi(2)=0.992903283056 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Emax=9.70913435997 nrkept=1024 nrkeptmult=1024 
Total energy=-2.47374244809  Egs=-1.58037403374
Kept: 1024 out of 1024, ratio=1
<SZd>=-7.38527321898e-05
Memory used: 17 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/67 (scale 0.36067376022224085) ***** 
[1] xi(3)=1.04266411315 zeta(4)=0
Stats: nr=49 min=1 max=400 sum3=230553856
Emax=10.2107510944 nrkept=3041 nrkeptmult=3041 
Total energy=-2.79981386665  Egs=-0.90406193775
Kept: 3041 out of 4096, ratio=0.742
<SZd>=-0.000104625235562
Memory used: 44 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/67 (scale 0.25503486164919736) ***** 
[1] xi(4)=1.01328740971 zeta(5)=0
Stats: nr=64 min=1 max=890 sum3=4292520476
Safeguard: keep additional 3 states
Emax=10.2083007966 nrkept=4672 nrkeptmult=4672 
Total energy=-3.19014829886  Egs=-1.53051402343
Kept: 4672 out of 12164, ratio=0.384
<SZd>=-0.000150185273828
Memory used: 80 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 6/67 (scale 0.18033688011112042) ***** 
[1] xi(5)=0.999788616061 zeta(6)=0
Stats: nr=77 min=1 max=1468 sum3=13977315586
Safeguard: keep additional 5 states
Emax=10.2375187984 nrkept=5286 nrkeptmult=5286 
Total energy=-3.34861460312  Egs=-0.878723776051
Kept: 5286 out of 18688, ratio=0.283
<SZd>=-0.000210441287198
Memory used: 116 MB  Time elapsed: 19 s

 ***** [NRG] Iteration 7/67 (scale 0.12751743082459868) ***** 
[1] xi(6)=0.998344002082 zeta(7)=0
Stats: nr=80 min=1 max=1482 sum3=20486955552
Safeguard: keep additional 5 states
Emax=10.020384773 nrkept=4746 nrkeptmult=4746 
Total energy=-3.53895479685  Egs=-1.49266019952
Kept: 4746 out of 21144, ratio=0.224
<SZd>=-0.000298931471219
Memory used: 116 MB  Time elapsed: 36 s

 ***** [NRG] Iteration 8/67 (scale 0.090168440055560212) ***** 
[1] xi(7)=0.998902918682 zeta(8)=0
Stats: nr=77 min=1 max=1514 sum3=15200412618
Safeguard: keep additional 2 states
Emax=10.000003643 nrkept=5004 nrkeptmult=5004 
Total energy=-3.62044498817  Egs=-0.903755141607
Kept: 5004 out of 18984, ratio=0.264
<SZd>=-0.000434897804195
Memory used: 130 MB  Time elapsed: 50 s

 ***** [NRG] Iteration 9/67 (scale 0.063758715412299341) ***** 
[1] xi(8)=0.999398551363 zeta(9)=0
Stats: nr=80 min=1 max=1411 sum3=17486360544
Safeguard: keep additional 3 states
Emax=10.0005621971 nrkept=4112 nrkeptmult=4112 
Total energy=-3.71367895465  Egs=-1.46229367835
Kept: 4112 out of 20016, ratio=0.205
<SZd>=-0.00061339268041
Memory used: 114 MB  Time elapsed: 66 s

 ***** [NRG] Iteration 10/67 (scale 0.045084220027780106) ***** 
[1] xi(9)=0.99968731683 zeta(10)=0
Stats: nr=77 min=1 max=1320 sum3=9984316544
Emax=10.0474091441 nrkept=4641 nrkeptmult=4641 
Total energy=-3.75566363755  Egs=-0.931250066507
Kept: 4641 out of 16448, ratio=0.282
<SZd>=-0.00086961326313
Memory used: 127 MB  Time elapsed: 76 s

 ***** [NRG] Iteration 11/67 (scale 0.031879357706149671) ***** 
[1] xi(10)=0.9998407918 zeta(11)=0
Stats: nr=76 min=1 max=1324 sum3=14394708342
Emax=10.0000088446 nrkept=3945 nrkeptmult=3945 
Total energy=-3.80121308538  Egs=-1.42880694936
Kept: 3945 out of 18564, ratio=0.213
<SZd>=-0.00130417691495
Memory used: 114 MB  Time elapsed: 89 s

 ***** [NRG] Iteration 12/67 (scale 0.022542110013890053) ***** 
[1] xi(11)=0.99991969299 zeta(12)=0
Stats: nr=77 min=1 max=1269 sum3=8848597338
Safeguard: keep additional 1 states
Emax=10.0000173713 nrkept=4320 nrkeptmult=4320 
Total energy=-3.82280935235  Egs=-0.958041059532
Kept: 4320 out of 15780, ratio=0.274
<SZd>=-0.00184334664986
Memory used: 119 MB  Time elapsed: 98 s

 ***** [NRG] Iteration 13/67 (scale 0.015939678853074835) ***** 
[1] xi(12)=0.999959672404 zeta(13)=0
Stats: nr=74 min=1 max=1241 sum3=11685770784
Emax=10.0019005828 nrkept=3797 nrkeptmult=3797 
Total energy=-3.8450573653  Egs=-1.39576293625
Kept: 3797 out of 17280, ratio=0.22
<SZd>=-0.00263229468661
Memory used: 115 MB  Time elapsed: 110 s

 ***** [NRG] Iteration 14/67 (scale 0.011271055006945026) ***** 
[1] xi(13)=0.999979792879 zeta(14)=0
Stats: nr=77 min=1 max=1222 sum3=7902449270
Safeguard: keep additional 3 states
Emax=10.0002943227 nrkept=3908 nrkeptmult=3908 
Total energy=-3.85611808674  Egs=-0.981338608607
Kept: 3908 out of 15188, ratio=0.257
<SZd>=-0.00390367907422
Memory used: 116 MB  Time elapsed: 118 s

 ***** [NRG] Iteration 15/67 (scale 0.0079698394265374176) ***** 
[1] xi(14)=0.999989885634 zeta(15)=0
Stats: nr=68 min=1 max=1143 sum3=8940336932
Safeguard: keep additional 4 states
Emax=10.0009686162 nrkept=3679 nrkeptmult=3679 
Total energy=-3.86701067704  Egs=-1.36672644404
Kept: 3679 out of 15632, ratio=0.235
<SZd>=-0.0056830540337
Memory used: 114 MB  Time elapsed: 127 s

 ***** [NRG] Iteration 16/67 (scale 0.0056355275034725132) ***** 
[1] xi(15)=0.999994940118 zeta(16)=0
Stats: nr=73 min=1 max=1195 sum3=7343343466
Safeguard: keep additional 2 states
Emax=10.0012083581 nrkept=3553 nrkeptmult=3553 
Total energy=-3.87264082083  Egs=-0.999044682671
Kept: 3553 out of 14716, ratio=0.241
<SZd>=-0.00846123709008
Memory used: 112 MB  Time elapsed: 135 s

 ***** [NRG] Iteration 17/67 (scale 0.0039849197132687088) ***** 
[1] xi(16)=0.999997469385 zeta(17)=0
Stats: nr=64 min=1 max=1045 sum3=6899597662
Safeguard: keep additional 3 states
Emax=10.0004418653 nrkept=3195 nrkeptmult=3195 
Total energy=-3.87799653839  Egs=-1.34399635315
Kept: 3195 out of 14212, ratio=0.225
<SZd>=-0.0125808808358
Memory used: 111 MB  Time elapsed: 142 s

 ***** [NRG] Iteration 18/67 (scale 0.0028177637517362566) ***** 
[1] xi(17)=0.999998734524 zeta(18)=0
Stats: nr=73 min=1 max=1057 sum3=4982574480
Safeguard: keep additional 1 states
Emax=10.0014809796 nrkept=3134 nrkeptmult=3134 
Total energy=-3.88084419536  Egs=-1.01060884494
Kept: 3134 out of 12780, ratio=0.245
<SZd>=-0.0185902120813
Memory used: 112 MB  Time elapsed: 148 s

 ***** [NRG] Iteration 19/67 (scale 0.0019924598566343544) ***** 
[1] xi(18)=0.99999936722 zeta(19)=0
Stats: nr=62 min=1 max=932 sum3=4820273528
Safeguard: keep additional 3 states
Emax=10.0016710228 nrkept=2799 nrkeptmult=2799 
Total energy=-3.88349214401  Egs=-1.32898469075
Kept: 2799 out of 12536, ratio=0.223
<SZd>=-0.0270922814349
Memory used: 110 MB  Time elapsed: 154 s

 ***** [NRG] Iteration 20/67 (scale 0.0014088818758681283) ***** 
[1] xi(19)=0.999999683599 zeta(20)=0
Stats: nr=73 min=1 max=928 sum3=3381045138
Safeguard: keep additional 1 states
Emax=10.0011257929 nrkept=2695 nrkeptmult=2695 
Total energy=-3.8849233925  Egs=-1.01587543744
Kept: 2695 out of 11196, ratio=0.241
<SZd>=-0.0387150439244
Memory used: 109 MB  Time elapsed: 158 s

 ***** [NRG] Iteration 21/67 (scale 0.0009962299283171772) ***** 
[1] xi(20)=0.999999841797 zeta(21)=0
Stats: nr=62 min=1 max=810 sum3=3106993828
Emax=10.0053709678 nrkept=2321 nrkeptmult=2321 
Total energy=-3.88624151253  Egs=-1.32310824813
Kept: 2321 out of 10780, ratio=0.215
<SZd>=-0.0542082587967
Memory used: 106 MB  Time elapsed: 162 s

 ***** [NRG] Iteration 22/67 (scale 0.00070444093793406415) ***** 
[1] xi(21)=0.999999920898 zeta(22)=0
Stats: nr=73 min=1 max=772 sum3=1938189182
Emax=10.0014289842 nrkept=2159 nrkeptmult=2159 
Total energy=-3.88695647272  Egs=-1.01493275783
Kept: 2159 out of 9284, ratio=0.233
<SZd>=-0.0747283199793
Memory used: 104 MB  Time elapsed: 164 s

 ***** [NRG] Iteration 23/67 (scale 0.0004981149641585886) ***** 
[1] xi(22)=0.999999960449 zeta(23)=0
Stats: nr=62 min=1 max=673 sum3=1644308414
Safeguard: keep additional 2 states
Emax=10.0010886235 nrkept=2024 nrkeptmult=2024 
Total energy=-3.88761720787  Egs=-1.3264712043
Kept: 2024 out of 8636, ratio=0.234
<SZd>=-0.101420395094
Memory used: 106 MB  Time elapsed: 166 s

 ***** [NRG] Iteration 24/67 (scale 0.00035222046896703208) ***** 
[1] xi(23)=0.999999980224 zeta(24)=0
Stats: nr=71 min=1 max=676 sum3=1296640364
Safeguard: keep additional 1 states
Emax=10.0041723498 nrkept=2020 nrkeptmult=2020 
Total energy=-3.88797285721  Egs=-1.00973499089
Kept: 2020 out of 8096, ratio=0.25
<SZd>=-0.13542716272
Memory used: 109 MB  Time elapsed: 168 s

 ***** [NRG] Iteration 25/67 (scale 0.0002490574820792943) ***** 
[1] xi(24)=0.999999990112 zeta(25)=0
Stats: nr=62 min=1 max=625 sum3=1331708212
Safeguard: keep additional 2 states
Emax=10.0013287607 nrkept=1819 nrkeptmult=1819 
Total energy=-3.88830547062  Egs=-1.33548853165
Kept: 1819 out of 8080, ratio=0.225
<SZd>=-0.177038314444
Memory used: 109 MB  Time elapsed: 170 s

 ***** [NRG] Iteration 26/67 (scale 0.00017611023448351604) ***** 
[1] xi(25)=0.999999995056 zeta(26)=0
Stats: nr=69 min=1 max=616 sum3=968321290
Safeguard: keep additional 3 states
Emax=10.0066392399 nrkept=1821 nrkeptmult=1821 
Total energy=-3.88848221432  Egs=-1.00359697744
Kept: 1821 out of 7276, ratio=0.25
<SZd>=-0.224037353522
Memory used: 112 MB  Time elapsed: 172 s

 ***** [NRG] Iteration 27/67 (scale 0.00012452874103964715) ***** 
[1] xi(26)=0.999999997528 zeta(27)=0
Stats: nr=61 min=1 max=582 sum3=1015194390
Safeguard: keep additional 3 states
Emax=10.001711823 nrkept=1604 nrkeptmult=1604 
Total energy=-3.88864967037  Egs=-1.34471814367
Kept: 1604 out of 7284, ratio=0.22
<SZd>=-0.271584675069
Memory used: 115 MB  Time elapsed: 174 s

 ***** [NRG] Iteration 28/67 (scale 8.8055117241758019e-05) ***** 
[1] xi(27)=0.999999998764 zeta(28)=0
Stats: nr=67 min=1 max=546 sum3=681138488
Emax=10.0073255952 nrkept=1537 nrkeptmult=1537 
Total energy=-3.88873761264  Egs=-0.998718424627
Kept: 1537 out of 6416, ratio=0.24
<SZd>=-0.312040153465
Memory used: 113 MB  Time elapsed: 175 s

 ***** [NRG] Iteration 29/67 (scale 6.2264370519823575e-05) ***** 
[1] xi(28)=0.999999999382 zeta(29)=0
Stats: nr=59 min=1 max=516 sum3=642098572
Safeguard: keep additional 3 states
Emax=10.0052533598 nrkept=1353 nrkeptmult=1353 
Total energy=-3.88882174549  Egs=-1.3512197934
Kept: 1353 out of 6148, ratio=0.22
<SZd>=-0.33972861901
Memory used: 112 MB  Time elapsed: 176 s

 ***** [NRG] Iteration 30/67 (scale 4.402755862087901e-05) ***** 
[1] xi(29)=0.999999999691 zeta(30)=0
Stats: nr=59 min=1 max=464 sum3=435097302
Safeguard: keep additional 1 states
Emax=10.0050895298 nrkept=1297 nrkeptmult=1297 
Total energy=-3.88886557798  Egs=-0.995569440489
Kept: 1297 out of 5412, ratio=0.24
<SZd>=-0.35433982602
Memory used: 117 MB  Time elapsed: 177 s

 ***** [NRG] Iteration 31/67 (scale 3.1132185259911788e-05) ***** 
[1] xi(30)=0.999999999846 zeta(31)=0
Stats: nr=57 min=1 max=457 sum3=402283822
Emax=10.0096184328 nrkept=1151 nrkeptmult=1151 
Total energy=-3.88890776218  Egs=-1.35500261271
Kept: 1151 out of 5188, ratio=0.222
<SZd>=-0.360530896658
Memory used: 117 MB  Time elapsed: 178 s

 ***** [NRG] Iteration 32/67 (scale 2.2013779310439505e-05) ***** 
[1] xi(31)=0.999999999923 zeta(32)=0
Stats: nr=58 min=1 max=393 sum3=281410544
Safeguard: keep additional 1 states
Emax=10.0001119128 nrkept=1153 nrkeptmult=1153 
Total energy=-3.88892963876  Egs=-0.993767951346
Kept: 1153 out of 4604, ratio=0.25
<SZd>=-0.3628581264
Memory used: 117 MB  Time elapsed: 178 s

 ***** [NRG] Iteration 33/67 (scale 1.5566092629955894e-05) ***** 
[1] xi(32)=0.999999999961 zeta(33)=0
Stats: nr=57 min=1 max=425 sum3=297934786
Emax=10.0074138997 nrkept=1086 nrkeptmult=1086 
Total energy=-3.88895076223  Egs=-1.35701799941
Kept: 1086 out of 4612, ratio=0.235
<SZd>=-0.363809738241
Memory used: 120 MB  Time elapsed: 179 s

 ***** [NRG] Iteration 34/67 (scale 1.1006889655219752e-05) ***** 
[1] xi(33)=0.999999999981 zeta(34)=0
Stats: nr=58 min=1 max=371 sum3=238660776
Safeguard: keep additional 1 states
Emax=10.004630111 nrkept=1116 nrkeptmult=1116 
Total energy=-3.88896168992  Egs=-0.992804188187
Kept: 1116 out of 4344, ratio=0.257
<SZd>=-0.364239820138
Memory used: 120 MB  Time elapsed: 179 s

 ***** [NRG] Iteration 35/67 (scale 7.7830463149779469e-06) ***** 
[1] xi(34)=0.99999999999 zeta(35)=0
Stats: nr=57 min=1 max=412 sum3=271853526
Safeguard: keep additional 1 states
Emax=10.0135740387 nrkept=1075 nrkeptmult=1075 
Total energy=-3.88897225971  Egs=-1.35805318287
Kept: 1075 out of 4464, ratio=0.241
<SZd>=-0.364446695885
Memory used: 120 MB  Time elapsed: 180 s

 ***** [NRG] Iteration 36/67 (scale 5.5034448276098762e-06) ***** 
[1] xi(35)=0.999999999995 zeta(36)=0
Stats: nr=58 min=1 max=371 sum3=234091648
Safeguard: keep additional 1 states
Emax=10.0472729407 nrkept=1116 nrkeptmult=1116 
Total energy=-3.88897772081  Egs=-0.992306014983
Kept: 1116 out of 4300, ratio=0.26
<SZd>=-0.364548186684
Memory used: 120 MB  Time elapsed: 181 s

 ***** [NRG] Iteration 37/67 (scale 3.8915231574889735e-06) ***** 
[1] xi(36)=0.999999999998 zeta(37)=0
Stats: nr=57 min=1 max=412 sum3=271853526
Safeguard: keep additional 3 states
Emax=10.0033311223 nrkept=1073 nrkeptmult=1073 
Total energy=-3.88898300774  Egs=-1.35857701231
Kept: 1073 out of 4464, ratio=0.24
<SZd>=-0.364598469599
Memory used: 122 MB  Time elapsed: 181 s

 ***** [NRG] Iteration 38/67 (scale 2.7517224138049381e-06) ***** 
[1] xi(37)=0.999999999999 zeta(38)=0
Stats: nr=58 min=1 max=371 sum3=232464242
Emax=10.0661469794 nrkept=1115 nrkeptmult=1115 
Total energy=-3.8889857376  Egs=-0.99205282323
Kept: 1115 out of 4292, ratio=0.26
<SZd>=-0.364623597731
Memory used: 122 MB  Time elapsed: 182 s

 ***** [NRG] Iteration 39/67 (scale 1.9457615787444867e-06) ***** 
[1] xi(38)=0.999999999999 zeta(39)=0
Stats: nr=57 min=1 max=413 sum3=272139644
Safeguard: keep additional 2 states
Emax=10.0037288731 nrkept=1063 nrkeptmult=1063 
Total energy=-3.88898838158  Egs=-1.35884039409
Kept: 1063 out of 4460, ratio=0.238
<SZd>=-0.364636045743
Memory used: 123 MB  Time elapsed: 182 s

 ***** [NRG] Iteration 40/67 (scale 1.3758612069024691e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=58 min=1 max=371 sum3=231767230
Emax=10.0717107269 nrkept=1115 nrkeptmult=1115 
Total energy=-3.88898974633  Egs=-0.991925200076
Kept: 1115 out of 4252, ratio=0.262
<SZd>=-0.364642402507
Memory used: 123 MB  Time elapsed: 183 s

 ***** [NRG] Iteration 41/67 (scale 9.7288078937224336e-07) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=57 min=1 max=413 sum3=272139644
Emax=10.0108912388 nrkept=1060 nrkeptmult=1060 
Total energy=-3.88899106845  Egs=-1.35897243948
Kept: 1060 out of 4460, ratio=0.238
<SZd>=-0.364645524214
Memory used: 124 MB  Time elapsed: 184 s

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

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

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

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

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

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

 ***** [NRG] Iteration 48/67 (scale 8.5991325431404316e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=58 min=1 max=371 sum3=231782012
Safeguard: keep additional 2 states
Emax=10.0759318791 nrkept=1117 nrkeptmult=1117 
Total energy=-3.88899350472  Egs=-0.991804931645
Kept: 1117 out of 4244, ratio=0.263
<SZd>=-0.364648710336
Memory used: 126 MB  Time elapsed: 187 s

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

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

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

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

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

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

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

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

 ***** [NRG] Iteration 57/67 (scale 3.8003155834853256e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190757309 nrkept=1061 nrkeptmult=1061 
Total energy=-3.88899374479  Egs=-1.35910419939
Kept: 1061 out of 4468, ratio=0.237
<SZd>=-0.364648741369
Memory used: 128 MB  Time elapsed: 192 s

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

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

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

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

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

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

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

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

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

 ***** [NRG] Iteration 67/67 (scale 1.1875986198391643e-10) ***** 
[1] xi(66)=1 zeta(67)=0
Stats: nr=57 min=1 max=413 sum3=273116836
Safeguard: keep additional 1 states
Emax=10.0190975755 nrkept=4468 nrkeptmult=4468 
Total energy=-3.88899375496  Egs=-1.35910470051
Kept: 4468 out of 4468, ratio=1
<SZd>=-0.36464874178
Memory used: 168 MB  Time elapsed: 198 s

** Iteration completed.

Spectrum: A_d-A_d-d FT
[A_d-A_d-d] pos=0.103477750261 neg=0.908102284999 sum= 1.01158003526
mu1=-0.00427997079676 mu2=0.00066301103522 mu3=-6.65144849313e-06 mu4=0.000226420127771

Spectrum: A_d-A_d-u FT
[A_d-A_d-u] pos=0.908102284999 neg=0.103477750261 sum= 1.01158003526
mu1=0.00427997079677 mu2=0.00066301103522 mu3=6.65144849315e-06 mu4=0.000226420127771


Total energy:             -3.888993754957126
TKW=0
Grand canonical Z=              6.963142981660885
Zft (last shell)=               6.963122372452432
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.5
>>  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.036753980447 
(-1 0) 2 states
0.0173892632005 0.0525489897136 
(-1 2) 2 states
0.0175034869175 0.0526427101507 
(0 -1) 1 states
0.0331842724671 
(0 1) 4 states
0 0.0332882751662 0.036753980447 0.0700421943691 
(0 3) 1 states
0.0333922166212 
(1 0) 2 states
0.0173892632005 0.0525489897136 
(1 2) 2 states
0.0175034869175 0.0526427101507 
(2 1) 1 states
0.036753980447 
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.0201394465967895) ***** 
[1] xi(0)=0.583569837945 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=2.34140558228 nrkept=64 nrkeptmult=64 
Kept: 64 out of 64, ratio=1
Memory used: 189 MB  Time elapsed: 223 s

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

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

 ***** [DMNRG] Iteration 4/67 (scale 0.36067376022224085) ***** 
[1] xi(3)=1.04266411315 zeta(4)=0
Stats: nr=49 min=1 max=400 sum3=230553856
Emax=10.2107510944 nrkept=3041 nrkeptmult=3041 
Kept: 3041 out of 4096, ratio=0.742
Memory used: 201 MB  Time elapsed: 223 s

 ***** [DMNRG] Iteration 5/67 (scale 0.25503486164919736) ***** 
[1] xi(4)=1.01328740971 zeta(5)=0
Stats: nr=64 min=1 max=890 sum3=4292520476
Safeguard: keep additional 3 states
Emax=10.2083007966 nrkept=4672 nrkeptmult=4672 
Kept: 4672 out of 12164, ratio=0.384
Memory used: 233 MB  Time elapsed: 225 s

 ***** [DMNRG] Iteration 6/67 (scale 0.18033688011112042) ***** 
[1] xi(5)=0.999788616061 zeta(6)=0
Stats: nr=77 min=1 max=1468 sum3=13977315586
Safeguard: keep additional 5 states
Emax=10.2375187984 nrkept=5286 nrkeptmult=5286 
Kept: 5286 out of 18688, ratio=0.283
Memory used: 257 MB  Time elapsed: 230 s

 ***** [DMNRG] Iteration 7/67 (scale 0.12751743082459868) ***** 
[1] xi(6)=0.998344002082 zeta(7)=0
Stats: nr=80 min=1 max=1482 sum3=20486955552
Safeguard: keep additional 5 states
Emax=10.020384773 nrkept=4746 nrkeptmult=4746 
Kept: 4746 out of 21144, ratio=0.224
Memory used: 242 MB  Time elapsed: 236 s

 ***** [DMNRG] Iteration 8/67 (scale 0.090168440055560212) ***** 
[1] xi(7)=0.998902918682 zeta(8)=0
Stats: nr=77 min=1 max=1514 sum3=15200412618
Safeguard: keep additional 2 states
Emax=10.000003643 nrkept=5004 nrkeptmult=5004 
Kept: 5004 out of 18984, ratio=0.264
Memory used: 249 MB  Time elapsed: 241 s

 ***** [DMNRG] Iteration 9/67 (scale 0.063758715412299341) ***** 
[1] xi(8)=0.999398551363 zeta(9)=0
Stats: nr=80 min=1 max=1411 sum3=17486360544
Safeguard: keep additional 3 states
Emax=10.0005621971 nrkept=4112 nrkeptmult=4112 
Kept: 4112 out of 20016, ratio=0.205
Memory used: 232 MB  Time elapsed: 246 s

 ***** [DMNRG] Iteration 10/67 (scale 0.045084220027780106) ***** 
[1] xi(9)=0.99968731683 zeta(10)=0
Stats: nr=77 min=1 max=1320 sum3=9984316544
Emax=10.0474091441 nrkept=4641 nrkeptmult=4641 
Kept: 4641 out of 16448, ratio=0.282
Memory used: 245 MB  Time elapsed: 249 s

 ***** [DMNRG] Iteration 11/67 (scale 0.031879357706149671) ***** 
[1] xi(10)=0.9998407918 zeta(11)=0
Stats: nr=76 min=1 max=1324 sum3=14394708342
Emax=10.0000088446 nrkept=3945 nrkeptmult=3945 
Kept: 3945 out of 18564, ratio=0.213
Memory used: 238 MB  Time elapsed: 253 s

 ***** [DMNRG] Iteration 12/67 (scale 0.022542110013890053) ***** 
[1] xi(11)=0.99991969299 zeta(12)=0
Stats: nr=77 min=1 max=1269 sum3=8848597338
Safeguard: keep additional 1 states
Emax=10.0000173713 nrkept=4320 nrkeptmult=4320 
Kept: 4320 out of 15780, ratio=0.274
Memory used: 242 MB  Time elapsed: 257 s

 ***** [DMNRG] Iteration 13/67 (scale 0.015939678853074835) ***** 
[1] xi(12)=0.999959672404 zeta(13)=0
Stats: nr=74 min=1 max=1241 sum3=11685770784
Emax=10.0019005828 nrkept=3797 nrkeptmult=3797 
Kept: 3797 out of 17280, ratio=0.22
Memory used: 234 MB  Time elapsed: 261 s

 ***** [DMNRG] Iteration 14/67 (scale 0.011271055006945026) ***** 
[1] xi(13)=0.999979792879 zeta(14)=0
Stats: nr=77 min=1 max=1222 sum3=7902449270
Safeguard: keep additional 3 states
Emax=10.0002943227 nrkept=3908 nrkeptmult=3908 
Kept: 3908 out of 15188, ratio=0.257
Memory used: 236 MB  Time elapsed: 264 s

 ***** [DMNRG] Iteration 15/67 (scale 0.0079698394265374176) ***** 
[1] xi(14)=0.999989885634 zeta(15)=0
Stats: nr=68 min=1 max=1143 sum3=8940336932
Safeguard: keep additional 4 states
Emax=10.0009686162 nrkept=3679 nrkeptmult=3679 
Kept: 3679 out of 15632, ratio=0.235
Memory used: 236 MB  Time elapsed: 267 s

 ***** [DMNRG] Iteration 16/67 (scale 0.0056355275034725132) ***** 
[1] xi(15)=0.999994940118 zeta(16)=0
Stats: nr=73 min=1 max=1195 sum3=7343343466
Safeguard: keep additional 2 states
Emax=10.0012083581 nrkept=3553 nrkeptmult=3553 
Kept: 3553 out of 14716, ratio=0.241
Memory used: 231 MB  Time elapsed: 269 s

 ***** [DMNRG] Iteration 17/67 (scale 0.0039849197132687088) ***** 
[1] xi(16)=0.999997469385 zeta(17)=0
Stats: nr=64 min=1 max=1045 sum3=6899597662
Safeguard: keep additional 3 states
Emax=10.0004418653 nrkept=3195 nrkeptmult=3195 
Kept: 3195 out of 14212, ratio=0.225
Memory used: 230 MB  Time elapsed: 272 s

 ***** [DMNRG] Iteration 18/67 (scale 0.0028177637517362566) ***** 
[1] xi(17)=0.999998734524 zeta(18)=0
Stats: nr=73 min=1 max=1057 sum3=4982574480
Safeguard: keep additional 1 states
Emax=10.0014809796 nrkept=3134 nrkeptmult=3134 
Kept: 3134 out of 12780, ratio=0.245
Memory used: 229 MB  Time elapsed: 274 s

 ***** [DMNRG] Iteration 19/67 (scale 0.0019924598566343544) ***** 
[1] xi(18)=0.99999936722 zeta(19)=0
Stats: nr=62 min=1 max=932 sum3=4820273528
Safeguard: keep additional 3 states
Emax=10.0016710228 nrkept=2799 nrkeptmult=2799 
Kept: 2799 out of 12536, ratio=0.223
Memory used: 228 MB  Time elapsed: 276 s

 ***** [DMNRG] Iteration 20/67 (scale 0.0014088818758681283) ***** 
[1] xi(19)=0.999999683599 zeta(20)=0
Stats: nr=73 min=1 max=928 sum3=3381045138
Safeguard: keep additional 1 states
Emax=10.0011257929 nrkept=2695 nrkeptmult=2695 
Kept: 2695 out of 11196, ratio=0.241
Memory used: 226 MB  Time elapsed: 278 s

 ***** [DMNRG] Iteration 21/67 (scale 0.0009962299283171772) ***** 
[1] xi(20)=0.999999841797 zeta(21)=0
Stats: nr=62 min=1 max=810 sum3=3106993828
Emax=10.0053709678 nrkept=2321 nrkeptmult=2321 
Kept: 2321 out of 10780, ratio=0.215
Memory used: 220 MB  Time elapsed: 279 s

 ***** [DMNRG] Iteration 22/67 (scale 0.00070444093793406415) ***** 
[1] xi(21)=0.999999920898 zeta(22)=0
Stats: nr=73 min=1 max=772 sum3=1938189182
Emax=10.0014289842 nrkept=2159 nrkeptmult=2159 
Kept: 2159 out of 9284, ratio=0.233
Memory used: 220 MB  Time elapsed: 280 s

 ***** [DMNRG] Iteration 23/67 (scale 0.0004981149641585886) ***** 
[1] xi(22)=0.999999960449 zeta(23)=0
Stats: nr=62 min=1 max=673 sum3=1644308414
Safeguard: keep additional 2 states
Emax=10.0010886235 nrkept=2024 nrkeptmult=2024 
Kept: 2024 out of 8636, ratio=0.234
Memory used: 215 MB  Time elapsed: 281 s

 ***** [DMNRG] Iteration 24/67 (scale 0.00035222046896703208) ***** 
[1] xi(23)=0.999999980224 zeta(24)=0
Stats: nr=71 min=1 max=676 sum3=1296640364
Safeguard: keep additional 1 states
Emax=10.0041723498 nrkept=2020 nrkeptmult=2020 
Kept: 2020 out of 8096, ratio=0.25
Memory used: 217 MB  Time elapsed: 282 s

 ***** [DMNRG] Iteration 25/67 (scale 0.0002490574820792943) ***** 
[1] xi(24)=0.999999990112 zeta(25)=0
Stats: nr=62 min=1 max=625 sum3=1331708212
Safeguard: keep additional 2 states
Emax=10.0013287607 nrkept=1819 nrkeptmult=1819 
Kept: 1819 out of 8080, ratio=0.225
Memory used: 215 MB  Time elapsed: 282 s

 ***** [DMNRG] Iteration 26/67 (scale 0.00017611023448351604) ***** 
[1] xi(25)=0.999999995056 zeta(26)=0
Stats: nr=69 min=1 max=616 sum3=968321290
Safeguard: keep additional 3 states
Emax=10.0066392399 nrkept=1821 nrkeptmult=1821 
Kept: 1821 out of 7276, ratio=0.25
Memory used: 216 MB  Time elapsed: 283 s

 ***** [DMNRG] Iteration 27/67 (scale 0.00012452874103964715) ***** 
[1] xi(26)=0.999999997528 zeta(27)=0
Stats: nr=61 min=1 max=582 sum3=1015194390
Safeguard: keep additional 3 states
Emax=10.001711823 nrkept=1604 nrkeptmult=1604 
Kept: 1604 out of 7284, ratio=0.22
Memory used: 215 MB  Time elapsed: 284 s

 ***** [DMNRG] Iteration 28/67 (scale 8.8055117241758019e-05) ***** 
[1] xi(27)=0.999999998764 zeta(28)=0
Stats: nr=67 min=1 max=546 sum3=681138488
Emax=10.0073255952 nrkept=1537 nrkeptmult=1537 
Kept: 1537 out of 6416, ratio=0.24
Memory used: 214 MB  Time elapsed: 284 s

 ***** [DMNRG] Iteration 29/67 (scale 6.2264370519823575e-05) ***** 
[1] xi(28)=0.999999999382 zeta(29)=0
Stats: nr=59 min=1 max=516 sum3=642098572
Safeguard: keep additional 3 states
Emax=10.0052533598 nrkept=1353 nrkeptmult=1353 
Kept: 1353 out of 6148, ratio=0.22
Memory used: 217 MB  Time elapsed: 285 s

 ***** [DMNRG] Iteration 30/67 (scale 4.402755862087901e-05) ***** 
[1] xi(29)=0.999999999691 zeta(30)=0
Stats: nr=59 min=1 max=464 sum3=435097302
Safeguard: keep additional 1 states
Emax=10.0050895298 nrkept=1297 nrkeptmult=1297 
Kept: 1297 out of 5412, ratio=0.24
Memory used: 217 MB  Time elapsed: 285 s

 ***** [DMNRG] Iteration 31/67 (scale 3.1132185259911788e-05) ***** 
[1] xi(30)=0.999999999846 zeta(31)=0
Stats: nr=57 min=1 max=457 sum3=402283822
Emax=10.0096184328 nrkept=1151 nrkeptmult=1151 
Kept: 1151 out of 5188, ratio=0.222
Memory used: 218 MB  Time elapsed: 285 s

 ***** [DMNRG] Iteration 32/67 (scale 2.2013779310439505e-05) ***** 
[1] xi(31)=0.999999999923 zeta(32)=0
Stats: nr=58 min=1 max=393 sum3=281410544
Safeguard: keep additional 1 states
Emax=10.0001119128 nrkept=1153 nrkeptmult=1153 
Kept: 1153 out of 4604, ratio=0.25
Memory used: 217 MB  Time elapsed: 286 s

 ***** [DMNRG] Iteration 33/67 (scale 1.5566092629955894e-05) ***** 
[1] xi(32)=0.999999999961 zeta(33)=0
Stats: nr=57 min=1 max=425 sum3=297934786
Emax=10.0074138997 nrkept=1086 nrkeptmult=1086 
Kept: 1086 out of 4612, ratio=0.235
Memory used: 218 MB  Time elapsed: 286 s

 ***** [DMNRG] Iteration 34/67 (scale 1.1006889655219752e-05) ***** 
[1] xi(33)=0.999999999981 zeta(34)=0
Stats: nr=58 min=1 max=371 sum3=238660776
Safeguard: keep additional 1 states
Emax=10.004630111 nrkept=1116 nrkeptmult=1116 
Kept: 1116 out of 4344, ratio=0.257
Memory used: 218 MB  Time elapsed: 286 s

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

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

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

 ***** [DMNRG] Iteration 38/67 (scale 2.7517224138049381e-06) ***** 
[1] xi(37)=0.999999999999 zeta(38)=0
Stats: nr=58 min=1 max=371 sum3=232464242
Emax=10.0661469794 nrkept=1115 nrkeptmult=1115 
Kept: 1115 out of 4292, ratio=0.26
Memory used: 218 MB  Time elapsed: 287 s

 ***** [DMNRG] Iteration 39/67 (scale 1.9457615787444867e-06) ***** 
[1] xi(38)=0.999999999999 zeta(39)=0
Stats: nr=57 min=1 max=413 sum3=272139644
Safeguard: keep additional 2 states
Emax=10.0037288731 nrkept=1063 nrkeptmult=1063 
Kept: 1063 out of 4460, ratio=0.238
Memory used: 218 MB  Time elapsed: 288 s

 ***** [DMNRG] Iteration 40/67 (scale 1.3758612069024691e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=58 min=1 max=371 sum3=231767230
Emax=10.0717107269 nrkept=1115 nrkeptmult=1115 
Kept: 1115 out of 4252, ratio=0.262
Memory used: 218 MB  Time elapsed: 288 s

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

** Iteration completed.

Spectrum: A_d-A_d-d DMNRG
[A_d-A_d-d] pos=0.135301936018 neg=0.864813926951 sum= 1.00011586297
mu1=-0.00379446937634 mu2=0.000662715844601 mu3=-5.05476245327e-06 mu4=0.000226374791343

Spectrum: A_d-A_d-d CFS
[A_d-A_d-d] pos=0.135351257048 neg=0.864648742951 sum= 0.999999999998
mu1=-0.00381945400461 mu2=0.000665488016349 mu3=-5.1271518497e-06 mu4=0.000226973348745

Spectrum: A_d-A_d-d FDM
[A_d-A_d-d] pos=0.135351257048 neg=0.86464874295 sum= 0.999999999998
mu1=-0.00381945400461 mu2=0.000665488016349 mu3=-5.1271518497e-06 mu4=0.000226973348745

Spectrum: A_d-A_d-u DMNRG
[A_d-A_d-u] pos=0.864813926951 neg=0.135301936018 sum= 1.00011586297
mu1=0.00379446937634 mu2=0.000662715844602 mu3=5.05476245328e-06 mu4=0.000226374791343

Spectrum: A_d-A_d-u CFS
[A_d-A_d-u] pos=0.864648742951 neg=0.135351257048 sum= 0.999999999998
mu1=0.00381945400462 mu2=0.000665488016349 mu3=5.1271518497e-06 mu4=0.000226973348745

Spectrum: A_d-A_d-u FDM
[A_d-A_d-u] pos=0.86464874295 neg=0.135351257048 sum= 0.999999999998
mu1=0.00381945400462 mu2=0.000665488016349 mu3=5.1271518497e-06 mu4=0.000226973348745

Master done. Terminating slave processes.
Master exiting.

Memory usage report (process 0)
===================
  after diag: 430308 kB
after recalc: 532232 kB
 after split: 488400 kB
  after trim: 263196 kB
 after trunc: 430788 kB
      diag_h: 399260 kB

Peak usage: 532232 kB

Timing report (process 0)
=============
         All: 296 s
          DM: 4 s
         FDM: 19 s
        diag: 175 s
    recalc d: 25 s
      spec d: 44 s
       Other: 26 s
