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

No MPI: single node calculation.

workdir=.
Tmin=1e-10 -> Nmax=67
TEMPERATURE T=1.29508547824e-10
Nmax=67
Tmin=1e-10
Ninit=0
Lambda=2
strategy=kept
lastall=0
lastalloverride=0
z=0.375
discretization=Z
band=flat
betabar=1
T=1.29508547824e-10
T_ratio=1
Tmin_ratio=-999
keep=10000
keepenergy=10
keepmin=0
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=100
discard_trim=1e-16
discard_immediately=1e-16
NN1=0
NN2even=1
NN2avg=0
NNtanh=0
goodE=2
savebins=0
broaden=1
broaden_max=2.5
broaden_min=-99
broaden_min_ratio=3
broaden_ratio=1.05
dumpenergies=0
saveall=0
dumpannotated=0
dumpscaled=1
dumpabs=0
dumpprecision=8
dumpgroups=1
grouptol=1e-06
dumpdiagonal=0
ops=A_d
specs=
specd=
spect=
specb=
reim=0
trsq=
specgt=A_d-A_d
gtp=0.7
speci1t=A_d-A_d
speci2t=A_d-A_d
speckubo=
kubop=0.7
kubor=4
kubob=2
cond=
specchit=
chitp=1
chitp_ratio=-999
dm=0
T0opt=-1
dmnrg=0
cfs=0
fdm=0
fdmexpv=0
fdmexpvn=0
mats=0
nromegan=100
width_td=12
width_custom=12
prec_td=6
prec_custom=6
prec_xy=6
safeguard=1e-05
safeguardmax=200
fixeps=1e-15
diag=dsyev
diagratio=1
dsyevrlimit=100
restart=1
restartfactor=2
checkdiag=0
checkrho=0
calc0=1
tdht=0
spin=2
tri=old
preccpp=2000
checksumrules=0
showEs=0
globalB=0
globalBx=0
polarized=0
pol2x2=0
mpi=0
saveF=0
trunc=
noimag=1
dumpsubspaces=0
done=0
resume=0
forcestop=-1
stopafter=
removefiles=1
store=one

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QS
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 		 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.375
>>  ops=	A_d
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=6
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||||> (z)

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

Number of invariant subspaces: 6
(-2 1) 1 states
0.0337033924096 
(-1 2) 2 states
0.0159983100249 0.0482303809588 
(0 1) 3 states
0 0.0337033924096 0.0642286909836 
(0 3) 1 states
0.030525298574 
(1 2) 2 states
0.0159983100249 0.0482303809588 
(2 1) 1 states
0.0337033924096 
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 GT
Spectrum: A_d-A_d I1T
Spectrum: A_d-A_d I2T

Before NRG iteration (N=-1)
GT[A_d-A_d]=0.226975746833
I1T[A_d-A_d]=-5.42101086243e-20
I2T[A_d-A_d]=0.000150169408761

 ***** [NRG] Iteration 1/67 (scale 1.1124699549106556) ***** 
[1] xi(0)=0.528127117018 zeta(1)=0
Stats: nr=10 min=1 max=8 sum3=1055
Emax=2.11893370802 nrkept=35 nrkeptmult=64 
Total energy=-1.17612504331  Egs=-1.01405031107
Kept: 35 out of 35, ratio=1
GT[A_d-A_d]=0.320955801434
I1T[A_d-A_d]=-9.71038570732e-18
I2T[A_d-A_d]=0.00018581755779
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/67 (scale 0.78663504898361736) ***** 
[1] xi(1)=0.716043705417 zeta(2)=0
Stats: nr=15 min=1 max=20 sum3=30066
Emax=4.22873602838 nrkept=126 nrkeptmult=256 
Total energy=-1.66073598641  Egs=-0.616055620364
Kept: 126 out of 126, ratio=1
GT[A_d-A_d]=0.453800164175
I1T[A_d-A_d]=9.11068638067e-18
I2T[A_d-A_d]=0.000240742063667
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/67 (scale 0.55623497745532779) ***** 
[1] xi(2)=1.0004392563 zeta(3)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=9.18404150142 nrkept=462 nrkeptmult=1024 
Total energy=-2.55174255875  Egs=-1.6018528292
Kept: 462 out of 462, ratio=1
GT[A_d-A_d]=0.641534288557
I1T[A_d-A_d]=4.76965787846e-16
I2T[A_d-A_d]=0.000287391719678
Memory used: 5 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/67 (scale 0.39331752449180868) ***** 
[1] xi(3)=1.10616421154 zeta(4)=0
Stats: nr=28 min=1 max=210 sum3=37226316
Safeguard: keep additional 9 states
Emax=10.0257297202 nrkept=1434 nrkeptmult=3470 
Total energy=-2.93976154956  Egs=-0.986528610215
Kept: 1434 out of 1716, ratio=0.836
GT[A_d-A_d]=0.90672941635
I1T[A_d-A_d]=7.33362369812e-16
I2T[A_d-A_d]=0.000325288315638
Memory used: 15 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/67 (scale 0.2781174887276639) ***** 
[1] xi(4)=1.0497443115 zeta(5)=0
Stats: nr=36 min=1 max=638 sum3=824058469
Safeguard: keep additional 13 states
Emax=10.0114089287 nrkept=1893 nrkeptmult=4702 
Total energy=-3.37404677998  Egs=-1.5615171574
Kept: 1893 out of 5389, ratio=0.351
GT[A_d-A_d]=1.28117210088
I1T[A_d-A_d]=-3.94382590118e-16
I2T[A_d-A_d]=0.000345521481125
Memory used: 24 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 6/67 (scale 0.19665876224590434) ***** 
[1] xi(5)=1.01133478261 zeta(6)=0
Stats: nr=42 min=1 max=721 sum3=1741815580
Safeguard: keep additional 1 states
Emax=10.0483939673 nrkept=2168 nrkeptmult=5380 
Total energy=-3.5485939692  Egs=-0.887563753712
Kept: 2168 out of 7126, ratio=0.304
GT[A_d-A_d]=1.80946300211
I1T[A_d-A_d]=-1.31476096896e-15
I2T[A_d-A_d]=0.000362199563992
Memory used: 27 MB  Time elapsed: 7 s

 ***** [NRG] Iteration 7/67 (scale 0.13905874436383195) ***** 
