From 5bdfa3f8bedb33e1a4e8f04a2de8942152eeaadc Mon Sep 17 00:00:00 2001 From: Fabian Joswig Date: Wed, 27 Oct 2021 18:16:59 +0100 Subject: [PATCH] DWF p_slash expression added to npr --- pyerrors/npr.py | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/pyerrors/npr.py b/pyerrors/npr.py index 7ff81cc0..6594f03f 100644 --- a/pyerrors/npr.py +++ b/pyerrors/npr.py @@ -85,6 +85,11 @@ def Zq(inv_prop, fermion='Wilson'): if fermion == 'Wilson': p_slash = -1j * (sin_mom[0] * gamma[0] + sin_mom[1] * gamma[1] + sin_mom[2] * gamma[2] + sin_mom[3] * gamma[3]) / np.sum(sin_mom ** 2) + elif fermion == 'DWF': + W = np.sum(1 - np.cos(2 * np.pi / L * mom)) + s2 = np.sum(sin_mom ** 2) + p_slash = -1j * (sin_mom[0] * gamma[0] + sin_mom[1] * gamma[1] + sin_mom[2] * gamma[2] + sin_mom[3] * gamma[3]) + p_slash /= 2 * (W - 1 + np.sqrt((1 - W) ** 2 + s2)) else: raise Exception("Fermion type '" + fermion + "' not implemented")