NRG Ljubljana 2.3.20 - (c) Rok Zitko <rok.zitko@ijs.si>
$Id: nrg.cc,v 1.31 2012/10/26 10:18:22 rokzitko Exp rokzitko $
Timestamp: Fri May 31 11:54:49 2013
Compiled on Jun  5 2013 at 18:59:39

No MPI: single node calculation.

workdir=.
Tmin=1e-10 -> Nmax=34
TEMPERATURE T=1e-09
Nmax=34
Tmin=1e-10
Ninit=0
Lambda=4
strategy=kept
lastall=0
lastalloverride=1
z=0.25
discretization=Z
band=asymode
betabar=1
T=1e-09
T_ratio=1
Tmin_ratio=-999
keep=8000
keepenergy=10
keepmin=500
keepblock=0
keepzrescale=0
log=
logall=0
logenumber=10
smooth=old
alpha=0.3
omega0=-1
omega0_ratio=1
brcut=0
loggauss_b=0.3
gauss_b=0.1
limitLL=0
eta=0
finite=0
bins=1000
discard_trim=1e-16
discard_immediately=1e-16
NN1=0
NN2even=1
NN2avg=0
NNtanh=0
goodE=2
savebins=1
broaden=0
broaden_max=4
broaden_min=1e-08
broaden_min_ratio=3
broaden_ratio=1.01
dumpenergies=0
saveall=0
dumpannotated=0
dumpscaled=1
dumpabs=0
dumpprecision=8
dumpgroups=1
grouptol=1e-06
dumpdiagonal=0
ops=A_d A_a self_d n_d n_d_ud Himp Hhyb Hpot n_a
specs=
specd=A_d-A_d self_d-A_d A_a-A_a
spect=
specb=
reim=0
trsq=
specgt=
gtp=0.7
speci1t=
speci2t=
speckubo=
kubop=0.7
kubor=4
kubob=2
cond=
specchit=
chitp=1
chitp_ratio=-999
dm=1
T0opt=-1
dmnrg=0
cfs=0
fdm=1
fdmexpv=1
fdmexpvn=0
mats=0
nromegan=100
width_td=12
width_custom=20
prec_td=6
prec_custom=14
prec_xy=12
safeguard=1e-05
safeguardmax=200
fixeps=1e-10
diag=dsyev
diagratio=1
dsyevrlimit=100
restart=1
restartfactor=2
checkdiag=0
checkrho=0
calc0=1
tdht=0
spin=2
tri=cpp
preccpp=2000
checksumrules=0
showEs=0
globalB=0
globalBx=0
polarized=0
pol2x2=0
mpi=0
saveF=0
trunc=
noimag=1
dumpsubspaces=0
done=1
resume=0
forcestop=-1
stopafter=
removefiles=1
store=one

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QS
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	../model.m	 Variant: 		 Channels: 	1
>>  U=	0	 Gamma=	-999	 Delta=	0.026090473124786528	 t=	0.05
>> 	 extraeps=	-0.15	 extraUf=	0.3
>>  Gamma^(1/2)=	0.05039694022504849
>>  Lambda=	4.	 BAND=	asymode
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d A_a self_d n_d n_d_ud Himp Hhyb Hpot n_a
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=2 subspaces=10
Reading <||||> (e)
Reading <||Himp||> (s)
Reading <||Hhyb||> (s)
Reading <||Hpot||> (s)
Reading <||n_d||> (s)
Reading <||n_d_ud||> (s)
Reading <||A_d||> (d)
Reading <||self_d||> (d)
Reading <||n_a||> (s)
Reading <||A_a||> (d)
Reading <||||> (z)
Reading <||||> (T)
Tridiagonalisation, ch=0. Using GMP version 5.0.5
mMAX=80
Using precision of 2000 digits.
  xi(0)=            0.1144305253662504 -->            0.07478096956776741
zeta(0)=           0.02593351375014982
  xi(1)=            0.6968772214949026 -->             0.9108260951574836
zeta(1)=          -0.03016320307855256
  xi(2)=            0.6327693341854034 -->              1.654072780726432
zeta(2)=          0.008770969097930048
  xi(3)=            0.2489581344918469 -->               1.30156393983131
zeta(3)=           0.01673998779235401
  xi(4)=            0.2132595788775466 -->              2.229860681252138
zeta(4)=          -0.02914734047892787
  xi(5)=           0.06828681255979968 -->              1.428025687536083
zeta(5)=           0.02399757856445574
  xi(6)=           0.03089896454985272 -->              1.292329029322201
zeta(6)=          -0.01068137250665488
  xi(7)=            0.0110134431065176 -->             0.9212601423182004
zeta(7)=         -0.002988320073562747
  xi(8)=          0.004453076685064295 -->              0.744988106073463
zeta(8)=         -0.004265640040074794
  xi(9)=          0.003634123832871339 -->              1.215958863034115
zeta(9)=          0.001464833932032978
  xi(10)=           0.00140545454559611 -->             0.9405155079478427
zeta(10)=          0.000984526360301138
  xi(11)=         0.0008516751217631156 -->              1.139864198755479
zeta(11)=         0.0003435472093770731
  xi(12)=         0.0003877558346785009 -->              1.037928624458531
zeta(12)=         0.0002167551274855588
  xi(13)=         0.0001932630827652651 -->              1.034637097437795
zeta(13)=         5.156187371365599e-05
  xi(14)=         9.495632386217794e-05 -->              1.016700488250563
zeta(14)=         1.949353390157129e-05
  xi(15)=         4.692148238509748e-05 -->              1.004779715769309
zeta(15)=         5.007939810954625e-06
  xi(16)=         2.336003193250208e-05 -->              1.000466526307198
zeta(16)=         1.289928142932856e-06
  xi(17)=         1.166493336921574e-05 -->             0.9991746073999636
zeta(17)=         3.401021049833176e-07
  xi(18)=         5.830222263286773e-06 -->             0.9987901099114028
zeta(18)=         9.191570050979866e-08
  xi(19)=         2.914844114425201e-06 -->               0.99869862312251