[1] xi(6)=1.00233827108 zeta(7)=0
Stats: nr=44 min=1 max=933 sum3=2698661041
Safeguard: keep additional 2 states
Emax=10.0154472091 nrkept=1921 nrkeptmult=4746 
Total energy=-3.75692462425  Egs=-1.49814854152
Kept: 1921 out of 8173, ratio=0.235
GT[A_d-A_d]=2.55398426542
I1T[A_d-A_d]=-2.24441571267e-15
I2T[A_d-A_d]=0.000377973045044
Memory used: 29 MB  Time elapsed: 13 s

 ***** [NRG] Iteration 8/67 (scale 0.09832938112295217) ***** 
[1] xi(7)=1.00063056638 zeta(8)=0
Stats: nr=42 min=1 max=745 sum3=1897784946
Safeguard: keep additional 2 states
Emax=10.0498064528 nrkept=2121 nrkeptmult=5273 
Total energy=-3.84581501183  Egs=-0.904006376994
Kept: 2121 out of 7230, ratio=0.293
GT[A_d-A_d]=3.60123858996
I1T[A_d-A_d]=-4.50285526182e-15
I2T[A_d-A_d]=0.000403299288505
Memory used: 29 MB  Time elapsed: 18 s

 ***** [NRG] Iteration 9/67 (scale 0.069529372181915974) ***** 
[1] xi(8)=1.00021749239 zeta(9)=0
Stats: nr=44 min=1 max=911 sum3=2522727317
Safeguard: keep additional 1 states
Emax=10.0025168425 nrkept=1680 nrkeptmult=4174 
Total energy=-3.94767202438  Egs=-1.46494940705
Kept: 1680 out of 7997, ratio=0.21
GT[A_d-A_d]=5.06998425302
I1T[A_d-A_d]=-6.90233526542e-15
I2T[A_d-A_d]=0.000434345723768
Memory used: 27 MB  Time elapsed: 24 s

 ***** [NRG] Iteration 10/67 (scale 0.049164690561476085) ***** 
[1] xi(9)=1.00008734701 zeta(10)=0
Stats: nr=42 min=1 max=660 sum3=1300854274
Emax=10.0325543922 nrkept=1873 nrkeptmult=4643 
Total energy=-3.99345663599  Egs=-0.931249868189
Kept: 1873 out of 6346, ratio=0.295
GT[A_d-A_d]=7.11940166554
I1T[A_d-A_d]=-1.06013388217e-14
I2T[A_d-A_d]=0.000479797500502
Memory used: 29 MB  Time elapsed: 28 s

 ***** [NRG] Iteration 11/67 (scale 0.034764686090957987) ***** 
[1] xi(10)=1.00003860055 zeta(11)=0
Stats: nr=42 min=1 max=815 sum3=1780810507
Safeguard: keep additional 7 states
Emax=10.0037835727 nrkept=1638 nrkeptmult=4060 
Total energy=-4.04314149458  Egs=-1.42917610303
Kept: 1638 out of 7057, ratio=0.232
GT[A_d-A_d]=9.95538384279
I1T[A_d-A_d]=-1.42455239656e-14
I2T[A_d-A_d]=0.00053611479821
Memory used: 28 MB  Time elapsed: 32 s

 ***** [NRG] Iteration 12/67 (scale 0.024582345280738042) ***** 
[1] xi(11)=1.00001806211 zeta(12)=0
Stats: nr=42 min=1 max=640 sum3=1188522258
Safeguard: keep additional 1 states
Emax=10.0243818761 nrkept=1792 nrkeptmult=4448 
Total energy=-4.06675613444  Egs=-0.960634129719
Kept: 1792 out of 6180, ratio=0.29
GT[A_d-A_d]=13.8211131044
I1T[A_d-A_d]=-1.70350530511e-14
I2T[A_d-A_d]=0.000610691890257
Memory used: 30 MB  Time elapsed: 36 s

 ***** [NRG] Iteration 13/67 (scale 0.017382343045478994) ***** 
