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 20:34:53

No MPI: single node calculation.

workdir=.
Tmin=1e-10 -> Nmax=67
TEMPERATURE T=1.18759861984e-10
Nmax=67
Tmin=1e-10
Ninit=0
Lambda=2
strategy=kept
lastall=0
lastalloverride=0
z=0.5
discretization=Z
band=flat
betabar=1
T=1.18759861984e-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=n_d^2 n_a^2 SaSd
specs=
specd=
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=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: 	../dqd.m	 Variant: 		 Channels: 	1
>>  U=	0.01	 Gamma=	0.0006	 Delta=	0	 t=	0.
>> 
>>  Gamma^(1/2)=	0.019544100476116797
>>  Lambda=	2.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	0.5
>>  ops=	n_d^2 n_a^2 SaSd
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=2 subspaces=10
Reading <||||> (e)
Reading <||n_d^2||> (s)
Reading <||n_a^2||> (s)
Reading <||SaSd||> (s)
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.5
# keep=10000 keepenergy=10 Nmax=67 channels=1 dots=2 betabar=1

Number of invariant subspaces: 10
(-3 1) 1 states
0.042770888302 
(-2 2) 3 states
0.0218015472497 0.0393051523992 0.0602744934515 
(-1 1) 6 states
0.00172225730669 0.0218015472497 0.0375828950925 0.042770888302 0.0602744934515 0.0786103047985 
(-1 3) 3 states
0.0183358113469 0.0393051523992 0.0568087575487 
(0 2) 8 states
0 0.0183358113469 0.0218015472497 0.0393051523992 0.0410274097059 0.0568087575487 0.0602744934515 0.0768880474918 
(0 4) 1 states
0.0358394164964 
(1 1) 6 states
0.00172225730669 0.0218015472497 0.0375828950925 0.042770888302 0.0602744934515 0.0786103047985 
(1 3) 3 states
0.0183358113469 0.0393051523992 0.0568087575487 
(2 2) 3 states
0.0218015472497 0.0393051523992 0.0602744934515 
(3 1) 1 states
0.042770888302 
Number of states (multiplicity taken into account): 64

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

Computing the following spectra:

Before NRG iteration (N=-1)
<SaSd>=3.65134923114e-13
<n_a^2>=1.49913361989
<n_d^2>=1.49913361989

 ***** [NRG] Iteration 1/67 (scale 1.0201394465967895) ***** 
[1] xi(0)=0.583569837945 zeta(1)=0
Stats: nr=15 min=1 max=20 sum3=30066
Emax=2.34657550001 nrkept=126 nrkeptmult=256 
Total energy=-1.19191711599  Egs=-1.11770187041
Kept: 126 out of 126, ratio=1
<SaSd>=1.9309117284e-12
<n_a^2>=1.49877482473
<n_d^2>=1.49877482473
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/67 (scale 0.72134752044448169) ***** 
[1] xi(1)=0.776754544846 zeta(2)=0
Stats: nr=21 min=1 max=75 sum3=1007202
Emax=4.64390358998 nrkept=462 nrkeptmult=1024 
Total energy=-1.66993416991  Egs=-0.662672346365
Kept: 462 out of 462, ratio=1
<SaSd>=1.03113449695e-11
<n_a^2>=1.49826753735
<n_d^2>=1.49826753735
Memory used: 5 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/67 (scale 0.51006972329839473) ***** 
[1] xi(2)=0.992903283056 zeta(3)=0
Stats: nr=28 min=1 max=210 sum3=37226316
Emax=9.71979427005 nrkept=1716 nrkeptmult=4096 
Total energy=-2.47388638692  Egs=-1.57616141537
Kept: 1716 out of 1716, ratio=1
<SaSd>=5.23882170675e-11
<n_a^2>=1.49755044224
<n_d^2>=1.49755044224
Memory used: 18 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/67 (scale 0.36067376022224085) ***** 
[1] xi(3)=1.04266411315 zeta(4)=0
Stats: nr=36 min=1 max=784 sum3=1484039727
Emax=10.2109807497 nrkept=4726 nrkeptmult=12162 
Total energy=-2.80129992769  Egs=-0.907783090655
Kept: 4726 out of 6435, ratio=0.734
<SaSd>=4.03691248515e-08
<n_a^2>=1.49653701466
<n_d^2>=1.49653701466
Memory used: 75 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 5/67 (scale 0.25503486164919736) ***** 
[1] xi(4)=1.01328740971 zeta(5)=0
Stats: nr=45 min=1 max=1707 sum3=10918885202
Safeguard: keep additional 3 states
Emax=10.2094147111 nrkept=7080 nrkeptmult=18696 
Total energy=-3.19037615236  Egs=-1.52558055064
Kept: 7080 out of 17876, ratio=0.396
<SaSd>=1.41992597567e-06
<n_a^2>=1.49510277959
<n_d^2>=1.49510277959
Memory used: 147 MB  Time elapsed: 79 s

 ***** [NRG] Iteration 6/67 (scale 0.18033688011112042) ***** 
[1] xi(5)=0.999788616061 zeta(6)=0
Stats: nr=52 min=1 max=2690 sum3=26872266871
Emax=10.0002202425 nrkept=7981 nrkeptmult=21088 
Total energy=-3.34956627607  Egs=-0.882737483408
Kept: 7981 out of 26847, ratio=0.297
<SaSd>=3.52912594462e-07
<n_a^2>=1.49308346839
<n_d^2>=1.49308346839
Memory used: 180 MB  Time elapsed: 316 s

 ***** [NRG] Iteration 7/67 (scale 0.12751743082459868) ***** 