zeta(19)=         2.422174883967601e-08
  xi(20)=         1.457338979921562e-06 -->             0.9986416943997959
zeta(20)=         6.633585504852539e-09
  xi(21)=         7.287106693225844e-07 -->             0.9986981307240993
zeta(21)=         1.637265467619852e-09
  xi(22)=         3.643168019875192e-07 -->              0.998592512648534
zeta(22)=         4.778007196524093e-10
  xi(23)=          1.82203658695768e-07 -->              0.998840615410299
zeta(23)=         8.466029765198793e-11
  xi(24)=         9.107096603602481e-08 -->             0.9985022300053988
zeta(24)=         3.994769882475309e-11
  xi(25)=         4.556017673974459e-08 -->             0.9990437140214347
zeta(25)=        -7.972177935061105e-12
  xi(26)=         2.282088912064963e-08 -->              1.000833072031361
zeta(26)=         7.416818421601249e-12
  xi(27)=         1.144310347902419e-08 -->              1.003697651562711
zeta(27)=        -3.171625254639905e-12
  xi(28)=         5.713943895678882e-09 -->              1.002363052954414
zeta(28)=        -6.692577148040026e-12
  xi(29)=         2.853615186698242e-09 -->              1.001185339834733
zeta(29)=         1.825967074873624e-12
  xi(30)=         1.427053003085628e-09 -->              1.001357542892517
zeta(30)=         1.175610211991365e-12
  xi(31)=         7.135253473237339e-10 -->              1.001355923070437
zeta(31)=         8.062202243270936e-15
  xi(32)=          3.56785559225573e-10 -->              1.001420157914667
zeta(32)=         7.480784218830556e-14
  xi(33)=         1.783927041638332e-10 -->               1.00141973437677
zeta(33)=        -5.423227387090856e-14
  xi(34)=         8.919780348115193e-11 -->               1.00143602943596
zeta(34)=         4.676774521902096e-15

# $Id: nrg.cc,v 1.31 2012/10/26 10:18:22 rokzitko Exp rokzitko $ QS
# disk=Z band=asymode Lambda=4 z=0.25
# keep=8000 keepenergy=10 Nmax=34 channels=1 dots=2 betabar=1

Number of invariant subspaces: 10
(-3 1) 1 states
0.0696194483634 
(-2 2) 3 states
0.015987812574 0.0644139790302 0.0964426252426 
(-1 1) 6 states
0.00127081030268 0.0345081403727 0.0584025549611 0.077040431865 0.0930117577717 0.12345309628 
(-1 3) 3 states
0.0107823432408 0.0428109894532 0.0912371559094 
(0 2) 8 states
0 0.0249244319275 0.0320938274897 0.0640974610585 0.0697081987579 0.0854003759797 0.101570962928 0.119100531513 
(0 4) 1 states
0.03760552012 
(1 1) 6 states
0.0167620859594 0.0506186030213 0.0614401567319 0.0898792857832 0.111408618135 0.125573957642 
(1 3) 3 states
0.0272408044948 0.0584904671109 0.0930971748575 
(2 2) 3 states
0.0481257514858 0.0827324592324 0.113982121848 
(3 1) 1 states
0.103617406223 
Number of states (multiplicity taken into account): 64

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

Computing the following spectra:

Before NRG iteration (N=-1)
<Hhyb>=-0.0526145769711
<Himp>=0.025979273892
<Hpot>=0
<n_a>=1.00000009247
<n_d>=0.995737937282
<n_d_ud>=0.24787350993

 ***** [NRG] Iteration 1/34 (scale 1.5302091698951845) ***** 
[1] xi(0)=0.0747809695678 zeta(1)=-0.0301632030786
Stats: nr=15 min=1 max=20 sum3=30066
Emax=0.497184529091 nrkept=126 nrkeptmult=256 
Total energy=-0.426768632707  Egs=-0.122709029637
Kept: 126 out of 126, ratio=1
<Hhyb>=-0.0563139164999
<Himp>=0.0258681534239
<Hpot>=0
<n_a>=1.0000007206
<n_d>=0.991478893471
<n_d_ud>=0.245757598876
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/34 (scale 0.76510458494759226) ***** 
[1] xi(1)=0.910826095157 zeta(2)=0.00877096909793
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=4.19278962132 nrkept=462 nrkeptmult=1024 
Total energy=-1.65080740555  Egs=-1.59983196667
Kept: 462 out of 462, ratio=1
<Hhyb>=-0.0636750265678
<Himp>=0.0256464614036
<Hpot>=0
<n_a>=1.00000545195
<n_d>=0.982981844787
<n_d_ud>=0.241563321252
Memory used: 7 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/34 (scale 0.38255229247379613) ***** 
[1] xi(2)=1.65407278073 zeta(3)=0.0167399877924
Stats: nr=28 min=1 max=210 sum3=37226316
Emax=10.0915024343 nrkept=1595 nrkeptmult=3850 
Total energy=-2.24289538978  Egs=-1.54773084849
Kept: 1595 out of 1716, ratio=0.929
<Hhyb>=-0.0781064620371
<Himp>=0.0252073894645
<Hpot>=0
<n_a>=1.0000385136
<n_d>=0.966153022369
<n_d_ud>=0.233362703411
Memory used: 30 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/34 (scale 0.19127614623689806) ***** 
[1] xi(3)=1.30156393983 zeta(4)=-0.0291473404789
Stats: nr=36 min=1 max=714 sum3=1151440282
Emax=10.0030831045 nrkept=1961 nrkeptmult=4935 
Total energy=-2.55133589096  Egs=-1.61254033633
Kept: 1961 out of 5992, ratio=0.327
<Hhyb>=-0.104886875565
<Himp>=0.0243615696617
<Hpot>=0
<n_a>=1.00023175068
<n_d>=0.933734300072
<n_d_ud>=0.217959090639
Memory used: 48 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 5/34 (scale 0.095638073118449032) ***** 
[1] xi(4)=2.22986068125 zeta(5)=0.0239975785645
Stats: nr=42 min=1 max=746 sum3=1939801406
Emax=10.0118719075 nrkept=1093 nrkeptmult=2656 
Total energy=-2.80486398141  Egs=-2.65091173613
Kept: 1093 out of 7406, ratio=0.148
<Hhyb>=-0.146195864863
<Himp>=0.0228808389141
<Hpot>=0
<n_a>=1.00080535907
<n_d>=0.876980605323
<n_d_ud>=0.192153919504
Memory used: 25 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 6/34 (scale 0.047819036559224516) ***** 
[1] xi(5)=1.42802568754 zeta(6)=-0.0106813725067
Stats: nr=36 min=1 max=491 sum3=372016035
Emax=10.0022079988 nrkept=676 nrkeptmult=1635 
Total energy=-2.86614072451  Egs=-1.28142989706
Kept: 676 out of 4107, ratio=0.165
<Hhyb>=-0.186895385541
<Himp>=0.0208902363408
<Hpot>=0
<n_a>=1.0002997895
<n_d>=0.800684458305
<n_d_ud>=0.158788377696
Memory used: 25 MB  Time elapsed: 9 s

 ***** [NRG] Iteration 7/34 (scale 0.023909518279612258) ***** 