[1] xi(12)=1.00000872507 zeta(13)=0
Stats: nr=42 min=1 max=778 sum3=1552502967
Emax=10.0009506982 nrkept=1551 nrkeptmult=3816 
Total energy=-4.09095191967  Egs=-1.39197490021
Kept: 1551 out of 6753, ratio=0.23
GT[A_d-A_d]=18.959743346
I1T[A_d-A_d]=-1.91637888017e-14
I2T[A_d-A_d]=0.000695855916194
Memory used: 28 MB  Time elapsed: 40 s

 ***** [NRG] Iteration 14/67 (scale 0.012291172640369021) ***** 
[1] xi(13)=1.00000428647 zeta(14)=0
Stats: nr=42 min=1 max=608 sum3=1015204456
Safeguard: keep additional 3 states
Emax=10.0024835866 nrkept=1622 nrkeptmult=3996 
Total energy=-4.10311526564  Egs=-0.989600123859
Kept: 1622 out of 5842, ratio=0.278
GT[A_d-A_d]=25.4838767191
I1T[A_d-A_d]=-1.96081479801e-14
I2T[A_d-A_d]=0.000787744294864
Memory used: 28 MB  Time elapsed: 43 s

 ***** [NRG] Iteration 15/67 (scale 0.0086911715227394968) ***** 
[1] xi(14)=1.00000212427 zeta(15)=0
Stats: nr=40 min=1 max=712 sum3=1175944449
Emax=10.0235597405 nrkept=1539 nrkeptmult=3784 
Total energy=-4.11489360055  Egs=-1.35520681909
Kept: 1539 out of 6105, ratio=0.252
GT[A_d-A_d]=33.1486926557
I1T[A_d-A_d]=-1.73483053366e-14
I2T[A_d-A_d]=0.000859021419525
Memory used: 29 MB  Time elapsed: 46 s

 ***** [NRG] Iteration 16/67 (scale 0.0061455863201845106) ***** 
[1] xi(15)=1.0000010574 zeta(16)=0
Stats: nr=42 min=1 max=602 sum3=986928322
Safeguard: keep additional 2 states
Emax=10.0042970752 nrkept=1481 nrkeptmult=3599 
Total energy=-4.12114709864  Egs=-1.01755923086
Kept: 1481 out of 5794, ratio=0.256
GT[A_d-A_d]=40.9874000243
I1T[A_d-A_d]=-9.28673979647e-15
I2T[A_d-A_d]=0.000876884404331
Memory used: 29 MB  Time elapsed: 49 s

 ***** [NRG] Iteration 17/67 (scale 0.0043455857613697484) ***** 
[1] xi(16)=1.00000052752 zeta(17)=0
Stats: nr=36 min=1 max=663 sum3=929268912
Safeguard: keep additional 3 states
Emax=10.004436577 nrkept=1361 nrkeptmult=3360 
Total energy=-4.12687745613  Egs=-1.31866169678
Kept: 1361 out of 5574, ratio=0.244
GT[A_d-A_d]=47.3009863582
I1T[A_d-A_d]=4.54180612153e-15
I2T[A_d-A_d]=0.000798516277845
Memory used: 28 MB  Time elapsed: 52 s

 ***** [NRG] Iteration 18/67 (scale 0.0030727931600922553) ***** 
[1] xi(17)=1.00000026346 zeta(18)=0
Stats: nr=40 min=1 max=540 sum3=707668730
Safeguard: keep additional 2 states
Emax=10.0203767056 nrkept=1306 nrkeptmult=3196 
Total energy=-4.13009405945  Egs=-1.04680112013
Kept: 1306 out of 5120, ratio=0.255
GT[A_d-A_d]=50.4847813197
I1T[A_d-A_d]=1.95305426494e-14
I2T[A_d-A_d]=0.000620702309098
Memory used: 28 MB  Time elapsed: 54 s

 ***** [NRG] Iteration 19/67 (scale 0.0021727928806848742) ***** 
[1] xi(18)=1.00000013166 zeta(19)=0
Stats: nr=34 min=1 max=592 sum3=656857738
Safeguard: keep additional 1 states
Emax=10.0234937293 nrkept=1178 nrkeptmult=2914 
Total energy=-4.13287477787  Egs=-1.2797899167
Kept: 1178 out of 4918, ratio=0.24
GT[A_d-A_d]=50.7987274648
I1T[A_d-A_d]=2.72753461569e-14
I2T[A_d-A_d]=0.000397151345278
Memory used: 27 MB  Time elapsed: 56 s

 ***** [NRG] Iteration 20/67 (scale 0.0015363965800461277) ***** 
[1] xi(19)=1.00000006581 zeta(20)=0
Stats: nr=40 min=1 max=471 sum3=471392010
Safeguard: keep additional 1 states
Emax=10.0118777693 nrkept=1110 nrkeptmult=2756 
Total energy=-4.1345358588  Egs=-1.0811537497
Kept: 1110 out of 4440, ratio=0.25
GT[A_d-A_d]=50.6843818193
I1T[A_d-A_d]=2.32404546318e-14
I2T[A_d-A_d]=0.000210563724013
Memory used: 29 MB  Time elapsed: 58 s

 ***** [NRG] Iteration 21/67 (scale 0.0010863964403424371) ***** 
