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=42
TEMPERATURE T=1e-10
Nmax=42
Tmin=1e-10
Ninit=0
Lambda=3
strategy=kept
lastall=1
lastalloverride=0
z=1
discretization=Z
band=flat
betabar=1
T=1e-10
T_ratio=1
Tmin_ratio=-999
keep=10000
keepenergy=10
keepmin=0
keepblock=0
keepzrescale=0
log=
logall=0
logenumber=10
smooth=new
alpha=0.5
omega0=1e-99
omega0_ratio=1
brcut=0
loggauss_b=0.3
gauss_b=0.1
limitLL=0
eta=0
finite=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=0.1
broaden_min=1e-08
broaden_min_ratio=3
broaden_ratio=1.02
dumpenergies=0
saveall=0
dumpannotated=0
dumpscaled=1
dumpabs=0
dumpprecision=8
dumpgroups=1
grouptol=1e-06
dumpdiagonal=0
ops=A_d SZd
specs=
specd=A_d-A_d-u A_d-A_d-d
spect=
specb=
reim=0
trsq=
specgt=
gtp=0.7
speci1t=
speci2t=
speckubo=
kubop=0.7
kubor=4
kubob=2
cond=
specchit=
chitp=1
chitp_ratio=-999
dm=1
T0opt=-1
dmnrg=0
cfs=0
fdm=1
fdmexpv=0
fdmexpvn=0
mats=0
nromegan=100
width_td=12
width_custom=12
prec_td=6
prec_custom=6
prec_xy=6
safeguard=1e-05
safeguardmax=200
fixeps=1e-15
diag=dsyev
diagratio=1
dsyevrlimit=100
restart=1
restartfactor=2
checkdiag=0
checkrho=0
calc0=1
tdht=0
spin=2
tri=old
preccpp=2000
checksumrules=0
showEs=0
globalB=0
globalBx=0
polarized=0
pol2x2=0
mpi=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: QSZ
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 	MAGFIELD	 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 	 extraB=	0.0001
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	3.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	1.
>>  ops=	A_d SZd
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=9
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||SZd||> (s)
Reading <||||> (z)

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

Number of invariant subspaces: 9
(-2 1) 1 states
0.0504489594881 
(-1 0) 2 states
0.0239209588286 0.0721722580855 
(-1 2) 2 states
0.023973220682 0.0722151388472 
(0 -1) 1 states
0.045644257426 
(0 1) 4 states
0 0.0456918334039 0.0504489594881 0.0961407835514 
(0 3) 1 states
0.0457394000411 
(1 0) 2 states
0.0239209588286 0.0721722580855 
(1 2) 2 states
0.023973220682 0.0722151388472 
(2 1) 1 states
0.0504489594881 
Number of states (multiplicity taken into account): 16

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

Computing the following spectra:

Before NRG iteration (N=-1)
<SZd>=-1.19199671072e-05

 ***** [NRG] Iteration 1/42 (scale 0.60682615108455829) ***** 
[1] xi(0)=0.832050294338 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=3.34063535785 nrkept=64 nrkeptmult=64 
Total energy=-1.01111738712  Egs=-1.58709841131
Kept: 64 out of 64, ratio=1
<SZd>=-2.06782478412e-05
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 2/42 (scale 0.35035124167997428) ***** 
[1] xi(1)=0.907691230187 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=6.97117080502 nrkept=256 nrkeptmult=256 
Total energy=-1.21867918236  Egs=-0.59243915975
Kept: 256 out of 256, ratio=1
<SZd>=-3.59068266186e-05
Memory used: 4 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 3/42 (scale 0.20227538369485276) ***** 
[1] xi(2)=0.965171190957 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Safeguard: keep additional 1 states
Emax=10.6127113378 nrkept=786 nrkeptmult=786 
Total energy=-1.53995521835  Egs=-1.58831010536
Kept: 786 out of 1024, ratio=0.768
<SZd>=-6.36229015219e-05
Memory used: 7 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 4/42 (scale 0.11678374722665809) ***** 
[1] xi(3)=0.987905295336 zeta(4)=0
Stats: nr=49 min=1 max=300 sum3=99197412
Safeguard: keep additional 15 states
Emax=10.2988760781 nrkept=976 nrkeptmult=976 
Total energy=-1.62444624004  Egs=-0.723482707975
Kept: 976 out of 3144, ratio=0.31
<SZd>=-0.000108965954893
Memory used: 13 MB  Time elapsed: 0 s

 ***** [NRG] Iteration 5/42 (scale 0.067425127898284254) ***** 