[1] xi(6)=0.998344002082 zeta(7)=0
Stats: nr=55 min=1 max=2842 sum3=30658380056
Safeguard: keep additional 1 states
Emax=10.0141025434 nrkept=6948 nrkeptmult=18324 
Total energy=-3.53926198629  Egs=-1.48760611778
Kept: 6948 out of 30290, ratio=0.229
<SaSd>=1.74821944404e-08
<n_a^2>=1.49013883382
<n_d^2>=1.49013883382
Memory used: 148 MB  Time elapsed: 625 s

 ***** [NRG] Iteration 8/67 (scale 0.090168440055560212) ***** 
[1] xi(7)=0.998902918682 zeta(8)=0
Stats: nr=52 min=1 max=2705 sum3=26202719572
Safeguard: keep additional 1 states
Emax=10.0090846566 nrkept=7299 nrkeptmult=19210 
Total energy=-3.62108279984  Egs=-0.907421859578
Kept: 7299 out of 26394, ratio=0.277
<SaSd>=-4.03543087283e-07
<n_a^2>=1.4861387781
<n_d^2>=1.4861387781
Memory used: 158 MB  Time elapsed: 856 s

 ***** [NRG] Iteration 9/67 (scale 0.063758715412299341) ***** 
[1] xi(8)=0.999398551363 zeta(9)=0
Stats: nr=53 min=1 max=2618 sum3=28302445696
Safeguard: keep additional 2 states
Emax=10.0006005757 nrkept=6013 nrkeptmult=15763 
Total energy=-3.71405733579  Egs=-1.45822473591
Kept: 6013 out of 27674, ratio=0.217
<SaSd>=1.80436828518e-07
<n_a^2>=1.48047639395
<n_d^2>=1.48047639395
Memory used: 119 MB  Time elapsed: 1113 s

 ***** [NRG] Iteration 10/67 (scale 0.045084220027780106) ***** 
[1] xi(9)=0.99968731683 zeta(10)=0
Stats: nr=52 min=1 max=2334 sum3=30058121709
Safeguard: keep additional 1 states
Emax=10.0047605328 nrkept=6651 nrkeptmult=17446 
Total energy=-3.75614308084  Egs=-0.933491696872
Kept: 6651 out of 22771, ratio=0.292
<SaSd>=1.34676051104e-05
<n_a^2>=1.47251355089
<n_d^2>=1.47251355089
Memory used: 138 MB  Time elapsed: 1269 s

 ***** [NRG] Iteration 11/67 (scale 0.031879357706149671) ***** 
[1] xi(10)=0.9998407918 zeta(11)=0
Stats: nr=53 min=1 max=2398 sum3=30212532218
Safeguard: keep additional 1 states
Emax=10.0004410325 nrkept=5846 nrkeptmult=15276 
Total energy=-3.80164995477  Egs=-1.42747147992
Kept: 5846 out of 25210, ratio=0.232
<SaSd>=2.26752742523e-06
<n_a^2>=1.4615594792
<n_d^2>=1.4615594792
Memory used: 116 MB  Time elapsed: 1488 s

 ***** [NRG] Iteration 12/67 (scale 0.022542110013890053) ***** 
[1] xi(11)=0.99991969299 zeta(12)=0
Stats: nr=52 min=1 max=2277 sum3=26438797643
Safeguard: keep additional 2 states
Emax=10.0000699998 nrkept=5972 nrkeptmult=15574 
Total energy=-3.82323502635  Egs=-0.95754441663
Kept: 5972 out of 22143, ratio=0.27
<SaSd>=-6.87289288658e-06
<n_a^2>=1.4460575201
<n_d^2>=1.4460575201
Memory used: 116 MB  Time elapsed: 1641 s

 ***** [NRG] Iteration 13/67 (scale 0.015939678853074835) ***** 
[1] xi(12)=0.999959672404 zeta(13)=0
Stats: nr=51 min=1 max=2169 sum3=22224000922
Emax=10.0011178269 nrkept=5779 nrkeptmult=15125 
Total energy=-3.84553569625  Egs=-1.39906644957
Kept: 5779 out of 22630, ratio=0.255
<SaSd>=1.28154695259e-05
<n_a^2>=1.42507745197
<n_d^2>=1.42507745197
Memory used: 113 MB  Time elapsed: 1812 s

 ***** [NRG] Iteration 14/67 (scale 0.011271055006945026) ***** 
[1] xi(13)=0.999979792879 zeta(14)=0
Stats: nr=52 min=1 max=2252 sum3=25043494154
Safeguard: keep additional 1 states
Emax=10.0005412233 nrkept=5649 nrkeptmult=14724 
Total energy=-3.85655135007  Egs=-0.977340082015
Kept: 5649 out of 21906, ratio=0.258
<SaSd>=2.89264805823e-05
<n_a^2>=1.39639605533
<n_d^2>=1.39639605533
Memory used: 111 MB  Time elapsed: 1966 s

 ***** [NRG] Iteration 15/67 (scale 0.0079698394265374176) ***** 
[1] xi(14)=0.999989885634 zeta(15)=0
Stats: nr=49 min=1 max=2054 sum3=23481696368
Safeguard: keep additional 2 states
Emax=10.0009873329 nrkept=4989 nrkeptmult=13093 
Total energy=-3.86751261975  Egs=-1.37534385471
Kept: 4989 out of 21402, ratio=0.233
<SaSd>=0.000105805634552
<n_a^2>=1.35861487487
<n_d^2>=1.35861487487
Memory used: 97 MB  Time elapsed: 2115 s

 ***** [NRG] Iteration 16/67 (scale 0.0056355275034725132) ***** 
[1] xi(15)=0.999994940118 zeta(16)=0
Stats: nr=50 min=1 max=1984 sum3=15543748672
Safeguard: keep additional 2 states
Emax=10.0000124971 nrkept=4724 nrkeptmult=12358 
Total energy=-3.87310710377  Egs=-0.99271701085
Kept: 4724 out of 18922, ratio=0.25
<SaSd>=0.000299887075436
<n_a^2>=1.31110988104
<n_d^2>=1.31110988104
Memory used: 97 MB  Time elapsed: 2225 s

 ***** [NRG] Iteration 17/67 (scale 0.0039849197132687088) ***** 
[1] xi(16)=0.999997469385 zeta(17)=0
Stats: nr=43 min=1 max=1750 sum3=11608477530
Safeguard: keep additional 2 states
Emax=10.0002829631 nrkept=4085 nrkeptmult=10697 
Total energy=-3.8785084568  Egs=-1.35544839469
Kept: 4085 out of 17910, ratio=0.228
<SaSd>=0.00075238595987
<n_a^2>=1.25516145966
<n_d^2>=1.25516145966
Memory used: 84 MB  Time elapsed: 2318 s

 ***** [NRG] Iteration 18/67 (scale 0.0028177637517362566) ***** 
[1] xi(17)=0.999998734524 zeta(18)=0
Stats: nr=50 min=1 max=1656 sum3=11924320549
Safeguard: keep additional 2 states
Emax=10.0009959858 nrkept=3608 nrkeptmult=9514 
Total energy=-3.88134416274  Egs=-1.00636752859
Kept: 3608 out of 15485, ratio=0.233
<SaSd>=0.00186355428013
<n_a^2>=1.19662708817
<n_d^2>=1.19662708817
Memory used: 73 MB  Time elapsed: 2388 s

 ***** [NRG] Iteration 19/67 (scale 0.0019924598566343544) ***** 
[1] xi(18)=0.99999936722 zeta(19)=0
Stats: nr=43 min=1 max=1373 sum3=11461927684
Safeguard: keep additional 1 states
Emax=10.0006055673 nrkept=2992 nrkeptmult=7984 
Total energy=-3.88400679197  Egs=-1.3363527598
Kept: 2992 out of 13690, ratio=0.219
<SaSd>=0.00401769569988
<n_a^2>=1.14483058942
<n_d^2>=1.14483058942
Memory used: 61 MB  Time elapsed: 2437 s

 ***** [NRG] Iteration 20/67 (scale 0.0014088818758681283) ***** 
[1] xi(19)=0.999999683599 zeta(20)=0
Stats: nr=48 min=1 max=1223 sum3=6534353258
Safeguard: keep additional 2 states
Emax=10.0003300641 nrkept=2473 nrkeptmult=6776 
Total energy=-3.88544606078  Egs=-1.02156812058
Kept: 2473 out of 11360, ratio=0.218
<SaSd>=0.00790490773322
<n_a^2>=1.10859932672
<n_d^2>=1.10859932672
Memory used: 51 MB  Time elapsed: 2465 s

 ***** [NRG] Iteration 21/67 (scale 0.0009962299283171772) ***** 
[1] xi(20)=0.999999841797 zeta(21)=0
Stats: nr=41 min=1 max=981 sum3=3812080826
Safeguard: keep additional 1 states
Emax=10.0004148587 nrkept=1981 nrkeptmult=5447 
Total energy=-3.88675644799  Egs=-1.31534615855
Kept: 1981 out of 9422, ratio=0.21
<SaSd>=0.0137706415609
<n_a^2>=1.08952553395
<n_d^2>=1.08952553395
Memory used: 43 MB  Time elapsed: 2482 s

 ***** [NRG] Iteration 22/67 (scale 0.00070444093793406415) ***** 
[1] xi(21)=0.999999920898 zeta(22)=0
Stats: nr=46 min=1 max=818 sum3=1995059482
Emax=10.0005434666 nrkept=1794 nrkeptmult=5010 
Total energy=-3.88748919826  Egs=-1.04018695102
Kept: 1794 out of 7540, ratio=0.238
<SaSd>=0.0220623115708
<n_a^2>=1.08207586737
<n_d^2>=1.08207586737
Memory used: 40 MB  Time elapsed: 2492 s

 ***** [NRG] Iteration 23/67 (scale 0.0004981149641585886) ***** 
[1] xi(22)=0.999999960449 zeta(23)=0
Stats: nr=41 min=1 max=733 sum3=1504612438
Safeguard: keep additional 2 states
Emax=10.0025342056 nrkept=1559 nrkeptmult=4417 
Total energy=-3.88813229694  Egs=-1.29106476476
Kept: 1559 out of 6856, ratio=0.227
<SaSd>=0.0336279493207
<n_a^2>=1.07980810274
<n_d^2>=1.07980810274
Memory used: 37 MB  Time elapsed: 2500 s

 ***** [NRG] Iteration 24/67 (scale 0.00035222046896703208) ***** 