[1] xi(20)=1.0000000329 zeta(21)=0
Stats: nr=34 min=1 max=505 sum3=407812833
Safeguard: keep additional 2 states
Emax=10.0038977216 nrkept=973 nrkeptmult=2444 
Total energy=-4.1358778607  Egs=-1.23527825889
Kept: 973 out of 4185, ratio=0.232
GT[A_d-A_d]=52.6452635668
I1T[A_d-A_d]=1.39972625859e-14
I2T[A_d-A_d]=9.95366395811e-05
Memory used: 26 MB  Time elapsed: 59 s

 ***** [NRG] Iteration 22/67 (scale 0.00076819829002306383) ***** 
[1] xi(21)=1.00000001645 zeta(22)=0
Stats: nr=40 min=1 max=397 sum3=276136656
Safeguard: keep additional 3 states
Emax=10.0025869864 nrkept=907 nrkeptmult=2253 
Total energy=-4.13674133807  Egs=-1.12402928231
Kept: 907 out of 3672, ratio=0.247
GT[A_d-A_d]=57.6428446835
I1T[A_d-A_d]=6.67014210233e-15
I2T[A_d-A_d]=4.72355170437e-05
Memory used: 29 MB  Time elapsed: 60 s

 ***** [NRG] Iteration 23/67 (scale 0.00054319822017121855) ***** 
[1] xi(22)=1.00000000822 zeta(23)=0
Stats: nr=34 min=1 max=421 sum3=232310447
Safeguard: keep additional 2 states
Emax=10.0009806187 nrkept=876 nrkeptmult=2234 
Total energy=-4.13738372003  Egs=-1.18259216556
Kept: 876 out of 3425, ratio=0.256
GT[A_d-A_d]=65.2983424779
I1T[A_d-A_d]=2.86937085115e-15
I2T[A_d-A_d]=2.43125785361e-05
Memory used: 29 MB  Time elapsed: 61 s

 ***** [NRG] Iteration 24/67 (scale 0.00038409914501153191) ***** 
[1] xi(23)=1.00000000411 zeta(24)=0
Stats: nr=40 min=1 max=358 sum3=198470116
Safeguard: keep additional 2 states
Emax=10.0255241367 nrkept=829 nrkeptmult=2051 
Total energy=-4.13783612633  Egs=-1.17783728097
Kept: 829 out of 3316, ratio=0.25
GT[A_d-A_d]=75.5421361184
I1T[A_d-A_d]=1.32199196846e-15
I2T[A_d-A_d]=1.33257149975e-05
Memory used: 28 MB  Time elapsed: 62 s

 ***** [NRG] Iteration 25/67 (scale 0.00027159911008560927) ***** 
[1] xi(24)=1.00000000206 zeta(25)=0
Stats: nr=34 min=1 max=385 sum3=175545704
Safeguard: keep additional 3 states
Emax=10.0045978413 nrkept=820 nrkeptmult=2078 
Total energy=-4.13814053354  Egs=-1.12079606224
Kept: 820 out of 3122, ratio=0.263
GT[A_d-A_d]=88.2105532032
I1T[A_d-A_d]=6.45061398596e-16
I2T[A_d-A_d]=7.60015042156e-06
Memory used: 29 MB  Time elapsed: 63 s

 ***** [NRG] Iteration 26/67 (scale 0.00019204957250576596) ***** 
[1] xi(25)=1.00000000103 zeta(26)=0
Stats: nr=38 min=1 max=335 sum3=163906502
Safeguard: keep additional 3 states
Emax=10.0074604307 nrkept=783 nrkeptmult=1957 
Total energy=-4.13837938251  Egs=-1.24368395626
Kept: 783 out of 3098, ratio=0.253
GT[A_d-A_d]=103.71433958
I1T[A_d-A_d]=1.87082195061e-16
I2T[A_d-A_d]=4.41645428602e-06
Memory used: 29 MB  Time elapsed: 63 s

 ***** [NRG] Iteration 27/67 (scale 0.00013579955504280464) ***** 
[1] xi(26)=1.00000000051 zeta(27)=0
Stats: nr=34 min=1 max=359 sum3=144739908
Emax=10.0015804894 nrkept=798 nrkeptmult=2024 
Total energy=-4.13852219469  Egs=-1.05163950859
Kept: 798 out of 2952, ratio=0.27
GT[A_d-A_d]=122.105661972
I1T[A_d-A_d]=-7.03116027641e-17
I2T[A_d-A_d]=2.59680873209e-06
Memory used: 29 MB  Time elapsed: 64 s

 ***** [NRG] Iteration 28/67 (scale 9.6024786252882978e-05) ***** 
[1] xi(27)=1.00000000026 zeta(28)=0
Stats: nr=38 min=1 max=327 sum3=150905436
Safeguard: keep additional 2 states
Emax=10.0411184273 nrkept=763 nrkeptmult=1925 
Total energy=-4.13864893419  Egs=-1.31986235163
Kept: 763 out of 3012, ratio=0.253
GT[A_d-A_d]=143.519494563
I1T[A_d-A_d]=-2.10372530617e-16
I2T[A_d-A_d]=1.54105469967e-06
Memory used: 29 MB  Time elapsed: 65 s

 ***** [NRG] Iteration 29/67 (scale 6.7899777521402319e-05) ***** 
[1] xi(28)=1.00000000013 zeta(29)=0
Stats: nr=34 min=1 max=349 sum3=134658679
Safeguard: keep additional 1 states
Emax=10.0003360016 nrkept=768 nrkeptmult=1922 
Total energy=-4.13871554383  Egs=-0.980999430758
Kept: 768 out of 2881, ratio=0.267
GT[A_d-A_d]=167.786793106
I1T[A_d-A_d]=-2.45372626002e-16
I2T[A_d-A_d]=9.13473958049e-07
Memory used: 29 MB  Time elapsed: 65 s

 ***** [NRG] Iteration 30/67 (scale 4.8012393126441489e-05) ***** 
[1] xi(29)=1.00000000006 zeta(30)=0
Stats: nr=38 min=1 max=317 sum3=138601396
Safeguard: keep additional 3 states
Emax=10.0174694642 nrkept=721 nrkeptmult=1795 
Total energy=-4.13878274002  Egs=-1.39955932591
Kept: 721 out of 2896, ratio=0.249
GT[A_d-A_d]=193.813019862
I1T[A_d-A_d]=-2.41029303264e-16
I2T[A_d-A_d]=5.44325504016e-07
Memory used: 30 MB  Time elapsed: 66 s

 ***** [NRG] Iteration 31/67 (scale 3.3949888760701159e-05) ***** 
[1] xi(30)=1.00000000003 zeta(31)=0
Stats: nr=33 min=1 max=332 sum3=115086630
Safeguard: keep additional 2 states
Emax=10.0003683685 nrkept=725 nrkeptmult=1794 
Total energy=-4.13881391174  Egs=-0.918168554055
Kept: 725 out of 2718, ratio=0.267
GT[A_d-A_d]=220.878146585
I1T[A_d-A_d]=-2.02345331724e-16
I2T[A_d-A_d]=3.19571515598e-07
Memory used: 31 MB  Time elapsed: 66 s

 ***** [NRG] Iteration 32/67 (scale 2.4006196563220745e-05) ***** 
[1] xi(31)=1.00000000002 zeta(32)=0
Stats: nr=36 min=1 max=303 sum3=121137592
Safeguard: keep additional 2 states
Emax=10.0044571199 nrkept=656 nrkeptmult=1606 
Total energy=-4.13884923281  Egs=-1.47133134204
Kept: 656 out of 2734, ratio=0.24
GT[A_d-A_d]=246.019475387
I1T[A_d-A_d]=-1.73511886831e-16
I2T[A_d-A_d]=1.87081755385e-07
Memory used: 31 MB  Time elapsed: 67 s

 ***** [NRG] Iteration 33/67 (scale 1.697494438035058e-05) ***** 
[1] xi(32)=1.00000000001 zeta(33)=0
Stats: nr=33 min=1 max=311 sum3=91356757
Safeguard: keep additional 2 states
Emax=10.0189232103 nrkept=688 nrkeptmult=1688 
Total energy=-4.13886401371  Egs=-0.870748081799
Kept: 688 out of 2473, ratio=0.278
GT[A_d-A_d]=268.925101242
I1T[A_d-A_d]=-1.32871656303e-16
I2T[A_d-A_d]=1.06171077554e-07
Memory used: 31 MB  Time elapsed: 67 s

 ***** [NRG] Iteration 34/67 (scale 1.2003098281610372e-05) ***** 
[1] xi(33)=1 zeta(34)=0
Stats: nr=36 min=1 max=286 sum3=104077260
Safeguard: keep additional 3 states
Emax=10.0020289857 nrkept=602 nrkeptmult=1444 
Total energy=-4.13888232503  Egs=-1.52554943702
Kept: 602 out of 2592, ratio=0.232
GT[A_d-A_d]=285.957126244
I1T[A_d-A_d]=-7.5293681815e-17
I2T[A_d-A_d]=5.95509105956e-08
Memory used: 31 MB  Time elapsed: 68 s

 ***** [NRG] Iteration 35/67 (scale 8.4874721901752898e-06) ***** 
[1] xi(34)=1 zeta(35)=0
Stats: nr=33 min=1 max=291 sum3=72419340
Safeguard: keep additional 2 states
Emax=10.0021340061 nrkept=612 nrkeptmult=1466 
Total energy=-4.13888945412  Egs=-0.839953883729
Kept: 612 out of 2262, ratio=0.271
GT[A_d-A_d]=299.589953588
I1T[A_d-A_d]=-2.02131409583e-17
I2T[A_d-A_d]=3.20534836419e-08
Memory used: 32 MB  Time elapsed: 68 s

 ***** [NRG] Iteration 36/67 (scale 6.0015491408051861e-06) ***** 
[1] xi(35)=1 zeta(36)=0
Stats: nr=36 min=1 max=260 sum3=75905034
Safeguard: keep additional 2 states
Emax=10.0058420339 nrkept=540 nrkeptmult=1278 
Total energy=-4.13889881968  Egs=-1.56052366949
Kept: 540 out of 2298, ratio=0.235
GT[A_d-A_d]=306.869393258
I1T[A_d-A_d]=6.6058494826e-18
I2T[A_d-A_d]=1.70864708884e-08
Memory used: 32 MB  Time elapsed: 68 s

 ***** [NRG] Iteration 37/67 (scale 4.2437360950876449e-06) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=33 min=1 max=263 sum3=53120001
Safeguard: keep additional 1 states
Emax=10.0287783248 nrkept=569 nrkeptmult=1328 
Total energy=-4.13890230819  Egs=-0.822036946394
Kept: 569 out of 2025, ratio=0.281
GT[A_d-A_d]=313.008596551
I1T[A_d-A_d]=1.78263394249e-17
I2T[A_d-A_d]=8.77865435248e-09
Memory used: 32 MB  Time elapsed: 69 s

 ***** [NRG] Iteration 38/67 (scale 3.0007745704025931e-06) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=34 min=1 max=245 sum3=62083694
Safeguard: keep additional 3 states
Emax=10.0217178284 nrkept=509 nrkeptmult=1181 
Total energy=-4.1389070515  Egs=-1.58069477345
Kept: 509 out of 2132, ratio=0.239
GT[A_d-A_d]=314.537917722
I1T[A_d-A_d]=1.9688919381e-17
I2T[A_d-A_d]=4.52190073295e-09
Memory used: 32 MB  Time elapsed: 69 s

 ***** [NRG] Iteration 39/67 (scale 2.1218680475438225e-06) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=33 min=1 max=249 sum3=44688876
Emax=10.3903501477 nrkept=567 nrkeptmult=1322 
Total energy=-4.1389087751  Egs=-0.812304784823
Kept: 567 out of 1902, ratio=0.298
GT[A_d-A_d]=317.37798101
I1T[A_d-A_d]=1.90419210397e-17
I2T[A_d-A_d]=2.26685528969e-09
Memory used: 32 MB  Time elapsed: 69 s

 ***** [NRG] Iteration 40/67 (scale 1.5003872852012965e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.1648320005 nrkept=499 nrkeptmult=1149 
Total energy=-4.13891116307  Egs=-1.5915721116
Kept: 499 out of 2124, ratio=0.235
GT[A_d-A_d]=316.752080748
I1T[A_d-A_d]=1.57101336921e-17
I2T[A_d-A_d]=1.15105537624e-09
Memory used: 32 MB  Time elapsed: 70 s

 ***** [NRG] Iteration 41/67 (scale 1.0609340237719112e-06) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.6535833464 nrkept=567 nrkeptmult=1322 
Total energy=-4.13891201948  Egs=-0.807221685895
Kept: 567 out of 1862, ratio=0.305
GT[A_d-A_d]=318.598296449
I1T[A_d-A_d]=1.1182734524e-17
I2T[A_d-A_d]=5.7198738692e-10
Memory used: 32 MB  Time elapsed: 70 s

 ***** [NRG] Iteration 42/67 (scale 7.5019364260064827e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Safeguard: keep additional 1 states
Emax=10.3807619665 nrkept=499 nrkeptmult=1149 
Total energy=-4.13891321771  Egs=-1.59722623563
Kept: 499 out of 2124, ratio=0.235
GT[A_d-A_d]=317.322093541
I1T[A_d-A_d]=6.82799803848e-18
I2T[A_d-A_d]=2.89268992506e-10
Memory used: 32 MB  Time elapsed: 70 s

 ***** [NRG] Iteration 43/67 (scale 5.3046701188595561e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=33 min=1 max=245 sum3=42367736
Emax=10.7458471319 nrkept=567 nrkeptmult=1322 
Total energy=-4.13891364454  Egs=-0.804622399568
Kept: 567 out of 1862, ratio=0.305
GT[A_d-A_d]=318.925228102
I1T[A_d-A_d]=2.86317215981e-18
I2T[A_d-A_d]=1.43340350117e-10
Memory used: 32 MB  Time elapsed: 70 s

 ***** [NRG] Iteration 44/67 (scale 3.7509682130032413e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.4988723747 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891424473  Egs=-1.60010951565
Kept: 498 out of 2124, ratio=0.234
GT[A_d-A_d]=317.459390941
I1T[A_d-A_d]=-4.98037692e-20
I2T[A_d-A_d]=7.24172324896e-11
Memory used: 32 MB  Time elapsed: 71 s

 ***** [NRG] Iteration 45/67 (scale 2.6523350594297781e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.7936487764 nrkept=567 nrkeptmult=1322 
Total energy=-4.1389144578  Egs=-0.803307843106
Kept: 567 out of 1857, ratio=0.305
GT[A_d-A_d]=319.014634318
I1T[A_d-A_d]=-9.39623615916e-19
I2T[A_d-A_d]=3.58563083038e-11
Memory used: 33 MB  Time elapsed: 71 s

 ***** [NRG] Iteration 46/67 (scale 1.8754841065016207e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5123877372 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891475817  Egs=-1.60156551534
Kept: 498 out of 2124, ratio=0.234
GT[A_d-A_d]=317.49064596
I1T[A_d-A_d]=-9.4622767083e-19
I2T[A_d-A_d]=1.81120553819e-11
Memory used: 32 MB  Time elapsed: 71 s

 ***** [NRG] Iteration 47/67 (scale 1.326167529714889e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8178715531 nrkept=567 nrkeptmult=1322 
Total energy=-4.13891486461  Egs=-0.802646773477
Kept: 567 out of 1857, ratio=0.305
GT[A_d-A_d]=319.040288857
I1T[A_d-A_d]=-4.8754702558e-19
I2T[A_d-A_d]=8.96534105931e-12
Memory used: 33 MB  Time elapsed: 71 s

 ***** [NRG] Iteration 48/67 (scale 9.3774205325081033e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5192088906 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891501487  Egs=-1.60229714396
Kept: 498 out of 2124, ratio=0.234
GT[A_d-A_d]=317.496629476
I1T[A_d-A_d]=-1.94446613213e-19
I2T[A_d-A_d]=4.52863597297e-12
Memory used: 33 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 49/67 (scale 6.6308376485744452e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8296445148 nrkept=567 nrkeptmult=1322 
Total energy=-4.13891506807  Egs=-0.802315282717
Kept: 567 out of 1857, ratio=0.305
GT[A_d-A_d]=319.048318941
I1T[A_d-A_d]=-9.46598609636e-20
I2T[A_d-A_d]=2.2414045258e-12
Memory used: 33 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 50/67 (scale 4.6887102662540517e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=34 min=1 max=244 sum3=61218744
Emax=10.5226346405 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891514321  Egs=-1.60266387035
Kept: 498 out of 2124, ratio=0.234
GT[A_d-A_d]=317.497203116
I1T[A_d-A_d]=-4.78780991963e-20
I2T[A_d-A_d]=1.13221489777e-12
Memory used: 33 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 51/67 (scale 3.3154188242872226e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.833719195 nrkept=567 nrkeptmult=1324 
Total energy=-4.13891516981  Egs=-0.802149297327
Kept: 567 out of 1857, ratio=0.305
GT[A_d-A_d]=319.051206475
I1T[A_d-A_d]=-1.01607350928e-19
I2T[A_d-A_d]=5.60354242828e-13
Memory used: 34 MB  Time elapsed: 72 s

 ***** [NRG] Iteration 52/67 (scale 2.3443551331270258e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5243633897 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891520738  Egs=-1.60284746218
Kept: 498 out of 2124, ratio=0.234
GT[A_d-A_d]=317.496880388
I1T[A_d-A_d]=-1.10880766147e-19
I2T[A_d-A_d]=2.83059329834e-13
Memory used: 34 MB  Time elapsed: 73 s

 ***** [NRG] Iteration 53/67 (scale 1.6577094121436113e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8356839004 nrkept=567 nrkeptmult=1324 
Total energy=-4.13891522068  Egs=-0.802066244501
Kept: 567 out of 1857, ratio=0.305
GT[A_d-A_d]=319.052303534
I1T[A_d-A_d]=-1.24655155422e-19
I2T[A_d-A_d]=1.40088621865e-13
Memory used: 35 MB  Time elapsed: 73 s

 ***** [NRG] Iteration 54/67 (scale 1.1721775665635129e-08) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5252232913 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891523947  Egs=-1.60293931534
Kept: 498 out of 2124, ratio=0.234
GT[A_d-A_d]=317.496557712
I1T[A_d-A_d]=-1.30185151854e-19
I2T[A_d-A_d]=7.07654460394e-14
Memory used: 35 MB  Time elapsed: 73 s

 ***** [NRG] Iteration 55/67 (scale 8.2885470607180565e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8366652827 nrkept=567 nrkeptmult=1324 
Total energy=-4.13891524612  Egs=-0.80202470304
Kept: 567 out of 1857, ratio=0.305
GT[A_d-A_d]=319.052768162
I1T[A_d-A_d]=-1.13544607606e-19
I2T[A_d-A_d]=3.50221405498e-14
Memory used: 36 MB  Time elapsed: 73 s

 ***** [NRG] Iteration 56/67 (scale 5.8608878328175646e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5256535002 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891525551  Egs=-1.60298525623
Kept: 498 out of 2124, ratio=0.234
GT[A_d-A_d]=317.496349764
I1T[A_d-A_d]=-8.44588701246e-20
I2T[A_d-A_d]=1.76914323522e-14
Memory used: 35 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 57/67 (scale 4.1442735303590282e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8371557229 nrkept=567 nrkeptmult=1324 
Total energy=-4.13891525883  Egs=-0.802003928546
Kept: 567 out of 1857, ratio=0.305
GT[A_d-A_d]=319.052972806
I1T[A_d-A_d]=-5.08670387916e-20
I2T[A_d-A_d]=8.7555316788e-15
Memory used: 36 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 58/67 (scale 2.9304439164087823e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5258686692 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891526353  Egs=-1.60300823026
Kept: 498 out of 2124, ratio=0.234
GT[A_d-A_d]=317.496227437
I1T[A_d-A_d]=-2.47344693832e-20
I2T[A_d-A_d]=4.42286648397e-15
Memory used: 35 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 59/67 (scale 2.0721367651795141e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Emax=10.8374008807 nrkept=567 nrkeptmult=1324 
Total energy=-4.13891526519  Egs=-0.801993540358
Kept: 567 out of 1857, ratio=0.305
GT[A_d-A_d]=319.053061476
I1T[A_d-A_d]=-1.21292457821e-20
I2T[A_d-A_d]=2.18888234152e-15
Memory used: 36 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 60/67 (scale 1.4652219582043911e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=34 min=1 max=243 sum3=60833952
Emax=10.5259762698 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891526754  Egs=-1.60301971817
Kept: 498 out of 2124, ratio=0.234
GT[A_d-A_d]=317.496154983
I1T[A_d-A_d]=-7.98837216357e-21
I2T[A_d-A_d]=1.10571761843e-15
Memory used: 36 MB  Time elapsed: 74 s

 ***** [NRG] Iteration 61/67 (scale 1.0360683825897571e-09) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 2 states
Emax=10.8375271922 nrkept=569 nrkeptmult=1330 
Total energy=-4.13891526837  Egs=-0.801988346029
Kept: 569 out of 1857, ratio=0.306
GT[A_d-A_d]=319.053357088
I1T[A_d-A_d]=-5.66828099691e-21
I2T[A_d-A_d]=5.47221800133e-16
Memory used: 36 MB  Time elapsed: 75 s

 ***** [NRG] Iteration 62/67 (scale 7.3261097910219557e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=34 min=1 max=244 sum3=61697330
Emax=10.5260300742 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891526955  Egs=-1.60302546234
Kept: 498 out of 2132, ratio=0.234
GT[A_d-A_d]=317.496112225
I1T[A_d-A_d]=-2.4223028013e-21
I2T[A_d-A_d]=2.76429527784e-16
Memory used: 36 MB  Time elapsed: 75 s

 ***** [NRG] Iteration 63/67 (scale 5.1803419129487853e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8375999333 nrkept=575 nrkeptmult=1340 
Total energy=-4.13891526996  Egs=-0.801985748805
Kept: 575 out of 1857, ratio=0.31
GT[A_d-A_d]=319.053781611
I1T[A_d-A_d]=8.59827589561e-23
I2T[A_d-A_d]=1.36805978057e-16
Memory used: 36 MB  Time elapsed: 75 s

 ***** [NRG] Iteration 64/67 (scale 3.6630548955109779e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.5260425466 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891527055  Egs=-1.60302833449
Kept: 498 out of 2154, ratio=0.231
GT[A_d-A_d]=317.496087174
I1T[A_d-A_d]=5.613242891e-22
I2T[A_d-A_d]=6.91073976552e-17
Memory used: 37 MB  Time elapsed: 75 s

 ***** [NRG] Iteration 65/67 (scale 2.5901709564743927e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 8 states
Emax=10.8376170598 nrkept=575 nrkeptmult=1340 
Total energy=-4.13891527076  Egs=-0.801984450179
Kept: 575 out of 1857, ratio=0.31
GT[A_d-A_d]=319.053776506
I1T[A_d-A_d]=-2.72179378086e-22
I2T[A_d-A_d]=3.42014915558e-17
Memory used: 37 MB  Time elapsed: 76 s

 ***** [NRG] Iteration 66/67 (scale 1.8315274477554889e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=34 min=1 max=248 sum3=64262334
Emax=10.526055999 nrkept=498 nrkeptmult=1142 
Total energy=-4.13891527105  Egs=-1.60302977058
Kept: 498 out of 2154, ratio=0.231
GT[A_d-A_d]=317.49606222
I1T[A_d-A_d]=-8.83840919413e-22
I2T[A_d-A_d]=1.7276850475e-17
Memory used: 37 MB  Time elapsed: 76 s

 ***** [NRG] Iteration 67/67 (scale 1.2950854782371963e-10) ***** 
[1] xi(66)=1 zeta(67)=0
Stats: nr=33 min=1 max=246 sum3=42402345
Safeguard: keep additional 12 states
Emax=10.8376394451 nrkept=579 nrkeptmult=1344 
Total energy=-4.13891527116  Egs=-0.801983800862
Kept: 579 out of 1857, ratio=0.312
GT[A_d-A_d]=319.05394871
I1T[A_d-A_d]=-9.2695740592e-22
I2T[A_d-A_d]=8.55038582003e-18
Memory used: 38 MB  Time elapsed: 76 s

** Iteration completed.

Spectrum: A_d-A_d GT

Spectrum: A_d-A_d I1T

Spectrum: A_d-A_d I2T


Total energy:             -4.138915271155344
TKW=2.15937843791e-05
Grand canonical Z=              13.43219079983303
Zft (last shell)=                13.4183239911852

Memory usage report
===================
  after diag: 59088 kB
after recalc: 39144 kB
 after split: 86044 kB
  after trim: 39144 kB
 after trunc: 59240 kB
      diag_h: 59088 kB

Peak usage: 86044 kB

Timing report
=============
         All: 76 s
        diag: 63 s
    recalc d: 1 s
     spec gt: 1 s
    spec i1t: 1 s
    spec i2t: 1 s
       Other: 7 s