[1] xi(4)=0.995912883703 zeta(5)=0
Stats: nr=60 min=1 max=308 sum3=168775504
Safeguard: keep additional 1 states
Emax=10.6236055056 nrkept=866 nrkeptmult=866 
Total energy=-1.72842051254  Egs=-1.54207008191
Kept: 866 out of 3904, ratio=0.222
<SZd>=-0.000192073916322
Memory used: 17 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 6/42 (scale 0.038927915742219356) ***** 
[1] xi(5)=0.998631389726 zeta(6)=0
Stats: nr=49 min=1 max=332 sum3=134201072
Safeguard: keep additional 3 states
Emax=10.322446913 nrkept=964 nrkeptmult=964 
Total energy=-1.75904030675  Egs=-0.786576769634
Kept: 964 out of 3464, ratio=0.278
<SZd>=-0.00033717954439
Memory used: 19 MB  Time elapsed: 1 s

 ***** [NRG] Iteration 7/42 (scale 0.022475042632761417) ***** 
[1] xi(6)=0.999543100929 zeta(7)=0
Stats: nr=60 min=1 max=304 sum3=160889638
Safeguard: keep additional 3 states
Emax=10.0138692035 nrkept=804 nrkeptmult=804 
Total energy=-1.79198654456  Egs=-1.46590323952
Kept: 804 out of 3856, ratio=0.209
<SZd>=-0.000602661379269
Memory used: 19 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 8/42 (scale 0.01297597191407312) ***** 
[1] xi(7)=0.999847622921 zeta(8)=0
Stats: nr=49 min=1 max=310 sum3=108932430
Emax=10.4023070595 nrkept=961 nrkeptmult=961 
Total energy=-1.80293969248  Egs=-0.84411002038
Kept: 961 out of 3216, ratio=0.299
<SZd>=-0.00109539975626
Memory used: 21 MB  Time elapsed: 2 s

 ***** [NRG] Iteration 9/42 (scale 0.0074916808775871389) ***** 
[1] xi(8)=0.999949199038 zeta(9)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Emax=10.0075273714 nrkept=765 nrkeptmult=765 
Total energy=-1.81331632123  Egs=-1.38508686192
Kept: 765 out of 3844, ratio=0.199
<SZd>=-0.00203218233418
Memory used: 22 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 10/42 (scale 0.0043253239713577068) ***** 
[1] xi(9)=0.99998306539 zeta(10)=0
Stats: nr=49 min=1 max=291 sum3=90963678
Safeguard: keep additional 1 states
Emax=10.0035616952 nrkept=810 nrkeptmult=810 
Total energy=-1.81718207525  Egs=-0.893749010814
Kept: 810 out of 3060, ratio=0.265
<SZd>=-0.00381112179452
Memory used: 23 MB  Time elapsed: 3 s

 ***** [NRG] Iteration 11/42 (scale 0.0024972269591957133) ***** 
[1] xi(10)=0.999994355024 zeta(11)=0
Stats: nr=56 min=1 max=266 sum3=105194994
Emax=10.0009295219 nrkept=649 nrkeptmult=649 
Total energy=-1.82045603495  Egs=-1.31103810253
Kept: 649 out of 3240, ratio=0.2
<SZd>=-0.00707375360271
Memory used: 22 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 12/42 (scale 0.0014417746571192356) ***** 
[1] xi(11)=0.999998118329 zeta(12)=0
Stats: nr=49 min=1 max=250 sum3=57164242
Safeguard: keep additional 2 states
Emax=10.0588555575 nrkept=635 nrkeptmult=635 
Total energy=-1.82180897873  Egs=-0.938387825329
Kept: 635 out of 2596, ratio=0.245
<SZd>=-0.0126119759419
Memory used: 23 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 13/42 (scale 0.00083240898639857098) ***** 
[1] xi(12)=0.999999372775 zeta(13)=0
Stats: nr=56 min=1 max=211 sum3=51921140
Safeguard: keep additional 2 states
Emax=10.078518166 nrkept=509 nrkeptmult=509 
Total energy=-1.82284399349  Egs=-1.24339690115
Kept: 509 out of 2540, ratio=0.2
<SZd>=-0.0214745978356
Memory used: 25 MB  Time elapsed: 4 s

 ***** [NRG] Iteration 14/42 (scale 0.00048059155237307853) ***** 
