Commit efbafffb by Fakher F. Assaad

Updated Prog Analysis Libraries to ALF-1.0

1 parent bc2f74fe
Showing with 534 additions and 512 deletions
......@@ -50,7 +50,7 @@
Integer :: Nunit, Norb, ierr
Integer :: no, no1, n, n1,m, nbins, n_skip, nb, N_rebin, N_cov
Integer :: no, no1, n, n1,m, nbins, n_skip, nb, N_rebin, N_cov, N_Back
real (Kind=Kind(0.d0)):: X, Y
Complex (Kind=Kind(0.d0)), allocatable :: Phase(:)
Type (Mat_C), allocatable :: Bins (:,:), Bins_R(:,:)
......@@ -66,10 +66,12 @@
NAMELIST /VAR_lattice/ L1, L2, Lattice_type, Model
NAMELIST /VAR_errors/ n_skip, N_rebin, N_Cov
NAMELIST /VAR_errors/ n_skip, N_rebin, N_Cov, N_Back
N_Back = 1
OPEN(UNIT=5,FILE='parameters',STATUS='old',ACTION='read',IOSTAT=ierr)
IF (ierr /= 0) THEN
WRITE(*,*) 'unable to open <parameters>',ierr
......@@ -147,13 +149,15 @@
bins (n,nb)%el = cmplx(0.d0,0.d0,kind(0.d0))
Enddo
Enddo
Bins0 = cmplx(0.d0,0.d0,kind(0.d0))
Open ( Unit=10, File="ineq", status="unknown" )
do nb = 1, nbins + n_skip
if (nb > n_skip ) then
Read(10,*,End=10) X,no,no1
Phase(nb-n_skip) = cmplx(X,0.d0,kind(0.d0))
Do no = 1,Norb
Read(10,*) Bins0(nb-n_skip,no)
Read(10,*) Z
if (N_Back == 1 ) Bins0(nb-n_skip,no) = Z
enddo
do n = 1,Nunit
Read(10,*) Xk_p(1), Xk_p(2)
......
......@@ -48,23 +48,23 @@
REAL (Kind=Kind(0.d0)), DIMENSION(:), ALLOCATABLE :: EN, SIGN
REAL (Kind=Kind(0.d0)) :: XM, XERR, X
Complex (Kind=Kind(0.d0)) Z1,Z2,Z3,Z4,Z5
! Complex (Kind=Kind(0.d0)) Z1,Z2,Z3,Z4,Z5
Complex (Kind=Kind(0.d0)), Allocatable :: Tmp(:)
Integer :: Nobs
Integer :: NST, NS, NS1, NS2, NSTEP, NC, NP, Nbins, NP_EFF, ISEED, I, IOBS
Integer :: N,N1, NBIN
Integer :: Nbins, Nbins_eff, I, IOBS, N_Back
Integer :: N,N1 !, NBIN
Integer :: n_skip, N_rebin, N_Cov, ierr
NAMELIST /VAR_errors/ n_skip, N_rebin, N_Cov
NAMELIST /VAR_errors/ n_skip, N_rebin, N_Cov, N_Back
OPEN(UNIT=5,FILE='parameters',STATUS='old',ACTION='read',IOSTAT=ierr)
IF (ierr /= 0) THEN
WRITE(*,*) 'unable to open <parameters>',ierr
STOP
END IF
READ(5,NML=VAR_errors)
CLOSE(5)
OPEN(UNIT=5,FILE='parameters',STATUS='old',ACTION='read',IOSTAT=ierr)
IF (ierr /= 0) THEN
WRITE(*,*) 'unable to open <parameters>',ierr
STOP
END IF
READ(5,NML=VAR_errors)
CLOSE(5)
! Count the number of bins
......@@ -72,65 +72,53 @@
Read(10,*) NOBS
allocate (Tmp(NOBS) )
rewind(10)
nbins = 0
Nbins = 0
do
read(10,*,End=10) N, (Tmp(I), I=1,size(Tmp,1)-1), X
Tmp(NOBS) = cmplx(X,0.d0,kind(0.d0))
nbins = nbins + 1
enddo
10 continue
Write(6,*) "# of bins: ", Nbins
Close(10)
NP = NBINS
ALLOCATE(OBS(NP,NOBS))
read(10,*,End=10) N, (Tmp(I), I=1,size(Tmp,1)-1), X
Tmp(NOBS) = cmplx(X,0.d0,kind(0.d0))
Nbins = Nbins + 1
enddo
10 continue
Write(6,*) "# of bins: ", Nbins
Close(10)
ALLOCATE(OBS(Nbins,NOBS))
NST = N_skip; NS1 = N_Rebin; NS2 = N_Rebin; NSTEP = 1
OPEN (UNIT=20, FILE='Var_scal', STATUS='old')
NC = 0
DO N = 1,NP
IF (N.GE.NST) THEN
NC = NC + 1
READ(20,*) N1, (Tmp(I), I=1,size(Tmp,1)-1),X
Tmp(NOBS) = cmplx(X,0.d0,kind(0.d0))
OBS(NC,:) = dble(Tmp(:))
ELSE
READ(20,*) N1, (Tmp(I), I=1,size(Tmp,1)-1),X
ENDIF
ENDDO
CLOSE(20)
2100 FORMAT(I6,2X,F16.8)
OPEN (UNIT=20, FILE='Var_scal', STATUS='old')
Nbins_eff = 0
DO N = 1,Nbins
IF (N > N_skip) THEN
Nbins_eff = Nbins_eff + 1
READ(20,*) N1, (Tmp(I), I=1,size(Tmp,1)-1),X
Tmp(NOBS) = cmplx(X,0.d0,kind(0.d0))
OBS(Nbins_eff,:) = dble(Tmp(:))
ELSE
READ(20,*) N1, (Tmp(I), I=1,size(Tmp,1)-1),X
ENDIF
ENDDO
CLOSE(20)
2100 FORMAT(I6,2X,F16.8)
OPEN (UNIT=21, FILE='Var_scalJ', STATUS='unknown')
WRITE(21,*) 'Effective number of bins, and bins: ', NC, NP
NP_EFF = NC
ALLOCATE (EN(NP_EFF), SIGN(NP_EFF))
DO IOBS = 1,NOBS
WRITE(21,*)
DO I = 1,NP_EFF
EN (I) = Real(OBS(I,IOBS), kind(0.d0))
SIGN(I) = Real(OBS(I,NOBS), kind(0.d0))
ENDDO
DO NBIN = NS1, NS2, NSTEP
if (NBIN.gt.0) then
IF (IOBS.EQ.NOBS ) then
CALL ERRCALCJ(EN,XM,XERR,NBIN)
else
CALL ERRCALCJ(EN,SIGN,XM,XERR,NBIN)
endif
WRITE(21,2001) IOBS, XM, XERR
! NBOOT = 40
! CALL BOOTSTRAP( EN,XM_BS,XERR_BS,NBOOT,ISEED)
! WRITE(21,2001) IOBS, XM_BS, XERR_BS
! IF (IOBS == 4) Write(22,"(F14.7,2x,F14.7)") XM/dble(L*L), XERR/dble(L*L)
endif
ENDDO
ENDDO
CLOSE(21)
2001 FORMAT('OBS : ', I4,4x,F12.6,2X, F12.6)
OPEN (UNIT=21, FILE='Var_scalJ', STATUS='unknown')
WRITE(21,*) 'Effective number of bins, and bins: ', Nbins_eff, Nbins
ALLOCATE (EN(Nbins_eff), SIGN(Nbins_eff))
DO IOBS = 1,NOBS
WRITE(21,*)
DO I = 1,Nbins_eff
EN (I) = Real(OBS(I,IOBS), kind(0.d0))
SIGN(I) = Real(OBS(I,NOBS), kind(0.d0))
ENDDO
IF (IOBS.EQ.NOBS ) then
CALL ERRCALCJ(EN, XM,XERR,N_Rebin)
else
CALL ERRCALCJ(EN,SIGN,XM,XERR,N_Rebin)
endif
WRITE(21,2001) IOBS, XM, XERR
ENDDO
CLOSE(21)
2001 FORMAT('OBS : ', I4,4x,F12.6,2X, F12.6)
DEALLOCATE (EN,SIGN,OBS)
DEALLOCATE (EN,SIGN,OBS)
END Program Cov_vec
END Program Cov_vec
......@@ -47,7 +47,7 @@
Integer :: Nunit, Norb
Integer :: no, no1, n, nbins, n_skip, nb, NT, NT1, Lt, N_rebin, N_cov, ierr
Integer :: no, no1, n, nbins, n_skip, nb, NT, NT1, Lt, N_rebin, N_cov, ierr, N_Back
real (Kind=Kind(0.d0)):: X, Y, dtau, X_diag
Complex (Kind=Kind(0.d0)), allocatable :: Xmean(:), Xcov(:,:)
Complex (Kind=Kind(0.d0)) :: Z, Z_diag
......@@ -59,10 +59,12 @@
Real (Kind=Kind(0.d0)), allocatable :: Xk_p(:,:)
Character (len=64) :: File_out
NAMELIST /VAR_errors/ n_skip, N_rebin, N_Cov
NAMELIST /VAR_errors/ n_skip, N_rebin, N_Cov, N_Back
N_Back = 1
OPEN(UNIT=5,FILE='parameters',STATUS='old',ACTION='read',IOSTAT=ierr)
IF (ierr /= 0) THEN
WRITE(*,*) 'unable to open <parameters>',ierr
......@@ -101,6 +103,7 @@
! Allocate space
Allocate ( bins(Nunit,Lt,Nbins), Phase(Nbins), Xk_p(2,Nunit), V_help(lt,Nbins), bins0(Nbins,Norb))
Allocate (Xmean(Lt), Xcov(Lt,Lt))
bins = 0.d0
bins0 = cmplx(0.d0,0.d0,Kind(0.d0))
......@@ -110,7 +113,8 @@
Read(10,*,End=10) Phase(nb-n_skip),no,no1,n, X
Z_diag = cmplx(0.d0,0.d0,kind(0.d0))
Do no = 1,Norb
Read(10,*) bins0(nb-n_skip,no)
Read(10,*) Z
If ( N_back == 1 ) bins0(nb-n_skip,no) = Z
Z_diag = Z_diag + bins0(nb-n_skip,no)*bins0(nb-n_skip,no)
Enddo
do n = 1,Nunit
......
-2.094395 -1.209200
1 1 0.24568282 0.00110286 0.00000000 0.00000000
1 2 -0.09259233 0.00045271 -0.05393494 0.00048046
2 1 -0.09259233 0.00045271 0.05393494 0.00048046
2 2 0.24520403 0.00158845 0.00000000 0.00000000
1 1 0.24505594 0.00039439 0.00000000 0.00000000
1 2 -0.09288481 0.00028114 -0.05357951 0.00018946
2 1 -0.09288481 0.00028114 0.05357951 0.00018946
2 2 0.24504256 0.00030479 0.00000000 0.00000000
-2.094395 1.209200
1 1 0.24563517 0.00102974 -0.00000000 0.00000000
1 2 -0.09326569 0.00069485 0.05367292 0.00021059
2 1 -0.09326569 0.00069485 -0.05367292 0.00021059
2 2 0.24553229 0.00119628 -0.00000000 0.00000000
1 1 0.24512725 0.00041733 -0.00000000 0.00000000
1 2 -0.09259853 0.00034718 0.05391398 0.00026912
2 1 -0.09259853 0.00034718 -0.05391398 0.00026912
2 2 0.24511178 0.00026040 -0.00000000 0.00000000
0.000000 -2.418399
1 1 0.24558600 0.00097753 0.00000000 0.00000000
1 2 0.00021017 0.00067370 -0.10756730 0.00056005
2 1 0.00021017 0.00067370 0.10756730 0.00056005
2 2 0.24565525 0.00104289 0.00000000 0.00000000
1 1 0.24511123 0.00052525 0.00000000 0.00000000
1 2 0.00033002 0.00059236 -0.10761237 0.00046790
2 1 0.00033002 0.00059236 0.10761237 0.00046790
2 2 0.24497686 0.00024583 0.00000000 0.00000000
0.000000 0.000000
1 1 0.20739435 0.00053608 0.00000000 0.00000000
1 2 -0.20717016 0.00066048 0.00000000 0.00000000
2 1 -0.20717016 0.00066048 0.00000000 0.00000000
2 2 0.20743329 0.00061804 0.00000000 0.00000000
1 1 0.20711030 0.00025534 0.00000000 0.00000000
1 2 -0.20676033 0.00024901 0.00000000 0.00000000
2 1 -0.20676033 0.00024901 0.00000000 0.00000000
2 2 0.20709055 0.00023829 0.00000000 0.00000000
0.000000 2.418399
1 1 0.24558600 0.00097753 -0.00000000 0.00000000
1 2 0.00021017 0.00067370 0.10756730 0.00056005
2 1 0.00021017 0.00067370 -0.10756730 0.00056005
2 2 0.24565525 0.00104289 -0.00000000 0.00000000
1 1 0.24511123 0.00052525 -0.00000000 0.00000000
1 2 0.00033002 0.00059236 0.10761237 0.00046790
2 1 0.00033002 0.00059236 -0.10761237 0.00046790
2 2 0.24497686 0.00024583 -0.00000000 0.00000000
2.094395 -1.209200
1 1 0.24563517 0.00102974 0.00000000 0.00000000
1 2 -0.09326569 0.00069485 -0.05367292 0.00021059
2 1 -0.09326569 0.00069485 0.05367292 0.00021059
2 2 0.24553229 0.00119628 0.00000000 0.00000000
1 1 0.24512725 0.00041733 0.00000000 0.00000000
1 2 -0.09259853 0.00034718 -0.05391398 0.00026912
2 1 -0.09259853 0.00034718 0.05391398 0.00026912
2 2 0.24511178 0.00026040 0.00000000 0.00000000
2.094395 1.209200
1 1 0.24568282 0.00110286 -0.00000000 0.00000000
1 2 -0.09259233 0.00045271 0.05393494 0.00048046
2 1 -0.09259233 0.00045271 -0.05393494 0.00048046
2 2 0.24520403 0.00158845 -0.00000000 0.00000000
1 1 0.24505594 0.00039439 -0.00000000 0.00000000
1 2 -0.09288481 0.00028114 0.05357951 0.00018946
2 1 -0.09288481 0.00028114 -0.05357951 0.00018946
2 2 0.24504256 0.00030479 -0.00000000 0.00000000
2.094395 3.627599
1 1 0.25202851 0.00108825 0.00000000 0.00000000
1 2 0.00018585 0.00052064 -0.00044436 0.00045577
2 1 0.00018585 0.00052064 0.00044436 0.00045577
2 2 0.25187526 0.00118423 0.00000000 0.00000000
1 1 0.25128113 0.00044155 0.00000000 0.00000000
1 2 0.00032748 0.00048702 0.00015958 0.00018628
2 1 0.00032748 0.00048702 -0.00015958 0.00018628
2 2 0.25112838 0.00027259 0.00000000 0.00000000
4.188790 0.000000
1 1 0.25202851 0.00108825 0.00000000 0.00000000
1 2 0.00018585 0.00052064 0.00044436 0.00045577
2 1 0.00018585 0.00052064 -0.00044436 0.00045577
2 2 0.25187526 0.00118423 0.00000000 0.00000000
1 1 0.25128113 0.00044155 0.00000000 0.00000000
1 2 0.00032748 0.00048702 -0.00015958 0.00018628
2 1 0.00032748 0.00048702 0.00015958 0.00018628
2 2 0.25112838 0.00027259 0.00000000 0.00000000
-1.000000 0.000000
1 1 -0.00497557 0.00009838 -0.00000000 0.00000000
1 2 -0.00237705 0.00002613 -0.00000000 0.00000000
2 1 -0.00230687 0.00005006 -0.00000000 0.00000000
2 2 -0.00487420 0.00014059 -0.00000000 0.00000000
1 1 -0.00490351 0.00005443 -0.00000000 0.00000000
1 2 -0.00236081 0.00003342 -0.00000000 0.00000000
2 1 -0.00229350 0.00005362 -0.00000000 0.00000000
2 2 -0.00491538 0.00001911 -0.00000000 0.00000000
-0.500000 0.866025
1 1 -0.00494330 0.00004269 -0.00000000 0.00000000
1 2 -0.06422131 0.00031014 -0.00000000 0.00000000
2 1 -0.00233095 0.00003103 -0.00000000 0.00000000
2 2 -0.00502461 0.00008337 -0.00000000 0.00000000
1 1 -0.00492194 0.00004018 -0.00000000 0.00000000
1 2 -0.06451524 0.00029840 -0.00000000 0.00000000
2 1 -0.00228230 0.00002736 -0.00000000 0.00000000
2 2 -0.00489348 0.00002183 -0.00000000 0.00000000
0.000000 1.732051
1 1 -0.00282807 0.00007084 0.00000000 0.00000000
1 2 -0.00235787 0.00004124 0.00000000 0.00000000
2 1 -0.00234226 0.00008886 0.00000000 0.00000000
2 2 -0.00280086 0.00010439 0.00000000 0.00000000
1 1 -0.00284687 0.00005022 0.00000000 0.00000000
1 2 -0.00230513 0.00001011 0.00000000 0.00000000
2 1 -0.00235089 0.00006509 0.00000000 0.00000000
2 2 -0.00286488 0.00000827 0.00000000 0.00000000
-0.500000 -0.866025
1 1 -0.00495918 0.00006878 -0.00000000 0.00000000
1 2 -0.00233394 0.00005916 -0.00000000 0.00000000
2 1 -0.06466722 0.00055844 -0.00000000 0.00000000
2 2 -0.00491519 0.00011044 -0.00000000 0.00000000
1 1 -0.00489817 0.00004468 -0.00000000 0.00000000
1 2 -0.00231266 0.00004615 -0.00000000 0.00000000
2 1 -0.06429402 0.00024020 -0.00000000 0.00000000
2 2 -0.00487041 0.00001864 -0.00000000 0.00000000
0.000000 0.000000
1 1 0.24280659 0.00098533 0.00000000 0.00000000
1 2 -0.06423268 0.00033782 0.00000000 0.00000000
2 1 -0.06423268 0.00033782 0.00000000 0.00000000
2 2 0.24266300 0.00117256 0.00000000 0.00000000
1 1 0.24225127 0.00040860 0.00000000 0.00000000
1 2 -0.06404578 0.00032313 0.00000000 0.00000000
2 1 -0.06404578 0.00032313 0.00000000 0.00000000
2 2 0.24217886 0.00026353 0.00000000 0.00000000
0.500000 0.866025
1 1 -0.00495918 0.00006878 -0.00000000 0.00000000
1 2 -0.06466722 0.00055844 -0.00000000 0.00000000
2 1 -0.00233394 0.00005916 -0.00000000 0.00000000
2 2 -0.00491519 0.00011044 -0.00000000 0.00000000
1 1 -0.00489817 0.00004468 -0.00000000 0.00000000
1 2 -0.06429402 0.00024020 -0.00000000 0.00000000
2 1 -0.00231266 0.00004615 -0.00000000 0.00000000
2 2 -0.00487041 0.00001864 -0.00000000 0.00000000
0.500000 -0.866025
1 1 -0.00494330 0.00004269 -0.00000000 0.00000000
1 2 -0.00233095 0.00003103 -0.00000000 0.00000000
2 1 -0.06422131 0.00031014 -0.00000000 0.00000000
2 2 -0.00502461 0.00008337 -0.00000000 0.00000000
1 1 -0.00492194 0.00004018 -0.00000000 0.00000000
1 2 -0.00228230 0.00002736 -0.00000000 0.00000000
2 1 -0.06451524 0.00029840 -0.00000000 0.00000000
2 2 -0.00489348 0.00002183 -0.00000000 0.00000000
1.000000 0.000000
1 1 -0.00497557 0.00009838 -0.00000000 0.00000000
1 2 -0.00230687 0.00005006 -0.00000000 0.00000000
2 1 -0.00237705 0.00002613 -0.00000000 0.00000000
2 2 -0.00487420 0.00014059 -0.00000000 0.00000000
1 1 -0.00490351 0.00005443 -0.00000000 0.00000000
1 2 -0.00229350 0.00005362 -0.00000000 0.00000000
2 1 -0.00236081 0.00003342 -0.00000000 0.00000000
2 2 -0.00491538 0.00001911 -0.00000000 0.00000000
1.500000 0.866025
1 1 -0.00282807 0.00007084 0.00000000 0.00000000
1 2 -0.00234226 0.00008886 0.00000000 0.00000000
2 1 -0.00235787 0.00004124 0.00000000 0.00000000
2 2 -0.00280086 0.00010439 0.00000000 0.00000000
1 1 -0.00284687 0.00005022 0.00000000 0.00000000
1 2 -0.00235089 0.00006509 0.00000000 0.00000000
2 1 -0.00230513 0.00001011 0.00000000 0.00000000
2 2 -0.00286488 0.00000827 0.00000000 0.00000000
This diff could not be displayed because it is too large.
2 ( -14.157725232625806 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.137815778863374 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.191926052350734 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.132634361342275 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.168006494084787 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.253161036462783 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.193584748486801 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.204304216839628 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.209642735886487 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.172620276779236 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.184070269041591 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.178898531144977 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.156845673760287 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.178894467470176 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.144145888120262 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.122879179406118 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.154783062201250 , 0.0000000000000000 ) 1.0000000000000000
2 ( -14.160323969985406 , 0.0000000000000000 ) 1.0000000000000000
Effective number of bins, and bins: 8 8
Effective number of bins, and bins: 9 10
OBS : 1 -14.179895 0.017149
OBS : 1 -14.161642 0.009650
OBS : 2 1.000000 0.000000
-2.094395 -1.209200
1 1 1.00000000 0.00000000 -0.00098062 0.00150455
1 2 0.72680964 0.00166193 0.42242686 0.00118591
2 1 0.72680964 0.00166193 -0.42242686 0.00118591
2 2 1.00000000 0.00000004 0.00154110 0.00096284
1 1 1.00000000 0.00000000 -0.00127207 0.00091531
1 2 0.72627353 0.00047739 0.41994328 0.00163322
2 1 0.72627353 0.00047739 -0.41994328 0.00163322
2 2 1.00000000 0.00000000 0.00052974 0.00078660
-2.094395 1.209200
1 1 1.00000000 0.00000000 0.00111106 0.00036854
1 2 0.72816073 0.00339911 -0.41936437 0.00099926
2 1 0.72816073 0.00339911 0.41936437 0.00099926
2 2 1.00000000 0.00000003 0.00034859 0.00039618
1 1 1.00000000 0.00000002 0.00185057 0.00095707
1 2 0.72643406 0.00054052 -0.42087897 0.00196433
2 1 0.72643406 0.00054052 0.42087897 0.00196433
2 2 1.00000000 0.00000003 -0.00012084 0.00039918
0.000000 -2.418399
1 1 1.00000000 0.00000000 -0.00009931 0.00171276
1 2 -0.00115008 0.00078257 0.84085773 0.00186779
2 1 -0.00115008 0.00078257 -0.84085773 0.00186779
2 2 1.00000000 0.00000000 -0.00080136 0.00090838
1 1 1.00000000 0.00000000 -0.00180229 0.00139738
1 2 -0.00143245 0.00186284 0.83876996 0.00102821
2 1 -0.00143245 0.00186284 -0.83876996 0.00102821
2 2 1.00000000 0.00000000 0.00027805 0.00081248
0.000000 0.000000
1 1 1.00000000 0.00000003 0.00000000 0.00000000
1 2 0.90774471 0.00134678 0.00000000 0.00000000
2 1 0.90774471 0.00134678 0.00000000 0.00000000
2 2 1.00000000 0.00000003 0.00000000 0.00000000
1 1 1.00000000 0.00000000 0.00000000 0.00000000
1 2 0.90687940 0.00026404 0.00000000 0.00000000
2 1 0.90687940 0.00026404 0.00000000 0.00000000
2 2 1.00000000 0.00000000 0.00000000 0.00000000
0.000000 2.418399
1 1 1.00000000 0.00000000 0.00009931 0.00171276
1 2 -0.00115008 0.00078257 -0.84085773 0.00186779
2 1 -0.00115008 0.00078257 0.84085773 0.00186779
2 2 1.00000000 0.00000000 0.00080136 0.00090838
1 1 1.00000000 0.00000000 0.00180229 0.00139738
1 2 -0.00143245 0.00186284 -0.83876996 0.00102821
2 1 -0.00143245 0.00186284 0.83876996 0.00102821
2 2 1.00000000 0.00000000 -0.00027805 0.00081248
2.094395 -1.209200
1 1 1.00000000 0.00000000 -0.00111106 0.00036854
1 2 0.72816073 0.00339911 0.41936437 0.00099926
2 1 0.72816073 0.00339911 -0.41936437 0.00099926
2 2 1.00000000 0.00000003 -0.00034859 0.00039618
1 1 1.00000000 0.00000002 -0.00185057 0.00095707
1 2 0.72643406 0.00054052 0.42087897 0.00196433
2 1 0.72643406 0.00054052 -0.42087897 0.00196433
2 2 1.00000000 0.00000003 0.00012084 0.00039918
2.094395 1.209200
1 1 1.00000000 0.00000000 0.00098062 0.00150455
1 2 0.72680964 0.00166193 -0.42242686 0.00118591
2 1 0.72680964 0.00166193 0.42242686 0.00118591
2 2 1.00000000 0.00000004 -0.00154110 0.00096284
1 1 1.00000000 0.00000000 0.00127207 0.00091531
1 2 0.72627353 0.00047739 -0.41994328 0.00163322
2 1 0.72627353 0.00047739 0.41994328 0.00163322
2 2 1.00000000 0.00000000 -0.00052974 0.00078660
2.094395 3.627599
1 1 1.00000000 0.00000000 0.00813331 0.00688072
1 2 -0.00217098 0.00933782 0.00262890 0.00518360
2 1 -0.00217098 0.00933782 -0.00262890 0.00518360
2 2 1.00000000 0.00000000 -0.00449365 0.00165944
1 1 1.00000000 0.00000000 0.00107634 0.00323507
1 2 -0.00340249 0.00563291 -0.00178117 0.00266545
2 1 -0.00340249 0.00563291 0.00178117 0.00266545
2 2 1.00000000 0.00000000 -0.00397774 0.00211691
4.188790 0.000000
1 1 1.00000000 0.00000000 -0.00813331 0.00688072
1 2 -0.00217098 0.00933782 -0.00262890 0.00518360
2 1 -0.00217098 0.00933782 0.00262890 0.00518360
2 2 1.00000000 0.00000000 0.00449365 0.00165944
1 1 1.00000000 0.00000000 -0.00107634 0.00323507
1 2 -0.00340249 0.00563291 0.00178117 0.00266545
2 1 -0.00340249 0.00563291 -0.00178117 0.00266545
2 2 1.00000000 0.00000000 0.00397774 0.00211691
-1.000000 0.000000
1 1 0.00154015 0.00166259 -0.00000000 0.00000000
1 2 -0.06090065 0.00140706 0.00000000 0.00000000
2 1 -0.06073376 0.00105812 0.00000000 0.00000000
2 2 -0.00122847 0.00057624 -0.00000000 0.00000000
1 1 0.00009581 0.00083417 -0.00000000 0.00000000
1 2 -0.06075056 0.00037534 0.00000000 0.00000000
2 1 -0.06042513 0.00057975 0.00000000 0.00000000
2 2 -0.00084421 0.00028027 -0.00000000 0.00000000
-0.500000 0.866025
1 1 -0.00179819 0.00111215 0.00000000 0.00000000
1 2 0.42385996 0.00111717 -0.00000000 0.00000000
2 1 -0.06018789 0.00162760 -0.00000000 0.00000000
2 2 0.00064350 0.00021914 0.00000000 0.00000000
1 1 -0.00091014 0.00069719 0.00000000 0.00000000
1 2 0.42474315 0.00169509 -0.00000000 0.00000000
2 1 -0.06078152 0.00077933 -0.00000000 0.00000000
2 2 0.00084228 0.00059778 0.00000000 0.00000000
0.000000 1.732051
1 1 -0.00038343 0.00051707 -0.00000000 0.00000000
1 2 -0.06097785 0.00264042 -0.00000000 0.00000000
2 1 -0.06133716 0.00211018 0.00000000 0.00000000
2 2 0.00038372 0.00027437 -0.00000000 0.00000000
1 1 -0.00025410 0.00018802 -0.00000000 0.00000000
1 2 -0.06084956 0.00065091 -0.00000000 0.00000000
2 1 -0.06163949 0.00136384 0.00000000 0.00000000
2 2 0.00007169 0.00028834 -0.00000000 0.00000000
-0.500000 -0.866025
1 1 -0.00135743 0.00125881 -0.00000000 0.00000000
1 2 -0.06133877 0.00057514 -0.00000000 0.00000000
2 1 0.42591156 0.00274497 -0.00000000 0.00000000
2 2 0.00072244 0.00025050 -0.00000000 0.00000000
1 1 0.00038452 0.00059784 -0.00000000 0.00000000
1 2 -0.05986236 0.00070299 -0.00000000 0.00000000
2 1 0.42393101 0.00118249 -0.00000000 0.00000000
2 2 0.00061006 0.00037241 -0.00000000 0.00000000
0.000000 0.000000
1 1 1.00000000 0.00000002 -0.00000000 0.00000000
1 2 0.42344926 0.00197215 -0.00000000 0.00000000
2 1 0.42344926 0.00197215 0.00000000 0.00000000
2 2 1.00000000 0.00000003 -0.00000000 0.00000000
1 1 1.00000000 0.00000000 0.00000000 0.00000000
1 2 0.42251386 0.00175373 -0.00000000 0.00000000
2 1 0.42251386 0.00175373 0.00000000 0.00000000
2 2 1.00000000 0.00000000 -0.00000000 0.00000000
0.500000 0.866025
1 1 0.00135743 0.00125881 0.00000000 0.00000000
1 2 0.42591156 0.00274497 0.00000000 0.00000000
2 1 -0.06133877 0.00057514 0.00000000 0.00000000
2 2 -0.00072244 0.00025050 0.00000000 0.00000000
1 1 -0.00038452 0.00059784 0.00000000 0.00000000
1 2 0.42393101 0.00118249 0.00000000 0.00000000
2 1 -0.05986236 0.00070299 0.00000000 0.00000000
2 2 -0.00061006 0.00037241 0.00000000 0.00000000
0.500000 -0.866025
1 1 0.00179819 0.00111215 -0.00000000 0.00000000
1 2 -0.06018789 0.00162760 0.00000000 0.00000000
2 1 0.42385996 0.00111717 0.00000000 0.00000000
2 2 -0.00064350 0.00021914 -0.00000000 0.00000000
1 1 0.00091014 0.00069719 -0.00000000 0.00000000
1 2 -0.06078152 0.00077933 0.00000000 0.00000000
2 1 0.42474315 0.00169509 0.00000000 0.00000000
2 2 -0.00084228 0.00059778 -0.00000000 0.00000000
1.000000 0.000000
1 1 -0.00154015 0.00166259 0.00000000 0.00000000
1 2 -0.06073376 0.00105812 -0.00000000 0.00000000
2 1 -0.06090065 0.00140706 -0.00000000 0.00000000
2 2 0.00122847 0.00057624 0.00000000 0.00000000
1 1 -0.00009581 0.00083417 0.00000000 0.00000000
1 2 -0.06042513 0.00057975 -0.00000000 0.00000000
2 1 -0.06075056 0.00037534 -0.00000000 0.00000000
2 2 0.00084421 0.00028027 0.00000000 0.00000000
1.500000 0.866025
1 1 0.00038343 0.00051707 -0.00000000 0.00000000
1 2 -0.06133716 0.00211018 -0.00000000 0.00000000
2 1 -0.06097785 0.00264042 0.00000000 0.00000000
2 2 -0.00038372 0.00027437 -0.00000000 0.00000000
1 1 0.00025410 0.00018802 -0.00000000 0.00000000
1 2 -0.06163949 0.00136384 -0.00000000 0.00000000
2 1 -0.06084956 0.00065091 0.00000000 0.00000000
2 2 -0.00007169 0.00028834 -0.00000000 0.00000000
This diff could not be displayed because it is too large.
2 ( -22.874856040464330 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.945061162244045 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.977849218799857 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.841062674072216 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.834637764599300 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.971523087238108 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.937710348912237 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.915004541949468 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.977478818931498 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.917479697950039 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.901644286641201 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.823374225681928 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.923820039274261 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.944946167575079 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.864521846724717 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.863983057715902 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.811305649134695 , 0.0000000000000000 ) 1.0000000000000000
2 ( -22.938851890058952 , 0.0000000000000000 ) 1.0000000000000000
Effective number of bins, and bins: 8 8
Effective number of bins, and bins: 9 10
OBS : 1 -22.912213 0.005736
OBS : 1 -22.881384 0.019186
OBS : 2 1.000000 0.000000
2 ( 17.999999999999968 , 0.0000000000000000 ) 1.0000000000000000
2 ( 17.999999999999787 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000142 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000025 , 0.0000000000000000 ) 1.0000000000000000
2 ( 17.999999999999613 , 0.0000000000000000 ) 1.0000000000000000
2 ( 17.999999999999872 , 0.0000000000000000 ) 1.0000000000000000
2 ( 17.999999999993339 , 0.0000000000000000 ) 1.0000000000000000
2 ( 17.999999999999858 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000000 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000018 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000039 , 0.0000000000000000 ) 1.0000000000000000
2 ( 17.999999999999989 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000014 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000014 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000014 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000011 , 0.0000000000000000 ) 1.0000000000000000
2 ( 18.000000000000011 , 0.0000000000000000 ) 1.0000000000000000
2 ( 17.999999999999982 , 0.0000000000000000 ) 1.0000000000000000
Effective number of bins, and bins: 8 8
Effective number of bins, and bins: 9 10
OBS : 1 18.000000 0.000000
......
2 ( 8.7171308078382292 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.8072453833804705 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7859231664490220 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7084283127301543 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.6666312705144648 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7183620507757400 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7441256004251997 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7107003251103254 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7678360830449993 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7448594211708421 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7175740176001533 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.6444756945372578 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7669743655140984 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7660517001049563 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7203759586044658 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7411038783095893 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.6565225869337237 , 0.0000000000000000 ) 1.0000000000000000
2 ( 8.7785279200731505 , 0.0000000000000000 ) 1.0000000000000000
Effective number of bins, and bins: 8 8
Effective number of bins, and bins: 9 10
OBS : 1 8.732318 0.017393
OBS : 1 8.719742 0.012097
OBS : 2 1.000000 0.000000
-2.094395 -1.209200
1 1 0.66474319 0.00899187 0.00000000 0.00000000
1 2 -0.28640589 0.00236890 -0.17808202 0.00557894
2 1 -0.28640589 0.00236890 0.17808202 0.00557894
2 2 0.66632806 0.00427310 0.00000000 0.00000000
1 1 0.67021820 0.00724193 0.00000000 0.00000000
1 2 -0.28497210 0.00629375 -0.17614094 0.00945168
2 1 -0.28497210 0.00629375 0.17614094 0.00945168
2 2 0.65019216 0.00829285 0.00000000 0.00000000
-2.094395 1.209200
1 1 0.66733232 0.00784675 -0.00000000 0.00000000
1 2 -0.29884303 0.00944811 0.17879555 0.00363852
2 1 -0.29884303 0.00944811 -0.17879555 0.00363852
2 2 0.66528493 0.00852377 -0.00000000 0.00000000
1 1 0.65333887 0.00692719 -0.00000000 0.00000000
1 2 -0.28043028 0.00434339 0.17421023 0.00267363
2 1 -0.28043028 0.00434339 -0.17421023 0.00267363
2 2 0.66886921 0.00734977 -0.00000000 0.00000000
0.000000 -2.418399
1 1 0.66628824 0.00235470 0.00000000 0.00000000
1 2 0.01185691 0.00919473 -0.32905521 0.00932207
2 1 0.01185691 0.00919473 0.32905521 0.00932207
2 2 0.66230644 0.01757582 0.00000000 0.00000000
1 1 0.65887748 0.00849640 0.00000000 0.00000000
1 2 0.00279780 0.00472879 -0.32732682 0.00587355
2 1 0.00279780 0.00472879 0.32732682 0.00587355
2 2 0.65945228 0.00443599 0.00000000 0.00000000
0.000000 0.000000
1 1 1.56994810 0.01845290 0.00000000 0.00000000
1 2 -1.53724864 0.01055898 0.00000000 0.00000000
2 1 -1.53724864 0.01055898 0.00000000 0.00000000
2 2 1.59007643 0.00478747 0.00000000 0.00000000
1 1 1.59990121 0.01598345 0.00000000 0.00000000
1 2 -1.55970816 0.02115567 0.00000000 0.00000000
2 1 -1.55970816 0.02115567 0.00000000 0.00000000
2 2 1.62088207 0.02421859 0.00000000 0.00000000
0.000000 2.418399
1 1 0.66628824 0.00235470 -0.00000000 0.00000000
1 2 0.01185691 0.00919473 0.32905521 0.00932207
2 1 0.01185691 0.00919473 -0.32905521 0.00932207
2 2 0.66230644 0.01757582 -0.00000000 0.00000000
1 1 0.65887748 0.00849640 -0.00000000 0.00000000
1 2 0.00279780 0.00472879 0.32732682 0.00587355
2 1 0.00279780 0.00472879 -0.32732682 0.00587355
2 2 0.65945228 0.00443599 -0.00000000 0.00000000
2.094395 -1.209200
1 1 0.66733232 0.00784675 0.00000000 0.00000000
1 2 -0.29884303 0.00944811 -0.17879555 0.00363852
2 1 -0.29884303 0.00944811 0.17879555 0.00363852
2 2 0.66528493 0.00852377 0.00000000 0.00000000
1 1 0.65333887 0.00692719 0.00000000 0.00000000
1 2 -0.28043028 0.00434339 -0.17421023 0.00267363
2 1 -0.28043028 0.00434339 0.17421023 0.00267363
2 2 0.66886921 0.00734977 0.00000000 0.00000000
2.094395 1.209200
1 1 0.66474319 0.00899187 -0.00000000 0.00000000
1 2 -0.28640589 0.00236890 0.17808202 0.00557894
2 1 -0.28640589 0.00236890 -0.17808202 0.00557894
2 2 0.66632806 0.00427310 -0.00000000 0.00000000
1 1 0.67021820 0.00724193 -0.00000000 0.00000000
1 2 -0.28497210 0.00629375 0.17614094 0.00945168
2 1 -0.28497210 0.00629375 -0.17614094 0.00945168
2 2 0.65019216 0.00829285 -0.00000000 0.00000000
2.094395 3.627599
1 1 0.62403253 0.00838723 0.00000000 0.00000000
1 2 -0.00165490 0.00334731 0.00078752 0.00736777
2 1 -0.00165490 0.00334731 -0.00078752 0.00736777
2 2 0.61905887 0.00350114 0.00000000 0.00000000
1 1 0.62748416 0.00801712 0.00000000 0.00000000
1 2 -0.00761816 0.00381029 -0.00200590 0.00806226
2 1 -0.00761816 0.00381029 0.00200590 0.00806226
2 2 0.62124046 0.00442354 0.00000000 0.00000000
4.188790 0.000000
1 1 0.62403253 0.00838723 0.00000000 0.00000000
1 2 -0.00165490 0.00334731 -0.00078752 0.00736777
2 1 -0.00165490 0.00334731 0.00078752 0.00736777
2 2 0.61905887 0.00350114 0.00000000 0.00000000
1 1 0.62748416 0.00801712 0.00000000 0.00000000
1 2 -0.00761816 0.00381029 0.00200590 0.00806226
2 1 -0.00761816 0.00381029 -0.00200590 0.00806226
2 2 0.62124046 0.00442354 0.00000000 0.00000000
-1.000000 0.000000
1 1 0.10515739 0.00373910 -0.00000000 0.00000000
1 2 -0.10294476 0.00084514 -0.00000000 0.00000000
2 1 -0.10297324 0.00248061 0.00000000 0.00000000
2 2 0.10711305 0.00516339 -0.00000000 0.00000000
1 1 0.10740166 0.00028833 -0.00000000 0.00000000
1 2 -0.10902470 0.00430658 -0.00000000 0.00000000
2 1 -0.10899575 0.00356211 0.00000000 0.00000000
2 2 0.11105387 0.00451528 -0.00000000 0.00000000
-0.500000 0.866025
1 1 0.10464238 0.00416285 0.00000000 0.00000000
1 2 -0.30026742 0.00241615 0.00000000 0.00000000
2 1 -0.10449245 0.00334322 0.00000000 0.00000000
2 2 0.10845359 0.00212339 -0.00000000 0.00000000
1 1 0.11118190 0.00508196 -0.00000000 0.00000000
1 2 -0.30106834 0.00336927 0.00000000 0.00000000