From 535786a6d41d9966f2f1dbec715d74310879af44 Mon Sep 17 00:00:00 2001 From: Fabian Joswig Date: Tue, 25 Jan 2022 10:49:59 +0000 Subject: [PATCH] fix: hadrons npr input routines can now handle floating point valued momenta --- pyerrors/input/hadrons.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/pyerrors/input/hadrons.py b/pyerrors/input/hadrons.py index 92e4bc40..aa9d487b 100644 --- a/pyerrors/input/hadrons.py +++ b/pyerrors/input/hadrons.py @@ -152,7 +152,7 @@ def read_ExternalLeg_hd5(path, filestem, ens_id, idl=None): raw_data = file['ExternalLeg/corr'][0][0].view('complex') corr_data.append(raw_data) if mom is None: - mom = np.array(str(file['ExternalLeg/info'].attrs['pIn'])[3:-2].strip().split(' '), dtype=int) + mom = np.array(str(file['ExternalLeg/info'].attrs['pIn'])[3:-2].strip().split(' '), dtype=float) file.close() corr_data = np.array(corr_data) @@ -197,9 +197,9 @@ def read_Bilinear_hd5(path, filestem, ens_id, idl=None): raw_data = file['Bilinear/Bilinear_' + str(i) + '/corr'][0][0].view('complex') corr_data[name].append(raw_data) if mom_in is None: - mom_in = np.array(str(file['Bilinear/Bilinear_' + str(i) + '/info'].attrs['pIn'])[3:-2].strip().split(' '), dtype=int) + mom_in = np.array(str(file['Bilinear/Bilinear_' + str(i) + '/info'].attrs['pIn'])[3:-2].strip().split(' '), dtype=float) if mom_out is None: - mom_out = np.array(str(file['Bilinear/Bilinear_' + str(i) + '/info'].attrs['pOut'])[3:-2].strip().split(' '), dtype=int) + mom_out = np.array(str(file['Bilinear/Bilinear_' + str(i) + '/info'].attrs['pOut'])[3:-2].strip().split(' '), dtype=float) file.close() @@ -262,9 +262,9 @@ def read_Fourquark_hd5(path, filestem, ens_id, idl=None, vertices=["VA", "AV"]): raw_data = file[tree + str(i) + '/corr'][0][0].view('complex') corr_data[name].append(raw_data) if mom_in is None: - mom_in = np.array(str(file[tree + str(i) + '/info'].attrs['pIn'])[3:-2].strip().split(' '), dtype=int) + mom_in = np.array(str(file[tree + str(i) + '/info'].attrs['pIn'])[3:-2].strip().split(' '), dtype=float) if mom_out is None: - mom_out = np.array(str(file[tree + str(i) + '/info'].attrs['pOut'])[3:-2].strip().split(' '), dtype=int) + mom_out = np.array(str(file[tree + str(i) + '/info'].attrs['pOut'])[3:-2].strip().split(' '), dtype=float) file.close()