[1] xi(23)=0.999999980224 zeta(24)=0
Stats: nr=46 min=1 max=628 sum3=946151505
Safeguard: keep additional 1 states
Emax=10.0080205804 nrkept=1481 nrkeptmult=4192 
Total energy=-3.88850666485  Egs=-1.06287948034
Kept: 1481 out of 5955, ratio=0.249
<SaSd>=0.0491600044482
<n_a^2>=1.07934224443
<n_d^2>=1.07934224443
Memory used: 37 MB  Time elapsed: 2504 s

 ***** [NRG] Iteration 25/67 (scale 0.0002490574820792943) ***** 
[1] xi(24)=0.999999990112 zeta(25)=0
Stats: nr=41 min=1 max=608 sum3=846966154
Emax=10.0011792405 nrkept=1455 nrkeptmult=4181 
Total energy=-3.8888212263  Egs=-1.26300743308
Kept: 1455 out of 5668, ratio=0.257
<SaSd>=0.0692483878382
<n_a^2>=1.07952252729
<n_d^2>=1.07952252729
Memory used: 38 MB  Time elapsed: 2509 s

 ***** [NRG] Iteration 26/67 (scale 0.00017611023448351604) ***** 
[1] xi(25)=0.999999995056 zeta(26)=0
Stats: nr=46 min=1 max=586 sum3=772041240
Safeguard: keep additional 1 states
Emax=10.0033228665 nrkept=1427 nrkeptmult=4034 
Total energy=-3.88901317154  Egs=-1.08991532454
Kept: 1427 out of 5574, ratio=0.256
<SaSd>=0.0944356568339
<n_a^2>=1.0800072665
<n_d^2>=1.0800072665
Memory used: 37 MB  Time elapsed: 2514 s

 ***** [NRG] Iteration 27/67 (scale 0.00012452874103964715) ***** 
[1] xi(26)=0.999999997528 zeta(27)=0
Stats: nr=41 min=1 max=584 sum3=755728920
Safeguard: keep additional 2 states
Emax=10.0115790441 nrkept=1350 nrkeptmult=3896 
Total energy=-3.88916646587  Egs=-1.23099554663
Kept: 1350 out of 5460, ratio=0.247
<SaSd>=0.123629801514
<n_a^2>=1.08063715078
<n_d^2>=1.08063715078
Memory used: 37 MB  Time elapsed: 2517 s

 ***** [NRG] Iteration 28/67 (scale 8.8055117241758019e-05) ***** 
[1] xi(27)=0.999999998764 zeta(28)=0
Stats: nr=44 min=1 max=545 sum3=626984916
Safeguard: keep additional 2 states
Emax=10.0011814442 nrkept=1268 nrkeptmult=3656 
Total energy=-3.88926522207  Egs=-1.1215270569
Kept: 1268 out of 5172, ratio=0.245
<SaSd>=0.15391704378
<n_a^2>=1.08131490856
<n_d^2>=1.08131490856
Memory used: 38 MB  Time elapsed: 2520 s

 ***** [NRG] Iteration 29/67 (scale 6.2264370519823575e-05) ***** 
[1] xi(28)=0.999999999382 zeta(29)=0
Stats: nr=39 min=1 max=520 sum3=525164968
Safeguard: keep additional 2 states
Emax=10.00993117 nrkept=1216 nrkeptmult=3568 
Total energy=-3.88933963098  Egs=-1.19504809957
Kept: 1216 out of 4864, ratio=0.25
<SaSd>=0.18008099269
<n_a^2>=1.08193276515
<n_d^2>=1.08193276515
Memory used: 37 MB  Time elapsed: 2523 s

 ***** [NRG] Iteration 30/67 (scale 4.402755862087901e-05) ***** 
[1] xi(29)=0.999999999691 zeta(30)=0
Stats: nr=44 min=1 max=483 sum3=447801360
Safeguard: keep additional 1 states
Emax=10.0080877726 nrkept=1142 nrkeptmult=3350 
Total energy=-3.88939060868  Egs=-1.15785887923
Kept: 1142 out of 4668, ratio=0.245
<SaSd>=0.197982140522
<n_a^2>=1.08241503257
<n_d^2>=1.08241503257
Memory used: 38 MB  Time elapsed: 2526 s

 ***** [NRG] Iteration 31/67 (scale 3.1132185259911788e-05) ***** 
[1] xi(30)=0.999999999846 zeta(31)=0
Stats: nr=39 min=1 max=466 sum3=375454518
Safeguard: keep additional 2 states
Emax=10.0110865449 nrkept=1089 nrkeptmult=3219 
Total energy=-3.88942658229  Egs=-1.15551170301
Kept: 1089 out of 4386, ratio=0.248
<SaSd>=0.206489267056
<n_a^2>=1.08273261537
<n_d^2>=1.08273261537
Memory used: 37 MB  Time elapsed: 2528 s

 ***** [NRG] Iteration 32/67 (scale 2.2013779310439505e-05) ***** 