[1] xi(6)=1.29232902932 zeta(7)=-0.00298832007356
Stats: nr=37 min=1 max=293 sum3=99347966
Emax=11.7322871381 nrkept=500 nrkeptmult=1251 
Total energy=-2.90350944909  Egs=-1.56292252096
Kept: 500 out of 2540, ratio=0.197
<Hhyb>=-0.209388977972
<Himp>=0.0193201547205
<Hpot>=0
<n_a>=0.996935754021
<n_d>=0.740506108422
<n_d_ud>=0.129401312837
Memory used: 26 MB  Time elapsed: 9 s

 ***** [NRG] Iteration 8/34 (scale 0.011954759139806129) ***** 
[1] xi(7)=0.921260142318 zeta(8)=-0.00426564004007
Stats: nr=33 min=1 max=235 sum3=39962035
Emax=14.2796160777 nrkept=500 nrkeptmult=1266 
Total energy=-2.91259285204  Egs=-0.759814801807
Kept: 500 out of 1885, ratio=0.265
<Hhyb>=-0.21893523384
<Himp>=0.0189039429654
<Hpot>=0
<n_a>=0.995877204504
<n_d>=0.724553474941
<n_d_ud>=0.113126887499
Memory used: 26 MB  Time elapsed: 9 s

 ***** [NRG] Iteration 9/34 (scale 0.0059773795699030645) ***** 
[1] xi(8)=0.744988106073 zeta(9)=0.00146483393203
Stats: nr=34 min=1 max=216 sum3=39459782
Emax=17.7279797771 nrkept=500 nrkeptmult=1238 
Total energy=-2.91889760126  Egs=-1.05476808766
Kept: 500 out of 1892, ratio=0.264
<Hhyb>=-0.221791579386
<Himp>=0.0192804698509
<Hpot>=0
<n_a>=0.998103365557
<n_d>=0.738985060129
<n_d_ud>=0.108545457387
Memory used: 26 MB  Time elapsed: 9 s

 ***** [NRG] Iteration 10/34 (scale 0.0029886897849515323) ***** 
[1] xi(9)=1.21595886303 zeta(10)=0.000984526360301
Stats: nr=32 min=1 max=236 sum3=41505866
Emax=21.1020306199 nrkept=500 nrkeptmult=1205 
Total energy=-2.92231263094  Egs=-1.14265110269
Kept: 500 out of 1886, ratio=0.265
<Hhyb>=-0.222497764905
<Himp>=0.0199157492179
<Hpot>=0
<n_a>=1.00162771528
<n_d>=0.76333415353
<n_d_ud>=0.110560039005
Memory used: 27 MB  Time elapsed: 10 s

 ***** [NRG] Iteration 11/34 (scale 0.0014943448924757661) ***** 
[1] xi(10)=0.940515507948 zeta(11)=0.000343547209377
Stats: nr=33 min=1 max=225 sum3=42319021
Emax=23.5001223623 nrkept=500 nrkeptmult=1221 
Total energy=-2.92415572858  Egs=-1.23338169698
Kept: 500 out of 1879, ratio=0.266
<Hhyb>=-0.22263429559
<Himp>=0.0203592330715
<Hpot>=0
<n_a>=1.00414382589
<n_d>=0.780332076544
<n_d_ud>=0.113420421416
Memory used: 26 MB  Time elapsed: 10 s

 ***** [NRG] Iteration 12/34 (scale 0.00074717244623788307) ***** 
[1] xi(11)=1.13986419876 zeta(12)=0.000216755127486
Stats: nr=32 min=1 max=238 sum3=40834017
Emax=25.1041667057 nrkept=500 nrkeptmult=1175 
Total energy=-2.9250002615  Egs=-1.13030523052
Kept: 500 out of 1881, ratio=0.266
<Hhyb>=-0.222655864696
<Himp>=0.0206098572712
<Hpot>=0
<n_a>=1.00560545475
<n_d>=0.789938042619
<n_d_ud>=0.11507052228
Memory used: 27 MB  Time elapsed: 10 s

 ***** [NRG] Iteration 13/34 (scale 0.00037358622311894153) ***** 
[1] xi(12)=1.03792862446 zeta(13)=5.15618737137e-05
Stats: nr=33 min=1 max=230 sum3=42551582
Emax=26.6338465679 nrkept=500 nrkeptmult=1176 
Total energy=-2.92545363263  Egs=-1.2135649045
Kept: 500 out of 1874, ratio=0.267
<Hhyb>=-0.222652804216
<Himp>=0.0207237481232
<Hpot>=0
<n_a>=1.0062857987
<n_d>=0.794303270166
<n_d_ud>=0.115801855094
Memory used: 27 MB  Time elapsed: 10 s

 ***** [NRG] Iteration 14/34 (scale 0.00018679311155947077) ***** 
[1] xi(13)=1.03463709744 zeta(14)=1.94935339016e-05
Stats: nr=32 min=1 max=245 sum3=43309706
Emax=27.4181001159 nrkept=500 nrkeptmult=1167 
Total energy=-2.92565436449  Egs=-1.07462132673
Kept: 500 out of 1874, ratio=0.267
<Hhyb>=-0.222650312787
<Himp>=0.020777492505
<Hpot>=0
<n_a>=1.00660603881
<n_d>=0.796363193784
<n_d_ud>=0.116120016497
Memory used: 28 MB  Time elapsed: 11 s

 ***** [NRG] Iteration 15/34 (scale 9.3396555779735383e-05) ***** 
[1] xi(14)=1.01670048825 zeta(15)=5.00793981095e-06
Stats: nr=33 min=1 max=232 sum3=42818201
Emax=27.9169622785 nrkept=500 nrkeptmult=1180 
Total energy=-2.92575590919  Egs=-1.08724250441
Kept: 500 out of 1871, ratio=0.267
<Hhyb>=-0.222649451435
<Himp>=0.020797689605
<Hpot>=0
<n_a>=1.00672557005
<n_d>=0.797137311598
<n_d_ud>=0.116230058213
Memory used: 28 MB  Time elapsed: 11 s

 ***** [NRG] Iteration 16/34 (scale 4.6698277889867692e-05) ***** 
[1] xi(15)=1.00477971577 zeta(16)=1.28992814293e-06
Stats: nr=32 min=1 max=242 sum3=43011074
Emax=30.1529136196 nrkept=500 nrkeptmult=1158 
Total energy=-2.92580596637  Egs=-1.0719276041
Kept: 500 out of 1874, ratio=0.267
<Hhyb>=-0.222649382368
<Himp>=0.0208040601707
<Hpot>=0
<n_a>=1.00676277878
<n_d>=0.7973814837
<n_d_ud>=0.116262035302
Memory used: 28 MB  Time elapsed: 11 s

 ***** [NRG] Iteration 17/34 (scale 2.3349138944933846e-05) ***** 
[1] xi(16)=1.00046652631 zeta(17)=3.40102104983e-07
Stats: nr=33 min=1 max=233 sum3=42404828
Emax=30.3361142298 nrkept=500 nrkeptmult=1166 
Total energy=-2.92583023516  Egs=-1.03938695577
Kept: 500 out of 1868, ratio=0.268
<Hhyb>=-0.222649432268
<Himp>=0.0208056451366
<Hpot>=0
<n_a>=1.00677189337
<n_d>=0.797442232539
<n_d_ud>=0.116269290421
Memory used: 28 MB  Time elapsed: 11 s

 ***** [NRG] Iteration 18/34 (scale 1.1674569472466923e-05) ***** 
[1] xi(17)=0.9991746074 zeta(18)=9.19157005098e-08
Stats: nr=33 min=1 max=238 sum3=42097666
Emax=32.0064905533 nrkept=500 nrkeptmult=1150 
Total energy=-2.92584282211  Egs=-1.07815159562
Kept: 500 out of 1870, ratio=0.267
<Hhyb>=-0.222649419171
<Himp>=0.0208061113894
<Hpot>=0
<n_a>=1.00677462951
<n_d>=0.797460103154
<n_d_ud>=0.116271668299
Memory used: 28 MB  Time elapsed: 12 s

 ***** [NRG] Iteration 19/34 (scale 5.8372847362334614e-06) ***** 
[1] xi(18)=0.998790109911 zeta(19)=2.42217488397e-08
Stats: nr=33 min=1 max=238 sum3=42169664
Emax=32.4072773313 nrkept=500 nrkeptmult=1148 
Total energy=-2.92584882313  Egs=-1.02804957242
Kept: 500 out of 1868, ratio=0.268
<Hhyb>=-0.222649431082
<Himp>=0.0208061942487
<Hpot>=0
<n_a>=1.00677508731
<n_d>=0.797463278999
<n_d_ud>=0.116271964258
Memory used: 29 MB  Time elapsed: 12 s

 ***** [NRG] Iteration 20/34 (scale 2.9186423681167307e-06) ***** 
[1] xi(19)=0.998698623123 zeta(20)=6.63358550485e-09
Stats: nr=33 min=1 max=236 sum3=42054289
Emax=32.3053167118 nrkept=500 nrkeptmult=1144 
Total energy=-2.92585197864  Egs=-1.08115509384
Kept: 500 out of 1867, ratio=0.268
<Hhyb>=-0.222649426783
<Himp>=0.0208062314069
<Hpot>=0
<n_a>=1.00677531184
<n_d>=0.797464703204
<n_d_ud>=0.116272182342
Memory used: 29 MB  Time elapsed: 12 s

 ***** [NRG] Iteration 21/34 (scale 1.4593211840583654e-06) ***** 
[1] xi(20)=0.9986416944 zeta(21)=1.63726546762e-09
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.7346637203 nrkept=500 nrkeptmult=1140 
Total energy=-2.9258534746  Egs=-1.02510577403
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.222649429281
<Himp>=0.02080623266
<Hpot>=0
<n_a>=1.00677531414
<n_d>=0.797464751234
<n_d_ud>=0.116272166309
Memory used: 29 MB  Time elapsed: 12 s

 ***** [NRG] Iteration 22/34 (scale 7.2966059202918268e-07) ***** 
[1] xi(21)=0.998698130724 zeta(22)=4.77800719652e-10
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.3964664823 nrkept=500 nrkeptmult=1144 
Total energy=-2.92585426422  Egs=-1.08218513617
Kept: 500 out of 1865, ratio=0.268
<Hhyb>=-0.222649428155
<Himp>=0.0208062369077
<Hpot>=0
<n_a>=1.00677534107
<n_d>=0.79746491404
<n_d_ud>=0.116272196851
Memory used: 29 MB  Time elapsed: 12 s

 ***** [NRG] Iteration 23/34 (scale 3.6483029601459134e-07) ***** 
[1] xi(22)=0.998592512649 zeta(23)=8.4660297652e-11
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.8349663324 nrkept=500 nrkeptmult=1140 
Total energy=-2.9258546379  Egs=-1.02424211061
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.222649428752
<Himp>=0.0208062359963
<Hpot>=0
<n_a>=1.00677533458
<n_d>=0.797464879106
<n_d_ud>=0.116272187152
Memory used: 29 MB  Time elapsed: 13 s

 ***** [NRG] Iteration 24/34 (scale 1.8241514800729567e-07) ***** 
