lattice/pyerrors
1
0
Fork 0
mirror of https://github.com/fjosw/pyerrors.git synced 2025-03-18 08:10:24 +01: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:49:11 +00:00
commit 53cfbecb59
2 changed files with 26 additions and 28 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

46
examples/06_gevp.ipynb
View file

File diff suppressed because one or more lines are too long

8
pyerrors/correlators.py
View file

@ -259,16 +259,16 @@ class Corr:
"Eigenvector" - Use the method described in arXiv:2004.10472 [hep-lat] to find the set of v(t) belonging to the state. "Eigenvector" - Use the method described in arXiv:2004.10472 [hep-lat] to find the set of v(t) belonging to the state.
The reference state is identified by its eigenvalue at t=ts The reference state is identified by its eigenvalue at t=ts
""" """
symmetric_corr = self.matrix_symmetric()
if sorted_list is None: if sorted_list is None:
if (ts is None): if (ts is None):
raise Exception("ts is required if sorted_list=None") raise Exception("ts is required if sorted_list=None")
if (self.content[t0] is None) or (self.content[ts] is None): if (self.content[t0] is None) or (self.content[ts] is None):
raise Exception("Corr not defined at t0/ts") raise Exception("Corr not defined at t0/ts")
G0, Gt = np.empty([self.N, self.N], dtype="double"), np.empty([self.N, self.N], dtype="double") G0, Gt = np.empty([self.N, self.N], dtype="double"), np.empty([self.N, self.N], dtype="double")
symmetric_corr = self.matrix_symmetric()
for i in range(self.N): for i in range(self.N):
for j in range(self.N): for j in range(self.N):
G0[i, j] = symmetric_corr.content[t0][i, j].value G0[i, j] = symmetric_corr[t0][i, j].value
Gt[i, j] = symmetric_corr[ts][i, j].value Gt[i, j] = symmetric_corr[ts][i, j].value
sp_vecs = _GEVP_solver(Gt, G0) sp_vecs = _GEVP_solver(Gt, G0)
@ -281,8 +281,8 @@ class Corr:
G0, Gt = np.empty([self.N, self.N], dtype="double"), np.empty([self.N, self.N], dtype="double") G0, Gt = np.empty([self.N, self.N], dtype="double"), np.empty([self.N, self.N], dtype="double")
for i in range(self.N): for i in range(self.N):
for j in range(self.N): for j in range(self.N):
G0[i, j] = self.content[t0][i, j].value G0[i, j] = symmetric_corr[t0][i, j].value
Gt[i, j] = self.content[t][i, j].value Gt[i, j] = symmetric_corr[t][i, j].value
sp_vecs = _GEVP_solver(Gt, G0) sp_vecs = _GEVP_solver(Gt, G0)
if sorted_list == "Eigenvalue": if sorted_list == "Eigenvalue":