Euclidean gamma matrix tests added

pyerrors/npr.py

@@ -134,10 +134,7 @@ def Zq(prop, fermion='Wilson'):
     sin_mom = np.sin(2 * np.pi / L * mom)
 
     if fermion == 'Wilson':
-        p_slash = -1j * (sin_mom[0] * gamma[0]
+        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)
-                   + sin_mom[1] * gamma[1]
-                   + sin_mom[2] * gamma[2]
-                   + sin_mom[3] * gamma[3]) / np.sum(sin_mom ** 2)
     else:
         raise Exception("Fermion type '" + fermion + "' not implemented")
 

tests/test_npr.py

@@ -0,0 +1,12 @@
+import numpy as np
+import pyerrors as pe
+import pytest
+
+np.random.seed(0)
+
+
+def test_gamma_matrices():
+    for matrix in pe.npr.gamma:
+        assert np.allclose(matrix @ matrix, np.identity(4))
+        assert np.allclose(matrix, matrix.T.conj())
+    assert np.allclose(pe.npr.gamma5, pe.npr.gamma[0] @ pe.npr.gamma[1] @ pe.npr.gamma[2] @ pe.npr.gamma[3])