[1] xi(23)=0.99884061541 zeta(24)=3.99476988248e-11
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.4218301218 nrkept=500 nrkeptmult=1144 
Total energy=-2.92585483541  Egs=-1.08275301994
Kept: 500 out of 1865, ratio=0.268
<Hhyb>=-0.222649428475
<Himp>=0.0208062367259
<Hpot>=0
<n_a>=1.00677533938
<n_d>=0.797464907072
<n_d_ud>=0.11627219314
Memory used: 29 MB  Time elapsed: 13 s

 ***** [NRG] Iteration 25/34 (scale 9.1207574003647835e-08) ***** 
[1] xi(24)=0.998502230005 zeta(25)=-7.97217793506e-12
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.8584853525 nrkept=500 nrkeptmult=1140 
Total energy=-2.92585492878  Egs=-1.02374344337
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.222649428629
<Himp>=0.0208062364066
<Hpot>=0
<n_a>=1.00677533721
<n_d>=0.797464894831
<n_d_ud>=0.116272190228
Memory used: 29 MB  Time elapsed: 13 s

 ***** [NRG] Iteration 26/34 (scale 4.5603787001823918e-08) ***** 
[1] xi(25)=0.999043714021 zeta(26)=7.4168184216e-12
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.4298043137 nrkept=500 nrkeptmult=1144 
Total energy=-2.92585497819  Egs=-1.08338809367
Kept: 500 out of 1865, ratio=0.268
<Hhyb>=-0.222649428566
<Himp>=0.0208062365551
<Hpot>=0
<n_a>=1.0067753382
<n_d>=0.797464900524
<n_d_ud>=0.116272191506
Memory used: 29 MB  Time elapsed: 13 s

 ***** [NRG] Iteration 27/34 (scale 2.2801893500911959e-08) ***** 
[1] xi(26)=1.00083307203 zeta(27)=-3.17162525464e-12
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.8673407534 nrkept=500 nrkeptmult=1140 
Total energy=-2.9258550016  Egs=-1.0268861456
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.222649428607
<Himp>=0.0208062364648
<Hpot>=0
<n_a>=1.0067753376
<n_d>=0.797464897062
<n_d_ud>=0.116272190706
Memory used: 29 MB  Time elapsed: 14 s

 ***** [NRG] Iteration 28/34 (scale 1.1400946750455979e-08) ***** 
[1] xi(27)=1.00369765156 zeta(28)=-6.69257714804e-12
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.4461867003 nrkept=500 nrkeptmult=1144 
Total energy=-2.92585501402  Egs=-1.0887822617
Kept: 500 out of 1865, ratio=0.268
<Hhyb>=-0.222649428589
<Himp>=0.0208062365068
<Hpot>=0
<n_a>=1.00677533788
<n_d>=0.797464898673
<n_d_ud>=0.116272191074
Memory used: 30 MB  Time elapsed: 14 s

 ***** [NRG] Iteration 29/34 (scale 5.7004733752279897e-09) ***** 
[1] xi(28)=1.00236305295 zeta(29)=1.82596707487e-12
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.8936122004 nrkept=500 nrkeptmult=1140 
Total energy=-2.92585501987  Egs=-1.02661257719
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.222649428598
<Himp>=0.0208062364859
<Hpot>=0
<n_a>=1.00677533774
<n_d>=0.797464897873
<n_d_ud>=0.11627219089
Memory used: 30 MB  Time elapsed: 14 s

 ***** [NRG] Iteration 30/34 (scale 2.8502366876139948e-09) ***** 
[1] xi(29)=1.00118533983 zeta(30)=1.17561021199e-12
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.4749320347 nrkept=500 nrkeptmult=1144 
Total energy=-2.92585502296  Egs=-1.0860559072
Kept: 500 out of 1865, ratio=0.268
<Hhyb>=-0.222649428597
<Himp>=0.0208062364876
<Hpot>=0
<n_a>=1.00677533775
<n_d>=0.797464897936
<n_d_ud>=0.116272190904
Memory used: 30 MB  Time elapsed: 14 s

 ***** [NRG] Iteration 31/34 (scale 1.4251183438069974e-09) ***** 
[1] xi(30)=1.00135754289 zeta(31)=8.06220224327e-15
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.9403661865 nrkept=500 nrkeptmult=1140 
Total energy=-2.92585502443  Egs=-1.02601871664
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.2226494286
<Himp>=0.0208062364822
<Hpot>=0
<n_a>=1.00677533771
<n_d>=0.797464897731
<n_d_ud>=0.116272190857
Memory used: 31 MB  Time elapsed: 14 s

 ***** [NRG] Iteration 32/34 (scale 7.1255917190349871e-10) ***** 
[1] xi(31)=1.00135592307 zeta(32)=7.48078421883e-14
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.5437372723 nrkept=500 nrkeptmult=1144 
Total energy=-2.9258550252  Egs=-1.08632554748
Kept: 500 out of 1865, ratio=0.268
<Hhyb>=-0.222649428599
<Himp>=0.0208062364843
<Hpot>=0
<n_a>=1.00677533773
<n_d>=0.79746489781
<n_d_ud>=0.116272190876
Memory used: 31 MB  Time elapsed: 15 s

 ***** [NRG] Iteration 33/34 (scale 3.5627958595174936e-10) ***** 
[1] xi(32)=1.00142015791 zeta(33)=-5.42322738709e-14
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.9788693712 nrkept=500 nrkeptmult=1140 
Total energy=-2.92585502557  Egs=-1.02610387253
Kept: 500 out of 1866, ratio=0.268
<Hhyb>=-0.222649428599
<Himp>=0.0208062364833
<Hpot>=0
<n_a>=1.00677533772
<n_d>=0.797464897772
<n_d_ud>=0.116272190867
Memory used: 31 MB  Time elapsed: 15 s

 ***** [NRG] Iteration 34/34 (scale 1.7813979297587468e-10) ***** 