[1] xi(13)=0.999999790925 zeta(14)=0
Stats: nr=47 min=1 max=200 sum3=28797164
Safeguard: keep additional 1 states
Emax=10.0168500688 nrkept=525 nrkeptmult=525 
Total energy=-1.82331632419  Egs=-0.982811074464
Kept: 525 out of 2036, ratio=0.258
<SZd>=-0.0354333621068
Memory used: 26 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 15/42 (scale 0.00027746966213285698) ***** 
[1] xi(14)=0.999999930308 zeta(15)=0
Stats: nr=56 min=1 max=174 sum3=27995250
Safeguard: keep additional 2 states
Emax=10.0154211494 nrkept=411 nrkeptmult=411 
Total energy=-1.82364568258  Egs=-1.18700684485
Kept: 411 out of 2100, ratio=0.196
<SZd>=-0.0570337552519
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 16/42 (scale 0.00016019718412435953) ***** 
[1] xi(15)=0.999999976769 zeta(16)=0
Stats: nr=47 min=1 max=164 sum3=15792954
Safeguard: keep additional 2 states
Emax=10.0306287105 nrkept=469 nrkeptmult=469 
Total energy=-1.82380892654  Egs=-1.01901893676
Kept: 469 out of 1644, ratio=0.285
<SZd>=-0.0891142753455
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 17/42 (scale 9.2489887377619005e-05) ***** 
[1] xi(16)=0.999999992256 zeta(17)=0
Stats: nr=52 min=1 max=162 sum3=21195436
Emax=10.0064958778 nrkept=387 nrkeptmult=387 
Total energy=-1.82391588789  Egs=-1.15646537213
Kept: 387 out of 1876, ratio=0.206
<SZd>=-0.132555130365
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 18/42 (scale 5.33990613747865e-05) ***** 
[1] xi(17)=0.999999997419 zeta(18)=0
Stats: nr=47 min=1 max=154 sum3=13021266
Safeguard: keep additional 1 states
Emax=10.0129995013 nrkept=398 nrkeptmult=398 
Total energy=-1.82397134579  Egs=-1.03855575792
Kept: 398 out of 1548, ratio=0.257
<SZd>=-0.183469996704
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 19/42 (scale 3.0829962459206333e-05) ***** 
[1] xi(18)=0.99999999914 zeta(19)=0
Stats: nr=50 min=1 max=139 sum3=13168796
Emax=10.0490886942 nrkept=332 nrkeptmult=332 
Total energy=-1.8240068562  Egs=-1.15181492367
Kept: 332 out of 1592, ratio=0.209
<SZd>=-0.229215489826
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 20/42 (scale 1.7799687124928835e-05) ***** 
[1] xi(19)=0.999999999713 zeta(20)=0
Stats: nr=46 min=1 max=135 sum3=8685326
Safeguard: keep additional 1 states
Emax=10.0019965573 nrkept=301 nrkeptmult=301 
Total energy=-1.82402545284  Egs=-1.04477332795
Kept: 301 out of 1328, ratio=0.227
<SZd>=-0.258182096723
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 21/42 (scale 1.0276654153068778e-05) ***** 
[1] xi(20)=0.999999999904 zeta(21)=0
Stats: nr=42 min=1 max=123 sum3=7086040
Safeguard: keep additional 1 states
Emax=10.0236810514 nrkept=262 nrkeptmult=262 
Total energy=-1.82403732776  Egs=-1.15552327795
Kept: 262 out of 1204, ratio=0.218
<SZd>=-0.270281798134
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 22/42 (scale 5.9332290416429446e-06) ***** 
[1] xi(21)=0.999999999968 zeta(22)=0
Stats: nr=42 min=1 max=110 sum3=4703098
Safeguard: keep additional 3 states
Emax=10.0334714762 nrkept=251 nrkeptmult=251 
Total energy=-1.82404353269  Egs=-1.04579450438
Kept: 251 out of 1048, ratio=0.24
<SZd>=-0.27437314703
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 23/42 (scale 3.4255513843562592e-06) ***** 
[1] xi(22)=0.999999999989 zeta(23)=0
Stats: nr=42 min=1 max=105 sum3=4212974
Safeguard: keep additional 1 states
Emax=10.0015498739 nrkept=237 nrkeptmult=237 
Total energy=-1.82404750001  Egs=-1.15815424818
Kept: 237 out of 1004, ratio=0.236
<SZd>=-0.275567379095
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 24/42 (scale 1.9777430138809818e-06) ***** 
[1] xi(23)=0.999999999996 zeta(24)=0
Stats: nr=42 min=1 max=95 sum3=3448026
Safeguard: keep additional 1 states
Emax=10.0159734498 nrkept=238 nrkeptmult=238 
Total energy=-1.82404956837  Egs=-1.04581825676
Kept: 238 out of 948, ratio=0.251
<SZd>=-0.276011155265
Memory used: 25 MB  Time elapsed: 5 s

 ***** [NRG] Iteration 25/42 (scale 1.1418504614520865e-06) ***** 