[1] xi(31)=0.999999999923 zeta(32)=0
Stats: nr=44 min=1 max=426 sum3=320771086
Safeguard: keep additional 1 states
Emax=10.0164688472 nrkept=1026 nrkeptmult=3078 
Total energy=-3.88945297098  Egs=-1.19873518631
Kept: 1026 out of 4186, ratio=0.245
<SaSd>=0.209027017793
<n_a^2>=1.08293221381
<n_d^2>=1.08293221381
Memory used: 38 MB  Time elapsed: 2529 s

 ***** [NRG] Iteration 33/67 (scale 1.5566092629955894e-05) ***** 
[1] xi(32)=0.999999999961 zeta(33)=0
Stats: nr=39 min=1 max=416 sum3=264153560
Emax=10.0333869149 nrkept=1008 nrkeptmult=3028 
Total energy=-3.88947030029  Egs=-1.11327311539
Kept: 1008 out of 3950, ratio=0.255
<SaSd>=0.209348443106
<n_a^2>=1.08307090628
<n_d^2>=1.08307090628
Memory used: 39 MB  Time elapsed: 2531 s

 ***** [NRG] Iteration 34/67 (scale 1.1006889655219752e-05) ***** 
[1] xi(33)=0.999999999981 zeta(34)=0
Stats: nr=44 min=1 max=382 sum3=250112229
Emax=10.0003428504 nrkept=969 nrkeptmult=2918 
Total energy=-3.88948398585  Egs=-1.24336259487
Kept: 969 out of 3885, ratio=0.249
<SaSd>=0.209306819058
<n_a^2>=1.08317952082
<n_d^2>=1.08317952082
Memory used: 39 MB  Time elapsed: 2532 s

 ***** [NRG] Iteration 35/67 (scale 7.7830463149779469e-06) ***** 
[1] xi(34)=0.99999999999 zeta(35)=0
Stats: nr=39 min=1 max=395 sum3=222465368
Safeguard: keep additional 3 states
Emax=10.0065289016 nrkept=1002 nrkeptmult=3000 
Total energy=-3.88949231298  Egs=-1.06990682422
Kept: 1002 out of 3734, ratio=0.268
<SaSd>=0.20926488491
<n_a^2>=1.08326690017
<n_d^2>=1.08326690017
Memory used: 38 MB  Time elapsed: 2533 s

 ***** [NRG] Iteration 36/67 (scale 5.5034448276098762e-06) ***** 
[1] xi(35)=0.999999999995 zeta(36)=0
Stats: nr=44 min=1 max=378 sum3=243490692
Safeguard: keep additional 2 states
Emax=10.004953535 nrkept=951 nrkeptmult=2834 
Total energy=-3.88949941297  Egs=-1.29009909167
Kept: 951 out of 3858, ratio=0.247
<SaSd>=0.209231476232
<n_a^2>=1.08333651514
<n_d^2>=1.08333651514
Memory used: 39 MB  Time elapsed: 2535 s

 ***** [NRG] Iteration 37/67 (scale 3.8915231574889735e-06) ***** 
[1] xi(36)=0.999999999998 zeta(37)=0
Stats: nr=38 min=1 max=388 sum3=213275672
Safeguard: keep additional 3 states
Emax=10.0141638415 nrkept=976 nrkeptmult=2866 
Total energy=-3.88950341182  Egs=-1.02758080714
Kept: 976 out of 3662, ratio=0.267
<SaSd>=0.209205400377
<n_a^2>=1.08339085059
<n_d^2>=1.08339085059
Memory used: 39 MB  Time elapsed: 2536 s

 ***** [NRG] Iteration 38/67 (scale 2.7517224138049381e-06) ***** 
[1] xi(37)=0.999999999999 zeta(38)=0
Stats: nr=42 min=1 max=379 sum3=232528533
Safeguard: keep additional 5 states
Emax=10.0027051975 nrkept=925 nrkeptmult=2704 
Total energy=-3.88950708952  Egs=-1.33650726649
Kept: 925 out of 3753, ratio=0.246
<SaSd>=0.209185376244
<n_a^2>=1.08343257589
<n_d^2>=1.08343257589
Memory used: 40 MB  Time elapsed: 2537 s

 ***** [NRG] Iteration 39/67 (scale 1.9457615787444867e-06) ***** 
[1] xi(38)=0.999999999999 zeta(39)=0
Stats: nr=40 min=1 max=379 sum3=198752862
Safeguard: keep additional 3 states
Emax=10.0012349076 nrkept=964 nrkeptmult=2814 
Total energy=-3.88950901312  Egs=-0.988607925141
Kept: 964 out of 3552, ratio=0.271
<SaSd>=0.209170393316
<n_a^2>=1.08346379661
<n_d^2>=1.08346379661
Memory used: 40 MB  Time elapsed: 2538 s

 ***** [NRG] Iteration 40/67 (scale 1.3758612069024691e-06) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=42 min=1 max=379 sum3=226139205
Safeguard: keep additional 1 states
Emax=10.009467048 nrkept=897 nrkeptmult=2606 
Total energy=-3.88951091161  Egs=-1.3798577964
Kept: 897 out of 3705, ratio=0.242
<SaSd>=0.209159383582
<n_a^2>=1.0834867382
<n_d^2>=1.0834867382
Memory used: 40 MB  Time elapsed: 2540 s

 ***** [NRG] Iteration 41/67 (scale 9.7288078937224336e-07) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=40 min=1 max=371 sum3=185209184
