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

feat: read_DistillationContraction_hd5 can now handle complex

Browse source 
normalization of triangle diagrams.
This commit is contained in:
Fabian Joswig 2022-12-01 14:23:37 +00:00
parent b4741edf10
commit fc3657be11
No known key found for this signature in database
1 changed files with 7 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

8
pyerrors/input/hadrons.py
View file

@ -172,9 +172,15 @@ def read_DistillationContraction_hd5(path, ens_id, diagrams=["direct"], idl=None
# "DistillationContraction/Metadata/DmfSuffix" contains info about different quarks, irrelevant in the SU(3) case. # "DistillationContraction/Metadata/DmfSuffix" contains info about different quarks, irrelevant in the SU(3) case.
for diagram in diagrams: for diagram in diagrams:
if diagram == "triangle" and "Identity" not in str(identifier):
part = "im"
else:
part = "re"
real_data = np.zeros(Nt) real_data = np.zeros(Nt)
for x0 in range(Nt): for x0 in range(Nt):
raw_data = h5file["DistillationContraction/Correlators/" + diagram + "/" + str(x0)][:]["re"].astype(np.double) raw_data = h5file["DistillationContraction/Correlators/" + diagram + "/" + str(x0)][:][part].astype(np.double)
real_data += np.roll(raw_data, -x0) real_data += np.roll(raw_data, -x0)
real_data /= Nt real_data /= Nt