[1] xi(24)=0.999999999999 zeta(25)=0
Stats: nr=41 min=1 max=99 sum3=3626542
Safeguard: keep additional 1 states
Emax=10.0780249476 nrkept=231 nrkeptmult=231 
Total energy=-1.82405089205  Egs=-1.15924111723
Kept: 231 out of 952, ratio=0.243
<SZd>=-0.276122361714
Memory used: 25 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 26/42 (scale 6.592476712936605e-07) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=40 min=1 max=93 sum3=3306120
Safeguard: keep additional 1 states
Emax=10.0546762932 nrkept=232 nrkeptmult=232 
Total energy=-1.82405158148  Egs=-1.04577454825
Kept: 232 out of 924, ratio=0.251
<SZd>=-0.276180025265
Memory used: 25 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 27/42 (scale 3.8061682048402881e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.108858851 nrkept=231 nrkeptmult=231 
Total energy=-1.82405202285  Egs=-1.15962880319
Kept: 231 out of 928, ratio=0.249
<SZd>=-0.27618689783
Memory used: 25 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 28/42 (scale 2.1974922376455349e-07) ***** 
[1] xi(27)=1 zeta(28)=0
Stats: nr=40 min=1 max=93 sum3=3306120
Safeguard: keep additional 1 states
Emax=10.0825890703 nrkept=232 nrkeptmult=232 
Total energy=-1.82405225265  Egs=-1.04575343922
Kept: 232 out of 924, ratio=0.251
<SZd>=-0.276196356374
Memory used: 25 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 29/42 (scale 1.2687227349467628e-07) ***** 
[1] xi(28)=1 zeta(29)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0026229622 nrkept=228 nrkeptmult=228 
Total energy=-1.82405239979  Egs=-1.15976091157
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276195345724
Memory used: 25 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 30/42 (scale 7.3249741254851168e-08) ***** 
[1] xi(29)=1 zeta(30)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0892566706 nrkept=232 nrkeptmult=232 
Total energy=-1.8240524764  Egs=-1.04574565294
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276197434338
Memory used: 25 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 31/42 (scale 4.2290757831558757e-08) ***** 
[1] xi(30)=1 zeta(31)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0069047426 nrkept=228 nrkeptmult=228 
Total energy=-1.82405252544  Egs=-1.15980526
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276196733509
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 32/42 (scale 2.4416580418283724e-08) ***** 
[1] xi(31)=1 zeta(32)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.091362489 nrkept=232 nrkeptmult=232 
Total energy=-1.82405255098  Egs=-1.0457429826
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276197310712
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 33/42 (scale 1.4096919277186252e-08) ***** 
[1] xi(32)=1 zeta(33)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0083430139 nrkept=228 nrkeptmult=228 
Total energy=-1.82405256733  Egs=-1.1598200674
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276197036927
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 34/42 (scale 8.1388601394279074e-09) ***** 
[1] xi(33)=1 zeta(34)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0920558999 nrkept=232 nrkeptmult=232 
Total energy=-1.82405257584  Egs=-1.04574209377
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276197216143
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 35/42 (scale 4.6989730923954177e-09) ***** 
[1] xi(34)=1 zeta(35)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0088258139 nrkept=228 nrkeptmult=228 
Total energy=-1.82405258129  Egs=-1.15982499574
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.276197120464
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 36/42 (scale 2.7129533798093026e-09) ***** 
[1] xi(35)=1 zeta(36)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0922861324 nrkept=232 nrkeptmult=232 
Total energy=-1.82405258413  Egs=-1.04574180729
Kept: 232 out of 912, ratio=0.254
<SZd>=-0.276197178734
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 37/42 (scale 1.5663243641318058e-09) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0089972888 nrkept=228 nrkeptmult=228 
Total energy=-1.82405258594  Egs=-1.1598266275
Kept: 228 out of 928, ratio=0.246
<SZd>=-0.27619714636
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 38/42 (scale 9.0431779326976753e-10) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923734014 nrkept=238 nrkeptmult=238 
Total energy=-1.82405258689  Egs=-1.04574172252
Kept: 238 out of 912, ratio=0.261
<SZd>=-0.276197165595
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 39/42 (scale 5.2210812137726861e-10) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=41 min=1 max=99 sum3=3556156
Safeguard: keep additional 1 states
Emax=10.0089577533 nrkept=228 nrkeptmult=228 
Total energy=-1.82405258749  Egs=-1.15982716004
Kept: 228 out of 952, ratio=0.239
<SZd>=-0.276197154776
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 40/42 (scale 3.0143926442325579e-10) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923886546 nrkept=238 nrkeptmult=238 
Total energy=-1.82405258781  Egs=-1.04574170505
Kept: 238 out of 912, ratio=0.261
<SZd>=-0.276197161166
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 41/42 (scale 1.7403604045908955e-10) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=41 min=1 max=99 sum3=3556156
Safeguard: keep additional 1 states
Emax=10.0089756078 nrkept=228 nrkeptmult=228 
Total energy=-1.82405258801  Egs=-1.15982732612
Kept: 228 out of 952, ratio=0.239
<SZd>=-0.276197157555
Memory used: 26 MB  Time elapsed: 6 s

 ***** [NRG] Iteration 42/42 (scale 1.0047975480775194e-10) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923971214 nrkept=912 nrkeptmult=912 
Total energy=-1.82405258812  Egs=-1.04574171005
Kept: 912 out of 912, ratio=1
<SZd>=-0.276197159818
Memory used: 26 MB  Time elapsed: 6 s

** Iteration completed.


Total energy:             -1.824052588116197
TKW=0.00016087348475
Grand canonical Z=               6.02967533461883
Zft (last shell)=                6.02968091547688
Calculating density matrices.
Calculating full density matrices.

>>  Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2012
SYMMETRY TYPE: QSZ
>>  $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
>>  Using sneg version 	1.233
>>  Model: 	SIAM	 Variant: 	MAGFIELD	 Channels: 	1
>>  U=	0.01	 Gamma=	0.001	 Delta=	0	 t=	0.
>> 	 extraB=	0.0001
>>  Gamma^(1/2)=	0.0252313252202016
>>  Lambda=	3.	 BAND=	flat
>>  DISCRETIZATION=	Z	 z=	1.
>>  ops=	A_d SZd
V7
>>  Number of channels, impurities, subspaces: 
channels=1 impurities=1 subspaces=9
Reading <||||> (e)
Reading <||A_d||> (d)
Reading <||SZd||> (s)
Reading <||||> (z)
Number of invariant subspaces: 9
(-2 1) 1 states
0.0504489594881 
(-1 0) 2 states
0.0239209588286 0.0721722580855 
(-1 2) 2 states
0.023973220682 0.0722151388472 
(0 -1) 1 states
0.045644257426 
(0 1) 4 states
0 0.0456918334039 0.0504489594881 0.0961407835514 
(0 3) 1 states
0.0457394000411 
(1 0) 2 states
0.0239209588286 0.0721722580855 
(1 2) 2 states
0.023973220682 0.0722151388472 
(2 1) 1 states
0.0504489594881 
Number of states (multiplicity taken into account): 16

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

Computing the following spectra:
Spectrum: A_d-A_d-d FDM
Spectrum: A_d-A_d-u FDM

 ***** [DMNRG] Iteration 1/42 (scale 0.60682615108455829) ***** 
[1] xi(0)=0.832050294338 zeta(1)=0
Stats: nr=16 min=1 max=9 sum3=3136
Emax=3.34063535785 nrkept=64 nrkeptmult=64 
Kept: 64 out of 64, ratio=1
Memory used: 30 MB  Time elapsed: 8 s

 ***** [DMNRG] Iteration 2/42 (scale 0.35035124167997428) ***** 
[1] xi(1)=0.907691230187 zeta(2)=0
Stats: nr=25 min=1 max=36 sum3=119716
Emax=6.97117080502 nrkept=256 nrkeptmult=256 
Kept: 256 out of 256, ratio=1
Memory used: 30 MB  Time elapsed: 8 s

 ***** [DMNRG] Iteration 3/42 (scale 0.20227538369485276) ***** 
[1] xi(2)=0.965171190957 zeta(3)=0
Stats: nr=36 min=1 max=100 sum3=5071504
Safeguard: keep additional 1 states
Emax=10.6127113378 nrkept=786 nrkeptmult=786 
Kept: 786 out of 1024, ratio=0.768
Memory used: 30 MB  Time elapsed: 8 s

 ***** [DMNRG] Iteration 4/42 (scale 0.11678374722665809) ***** 
[1] xi(3)=0.987905295336 zeta(4)=0
Stats: nr=49 min=1 max=300 sum3=99197412
Safeguard: keep additional 15 states
Emax=10.2988760781 nrkept=976 nrkeptmult=976 
Kept: 976 out of 3144, ratio=0.31
Memory used: 30 MB  Time elapsed: 8 s

 ***** [DMNRG] Iteration 5/42 (scale 0.067425127898284254) ***** 
[1] xi(4)=0.995912883703 zeta(5)=0
Stats: nr=60 min=1 max=308 sum3=168775504
Safeguard: keep additional 1 states
Emax=10.6236055056 nrkept=866 nrkeptmult=866 
Kept: 866 out of 3904, ratio=0.222
Memory used: 31 MB  Time elapsed: 8 s

 ***** [DMNRG] Iteration 6/42 (scale 0.038927915742219356) ***** 
[1] xi(5)=0.998631389726 zeta(6)=0
Stats: nr=49 min=1 max=332 sum3=134201072
Safeguard: keep additional 3 states
Emax=10.322446913 nrkept=964 nrkeptmult=964 
Kept: 964 out of 3464, ratio=0.278
Memory used: 32 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 7/42 (scale 0.022475042632761417) ***** 
[1] xi(6)=0.999543100929 zeta(7)=0
Stats: nr=60 min=1 max=304 sum3=160889638
Safeguard: keep additional 3 states
Emax=10.0138692035 nrkept=804 nrkeptmult=804 
Kept: 804 out of 3856, ratio=0.209
Memory used: 34 MB  Time elapsed: 9 s

 ***** [DMNRG] Iteration 8/42 (scale 0.01297597191407312) ***** 
[1] xi(7)=0.999847622921 zeta(8)=0
Stats: nr=49 min=1 max=310 sum3=108932430
Emax=10.4023070595 nrkept=961 nrkeptmult=961 
Kept: 961 out of 3216, ratio=0.299
Memory used: 34 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 9/42 (scale 0.0074916808775871389) ***** 
[1] xi(8)=0.999949199038 zeta(9)=0
Stats: nr=60 min=1 max=303 sum3=159977692
Emax=10.0075273714 nrkept=765 nrkeptmult=765 
Kept: 765 out of 3844, ratio=0.199
Memory used: 36 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 10/42 (scale 0.0043253239713577068) ***** 
[1] xi(9)=0.99998306539 zeta(10)=0
Stats: nr=49 min=1 max=291 sum3=90963678
Safeguard: keep additional 1 states
Emax=10.0035616952 nrkept=810 nrkeptmult=810 
Kept: 810 out of 3060, ratio=0.265
Memory used: 34 MB  Time elapsed: 10 s

 ***** [DMNRG] Iteration 11/42 (scale 0.0024972269591957133) ***** 
[1] xi(10)=0.999994355024 zeta(11)=0
Stats: nr=56 min=1 max=266 sum3=105194994
Emax=10.0009295219 nrkept=649 nrkeptmult=649 
Kept: 649 out of 3240, ratio=0.2
Memory used: 35 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 12/42 (scale 0.0014417746571192356) ***** 
[1] xi(11)=0.999998118329 zeta(12)=0
Stats: nr=49 min=1 max=250 sum3=57164242
Safeguard: keep additional 2 states
Emax=10.0588555575 nrkept=635 nrkeptmult=635 
Kept: 635 out of 2596, ratio=0.245
Memory used: 35 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 13/42 (scale 0.00083240898639857098) ***** 
[1] xi(12)=0.999999372775 zeta(13)=0
Stats: nr=56 min=1 max=211 sum3=51921140
Safeguard: keep additional 2 states
Emax=10.078518166 nrkept=509 nrkeptmult=509 
Kept: 509 out of 2540, ratio=0.2
Memory used: 34 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 14/42 (scale 0.00048059155237307853) ***** 
[1] xi(13)=0.999999790925 zeta(14)=0
Stats: nr=47 min=1 max=200 sum3=28797164
Safeguard: keep additional 1 states
Emax=10.0168500688 nrkept=525 nrkeptmult=525 
Kept: 525 out of 2036, ratio=0.258
Memory used: 33 MB  Time elapsed: 11 s

 ***** [DMNRG] Iteration 15/42 (scale 0.00027746966213285698) ***** 
[1] xi(14)=0.999999930308 zeta(15)=0
Stats: nr=56 min=1 max=174 sum3=27995250
Safeguard: keep additional 2 states
Emax=10.0154211494 nrkept=411 nrkeptmult=411 
Kept: 411 out of 2100, ratio=0.196
Memory used: 33 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 16/42 (scale 0.00016019718412435953) ***** 
[1] xi(15)=0.999999976769 zeta(16)=0
Stats: nr=47 min=1 max=164 sum3=15792954
Safeguard: keep additional 2 states
Emax=10.0306287105 nrkept=469 nrkeptmult=469 
Kept: 469 out of 1644, ratio=0.285
Memory used: 33 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 17/42 (scale 9.2489887377619005e-05) ***** 
[1] xi(16)=0.999999992256 zeta(17)=0
Stats: nr=52 min=1 max=162 sum3=21195436
Emax=10.0064958778 nrkept=387 nrkeptmult=387 
Kept: 387 out of 1876, ratio=0.206
Memory used: 34 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 18/42 (scale 5.33990613747865e-05) ***** 
[1] xi(17)=0.999999997419 zeta(18)=0
Stats: nr=47 min=1 max=154 sum3=13021266
Safeguard: keep additional 1 states
Emax=10.0129995013 nrkept=398 nrkeptmult=398 
Kept: 398 out of 1548, ratio=0.257
Memory used: 33 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 19/42 (scale 3.0829962459206333e-05) ***** 
[1] xi(18)=0.99999999914 zeta(19)=0
Stats: nr=50 min=1 max=139 sum3=13168796
Emax=10.0490886942 nrkept=332 nrkeptmult=332 
Kept: 332 out of 1592, ratio=0.209
Memory used: 34 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 20/42 (scale 1.7799687124928835e-05) ***** 
[1] xi(19)=0.999999999713 zeta(20)=0
Stats: nr=46 min=1 max=135 sum3=8685326
Safeguard: keep additional 1 states
Emax=10.0019965573 nrkept=301 nrkeptmult=301 
Kept: 301 out of 1328, ratio=0.227
Memory used: 34 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 21/42 (scale 1.0276654153068778e-05) ***** 
[1] xi(20)=0.999999999904 zeta(21)=0
Stats: nr=42 min=1 max=123 sum3=7086040
Safeguard: keep additional 1 states
Emax=10.0236810514 nrkept=262 nrkeptmult=262 
Kept: 262 out of 1204, ratio=0.218
Memory used: 34 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 22/42 (scale 5.9332290416429446e-06) ***** 
[1] xi(21)=0.999999999968 zeta(22)=0
Stats: nr=42 min=1 max=110 sum3=4703098
Safeguard: keep additional 3 states
Emax=10.0334714762 nrkept=251 nrkeptmult=251 
Kept: 251 out of 1048, ratio=0.24
Memory used: 34 MB  Time elapsed: 12 s

 ***** [DMNRG] Iteration 23/42 (scale 3.4255513843562592e-06) ***** 
[1] xi(22)=0.999999999989 zeta(23)=0
Stats: nr=42 min=1 max=105 sum3=4212974
Safeguard: keep additional 1 states
Emax=10.0015498739 nrkept=237 nrkeptmult=237 
Kept: 237 out of 1004, ratio=0.236
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 24/42 (scale 1.9777430138809818e-06) ***** 
[1] xi(23)=0.999999999996 zeta(24)=0
Stats: nr=42 min=1 max=95 sum3=3448026
Safeguard: keep additional 1 states
Emax=10.0159734498 nrkept=238 nrkeptmult=238 
Kept: 238 out of 948, ratio=0.251
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 25/42 (scale 1.1418504614520865e-06) ***** 
[1] xi(24)=0.999999999999 zeta(25)=0
Stats: nr=41 min=1 max=99 sum3=3626542
Safeguard: keep additional 1 states
Emax=10.0780249476 nrkept=231 nrkeptmult=231 
Kept: 231 out of 952, ratio=0.243
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 26/42 (scale 6.592476712936605e-07) ***** 
[1] xi(25)=1 zeta(26)=0
Stats: nr=40 min=1 max=93 sum3=3306120
Safeguard: keep additional 1 states
Emax=10.0546762932 nrkept=232 nrkeptmult=232 
Kept: 232 out of 924, ratio=0.251
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 27/42 (scale 3.8061682048402881e-07) ***** 
[1] xi(26)=1 zeta(27)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.108858851 nrkept=231 nrkeptmult=231 
Kept: 231 out of 928, ratio=0.249
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 28/42 (scale 2.1974922376455349e-07) ***** 
[1] xi(27)=1 zeta(28)=0
Stats: nr=40 min=1 max=93 sum3=3306120
Safeguard: keep additional 1 states
Emax=10.0825890703 nrkept=232 nrkeptmult=232 
Kept: 232 out of 924, ratio=0.251
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 29/42 (scale 1.2687227349467628e-07) ***** 
[1] xi(28)=1 zeta(29)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0026229622 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 30/42 (scale 7.3249741254851168e-08) ***** 
[1] xi(29)=1 zeta(30)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0892566706 nrkept=232 nrkeptmult=232 
Kept: 232 out of 912, ratio=0.254
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 31/42 (scale 4.2290757831558757e-08) ***** 
[1] xi(30)=1 zeta(31)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0069047426 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 32/42 (scale 2.4416580418283724e-08) ***** 
[1] xi(31)=1 zeta(32)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.091362489 nrkept=232 nrkeptmult=232 
Kept: 232 out of 912, ratio=0.254
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 33/42 (scale 1.4096919277186252e-08) ***** 
[1] xi(32)=1 zeta(33)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0083430139 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 34/42 (scale 8.1388601394279074e-09) ***** 
[1] xi(33)=1 zeta(34)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0920558999 nrkept=232 nrkeptmult=232 
Kept: 232 out of 912, ratio=0.254
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 35/42 (scale 4.6989730923954177e-09) ***** 
[1] xi(34)=1 zeta(35)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0088258139 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 36/42 (scale 2.7129533798093026e-09) ***** 
[1] xi(35)=1 zeta(36)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 1 states
Emax=10.0922861324 nrkept=232 nrkeptmult=232 
Kept: 232 out of 912, ratio=0.254
Memory used: 34 MB  Time elapsed: 13 s

 ***** [DMNRG] Iteration 37/42 (scale 1.5663243641318058e-09) ***** 
[1] xi(36)=1 zeta(37)=0
Stats: nr=41 min=1 max=95 sum3=3266548
Safeguard: keep additional 1 states
Emax=10.0089972888 nrkept=228 nrkeptmult=228 
Kept: 228 out of 928, ratio=0.246
Memory used: 34 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 38/42 (scale 9.0431779326976753e-10) ***** 
[1] xi(37)=1 zeta(38)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923734014 nrkept=238 nrkeptmult=238 
Kept: 238 out of 912, ratio=0.261
Memory used: 34 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 39/42 (scale 5.2210812137726861e-10) ***** 
[1] xi(38)=1 zeta(39)=0
Stats: nr=41 min=1 max=99 sum3=3556156
Safeguard: keep additional 1 states
Emax=10.0089577533 nrkept=228 nrkeptmult=228 
Kept: 228 out of 952, ratio=0.239
Memory used: 34 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 40/42 (scale 3.0143926442325579e-10) ***** 
[1] xi(39)=1 zeta(40)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923886546 nrkept=238 nrkeptmult=238 
Kept: 238 out of 912, ratio=0.261
Memory used: 34 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 41/42 (scale 1.7403604045908955e-10) ***** 
[1] xi(40)=1 zeta(41)=0
Stats: nr=41 min=1 max=99 sum3=3556156
Safeguard: keep additional 1 states
Emax=10.0089756078 nrkept=228 nrkeptmult=228 
Kept: 228 out of 952, ratio=0.239
Memory used: 34 MB  Time elapsed: 14 s

 ***** [DMNRG] Iteration 42/42 (scale 1.0047975480775194e-10) ***** 
[1] xi(41)=1 zeta(42)=0
Stats: nr=39 min=1 max=91 sum3=3180138
Safeguard: keep additional 7 states
Emax=10.0923971214 nrkept=912 nrkeptmult=912 
Kept: 912 out of 912, ratio=1
Memory used: 34 MB  Time elapsed: 14 s

** Iteration completed.

Spectrum: A_d-A_d-d FDM
[A_d-A_d-d] pos=0.223802846136 neg=0.776197153865 sum= 1
mu1=-0.00283692837452 mu2=0.000669610759636 mu3=-3.9371439046e-06 mu4=0.00016604219329

Spectrum: A_d-A_d-u FDM
[A_d-A_d-u] pos=0.776197153865 neg=0.223802846136 sum= 1
mu1=0.00283692837452 mu2=0.000669610759636 mu3=3.93714390459e-06 mu4=0.00016604219329


Memory usage report
===================
  after diag: 39048 kB
after recalc: 43780 kB
 after split: 40680 kB
  after trim: 37176 kB
 after trunc: 39048 kB
      diag_h: 27196 kB

Peak usage: 43780 kB

Timing report
=============
         All: 15 s
          DM: 0 s
         FDM: 0 s
     broaden: 0 s
        diag: 4 s
    recalc d: 1 s
    recalc f: 0 s
    recalc s: 0 s
      spec d: 4 s
       Other: 2 s
