In [27]:
# output = '/Runs/run4/output_Run3.txt'+repr(3)+'.txt' #link de harry potter, se puede cambiar el libro cambiando el link
output = 'Runs/Maj/run5/outputMajorana_Run.txt'
#A continuacion se presentan los archivos donde se escribiran los respectivos codigos y textos
#writeHuff = "descomp_huffman_" + book
#writeBin = "descomp_plano_" + book
from __future__ import division
%pylab inline
import string
import os
infile = open(os.path.abspath(output), 'r')
#load the full text by lines
text = infile.readlines()
#text
Populating the interactive namespace from numpy and matplotlib
WARNING: pylab import has clobbered these variables: ['text']
`%matplotlib` prevents importing * from pylab and numpy
In [28]:
#Nsites = int(text[4].split(' ')[-1].split('\n')[0])
#Nsites
text
Out[28]:
[' Got Model = 1chNupPdn_Majorana Got 1chNupPdn Majorana \n',
' - Model No: 7\n',
' Symmetry : OneChNupPdn - SymNo : 6\n',
' BandType : SquareWilson - BandNo : 0\n',
' NRG: Model : 1chNupPdn_Majorana - Symmetry : OneChNupPdn - Band Type: SquareWilson\n',
' Nsitesmax = 5\n',
' Ncutoff = 1000\n',
' U = 0.5\n',
' Gamma = 0.0282691\n',
' ed = -0.25\n',
' Lambda = 2.5\n',
' ChemPot = 0\n',
' UpdateBefCut = 0\n',
' calcdens = 0\n',
' Oliveira z = 1\n',
' Additional parameters : \n',
' Param 3 = 0\n',
' Param 4 = 1\n',
' Param 5 = 0\n',
' Param 6 = 0\n',
' Param 7 = 0\n',
' Param 8 = 0\n',
' Param 9 = 0\n',
' SetSingleSite begin... \n',
' Got OneChNupPdn \n',
' .. SetSingleSite done!\n',
' Fsq = 2 chi_m1 = 0.121207 eN(0) = 0\n',
' sqrt(ALambda)*chi_m1 = 0.125319\n',
'Zstring = _zEQ1\n',
' Majorana + QD effective model \n',
' Lambda= 2.5 HalfLambdaFactor= 0.7\n',
' U1~= 0.285714 U2~= 0 ed1~= -0.142857 ed2~= 0 hz~= 0\n',
' t1~= 0.571429 t2~= 0 phimag= 0 em~= 0 t12~= 0 t22~= 0 phimag2= 0 U+~= 0.285714 ed+~= -0.142857\n',
' c1_up: \n',
' Block (0) x Block (2) : 0.000000 -1.000000 x 1.000000 -1.000000 \n',
'Nsti = 0 to 3 ; Nstj = 8 to 15\n',
'(1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (1,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (1,0) (0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (1) x Block (3) : 0.000000 1.000000 x 1.000000 1.000000 \n',
'Nsti = 4 to 7 ; Nstj = 16 to 23\n',
'(1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (1,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (1,0) (0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (2) x Block (4) : 1.000000 -1.000000 x 2.000000 -1.000000 \n',
'Nsti = 8 to 15 ; Nstj = 24 to 27\n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(1,0) (0,0) (0,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) \n',
'(0,0) (0,0) (1,0) (0,0) \n',
'(0,0) (0,0) (0,0) (1,0) \n',
'\n',
' Block (3) x Block (5) : 1.000000 1.000000 x 2.000000 1.000000 \n',
'Nsti = 16 to 23 ; Nstj = 28 to 31\n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(1,0) (0,0) (0,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) \n',
'(0,0) (0,0) (1,0) (0,0) \n',
'(0,0) (0,0) (0,0) (1,0) \n',
'\n',
' c1_dn \n',
' Block (0) x Block (1) : 0.000000 -1.000000 x 0.000000 1.000000 \n',
'Nsti = 0 to 3 ; Nstj = 4 to 7\n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (1,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (1) x Block (0) : 0.000000 1.000000 x 0.000000 -1.000000 \n',
'Nsti = 4 to 7 ; Nstj = 0 to 3\n',
'(0,0) (0,0) (0,0) (1,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(1,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (2) x Block (3) : 1.000000 -1.000000 x 1.000000 1.000000 \n',
'Nsti = 8 to 15 ; Nstj = 16 to 23\n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (-1,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (-1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (1,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (1,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (3) x Block (2) : 1.000000 1.000000 x 1.000000 -1.000000 \n',
'Nsti = 16 to 23 ; Nstj = 8 to 15\n',
'(0,0) (0,0) (0,0) (-1,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(-1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (1,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (1,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (4) x Block (5) : 2.000000 -1.000000 x 2.000000 1.000000 \n',
'Nsti = 24 to 27 ; Nstj = 28 to 31\n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (-1,0) (0,0) \n',
'(0,0) (-1,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (5) x Block (4) : 2.000000 1.000000 x 2.000000 -1.000000 \n',
'Nsti = 28 to 31 ; Nstj = 24 to 27\n',
'(0,0) (0,0) (0,0) (-1,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(-1,0) (0,0) (0,0) (0,0) \n',
'\n',
' c2_dn \n',
' Block (0) x Block (1) : 0.000000 -1.000000 x 0.000000 1.000000 \n',
'Nsti = 0 to 3 ; Nstj = 4 to 7\n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (-1,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) \n',
'\n',
' Block (1) x Block (0) : 0.000000 1.000000 x 0.000000 -1.000000 \n',
'Nsti = 4 to 7 ; Nstj = 0 to 3\n',
'(0,0) (0,0) (1,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(-1,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (2) x Block (3) : 1.000000 -1.000000 x 1.000000 1.000000 \n',
'Nsti = 8 to 15 ; Nstj = 16 to 23\n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (-1,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (-1,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (1,0) (0,0) (0,0) \n',
'\n',
' Block (3) x Block (2) : 1.000000 1.000000 x 1.000000 -1.000000 \n',
'Nsti = 16 to 23 ; Nstj = 8 to 15\n',
'(0,0) (0,0) (-1,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (-1,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (1,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (4) x Block (5) : 2.000000 -1.000000 x 2.000000 1.000000 \n',
'Nsti = 24 to 27 ; Nstj = 28 to 31\n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (1,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (-1,0) (0,0) (0,0) \n',
'\n',
' Block (5) x Block (4) : 2.000000 1.000000 x 2.000000 -1.000000 \n',
'Nsti = 28 to 31 ; Nstj = 24 to 27\n',
'(0,0) (0,0) (1,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(-1,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Setting up 1-body observables...\n',
' Ndot2: \n',
' Block (0) x Block (0) : 0.000000 -1.000000 x 0.000000 -1.000000 \n',
'Nsti = 0 to 3 ; Nstj = 0 to 3\n',
'(1,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (1,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (1) x Block (1) : 0.000000 1.000000 x 0.000000 1.000000 \n',
'Nsti = 4 to 7 ; Nstj = 4 to 7\n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (1,0) \n',
'\n',
' Block (2) x Block (2) : 1.000000 -1.000000 x 1.000000 -1.000000 \n',
'Nsti = 8 to 15 ; Nstj = 8 to 15\n',
'(1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (1,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (2,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (1,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (2,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (1,0) \n',
'\n',
' Block (3) x Block (3) : 1.000000 1.000000 x 1.000000 1.000000 \n',
'Nsti = 16 to 23 ; Nstj = 16 to 23\n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (1,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (1,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (2,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (1,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (2,0) \n',
'\n',
' Block (4) x Block (4) : 2.000000 -1.000000 x 2.000000 -1.000000 \n',
'Nsti = 24 to 27 ; Nstj = 24 to 27\n',
'(2,0) (0,0) (0,0) (0,0) \n',
'(0,0) (1,0) (0,0) (0,0) \n',
'(0,0) (0,0) (2,0) (0,0) \n',
'(0,0) (0,0) (0,0) (1,0) \n',
'\n',
' Block (5) x Block (5) : 2.000000 1.000000 x 2.000000 1.000000 \n',
'Nsti = 28 to 31 ; Nstj = 28 to 31\n',
'(1,0) (0,0) (0,0) (0,0) \n',
'(0,0) (2,0) (0,0) (0,0) \n',
'(0,0) (0,0) (1,0) (0,0) \n',
'(0,0) (0,0) (0,0) (2,0) \n',
'\n',
' Sz1: \n',
' Block (0) x Block (0) : 0.000000 -1.000000 x 0.000000 -1.000000 \n',
'Nsti = 0 to 3 ; Nstj = 0 to 3\n',
'(-0.5,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (-0.5,0) \n',
'\n',
' Block (1) x Block (1) : 0.000000 1.000000 x 0.000000 1.000000 \n',
'Nsti = 4 to 7 ; Nstj = 4 to 7\n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (-0.5,0) (0,0) (0,0) \n',
'(0,0) (0,0) (-0.5,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (2) x Block (2) : 1.000000 -1.000000 x 1.000000 -1.000000 \n',
'Nsti = 8 to 15 ; Nstj = 8 to 15\n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0.5,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0.5,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (-0.5,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (-0.5,0) \n',
'\n',
' Block (3) x Block (3) : 1.000000 1.000000 x 1.000000 1.000000 \n',
'Nsti = 16 to 23 ; Nstj = 16 to 23\n',
'(0.5,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0.5,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (-0.5,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (-0.5,0) \n',
'\n',
' Block (4) x Block (4) : 2.000000 -1.000000 x 2.000000 -1.000000 \n',
'Nsti = 24 to 27 ; Nstj = 24 to 27\n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0.5,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0.5,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'\n',
' Block (5) x Block (5) : 2.000000 1.000000 x 2.000000 1.000000 \n',
'Nsti = 28 to 31 ; Nstj = 28 to 31\n',
'(0.5,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0.5,0) \n',
'\n',
' Setting up H_(N =-1 : ibl = 0 size = 4 (Nbls = 6)\n',
' Limits : 0 3\n',
' Setting up H_(N =-1 : ibl = 1 size = 4 (Nbls = 6)\n',
' Limits : 4 7\n',
' Setting up H_(N =-1 : ibl = 2 size = 8 (Nbls = 6)\n',
' Limits : 8 15\n',
' Setting up H_(N =-1 : ibl = 3 size = 8 (Nbls = 6)\n',
' Limits : 16 23\n',
' Setting up H_(N =-1 : ibl = 4 size = 4 (Nbls = 6)\n',
' Limits : 24 27\n',
' Setting up H_(N =-1 : ibl = 5 size = 4 (Nbls = 6)\n',
' Limits : 28 31\n',
' Done setting up HN. \n',
'Updating Aeig... \n',
' ******************* \n',
'Block no. 0 : 0.000000 -1.000000 \n',
' Beg: 0 End: 3\n',
'0 Eold = 0.000000e+00 \n',
'1 Eold = 1.428571e-01 \n',
'2 Eold = 1.428571e-01 \n',
'3 Eold = 0.000000e+00 \n',
' Block (0) x Block (0) : 0.000000 -1.000000 x 0.000000 -1.000000 \n',
'Nsti = 0 to 3 ; Nstj = 0 to 3\n',
'(0,0) (0,0) (-0.404061,-0) (0,0) \n',
'(0,-0) (0.142857,0) (0,0) (0.404061,0) \n',
'(-0.404061,0) (0,-0) (0.142857,0) (0,0) \n',
'(0,-0) (0.404061,-0) (0,-0) (0,0) \n',
'\n',
' Diagonalizing block 0 of 5\n',
'Upper triangular\n',
'Block no. 0 : 0.000000 -1.000000 \n',
' dEn : \n',
'0 : -3.388973e-01\n',
'1 : -3.388973e-01\n',
'2 : 4.817545e-01\n',
'3 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(3.384437e-01,0.000000e+00) (5.765270e-01,0.000000e+00) (2.838623e-01,0.000000e+00) (-6.873825e-01,0.000000e+00) \n',
'(6.873825e-01,0.000000e+00) (-2.838623e-01,0.000000e+00) (5.765270e-01,0.000000e+00) (3.384437e-01,0.000000e+00) \n',
'(-6.265361e-01,0.000000e+00) (1.703488e-01,0.000000e+00) (7.470074e-01,0.000000e+00) (1.428763e-01,0.000000e+00) \n',
'(1.428763e-01,0.000000e+00) (7.470074e-01,0.000000e+00) (-1.703488e-01,0.000000e+00) (6.265361e-01,0.000000e+00) \n',
' ******************* \n',
'Block no. 1 : 0.000000 1.000000 \n',
' Beg: 4 End: 7\n',
'4 Eold = 1.428571e-01 \n',
'5 Eold = 0.000000e+00 \n',
'6 Eold = 0.000000e+00 \n',
'7 Eold = 1.428571e-01 \n',
' Block (1) x Block (1) : 0.000000 1.000000 x 0.000000 1.000000 \n',
'Nsti = 4 to 7 ; Nstj = 4 to 7\n',
'(0.142857,0) (0,0) (0.404061,-0) (0,-0) \n',
'(0,-0) (0,0) (0,-0) (-0.404061,0) \n',
'(0.404061,0) (0,0) (0,0) (0,0) \n',
'(0,0) (-0.404061,-0) (0,-0) (0.142857,0) \n',
'\n',
' Diagonalizing block 1 of 5\n',
'Upper triangular\n',
'Block no. 1 : 0.000000 1.000000 \n',
' dEn : \n',
'4 : -3.388973e-01\n',
'5 : -3.388973e-01\n',
'6 : 4.817545e-01\n',
'7 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(-6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (7.661846e-01,-0.000000e+00) (0.000000e+00,-0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) (-1.701272e-16,0.000000e+00) (6.426206e-01,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (1.426904e-16,0.000000e+00) (7.661846e-01,0.000000e+00) \n',
'(-7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) \n',
' ******************* \n',
'Block no. 2 : 1.000000 -1.000000 \n',
' Beg: 8 End: 15\n',
'8 Eold = 1.428571e-01 \n',
'9 Eold = 0.000000e+00 \n',
'10 Eold = 0.000000e+00 \n',
'11 Eold = 1.428571e-01 \n',
'12 Eold = 0.000000e+00 \n',
'13 Eold = 1.428571e-01 \n',
'14 Eold = 0.000000e+00 \n',
'15 Eold = 0.000000e+00 \n',
' Block (2) x Block (2) : 1.000000 -1.000000 x 1.000000 -1.000000 \n',
'Nsti = 8 to 15 ; Nstj = 8 to 15\n',
'(0.142857,0) (0,0) (-0.404061,-0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,-0) (0,0) (0,0) (0.404061,0) (0,0) (-0,0) (0,0) (0,0) \n',
'(-0.404061,0) (0,-0) (0,0) (0,0) (0,0) (0,0) (-0,0) (0,0) \n',
'(0,-0) (0.404061,-0) (0,-0) (0.142857,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,-0) (0,-0) (0,-0) (0,-0) (0,0) (0,0) (0.404061,0) (0,0) \n',
'(0,-0) (-0,-0) (0,-0) (0,-0) (0,-0) (0.142857,0) (0,0) (-0.404061,-0) \n',
'(0,-0) (0,-0) (-0,-0) (0,-0) (0.404061,-0) (0,-0) (0,0) (-0,0) \n',
'(0,-0) (0,-0) (0,-0) (0,-0) (0,-0) (-0.404061,0) (-0,-0) (0,0) \n',
'\n',
' Diagonalizing block 2 of 5\n',
'Upper triangular\n',
'Block no. 2 : 1.000000 -1.000000 \n',
' dEn : \n',
'8 : -4.040610e-01\n',
'9 : -3.388973e-01\n',
'10 : -3.388973e-01\n',
'11 : -3.388973e-01\n',
'12 : 4.040610e-01\n',
'13 : 4.817545e-01\n',
'14 : 4.817545e-01\n',
'15 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-7.071068e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (7.071068e-01,-0.000000e+00) (0.000000e+00,-0.000000e+00) \n',
'(-6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (-7.661846e-01,0.000000e+00) (0.000000e+00,-0.000000e+00) (0.000000e+00,-0.000000e+00) (0.000000e+00,-0.000000e+00) (0.000000e+00,-0.000000e+00) (0.000000e+00,-0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (-7.661846e-01,0.000000e+00) (-1.701272e-16,0.000000e+00) (6.426206e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (1.426904e-16,0.000000e+00) (-7.661846e-01,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (-7.071068e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (-7.071068e-01,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) (-1.701272e-16,0.000000e+00) (-6.426206e-01,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (6.426206e-01,0.000000e+00) (1.426904e-16,0.000000e+00) (7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(-7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) \n',
' ******************* \n',
'Block no. 3 : 1.000000 1.000000 \n',
' Beg: 16 End: 23\n',
'16 Eold = 0.000000e+00 \n',
'17 Eold = 1.428571e-01 \n',
'18 Eold = 1.428571e-01 \n',
'19 Eold = 0.000000e+00 \n',
'20 Eold = 1.428571e-01 \n',
'21 Eold = 0.000000e+00 \n',
'22 Eold = 0.000000e+00 \n',
'23 Eold = 1.428571e-01 \n',
' Block (3) x Block (3) : 1.000000 1.000000 x 1.000000 1.000000 \n',
'Nsti = 16 to 23 ; Nstj = 16 to 23\n',
'(0,0) (0,0) (0.404061,-0) (0,-0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,-0) (0.142857,0) (0,-0) (-0.404061,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0.404061,0) (0,0) (0.142857,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (-0.404061,-0) (0,-0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,-0) (0,-0) (0,-0) (0,-0) (0.142857,0) (0,0) (-0.404061,0) (0,-0) \n',
'(0,-0) (0,-0) (0,-0) (0,-0) (0,-0) (0,0) (0,-0) (0.404061,-0) \n',
'(0,-0) (0,-0) (0,-0) (0,-0) (-0.404061,-0) (0,0) (0,0) (-0,0) \n',
'(0,-0) (0,-0) (0,-0) (0,-0) (0,0) (0.404061,0) (-0,-0) (0.142857,0) \n',
'\n',
' Diagonalizing block 3 of 5\n',
'Upper triangular\n',
'Block no. 3 : 1.000000 1.000000 \n',
' dEn : \n',
'16 : -3.388973e-01\n',
'17 : -3.388973e-01\n',
'18 : -3.388973e-01\n',
'19 : -3.388973e-01\n',
'20 : 4.817545e-01\n',
'21 : 4.817545e-01\n',
'22 : 4.817545e-01\n',
'23 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(-3.384437e-01,0.000000e+00) (-5.765270e-01,0.000000e+00) (2.838623e-01,0.000000e+00) (-6.873825e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(-6.873825e-01,0.000000e+00) (2.838623e-01,0.000000e+00) (5.765270e-01,0.000000e+00) (3.384437e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) \n',
'(-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (-7.661846e-01,0.000000e+00) (0.000000e+00,-0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (-7.661846e-01,0.000000e+00) (-1.701272e-16,0.000000e+00) (6.426206e-01,0.000000e+00) \n',
'(6.265361e-01,0.000000e+00) (-1.703488e-01,0.000000e+00) (7.470074e-01,0.000000e+00) (1.428763e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (6.426206e-01,0.000000e+00) (1.426904e-16,0.000000e+00) (7.661846e-01,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (-7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (6.426206e-01,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(-1.428763e-01,0.000000e+00) (-7.470074e-01,0.000000e+00) (-1.703488e-01,0.000000e+00) (6.265361e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
' ******************* \n',
'Block no. 4 : 2.000000 -1.000000 \n',
' Beg: 24 End: 27\n',
'24 Eold = 1.428571e-01 \n',
'25 Eold = 0.000000e+00 \n',
'26 Eold = 0.000000e+00 \n',
'27 Eold = 1.428571e-01 \n',
' Block (4) x Block (4) : 2.000000 -1.000000 x 2.000000 -1.000000 \n',
'Nsti = 24 to 27 ; Nstj = 24 to 27\n',
'(0.142857,0) (0,0) (0.404061,0) (0,0) \n',
'(0,-0) (0,0) (0,0) (-0.404061,-0) \n',
'(0.404061,-0) (0,-0) (0,0) (-0,0) \n',
'(0,-0) (-0.404061,0) (-0,-0) (0.142857,0) \n',
'\n',
' Diagonalizing block 4 of 5\n',
'Upper triangular\n',
'Block no. 4 : 2.000000 -1.000000 \n',
' dEn : \n',
'24 : -3.388973e-01\n',
'25 : -3.388973e-01\n',
'26 : 4.817545e-01\n',
'27 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(-6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (7.661846e-01,-0.000000e+00) (0.000000e+00,-0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) (-1.701272e-16,0.000000e+00) (6.426206e-01,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (1.426904e-16,0.000000e+00) (7.661846e-01,0.000000e+00) \n',
'(-7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) \n',
' ******************* \n',
'Block no. 5 : 2.000000 1.000000 \n',
' Beg: 28 End: 31\n',
'28 Eold = 0.000000e+00 \n',
'29 Eold = 1.428571e-01 \n',
'30 Eold = 1.428571e-01 \n',
'31 Eold = 0.000000e+00 \n',
' Block (5) x Block (5) : 2.000000 1.000000 x 2.000000 1.000000 \n',
'Nsti = 28 to 31 ; Nstj = 28 to 31\n',
'(0,0) (0,0) (-0.404061,0) (0,-0) \n',
'(0,-0) (0.142857,0) (0,-0) (0.404061,-0) \n',
'(-0.404061,-0) (0,0) (0.142857,0) (-0,0) \n',
'(0,0) (0.404061,0) (-0,-0) (0,0) \n',
'\n',
' Diagonalizing block 5 of 5\n',
'Upper triangular\n',
'Block no. 5 : 2.000000 1.000000 \n',
' dEn : \n',
'28 : -3.388973e-01\n',
'29 : -3.388973e-01\n',
'30 : 4.817545e-01\n',
'31 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(3.384437e-01,0.000000e+00) (5.765270e-01,0.000000e+00) (2.838623e-01,0.000000e+00) (-6.873825e-01,0.000000e+00) \n',
'(6.873825e-01,0.000000e+00) (-2.838623e-01,0.000000e+00) (5.765270e-01,0.000000e+00) (3.384437e-01,0.000000e+00) \n',
'(-6.265361e-01,0.000000e+00) (1.703488e-01,0.000000e+00) (7.470074e-01,0.000000e+00) (1.428763e-01,0.000000e+00) \n',
'(1.428763e-01,0.000000e+00) (7.470074e-01,0.000000e+00) (-1.703488e-01,0.000000e+00) (6.265361e-01,0.000000e+00) \n',
'Emin = -0.404061\n',
'En = 0.0651637 ---| 0 -1 >_N= -1\n',
'En = 0.0651637 ---| 0 -1 >_N= -1\n',
'En = 0.885815 ---| 0 -1 >_N= -1\n',
'En = 0.885815 ---| 0 -1 >_N= -1\n',
'En = 0.0651637 ---| 0 1 >_N= -1\n',
'En = 0.0651637 ---| 0 1 >_N= -1\n',
'En = 0.885815 ---| 0 1 >_N= -1\n',
'En = 0.885815 ---| 0 1 >_N= -1\n',
'En = 0 ---| 1 -1 >_N= -1\n',
'En = 0.0651637 ---| 1 -1 >_N= -1\n',
'En = 0.0651637 ---| 1 -1 >_N= -1\n',
'En = 0.0651637 ---| 1 -1 >_N= -1\n',
'En = 0.808122 ---| 1 -1 >_N= -1\n',
'En = 0.885815 ---| 1 -1 >_N= -1\n',
'En = 0.885815 ---| 1 -1 >_N= -1\n',
'En = 0.885815 ---| 1 -1 >_N= -1\n',
'En = 0.0651637 ---| 1 1 >_N= -1\n',
'En = 0.0651637 ---| 1 1 >_N= -1\n',
'En = 0.0651637 ---| 1 1 >_N= -1\n',
'En = 0.0651637 ---| 1 1 >_N= -1\n',
'En = 0.885815 ---| 1 1 >_N= -1\n',
'En = 0.885815 ---| 1 1 >_N= -1\n',
'En = 0.885815 ---| 1 1 >_N= -1\n',
'En = 0.885815 ---| 1 1 >_N= -1\n',
'En = 0.0651637 ---| 2 -1 >_N= -1\n',
'En = 0.0651637 ---| 2 -1 >_N= -1\n',
'En = 0.885815 ---| 2 -1 >_N= -1\n',
'En = 0.885815 ---| 2 -1 >_N= -1\n',
'En = 0.0651637 ---| 2 1 >_N= -1\n',
'En = 0.0651637 ---| 2 1 >_N= -1\n',
'En = 0.885815 ---| 2 1 >_N= -1\n',
'En = 0.885815 ---| 2 1 >_N= -1\n',
' \n',
' \n',
' \n',
' \n',
' \n',
' \n',
' \n',
' \n',
' \n',
' PRINTING MATRICES \n',
' Holaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa \n',
' Setting up H_(N =-1 : ibl = 0 size = 4 (Nbls = 6)\n',
' Limits : 0 3\n',
' Setting up H_(N =-1 : ibl = 1 size = 4 (Nbls = 6)\n',
' Limits : 4 7\n',
' Setting up H_(N =-1 : ibl = 2 size = 8 (Nbls = 6)\n',
' Limits : 8 15\n',
' Setting up H_(N =-1 : ibl = 3 size = 8 (Nbls = 6)\n',
' Limits : 16 23\n',
' Setting up H_(N =-1 : ibl = 4 size = 4 (Nbls = 6)\n',
' Limits : 24 27\n',
' Setting up H_(N =-1 : ibl = 5 size = 4 (Nbls = 6)\n',
' Limits : 28 31\n',
' Done setting up HN. \n',
'Updating Aeig... \n',
' ******************* \n',
'Block no. 0 : 0.000000 -1.000000 \n',
' Beg: 0 End: 3\n',
'0 Eold = 0.000000e+00 \n',
'1 Eold = 1.428571e-01 \n',
'2 Eold = 1.428571e-01 \n',
'3 Eold = 0.000000e+00 \n',
' Block (0) x Block (0) : 0.000000 -1.000000 x 0.000000 -1.000000 \n',
'Nsti = 0 to 3 ; Nstj = 0 to 3\n',
'(0,0) (0,0) (-0.404061,-0) (0,0) \n',
'(0,0) (0.142857,0) (0,0) (0.404061,0) \n',
'(-0.404061,0) (0,-0) (0.142857,0) (0,0) \n',
'(0,-0) (0.404061,-0) (0,0) (0,0) \n',
'\n',
' Diagonalizing block 0 of 5\n',
'Regular block\n',
'Block no. 0 : 0.000000 -1.000000 \n',
' dEn : \n',
'0 : -3.388973e-01\n',
'1 : -3.388973e-01\n',
'2 : 4.817545e-01\n',
'3 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(7.186642e-01,0.000000e+00) (-2.227934e-01,0.000000e+00) (6.027639e-01,0.000000e+00) (2.656324e-01,0.000000e+00) \n',
'(-2.656324e-01,0.000000e+00) (-6.027639e-01,0.000000e+00) (-2.227934e-01,-0.000000e+00) (7.186642e-01,-0.000000e+00) \n',
'(6.426206e-01,0.000000e+00) (1.785776e-16,0.000000e+00) (-7.661846e-01,0.000000e+00) (-8.375841e-17,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) (-5.182684e-17,0.000000e+00) (6.426206e-01,0.000000e+00) \n',
' ******************* \n',
'Block no. 1 : 0.000000 1.000000 \n',
' Beg: 4 End: 7\n',
'4 Eold = 1.428571e-01 \n',
'5 Eold = 0.000000e+00 \n',
'6 Eold = 0.000000e+00 \n',
'7 Eold = 1.428571e-01 \n',
' Block (1) x Block (1) : 0.000000 1.000000 x 0.000000 1.000000 \n',
'Nsti = 4 to 7 ; Nstj = 4 to 7\n',
'(0.142857,0) (0,0) (0.404061,-0) (0,-0) \n',
'(0,0) (0,0) (0,-0) (-0.404061,0) \n',
'(0.404061,0) (0,0) (0,0) (0,0) \n',
'(0,0) (-0.404061,-0) (0,0) (0.142857,0) \n',
'\n',
' Diagonalizing block 1 of 5\n',
'Regular block\n',
'Block no. 1 : 0.000000 1.000000 \n',
' dEn : \n',
'4 : -3.388973e-01\n',
'5 : -3.388973e-01\n',
'6 : 4.817545e-01\n',
'7 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(6.426206e-01,0.000000e+00) (-1.701272e-16,0.000000e+00) (-7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(-0.000000e+00,0.000000e+00) (-7.661846e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) \n',
'(-7.661846e-01,0.000000e+00) (-1.426904e-16,0.000000e+00) (-6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) \n',
' ******************* \n',
'Block no. 2 : 1.000000 -1.000000 \n',
' Beg: 8 End: 15\n',
'8 Eold = 1.428571e-01 \n',
'9 Eold = 0.000000e+00 \n',
'10 Eold = 0.000000e+00 \n',
'11 Eold = 1.428571e-01 \n',
'12 Eold = 0.000000e+00 \n',
'13 Eold = 1.428571e-01 \n',
'14 Eold = 0.000000e+00 \n',
'15 Eold = 0.000000e+00 \n',
' Block (2) x Block (2) : 1.000000 -1.000000 x 1.000000 -1.000000 \n',
'Nsti = 8 to 15 ; Nstj = 8 to 15\n',
'(0.142857,0) (0,0) (-0.404061,-0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0.404061,0) (0,0) (-0,0) (0,0) (0,0) \n',
'(-0.404061,0) (0,-0) (0,0) (0,0) (0,0) (0,0) (-0,0) (0,0) \n',
'(0,-0) (0.404061,-0) (0,0) (0.142857,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0.404061,0) (0,0) \n',
'(0,0) (-0,0) (0,0) (0,0) (0,0) (0.142857,0) (0,0) (-0.404061,-0) \n',
'(0,0) (0,0) (-0,0) (0,0) (0.404061,-0) (0,-0) (0,0) (-0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,-0) (-0.404061,0) (-0,0) (0,0) \n',
'\n',
' Diagonalizing block 2 of 5\n',
'Regular block\n',
'Block no. 2 : 1.000000 -1.000000 \n',
' dEn : \n',
'8 : -4.040610e-01\n',
'9 : -3.388973e-01\n',
'10 : -3.388973e-01\n',
'11 : -3.388973e-01\n',
'12 : 4.040610e-01\n',
'13 : 4.817545e-01\n',
'14 : 4.817545e-01\n',
'15 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (-7.071068e-01,0.000000e+00) (1.177569e-16,0.000000e+00) (7.071068e-01,0.000000e+00) (2.747662e-16,0.000000e+00) \n',
'(-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (2.977225e-16,0.000000e+00) (6.426206e-01,-0.000000e+00) (-2.497082e-16,0.000000e+00) (7.661846e-01,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) \n',
'(6.426206e-01,0.000000e+00) (-1.701272e-16,0.000000e+00) (7.661846e-01,-0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (7.071068e-01,0.000000e+00) (1.177569e-16,0.000000e+00) (7.071068e-01,0.000000e+00) (-2.747662e-16,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (2.497082e-16,0.000000e+00) (-7.661846e-01,0.000000e+00) (2.977225e-16,0.000000e+00) (6.426206e-01,0.000000e+00) \n',
'(-7.661846e-01,0.000000e+00) (-1.426904e-16,0.000000e+00) (6.426206e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (6.426206e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
' ******************* \n',
'Block no. 3 : 1.000000 1.000000 \n',
' Beg: 16 End: 23\n',
'16 Eold = 0.000000e+00 \n',
'17 Eold = 1.428571e-01 \n',
'18 Eold = 1.428571e-01 \n',
'19 Eold = 0.000000e+00 \n',
'20 Eold = 1.428571e-01 \n',
'21 Eold = 0.000000e+00 \n',
'22 Eold = 0.000000e+00 \n',
'23 Eold = 1.428571e-01 \n',
' Block (3) x Block (3) : 1.000000 1.000000 x 1.000000 1.000000 \n',
'Nsti = 16 to 23 ; Nstj = 16 to 23\n',
'(0,0) (0,0) (0.404061,-0) (0,-0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0.142857,0) (0,-0) (-0.404061,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0.404061,0) (0,0) (0.142857,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (-0.404061,-0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0.142857,0) (0,0) (-0.404061,0) (0,-0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0,0) (0,-0) (0.404061,-0) \n',
'(0,0) (0,0) (0,0) (0,0) (-0.404061,-0) (0,0) (0,0) (-0,0) \n',
'(0,0) (0,0) (0,0) (0,0) (0,0) (0.404061,0) (-0,0) (0.142857,0) \n',
'\n',
' Diagonalizing block 3 of 5\n',
'Regular block\n',
'Block no. 3 : 1.000000 1.000000 \n',
' dEn : \n',
'16 : -3.388973e-01\n',
'17 : -3.388973e-01\n',
'18 : -3.388973e-01\n',
'19 : -3.388973e-01\n',
'20 : 4.817545e-01\n',
'21 : 4.817545e-01\n',
'22 : 4.817545e-01\n',
'23 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(7.186642e-01,0.000000e+00) (-2.227934e-01,0.000000e+00) (-6.027639e-01,0.000000e+00) (-2.656324e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (-1.701272e-16,0.000000e+00) (-7.661846e-01,0.000000e+00) (-0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (-7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (6.426206e-01,0.000000e+00) \n',
'(-2.656324e-01,0.000000e+00) (-6.027639e-01,0.000000e+00) (2.227934e-01,0.000000e+00) (-7.186642e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(-6.426206e-01,0.000000e+00) (-1.785776e-16,0.000000e+00) (-7.661846e-01,0.000000e+00) (-8.375841e-17,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (-7.661846e-01,0.000000e+00) (1.426904e-16,0.000000e+00) (6.426206e-01,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (-7.661846e-01,0.000000e+00) (-5.182684e-17,0.000000e+00) (6.426206e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (0.000000e+00,0.000000e+00) (6.426206e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) \n',
' ******************* \n',
'Block no. 4 : 2.000000 -1.000000 \n',
' Beg: 24 End: 27\n',
'24 Eold = 1.428571e-01 \n',
'25 Eold = 0.000000e+00 \n',
'26 Eold = 0.000000e+00 \n',
'27 Eold = 1.428571e-01 \n',
' Block (4) x Block (4) : 2.000000 -1.000000 x 2.000000 -1.000000 \n',
'Nsti = 24 to 27 ; Nstj = 24 to 27\n',
'(0.142857,0) (0,0) (0.404061,0) (0,0) \n',
'(0,0) (0,0) (0,0) (-0.404061,-0) \n',
'(0.404061,-0) (0,-0) (0,0) (-0,0) \n',
'(0,-0) (-0.404061,0) (-0,0) (0.142857,0) \n',
'\n',
' Diagonalizing block 4 of 5\n',
'Regular block\n',
'Block no. 4 : 2.000000 -1.000000 \n',
' dEn : \n',
'24 : -3.388973e-01\n',
'25 : -3.388973e-01\n',
'26 : 4.817545e-01\n',
'27 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(-6.426206e-01,0.000000e+00) (-1.701272e-16,0.000000e+00) (7.661846e-01,-0.000000e+00) (-0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (6.426206e-01,0.000000e+00) \n',
'(-7.661846e-01,0.000000e+00) (1.426904e-16,0.000000e+00) (-6.426206e-01,0.000000e+00) (0.000000e+00,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (-6.426206e-01,0.000000e+00) (0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) \n',
' ******************* \n',
'Block no. 5 : 2.000000 1.000000 \n',
' Beg: 28 End: 31\n',
'28 Eold = 0.000000e+00 \n',
'29 Eold = 1.428571e-01 \n',
'30 Eold = 1.428571e-01 \n',
'31 Eold = 0.000000e+00 \n',
' Block (5) x Block (5) : 2.000000 1.000000 x 2.000000 1.000000 \n',
'Nsti = 28 to 31 ; Nstj = 28 to 31\n',
'(0,0) (0,0) (-0.404061,0) (0,-0) \n',
'(0,0) (0.142857,0) (0,-0) (0.404061,-0) \n',
'(-0.404061,-0) (0,0) (0.142857,0) (-0,0) \n',
'(0,0) (0.404061,0) (-0,0) (0,0) \n',
'\n',
' Diagonalizing block 5 of 5\n',
'Regular block\n',
'Block no. 5 : 2.000000 1.000000 \n',
' dEn : \n',
'28 : -3.388973e-01\n',
'29 : -3.388973e-01\n',
'30 : 4.817545e-01\n',
'31 : 4.817545e-01\n',
'Kept States: \n',
' dEigVec (rows): \n',
'(7.186642e-01,0.000000e+00) (2.227934e-01,0.000000e+00) (6.027639e-01,0.000000e+00) (-2.656324e-01,0.000000e+00) \n',
'(-2.656324e-01,0.000000e+00) (6.027639e-01,0.000000e+00) (-2.227934e-01,0.000000e+00) (-7.186642e-01,0.000000e+00) \n',
'(-6.426206e-01,0.000000e+00) (1.785776e-16,0.000000e+00) (7.661846e-01,0.000000e+00) (-8.375841e-17,0.000000e+00) \n',
'(0.000000e+00,0.000000e+00) (7.661846e-01,0.000000e+00) (5.182684e-17,0.000000e+00) (6.426206e-01,0.000000e+00) \n',
'Emin = -0.404061\n',
'\n',
' Holaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa \n',
' \n',
' \n',
' \n',
' Use parity: 1 1\n',
' NoCommon = 2\n',
' ParityQN : 0 1 \n',
' PosInOld : 0 1 \n',
' PosInSingleSite : 0 1 \n',
' MatchOld : 1 1 \n',
'No blocks = 8\n',
'No states = 128\n',
'\n',
' \n',
' \n',
' \n',
' \n',
' Helloooooooooooooo, everything ready to start \n',
' \n',
' \n',
' \n',
' Setting up H_(N =0 : ibl = 0 size = 8 (Nbls = 8)\n',
' Limits : 0 7\n',
' Setting up H_(N =0 : ibl = 1 size = 24 (Nbls = 8)\n',
' Limits : 8 31\n',
' Setting up H_(N =0 : ibl = 2 size = 8 (Nbls = 8)\n',
' Limits : 32 39\n',
' Setting up H_(N =0 : ibl = 3 size = 24 (Nbls = 8)\n',
' Limits : 40 63\n',
' Setting up H_(N =0 : ibl = 4 size = 24 (Nbls = 8)\n',
' Limits : 64 87\n',
' Setting up H_(N =0 : ibl = 5 size = 24 (Nbls = 8)\n',
' Limits : 88 111\n',
' Setting up H_(N =0 : ibl = 6 size = 8 (Nbls = 8)\n',
' Limits : 112 119\n',
' Setting up H_(N =0 : ibl = 7 size = 8 (Nbls = 8)\n',
' Limits : 120 127\n',
' Done setting up HN. \n',
'Updating Aeig... \n',
' Diagonalizing block 0 of 7\n',
'Upper triangular\n',
' Diagonalizing block 1 of 7\n',
'Upper triangular\n',
' Diagonalizing block 2 of 7\n',
'Upper triangular\n',
' Diagonalizing block 3 of 7\n',
'Upper triangular\n',
' Diagonalizing block 4 of 7\n',
'Upper triangular\n',
' Diagonalizing block 5 of 7\n',
'Upper triangular\n',
' Diagonalizing block 6 of 7\n',
'Upper triangular\n',
' Diagonalizing block 7 of 7\n',
'Upper triangular\n',
'Emin = -0.0618324\n',
'En = 0.117537 ---| 0 -1 >_N= 0\n',
'En = 0.117537 ---| 0 -1 >_N= 0\n',
'En = 0.117537 ---| 0 -1 >_N= 0\n',
'En = 0.117537 ---| 0 -1 >_N= 0\n',
'En = 0.957107 ---| 0 -1 >_N= 0\n',
'En = 0.957107 ---| 0 -1 >_N= 0\n',
'En = 0.957107 ---| 0 -1 >_N= 0\n',
'En = 0.957107 ---| 0 -1 >_N= 0\n',
'En = 0 ---| 1 -1 >_N= 0\n',
'En = 0.0383819 ---| 1 -1 >_N= 0\n',
'En = 0.0383819 ---| 1 -1 >_N= 0\n',
'En = 0.0383819 ---| 1 -1 >_N= 0\n',
'En = 0.117537 ---| 1 -1 >_N= 0\n',
'En = 0.117537 ---| 1 -1 >_N= 0\n',
'En = 0.117537 ---| 1 -1 >_N= 0\n',
'En = 0.117537 ---| 1 -1 >_N= 0\n',
'En = 0.188828 ---| 1 -1 >_N= 0\n',
'En = 0.21561 ---| 1 -1 >_N= 0\n',
'En = 0.21561 ---| 1 -1 >_N= 0\n',
'En = 0.21561 ---| 1 -1 >_N= 0\n',
'En = 0.812046 ---| 1 -1 >_N= 0\n',
'En = 0.859034 ---| 1 -1 >_N= 0\n',
'En = 0.859034 ---| 1 -1 >_N= 0\n',
'En = 0.859034 ---| 1 -1 >_N= 0\n',
'En = 0.957107 ---| 1 -1 >_N= 0\n',
'En = 0.957107 ---| 1 -1 >_N= 0\n',
'En = 0.957107 ---| 1 -1 >_N= 0\n',
'En = 0.957107 ---| 1 -1 >_N= 0\n',
'En = 1.00556 ---| 1 -1 >_N= 0\n',
'En = 1.03626 ---| 1 -1 >_N= 0\n',
'En = 1.03626 ---| 1 -1 >_N= 0\n',
'En = 1.03626 ---| 1 -1 >_N= 0\n',
'En = 0.126996 ---| 0 1 >_N= 0\n',
'En = 0.126996 ---| 0 1 >_N= 0\n',
'En = 0.126996 ---| 0 1 >_N= 0\n',
'En = 0.126996 ---| 0 1 >_N= 0\n',
'En = 0.947648 ---| 0 1 >_N= 0\n',
'En = 0.947648 ---| 0 1 >_N= 0\n',
'En = 0.947648 ---| 0 1 >_N= 0\n',
'En = 0.947648 ---| 0 1 >_N= 0\n',
'En = 0.0618324 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.126996 ---| 1 1 >_N= 0\n',
'En = 0.869954 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.947648 ---| 1 1 >_N= 0\n',
'En = 0.0618324 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.126996 ---| 2 -1 >_N= 0\n',
'En = 0.869954 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.947648 ---| 2 -1 >_N= 0\n',
'En = 0.0618324 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 2 1 >_N= 0\n',
'En = 0.869954 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.947648 ---| 2 1 >_N= 0\n',
'En = 0.126996 ---| 3 -1 >_N= 0\n',
'En = 0.126996 ---| 3 -1 >_N= 0\n',
'En = 0.126996 ---| 3 -1 >_N= 0\n',
'En = 0.126996 ---| 3 -1 >_N= 0\n',
'En = 0.947648 ---| 3 -1 >_N= 0\n',
'En = 0.947648 ---| 3 -1 >_N= 0\n',
'En = 0.947648 ---| 3 -1 >_N= 0\n',
'En = 0.947648 ---| 3 -1 >_N= 0\n',
'En = 0.126996 ---| 3 1 >_N= 0\n',
'En = 0.126996 ---| 3 1 >_N= 0\n',
'En = 0.126996 ---| 3 1 >_N= 0\n',
'En = 0.126996 ---| 3 1 >_N= 0\n',
'En = 0.947648 ---| 3 1 >_N= 0\n',
'En = 0.947648 ---| 3 1 >_N= 0\n',
'En = 0.947648 ---| 3 1 >_N= 0\n',
'En = 0.947648 ---| 3 1 >_N= 0\n',
'Block i : 0/7 (sz=8) Nshell=0\n',
' Block j : 0/7 (sz=8) Time to here 1e-05; last lap 1e-05\n',
' Updating Op : 3 of 5\n',
' Updating Op : 4 of 5\n',
' Updating Op : 5 of 5\n',
' Block j : 1/7 (sz=24) Time to here 0.000238; last lap 0.000216\n',
' Updating Op : 1 of 5\n',
' Block j : 2/7 (sz=8) Time to here 0.000579; last lap 0.000333\n',
' Updating Op : 2 of 5\n',
' Block j : 3/7 (sz=24) Time to here 0.000684; last lap 9.3e-05\n',
' Block j : 4/7 (sz=24) Time to here 0.000694; last lap 1e-06\n',
' Block j : 5/7 (sz=24) Time to here 0.000703; last lap 1e-06\n',
' Block j : 6/7 (sz=8) Time to here 0.000712; last lap 1e-06\n',
' Block j : 7/7 (sz=8) Time to here 0.00072; last lap 1e-06\n',
'Block i : 1/7 (sz=24) Nshell=0\n',
' Block j : 0/7 (sz=8) Time to here 0.000733; last lap 5e-06\n',
' Block j : 1/7 (sz=24) Time to here 0.000742; last lap 1e-06\n',
' Updating Op : 3 of 5\n',
' Updating Op : 4 of 5\n',
' Updating Op : 5 of 5\n',
' Block j : 2/7 (sz=8) Time to here 0.004779; last lap 0.00403\n',
' Block j : 3/7 (sz=24) Time to here 0.004804; last lap 2e-06\n',
' Updating Op : 2 of 5\n',
' Block j : 4/7 (sz=24) Time to here 0.006026; last lap 0.001215\n',
' Updating Op : 1 of 5\n',
' Block j : 5/7 (sz=24) Time to here 0.007261; last lap 0.001225\n',
' Block j : 6/7 (sz=8) Time to here 0.007274; last lap 1e-06\n',
' Block j : 7/7 (sz=8) Time to here 0.007283; last lap 1e-06\n',
'Block i : 2/7 (sz=8) Nshell=0\n',
' Block j : 0/7 (sz=8) Time to here 0.007295; last lap 6e-06\n',
' Updating Op : 2 of 5\n',
' Block j : 1/7 (sz=24) Time to here 0.007369; last lap 6.6e-05\n',
' Block j : 2/7 (sz=8) Time to here 0.007378; last lap 1e-06\n',
' Updating Op : 3 of 5\n',
' Updating Op : 4 of 5\n',
' Updating Op : 5 of 5\n',
' Block j : 3/7 (sz=24) Time to here 0.007588; last lap 0.000203\n',
' Updating Op : 1 of 5\n',
' Block j : 4/7 (sz=24) Time to here 0.007895; last lap 0.0003\n',
' Block j : 5/7 (sz=24) Time to here 0.007904; last lap 1e-06\n',
' Block j : 6/7 (sz=8) Time to here 0.007913; last lap 1e-06\n',
' Block j : 7/7 (sz=8) Time to here 0.007921; last lap 1e-06\n',
'Block i : 3/7 (sz=24) Nshell=0\n',
' Block j : 0/7 (sz=8) Time to here 0.007933; last lap 5e-06\n',
' Block j : 1/7 (sz=24) Time to here 0.007941; last lap 0\n',
' Updating Op : 2 of 5\n',
' Block j : 2/7 (sz=8) Time to here 0.009127; last lap 0.00118\n',
' Block j : 3/7 (sz=24) Time to here 0.00914; last lap 1e-06\n',
' Updating Op : 3 of 5\n',
' Updating Op : 4 of 5\n',
' Updating Op : 5 of 5\n',
' Block j : 4/7 (sz=24) Time to here 0.012905; last lap 0.003756\n',
' Block j : 5/7 (sz=24) Time to here 0.012927; last lap 3e-06\n',
' Updating Op : 1 of 5\n',
' Block j : 6/7 (sz=8) Time to here 0.014156; last lap 0.001222\n',
' Block j : 7/7 (sz=8) Time to here 0.014169; last lap 1e-06\n',
'Block i : 4/7 (sz=24) Nshell=0\n',
' Block j : 0/7 (sz=8) Time to here 0.014181; last lap 5e-06\n',
' Block j : 1/7 (sz=24) Time to here 0.014189; last lap 1e-06\n',
' Block j : 2/7 (sz=8) Time to here 0.014197; last lap 1e-06\n',
' Block j : 3/7 (sz=24) Time to here 0.014204; last lap 1e-06\n',
' Block j : 4/7 (sz=24) Time to here 0.014212; last lap 1e-06\n',
' Updating Op : 3 of 5\n',
' Updating Op : 4 of 5\n',
' Updating Op : 5 of 5\n',
' Block j : 5/7 (sz=24) Time to here 0.017911; last lap 0.003693\n',
' Updating Op : 2 of 5\n',
...]
In [29]:
Nsites = int(text[5].split(' ')[-1].split('\n')[0])
A = [] # Este sera el alfabeto
for x in text:
#print(list(x))
a=x.split(' ')
if x.split(' ')[0]=='En':
#En , |n,m> , N
Ens = [float(a[2]),int(a[4]),float(a[6]),int(a[-1].split('\n')[0])]
A.append(Ens)
A=array(A)
A
Out[29]:
array([[ 0.0651637, 0. , -1. , -1. ],
[ 0.0651637, 0. , -1. , -1. ],
[ 0.885815 , 0. , -1. , -1. ],
...,
[ 6.34987 , 7. , 1. , 5. ],
[ 6.69751 , 7. , 1. , 5. ],
[ 6.69751 , 7. , 1. , 5. ]])
In [35]:
numEst = 200
minEns = zeros((Nsites+1,numEst))
for n in range(0,Nsites+1):
An = A[A[:,3]==n]
An = sorted(An, key=lambda An_entry: An_entry[0]); AnMins = An[0:numEst]
for m in range(0,len(AnMins)):
minEns[n,m] = AnMins[m][0]
minEns
#plot(range(-1,Nsites+1),minEns,'.')
Out[35]:
array([[ 0. , 0.0383819, 0.0383819, ..., 0. , 0. ,
0. ],
[ 0. , 0.0752542, 0.0752542, ..., 0.50102 , 0.50102 ,
0.50102 ],
[ 0. , 0.0965742, 0.0965742, ..., 1.20435 , 1.20435 ,
1.20435 ],
[ 0. , 0.103238 , 0.103238 , ..., 0.853325 , 0.853325 ,
0.853325 ],
[ 0. , 0.243981 , 0.243981 , ..., 1.60256 , 1.60256 ,
1.60256 ],
[ 0. , 0.109371 , 0.109371 , ..., 1.67506 , 1.67506 ,
1.67506 ]])
In [36]:
print size(minEns[6::2,5])
size(range(Nsites+1)[4::2])
0
Out[36]:
1
In [37]:
print len(set(minEns[-1]))
set(minEns[-1])
45
Out[37]:
{0.0,
0.109371,
0.200492,
0.41869899999999999,
0.49082700000000001,
0.49760900000000002,
0.50492999999999999,
0.58371499999999998,
0.74963100000000005,
0.764513,
0.807056,
0.84488300000000005,
0.89525299999999997,
0.99843599999999999,
1.0041199999999999,
1.0135400000000001,
1.0641499999999999,
1.1049899999999999,
1.1254299999999999,
1.1304099999999999,
1.14764,
1.15385,
1.15795,
1.16083,
1.17076,
1.3247199999999999,
1.33328,
1.35789,
1.36849,
1.37141,
1.37649,
1.3882000000000001,
1.42306,
1.4268400000000001,
1.4276800000000001,
1.42875,
1.4297500000000001,
1.4339,
1.4509399999999999,
1.5934299999999999,
1.61802,
1.6193299999999999,
1.65273,
1.6742999999999999,
1.67506}
In [38]:
for m in range(0,numEst):
xlabel('Iteration Number(N)',fontsize=18)
ylabel('Energy Espectrum ',fontsize=18)
title('Espectrum Evolution of Odd Numbers',fontsize=20)
#label_size = 16
xticks(fontsize=16, rotation=0)
yticks(fontsize=16, rotation=0)
plot(range(Nsites+1)[1::2],minEns[1::2,m],'-');
In [26]:
print len(set(minEns[-2]))
set(minEns[-2])
1
Out[26]:
{0.0}
In [20]:
for m in range(0,numEst):
xlabel('Iteration Number(N)',fontsize=18)
ylabel('Energy Espectrum ',fontsize=18)
title('Espectrum Evolution of Even Numbers',fontsize=20)
#label_size = 16
xticks(fontsize=16, rotation=0)
yticks(fontsize=16, rotation=0)
plot(range(Nsites+1)[0::2],minEns[0:-1:2,m],'-');
In [16]:
for m in range(0,10):
plot(range(Nsites+1)[0:8:2],minEns[0:8:2,m],'-');
In [17]:
An = A[A[:,3]==0]
An = sorted(An, key=lambda An_entry: An_entry[0])
a=An[0:4]
a[2][0]
Out[17]:
0.106208
In [23]:
A= zeros((4, 3))
In [24]:
[2,3]
Out[24]:
[2, 3]
In [11]:
range(6)
Out[11]:
[0, 1, 2, 3, 4, 5]
In [12]:
len(minEns[::2,n])
Out[12]:
32
In [12]:
In [12]:
In [ ]:
Content source: cifu9502/nrgcode
Similar notebooks: