A method for calculating lattice quantum chromodynamics (QCD) that includes loading an aligned spinor and gauge matrix, and constructing a half-spinor for plus and minus from the aligned spinors. In a following step, the method may include calculating a half spinor minus to be transferred to a positive neighbor of the gauge matrix by multiplying the gauge matrix and half-spinor for minus. The half-spinor for plus (Hp) and the half spinor minus may be stored in a shared memory. In a following step, the method may continue with calculating a half spinor plus to be transferred to a negative neighbor of the gauge matrix by multiplying the gauge matrix and the half-spinor for plus. Thereafter, the sum of half-spinors may be calculated.