lattice/pyerrors
1
0
Fork 0
mirror of https://github.com/fjosw/pyerrors.git synced 2026-08-04 11:31:21 +02:00
Code Issues Projects Releases Packages Wiki Activity

correct docs, implement failesafe for therm. trajectories

Browse source
This commit is contained in:
Justus Kuhlmann 2026-05-20 09:54:49 +02:00
commit 52a2c67614
Signed by: jkuhl
GPG key ID: 00ED992DD79B85A6
1 changed files with 3 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

4
pyerrors/input/openQCD.py
View file

@ -250,7 +250,7 @@ def _extract_flowed_energy_density(path, prefix, dtr_read, xmin, spatial_extent,
dtr_read : int dtr_read : int
Determines how many trajectories should be skipped Determines how many trajectories should be skipped
when reading the ms.dat files. when reading the ms.dat files.
Corresponds to dtr_cnfg / dtr_ms (dncnfg / dnms) in the openQCD input file. Corresponds to dtr_ms (dnms) in the openQCD input file.
xmin : int xmin : int
First timeslice where the boundary First timeslice where the boundary
effects have sufficiently decayed. effects have sufficiently decayed.
@ -377,6 +377,8 @@ def _extract_flowed_energy_density(path, prefix, dtr_read, xmin, spatial_extent,
for current in range(0, len(item), tmax)]) for current in range(0, len(item), tmax)])
diffmeas = configlist[-1][-1] - configlist[-1][-2] diffmeas = configlist[-1][-1] - configlist[-1][-2]
if not np.isclose(configlist[-1][0]/diffmeas, int(configlist[-1][0]/diffmeas), 1e-12):
raise ValueError("Spacing of the first configuration after thermalization is irregular.")
configlist[-1] = [item // diffmeas for item in configlist[-1]] configlist[-1] = [item // diffmeas for item in configlist[-1]]
if kwargs.get('assume_thermalization', True) and configlist[-1][0] > 1: if kwargs.get('assume_thermalization', True) and configlist[-1][0] > 1:
warnings.warn('Assume thermalization and that the first measurement belongs to the first config.') warnings.warn('Assume thermalization and that the first measurement belongs to the first config.')