[1] xi(33)=1.00141973438 zeta(34)=4.6767745219e-15
Stats: nr=32 min=1 max=235 sum3=42515861
Safeguard: keep additional 3 states
Emax=32.501614643 nrkept=503 nrkeptmult=1152 
Total energy=-2.92585502576  Egs=-1.08630491541
Kept: 503 out of 1865, ratio=0.27
<Hhyb>=-0.222649428599
<Himp>=0.0208062364838
<Hpot>=0
<n_a>=1.00677533772
<n_d>=0.797464897792
<n_d_ud>=0.116272190871
Memory used: 32 MB  Time elapsed: 15 s

** Iteration completed.


Total energy:             -2.925855025759782
TKW=0.00102089882916
Grand canonical Z=              75.27182113152152
Zft (last shell)=               73.97850067891696
Calculating density matrices.
Calculating full density matrices.

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QS
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	../model.m	 Variant: 		 Channels: 	1
>>  U=	0	 Gamma=	-999	 Delta=	0.026090473124786528	 t=	0.05
>> 	 extraeps=	-0.15	 extraUf=	0.3
>>  Gamma^(1/2)=	0.05039694022504849
>>  Lambda=	4.	 BAND=	asymode
>>  DISCRETIZATION=	Z	 z=	0.25
>>  ops=	A_d A_a self_d n_d n_d_ud Himp Hhyb Hpot n_a
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=2 subspaces=10
Reading <||||> (e)
Reading <||Himp||> (s)
Reading <||Hhyb||> (s)
Reading <||Hpot||> (s)
Reading <||n_d||> (s)
Reading <||n_d_ud||> (s)
Reading <||A_d||> (d)
Reading <||self_d||> (d)
Reading <||n_a||> (s)
Reading <||A_a||> (d)
Reading <||||> (z)
Reading <||||> (T)
Number of invariant subspaces: 10
(-3 1) 1 states
0.0696194483634 
(-2 2) 3 states
0.015987812574 0.0644139790302 0.0964426252426 
(-1 1) 6 states
0.00127081030268 0.0345081403727 0.0584025549611 0.077040431865 0.0930117577717 0.12345309628 
(-1 3) 3 states
0.0107823432408 0.0428109894532 0.0912371559094 
(0 2) 8 states
0 0.0249244319275 0.0320938274897 0.0640974610585 0.0697081987579 0.0854003759797 0.101570962928 0.119100531513 
(0 4) 1 states
0.03760552012 
(1 1) 6 states
0.0167620859594 0.0506186030213 0.0614401567319 0.0898792857832 0.111408618135 0.125573957642 
(1 3) 3 states
0.0272408044948 0.0584904671109 0.0930971748575 
(2 2) 3 states
0.0481257514858 0.0827324592324 0.113982121848 
(3 1) 1 states
0.103617406223 
Number of states (multiplicity taken into account): 64

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

Computing the following spectra:
Spectrum: A_a-A_a FDM
Spectrum: A_d-A_d FDM
Spectrum: self_d-A_d FDM

 ***** [DMNRG] Iteration 1/34 (scale 1.5302091698951845) ***** 
[1] xi(0)=0.0747809695678 zeta(1)=-0.0301632030786
Stats: nr=15 min=1 max=20 sum3=30066
Emax=0.497184529091 nrkept=126 nrkeptmult=256 
<Hhyb>_fdm=-0.222649428592
<Himp>_fdm=0.0208062364995
<Hpot>_fdm=0
<n_a>_fdm=1.00677533783
<n_d>_fdm=0.797464898393
<n_d_ud>_fdm=0.116272191009
Kept: 126 out of 126, ratio=1
Memory used: 32 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 2/34 (scale 0.76510458494759226) ***** 
[1] xi(1)=0.910826095157 zeta(2)=0.00877096909793
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=4.19278962132 nrkept=462 nrkeptmult=1024 
Kept: 462 out of 462, ratio=1
Memory used: 33 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 3/34 (scale 0.38255229247379613) ***** 
[1] xi(2)=1.65407278073 zeta(3)=0.0167399877924
Stats: nr=28 min=1 max=210 sum3=37226316
Emax=10.0915024343 nrkept=1595 nrkeptmult=3850 
Kept: 1595 out of 1716, ratio=0.929
Memory used: 47 MB  Time elapsed: 17 s

 ***** [DMNRG] Iteration 4/34 (scale 0.19127614623689806) ***** 
[1] xi(3)=1.30156393983 zeta(4)=-0.0291473404789
Stats: nr=36 min=1 max=714 sum3=1151440282
Emax=10.0030831045 nrkept=1961 nrkeptmult=4935 
Kept: 1961 out of 5992, ratio=0.327
Memory used: 56 MB  Time elapsed: 20 s

 ***** [DMNRG] Iteration 5/34 (scale 0.095638073118449032) ***** 
[1] xi(4)=2.22986068125 zeta(5)=0.0239975785645
Stats: nr=42 min=1 max=746 sum3=1939801406
Emax=10.0118719075 nrkept=1093 nrkeptmult=2656 
Kept: 1093 out of 7406, ratio=0.148
Memory used: 42 MB  Time elapsed: 24 s

 ***** [DMNRG] Iteration 6/34 (scale 0.047819036559224516) ***** 
[1] xi(5)=1.42802568754 zeta(6)=-0.0106813725067
Stats: nr=36 min=1 max=491 sum3=372016035
Emax=10.0022079988 nrkept=676 nrkeptmult=1635 
Kept: 676 out of 4107, ratio=0.165
Memory used: 51 MB  Time elapsed: 26 s

 ***** [DMNRG] Iteration 7/34 (scale 0.023909518279612258) ***** 
[1] xi(6)=1.29232902932 zeta(7)=-0.00298832007356
Stats: nr=37 min=1 max=293 sum3=99347966
Emax=11.7322871381 nrkept=500 nrkeptmult=1251 
Kept: 500 out of 2540, ratio=0.197
Memory used: 49 MB  Time elapsed: 26 s

 ***** [DMNRG] Iteration 8/34 (scale 0.011954759139806129) ***** 
