latqcdtools.physics.diracFreespectra

_momenta_1d(L: int, antiperiodic: bool = False) -> numpy.ndarray:
'''
Return exactly L momentum modes.
n = -L//2, ..., L//2-1
p = 2π (n + shift)/L, shift=1/2 for antiperiodic
Works for even/odd L.
'''

class DiracOp(Lx, Ly, Lz, Lt, fermion='Wilson', bc_t='anti'):

class GammaMatrix():
'''The 4x4 gamma matrices used in Euclidean QFT.'''