Safeguard: keep additional 1 states
Emax=10.0042973432 nrkept=934 nrkeptmult=2696 
Total energy=-3.88951184051  Egs=-0.954790068014
Kept: 934 out of 3446, ratio=0.271
<SaSd>=0.209151562664
<n_a^2>=1.08350303509
<n_d^2>=1.08350303509
Memory used: 40 MB  Time elapsed: 2541 s

 ***** [NRG] Iteration 42/67 (scale 6.8793060345123453e-07) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=42 min=1 max=376 sum3=210859581
Safeguard: keep additional 3 states
Emax=10.0027847594 nrkept=870 nrkeptmult=2480 
Total energy=-3.88951281594  Egs=-1.41792705766
Kept: 870 out of 3585, ratio=0.243
<SaSd>=0.209146115302
<n_a^2>=1.08351438607
<n_d^2>=1.08351438607
Memory used: 40 MB  Time elapsed: 2542 s

 ***** [NRG] Iteration 43/67 (scale 4.8644039468612168e-07) ***** 
[1] xi(42)=1 zeta(43)=0
Stats: nr=40 min=1 max=361 sum3=171343234
Safeguard: keep additional 2 states
Emax=10.0165118452 nrkept=916 nrkeptmult=2596 
Total energy=-3.88951326685  Egs=-0.926949057124
Kept: 916 out of 3334, ratio=0.275
<SaSd>=0.209142459995
<n_a^2>=1.08352200284
<n_d^2>=1.08352200284
Memory used: 40 MB  Time elapsed: 2543 s

 ***** [NRG] Iteration 44/67 (scale 3.4396530172561726e-07) ***** 
[1] xi(43)=1 zeta(44)=0
Stats: nr=42 min=1 max=373 sum3=199243573
Emax=10.018010445 nrkept=853 nrkeptmult=2416 
Total energy=-3.88951376546  Egs=-1.44960024382
Kept: 853 out of 3505, ratio=0.243
<SaSd>=0.209140049943
<n_a^2>=1.08352702481
<n_d^2>=1.08352702481
Memory used: 40 MB  Time elapsed: 2544 s

 ***** [NRG] Iteration 45/67 (scale 2.4322019734306084e-07) ***** 
[1] xi(44)=1 zeta(45)=0
Stats: nr=40 min=1 max=355 sum3=163654654
Safeguard: keep additional 1 states
Emax=10.0127506569 nrkept=895 nrkeptmult=2491 
Total energy=-3.88951398555  Egs=-0.904920787358
Kept: 895 out of 3268, ratio=0.274
<SaSd>=0.209138513439
<n_a^2>=1.08353022651
<n_d^2>=1.08353022651
Memory used: 40 MB  Time elapsed: 2545 s

 ***** [NRG] Iteration 46/67 (scale 1.7198265086280863e-07) ***** 
[1] xi(45)=1 zeta(46)=0
Stats: nr=42 min=1 max=373 sum3=192533869
Safeguard: keep additional 3 states
Emax=10.014325589 nrkept=833 nrkeptmult=2326 
Total energy=-3.88951423922  Egs=-1.47491836411
Kept: 833 out of 3421, ratio=0.243
<SaSd>=0.209137549513
<n_a^2>=1.0835322351
<n_d^2>=1.0835322351
Memory used: 40 MB  Time elapsed: 2546 s

 ***** [NRG] Iteration 47/67 (scale 1.2161009867153042e-07) ***** 
[1] xi(46)=1 zeta(47)=0
Stats: nr=40 min=1 max=349 sum3=156303284
Emax=10.0203118416 nrkept=890 nrkeptmult=2434 
Total energy=-3.88951434719  Egs=-0.887907721147
Kept: 890 out of 3188, ratio=0.279
<SaSd>=0.209136961031
<n_a^2>=1.08353346136
<n_d^2>=1.08353346136
Memory used: 41 MB  Time elapsed: 2547 s

 ***** [NRG] Iteration 48/67 (scale 8.5991325431404316e-08) ***** 
[1] xi(47)=1 zeta(48)=0
Stats: nr=40 min=1 max=374 sum3=191756657
Safeguard: keep additional 1 states
Emax=10.0077485995 nrkept=812 nrkeptmult=2246 
Total energy=-3.88951447572  Egs=-1.49468306489
Kept: 812 out of 3395, ratio=0.239
<SaSd>=0.209136604761
<n_a^2>=1.08353420374
<n_d^2>=1.08353420374
Memory used: 41 MB  Time elapsed: 2548 s

 ***** [NRG] Iteration 49/67 (scale 6.080504933576521e-08) ***** 
[1] xi(48)=1 zeta(49)=0
Stats: nr=40 min=1 max=338 sum3=145360274
Emax=10.009053317 nrkept=893 nrkeptmult=2437 
Total energy=-3.88951452892  Egs=-0.874890559629
Kept: 893 out of 3104, ratio=0.288
<SaSd>=0.209136394005
<n_a^2>=1.0835346429
<n_d^2>=1.0835346429
Memory used: 41 MB  Time elapsed: 2549 s

 ***** [NRG] Iteration 50/67 (scale 4.2995662715702158e-08) ***** 
[1] xi(49)=1 zeta(50)=0
Stats: nr=40 min=1 max=375 sum3=193041116
Safeguard: keep additional 1 states
Emax=10.0363299468 nrkept=808 nrkeptmult=2230 
Total energy=-3.88951459384  Egs=-1.50998037868
Kept: 808 out of 3404, ratio=0.237
<SaSd>=0.209136268953
<n_a^2>=1.08353490348
<n_d^2>=1.08353490348
Memory used: 41 MB  Time elapsed: 2550 s

 ***** [NRG] Iteration 51/67 (scale 3.0402524667882605e-08) ***** 