[1] xi(7)=0.921260142318 zeta(8)=-0.00426564004007
Stats: nr=33 min=1 max=235 sum3=39962035
Emax=14.2796160777 nrkept=500 nrkeptmult=1266 
Kept: 500 out of 1885, ratio=0.265
Memory used: 49 MB  Time elapsed: 27 s

 ***** [DMNRG] Iteration 9/34 (scale 0.0059773795699030645) ***** 
[1] xi(8)=0.744988106073 zeta(9)=0.00146483393203
Stats: nr=34 min=1 max=216 sum3=39459782
Emax=17.7279797771 nrkept=500 nrkeptmult=1238 
Kept: 500 out of 1892, ratio=0.264
Memory used: 50 MB  Time elapsed: 27 s

 ***** [DMNRG] Iteration 10/34 (scale 0.0029886897849515323) ***** 
[1] xi(9)=1.21595886303 zeta(10)=0.000984526360301
Stats: nr=32 min=1 max=236 sum3=41505866
Emax=21.1020306199 nrkept=500 nrkeptmult=1205 
Kept: 500 out of 1886, ratio=0.265
Memory used: 52 MB  Time elapsed: 27 s

 ***** [DMNRG] Iteration 11/34 (scale 0.0014943448924757661) ***** 
[1] xi(10)=0.940515507948 zeta(11)=0.000343547209377
Stats: nr=33 min=1 max=225 sum3=42319021
Emax=23.5001223623 nrkept=500 nrkeptmult=1221 
Kept: 500 out of 1879, ratio=0.266
Memory used: 54 MB  Time elapsed: 28 s

 ***** [DMNRG] Iteration 12/34 (scale 0.00074717244623788307) ***** 
[1] xi(11)=1.13986419876 zeta(12)=0.000216755127486
Stats: nr=32 min=1 max=238 sum3=40834017
Emax=25.1041667057 nrkept=500 nrkeptmult=1175 
Kept: 500 out of 1881, ratio=0.266
Memory used: 55 MB  Time elapsed: 28 s

 ***** [DMNRG] Iteration 13/34 (scale 0.00037358622311894153) ***** 
[1] xi(12)=1.03792862446 zeta(13)=5.15618737137e-05
Stats: nr=33 min=1 max=230 sum3=42551582
Emax=26.6338465679 nrkept=500 nrkeptmult=1176 
Kept: 500 out of 1874, ratio=0.267
Memory used: 56 MB  Time elapsed: 29 s

 ***** [DMNRG] Iteration 14/34 (scale 0.00018679311155947077) ***** 
[1] xi(13)=1.03463709744 zeta(14)=1.94935339016e-05
Stats: nr=32 min=1 max=245 sum3=43309706
Emax=27.4181001159 nrkept=500 nrkeptmult=1167 
Kept: 500 out of 1874, ratio=0.267
Memory used: 58 MB  Time elapsed: 29 s

 ***** [DMNRG] Iteration 15/34 (scale 9.3396555779735383e-05) ***** 
[1] xi(14)=1.01670048825 zeta(15)=5.00793981095e-06
Stats: nr=33 min=1 max=232 sum3=42818201
Emax=27.9169622785 nrkept=500 nrkeptmult=1180 
Kept: 500 out of 1871, ratio=0.267
Memory used: 60 MB  Time elapsed: 29 s

 ***** [DMNRG] Iteration 16/34 (scale 4.6698277889867692e-05) ***** 
[1] xi(15)=1.00477971577 zeta(16)=1.28992814293e-06
Stats: nr=32 min=1 max=242 sum3=43011074
Emax=30.1529136196 nrkept=500 nrkeptmult=1158 
Kept: 500 out of 1874, ratio=0.267
Memory used: 61 MB  Time elapsed: 30 s

 ***** [DMNRG] Iteration 17/34 (scale 2.3349138944933846e-05) ***** 
[1] xi(16)=1.00046652631 zeta(17)=3.40102104983e-07
Stats: nr=33 min=1 max=233 sum3=42404828
Emax=30.3361142298 nrkept=500 nrkeptmult=1166 
Kept: 500 out of 1868, ratio=0.268
Memory used: 63 MB  Time elapsed: 30 s

 ***** [DMNRG] Iteration 18/34 (scale 1.1674569472466923e-05) ***** 
[1] xi(17)=0.9991746074 zeta(18)=9.19157005098e-08
Stats: nr=33 min=1 max=238 sum3=42097666
Emax=32.0064905533 nrkept=500 nrkeptmult=1150 
Kept: 500 out of 1870, ratio=0.267
Memory used: 65 MB  Time elapsed: 31 s

 ***** [DMNRG] Iteration 19/34 (scale 5.8372847362334614e-06) ***** 
[1] xi(18)=0.998790109911 zeta(19)=2.42217488397e-08
Stats: nr=33 min=1 max=238 sum3=42169664
Emax=32.4072773313 nrkept=500 nrkeptmult=1148 
Kept: 500 out of 1868, ratio=0.268
Memory used: 67 MB  Time elapsed: 31 s

 ***** [DMNRG] Iteration 20/34 (scale 2.9186423681167307e-06) ***** 
[1] xi(19)=0.998698623123 zeta(20)=6.63358550485e-09
Stats: nr=33 min=1 max=236 sum3=42054289
Emax=32.3053167118 nrkept=500 nrkeptmult=1144 
Kept: 500 out of 1867, ratio=0.268
Memory used: 68 MB  Time elapsed: 31 s

 ***** [DMNRG] Iteration 21/34 (scale 1.4593211840583654e-06) ***** 
[1] xi(20)=0.9986416944 zeta(21)=1.63726546762e-09
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.7346637203 nrkept=500 nrkeptmult=1140 
Kept: 500 out of 1866, ratio=0.268
Memory used: 70 MB  Time elapsed: 32 s

 ***** [DMNRG] Iteration 22/34 (scale 7.2966059202918268e-07) ***** 
[1] xi(21)=0.998698130724 zeta(22)=4.77800719652e-10
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.3964664823 nrkept=500 nrkeptmult=1144 
Kept: 500 out of 1865, ratio=0.268
Memory used: 71 MB  Time elapsed: 32 s

 ***** [DMNRG] Iteration 23/34 (scale 3.6483029601459134e-07) ***** 
