lattice/pyerrors
1
0
Fork 0
mirror of https://github.com/fjosw/pyerrors.git synced 2025-05-15 20:13:41 +02:00
Code Issues Projects Releases Packages Wiki Activity

Merge branch 'develop' into documentation

Browse source
This commit is contained in:
fjosw 2022-05-03 09:27:51 +00:00
commit b14474be35
2 changed files with 4 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

6
pyerrors/correlators.py
View file

@ -241,7 +241,7 @@ class Corr:
if self.N == 1: if self.N == 1:
raise Exception("Trying to symmetrize a correlator matrix, that already has N=1.") raise Exception("Trying to symmetrize a correlator matrix, that already has N=1.")
def GEVP(self, t0, ts=None, state=0, sorted_list=None): def GEVP(self, t0, ts=None, state=0, sorted_list="Eigenvalue"):
"""Solve the general eigenvalue problem on the current correlator """Solve the general eigenvalue problem on the current correlator
Parameters Parameters
@ -252,7 +252,7 @@ class Corr:
fixed time G(t_s)v= lambda G(t_0)v if return_list=False fixed time G(t_s)v= lambda G(t_0)v if return_list=False
If return_list=True and sorting=Eigenvector it gives a reference point for the sorting method. If return_list=True and sorting=Eigenvector it gives a reference point for the sorting method.
state : int state : int
The state one is interested in ordered by energy. The lowest state is zero. The state one is interested in, ordered by energy. The lowest state is zero.
sorted_list : string sorted_list : string
if this argument is set, a list of vectors (len=self.T) is returned. If it is left as None, only one vector is returned. if this argument is set, a list of vectors (len=self.T) is returned. If it is left as None, only one vector is returned.
"Eigenvalue" - The eigenvector is chosen according to which eigenvalue it belongs individually on every timeslice. "Eigenvalue" - The eigenvector is chosen according to which eigenvalue it belongs individually on every timeslice.
@ -1170,7 +1170,7 @@ class Corr:
evecs = [] evecs = []
for i in range(Ntrunc): for i in range(Ntrunc):
evecs.append(basematrix.GEVP(t0proj, tproj, state=i)) evecs.append(basematrix.GEVP(t0proj, tproj, state=i, sorted_list=None))
tmpmat = np.empty((Ntrunc, Ntrunc), dtype=object) tmpmat = np.empty((Ntrunc, Ntrunc), dtype=object)
rmat = [] rmat = []

2
setup.py
View file

@ -9,5 +9,5 @@ setup(name='pyerrors',
author_email='fabian.joswig@ed.ac.uk', author_email='fabian.joswig@ed.ac.uk',
packages=find_packages(), packages=find_packages(),
python_requires='>=3.6.0', python_requires='>=3.6.0',
install_requires=['numpy>=1.16', 'autograd @ git+https://github.com/HIPS/autograd.git', 'numdifftools', 'matplotlib>=3.3', 'scipy', 'iminuit>=2', 'h5py', 'lxml'] install_requires=['numpy>=1.16', 'autograd>=1.4', 'numdifftools', 'matplotlib>=3.3', 'scipy', 'iminuit>=2', 'h5py', 'lxml']
) )