[1] xi(50)=1 zeta(51)=0
Stats: nr=40 min=1 max=336 sum3=143251600
Emax=10.0523189743 nrkept=901 nrkeptmult=2449 
Total energy=-3.88951462014  Egs=-0.864905245754
Kept: 901 out of 3088, ratio=0.292
<SaSd>=0.209136196444
<n_a^2>=1.08353505457
<n_d^2>=1.08353505457
Memory used: 41 MB  Time elapsed: 2551 s

 ***** [NRG] Iteration 52/67 (scale 2.1497831357851079e-08) ***** 
[1] xi(51)=1 zeta(52)=0
Stats: nr=40 min=1 max=379 sum3=197528124
Safeguard: keep additional 1 states
Emax=10.0266138024 nrkept=796 nrkeptmult=2170 
Total energy=-3.88951465286  Egs=-1.52185756069
Kept: 796 out of 3432, ratio=0.232
<SaSd>=0.209136153818
<n_a^2>=1.08353514339
<n_d^2>=1.08353514339
Memory used: 42 MB  Time elapsed: 2552 s

 ***** [NRG] Iteration 53/67 (scale 1.5201262333941303e-08) ***** 
[1] xi(52)=1 zeta(53)=0
Stats: nr=40 min=1 max=334 sum3=140598232
Emax=10.0334966169 nrkept=901 nrkeptmult=2449 
Total energy=-3.88951466589  Egs=-0.857161940723
Kept: 901 out of 3040, ratio=0.296
<SaSd>=0.209136129503
<n_a^2>=1.08353519406
<n_d^2>=1.08353519406
Memory used: 41 MB  Time elapsed: 2553 s

 ***** [NRG] Iteration 54/67 (scale 1.0748915678925539e-08) ***** 
[1] xi(53)=1 zeta(54)=0
Stats: nr=40 min=1 max=379 sum3=197528124
Safeguard: keep additional 1 states
Emax=10.036596885 nrkept=796 nrkeptmult=2166 
Total energy=-3.88951468234  Egs=-1.53118174048
Kept: 796 out of 3432, ratio=0.232
<SaSd>=0.209136115201
<n_a^2>=1.08353522386
<n_d^2>=1.08353522386
Memory used: 42 MB  Time elapsed: 2554 s

 ***** [NRG] Iteration 55/67 (scale 7.6006311669706513e-09) ***** 
[1] xi(54)=1 zeta(55)=0
Stats: nr=40 min=1 max=336 sum3=141913786
Emax=10.0162726311 nrkept=901 nrkeptmult=2449 
Total energy=-3.88951468881  Egs=-0.851062009931
Kept: 901 out of 3040, ratio=0.296
<SaSd>=0.209136107086
<n_a^2>=1.08353524077
<n_d^2>=1.08353524077
Memory used: 41 MB  Time elapsed: 2555 s

 ***** [NRG] Iteration 56/67 (scale 5.3744578394627697e-09) ***** 
[1] xi(55)=1 zeta(56)=0
Stats: nr=40 min=1 max=379 sum3=197528124
Safeguard: keep additional 1 states
Emax=10.0308782171 nrkept=802 nrkeptmult=2172 
Total energy=-3.88951469708  Egs=-1.53861565116
Kept: 802 out of 3432, ratio=0.234
<SaSd>=0.209136102297
<n_a^2>=1.08353525075
<n_d^2>=1.08353525075
Memory used: 42 MB  Time elapsed: 2556 s

 ***** [NRG] Iteration 57/67 (scale 3.8003155834853256e-09) ***** 
[1] xi(56)=1 zeta(57)=0
Stats: nr=40 min=1 max=336 sum3=143905192
Safeguard: keep additional 1 states
Emax=10.0101187169 nrkept=903 nrkeptmult=2453 
Total energy=-3.8895147003  Egs=-0.846170181211
Kept: 903 out of 3058, ratio=0.295
<SaSd>=0.209136099581
<n_a^2>=1.08353525641
<n_d^2>=1.08353525641
Memory used: 42 MB  Time elapsed: 2557 s

 ***** [NRG] Iteration 58/67 (scale 2.6872289197313849e-09) ***** 
[1] xi(57)=1 zeta(58)=0
Stats: nr=40 min=1 max=381 sum3=199615376
Safeguard: keep additional 1 states
Emax=10.0183633622 nrkept=804 nrkeptmult=2174 
Total energy=-3.88951470445  Egs=-1.54464511911
Kept: 804 out of 3440, ratio=0.234
<SaSd>=0.209136097967
<n_a^2>=1.08353525977
<n_d^2>=1.08353525977
Memory used: 42 MB  Time elapsed: 2558 s

 ***** [NRG] Iteration 59/67 (scale 1.9001577917426628e-09) ***** 
[1] xi(58)=1 zeta(59)=0
Stats: nr=40 min=1 max=336 sum3=144834322
Emax=10.0110856975 nrkept=916 nrkeptmult=2478 
Total energy=-3.88951470605  Egs=-0.84217591977
Kept: 916 out of 3064, ratio=0.299
<SaSd>=0.209136097048
<n_a^2>=1.08353526169
<n_d^2>=1.08353526169
Memory used: 42 MB  Time elapsed: 2559 s

 ***** [NRG] Iteration 60/67 (scale 1.3436144598656924e-09) ***** 
[1] xi(59)=1 zeta(60)=0
Stats: nr=40 min=1 max=388 sum3=208879262
Safeguard: keep additional 1 states
Emax=10.0061306544 nrkept=804 nrkeptmult=2174 
Total energy=-3.88951470813  Egs=-1.54961972948
Kept: 804 out of 3488, ratio=0.231
<SaSd>=0.209136096498
<n_a^2>=1.08353526283
<n_d^2>=1.08353526283
Memory used: 43 MB  Time elapsed: 2560 s

 ***** [NRG] Iteration 61/67 (scale 9.5007889587133141e-10) ***** 
[1] xi(60)=1 zeta(61)=0
Stats: nr=40 min=1 max=336 sum3=144834322
Emax=10.014723154 nrkept=923 nrkeptmult=2491 
Total energy=-3.88951470893  Egs=-0.838858653104
Kept: 923 out of 3064, ratio=0.301
<SaSd>=0.209136096183
<n_a^2>=1.08353526349
<n_d^2>=1.08353526349
Memory used: 42 MB  Time elapsed: 2561 s

 ***** [NRG] Iteration 62/67 (scale 6.7180722993284622e-10) ***** 
[1] xi(61)=1 zeta(62)=0
Stats: nr=40 min=1 max=393 sum3=215138583
Emax=10.0389649719 nrkept=805 nrkeptmult=2176 
Total energy=-3.88951470997  Egs=-1.55379005761
Kept: 805 out of 3513, ratio=0.229
<SaSd>=0.209136095992
<n_a^2>=1.08353526389
<n_d^2>=1.08353526389
Memory used: 43 MB  Time elapsed: 2562 s

 ***** [NRG] Iteration 63/67 (scale 4.750394479356657e-10) ***** 
[1] xi(62)=1 zeta(63)=0
Stats: nr=40 min=1 max=337 sum3=145863974
Safeguard: keep additional 1 states
Emax=10.0001834972 nrkept=924 nrkeptmult=2492 
Total energy=-3.88951471037  Egs=-0.83606089792
Kept: 924 out of 3068, ratio=0.301
<SaSd>=0.209136095882
<n_a^2>=1.08353526412
<n_d^2>=1.08353526412
Memory used: 42 MB  Time elapsed: 2563 s

 ***** [NRG] Iteration 64/67 (scale 3.3590361496642311e-10) ***** 
[1] xi(63)=1 zeta(64)=0
Stats: nr=40 min=1 max=394 sum3=216028273
Emax=10.0312256197 nrkept=805 nrkeptmult=2176 
Total energy=-3.88951471089  Egs=-1.55733671256
Kept: 805 out of 3517, ratio=0.229
<SaSd>=0.209136095814
<n_a^2>=1.08353526426
<n_d^2>=1.08353526426
Memory used: 43 MB  Time elapsed: 2564 s

 ***** [NRG] Iteration 65/67 (scale 2.3751972396783285e-10) ***** 
[1] xi(64)=1 zeta(65)=0
Stats: nr=40 min=1 max=337 sum3=145863974
Safeguard: keep additional 2 states
Emax=10.0074717239 nrkept=932 nrkeptmult=2512 
Total energy=-3.88951471109  Egs=-0.833668942945
Kept: 932 out of 3068, ratio=0.304
<SaSd>=0.209136095775
<n_a^2>=1.08353526434
<n_d^2>=1.08353526434
Memory used: 42 MB  Time elapsed: 2565 s

 ***** [NRG] Iteration 66/67 (scale 1.6795180748321155e-10) ***** 
[1] xi(65)=1 zeta(66)=0
Stats: nr=40 min=1 max=398 sum3=222658263
Emax=10.0242025976 nrkept=805 nrkeptmult=2176 
Total energy=-3.88951471135  Egs=-1.56039137266
Kept: 805 out of 3549, ratio=0.227
<SaSd>=0.209136095751
<n_a^2>=1.08353526439
<n_d^2>=1.08353526439
Memory used: 43 MB  Time elapsed: 2566 s

 ***** [NRG] Iteration 67/67 (scale 1.1875986198391643e-10) ***** 
[1] xi(66)=1 zeta(67)=0
Stats: nr=40 min=1 max=337 sum3=145863974
Safeguard: keep additional 4 states
Emax=10.1033916348 nrkept=937 nrkeptmult=2517 
Total energy=-3.88951471145  Egs=-0.831599455334
Kept: 937 out of 3068, ratio=0.305
<SaSd>=0.209136095737
<n_a^2>=1.08353526442
<n_d^2>=1.08353526442
Memory used: 42 MB  Time elapsed: 2567 s

** Iteration completed.


Total energy:               -3.8895147114511
TKW=3.33130842244e-07
Grand canonical Z=              26.19759207209745
Zft (last shell)=               26.16431133025298

Memory usage report
===================
  after diag: 568880 kB
after recalc: 328116 kB
 after split: 921936 kB
  after trim: 184364 kB
 after trunc: 569620 kB
      diag_h: 572136 kB

Peak usage: 921936 kB

Timing report
=============
         All: 2567 s
        diag: 2449 s
    recalc s: 34 s
       Other: 83 s