[1] xi(22)=0.998592512649 zeta(23)=8.4660297652e-11
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.8349663324 nrkept=500 nrkeptmult=1140 
Kept: 500 out of 1866, ratio=0.268
Memory used: 73 MB  Time elapsed: 32 s

 ***** [DMNRG] Iteration 24/34 (scale 1.8241514800729567e-07) ***** 
[1] xi(23)=0.99884061541 zeta(24)=3.99476988248e-11
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.4218301218 nrkept=500 nrkeptmult=1144 
Kept: 500 out of 1865, ratio=0.268
Memory used: 75 MB  Time elapsed: 33 s

 ***** [DMNRG] Iteration 25/34 (scale 9.1207574003647835e-08) ***** 
[1] xi(24)=0.998502230005 zeta(25)=-7.97217793506e-12
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.8584853525 nrkept=500 nrkeptmult=1140 
Kept: 500 out of 1866, ratio=0.268
Memory used: 77 MB  Time elapsed: 33 s

 ***** [DMNRG] Iteration 26/34 (scale 4.5603787001823918e-08) ***** 
[1] xi(25)=0.999043714021 zeta(26)=7.4168184216e-12
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.4298043137 nrkept=500 nrkeptmult=1144 
Kept: 500 out of 1865, ratio=0.268
Memory used: 78 MB  Time elapsed: 33 s

 ***** [DMNRG] Iteration 27/34 (scale 2.2801893500911959e-08) ***** 
[1] xi(26)=1.00083307203 zeta(27)=-3.17162525464e-12
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.8673407534 nrkept=500 nrkeptmult=1140 
Kept: 500 out of 1866, ratio=0.268
Memory used: 80 MB  Time elapsed: 34 s

 ***** [DMNRG] Iteration 28/34 (scale 1.1400946750455979e-08) ***** 
[1] xi(27)=1.00369765156 zeta(28)=-6.69257714804e-12
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.4461867003 nrkept=500 nrkeptmult=1144 
Kept: 500 out of 1865, ratio=0.268
Memory used: 81 MB  Time elapsed: 34 s

 ***** [DMNRG] Iteration 29/34 (scale 5.7004733752279897e-09) ***** 
[1] xi(28)=1.00236305295 zeta(29)=1.82596707487e-12
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.8936122004 nrkept=500 nrkeptmult=1140 
Kept: 500 out of 1866, ratio=0.268
Memory used: 83 MB  Time elapsed: 35 s

 ***** [DMNRG] Iteration 30/34 (scale 2.8502366876139948e-09) ***** 
[1] xi(29)=1.00118533983 zeta(30)=1.17561021199e-12
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.4749320347 nrkept=500 nrkeptmult=1144 
Kept: 500 out of 1865, ratio=0.268
Memory used: 84 MB  Time elapsed: 35 s

 ***** [DMNRG] Iteration 31/34 (scale 1.4251183438069974e-09) ***** 
[1] xi(30)=1.00135754289 zeta(31)=8.06220224327e-15
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.9403661865 nrkept=500 nrkeptmult=1140 
Kept: 500 out of 1866, ratio=0.268
Memory used: 86 MB  Time elapsed: 35 s

 ***** [DMNRG] Iteration 32/34 (scale 7.1255917190349871e-10) ***** 
[1] xi(31)=1.00135592307 zeta(32)=7.48078421883e-14
Stats: nr=32 min=1 max=235 sum3=42515861
Emax=32.5437372723 nrkept=500 nrkeptmult=1144 
Kept: 500 out of 1865, ratio=0.268
Memory used: 88 MB  Time elapsed: 36 s

 ***** [DMNRG] Iteration 33/34 (scale 3.5627958595174936e-10) ***** 
[1] xi(32)=1.00142015791 zeta(33)=-5.42322738709e-14
Stats: nr=33 min=1 max=236 sum3=41601690
Emax=32.9788693712 nrkept=500 nrkeptmult=1140 
Kept: 500 out of 1866, ratio=0.268
Memory used: 90 MB  Time elapsed: 36 s

 ***** [DMNRG] Iteration 34/34 (scale 1.7813979297587468e-10) ***** 
[1] xi(33)=1.00141973438 zeta(34)=4.6767745219e-15
Stats: nr=32 min=1 max=235 sum3=42515861
Safeguard: keep additional 3 states
Emax=32.501614643 nrkept=503 nrkeptmult=1152 
Kept: 503 out of 1865, ratio=0.27
Memory used: 91 MB  Time elapsed: 37 s

** Iteration completed.

Spectrum: A_a-A_a FDM
[A_a-A_a] pos=0.496612332726 neg=0.503387667273 sum= 1
mu1=0.00104472857349 mu2=0.025036187937 mu3=0.000115561931374 mu4=0.000988996178605

Spectrum: A_d-A_d FDM
[A_d-A_d] pos=0.601267551259 neg=0.39873244874 sum= 1
mu1=0.0263079399171 mu2=0.00576815460237 mu3=0.00035664192472 mu4=0.000132255641193

Spectrum: self_d-A_d FDM
[self_d-A_d] pos=0.0291906850892 neg=-0.00310021196458 sum= 0.0260904731246
mu1=0.00318936887177 mu2=0.000218796163911 mu3=8.01291666679e-05 mu4=5.93726786446e-06


FDM expectation values
<Hhyb>=-0.222649428592
<Himp>=0.0208062364995
<Hpot>=0
<n_a>=1.00677533783
<n_d>=0.797464898393
<n_d_ud>=0.116272191009

Memory usage report
===================
  after diag: 96088 kB
after recalc: 175624 kB
 after split: 112880 kB
  after trim: 93756 kB
 after trunc: 96088 kB
      diag_h: 74712 kB

Peak usage: 175624 kB

Timing report
=============
         All: 37 s
         FDM: 1 s
        diag: 11 s
    recalc d: 4 s
    recalc s: 1 s
      spec d: 13 s
       Other: 5 s
