Innermost data sharing method of lattice quantum chromodynamics calculation

A non-transitory computer readable storage medium that includes a computer readable program for calculating lattice quantum chromodynamics (QCD), wherein the non-transitory computer readable program when executed on a computer causes the computer to perform the steps of calculating lattice quantum chromodynamics (QCD) that includes:

• loading an aligned spinor (W) array and gauge matrix (U) into a hardware processor;
• constructing a half-spinor for plus (Hp) and half-spinor for minus (Hn) from the aligned spinor (W) array;
• calculating a half spinor minus (Hn′) to be transferred to a positive neighbor of the gauge matrix by multiplying the gauge matrix (U) and half-spinor for minus (Hn) before storing in a shared memory;
• storing the half-spinor for plus (Hp) and the half spinor minus (Hn′) in a shared memory;
• calculating a half spinor plus (Hp′) to be transferred to a negative neighbor of the gauge matrix by multiplying the gauge matrix (U) and the half-spinor for plus (Hp);
• and
• calculating the sum of half-spinors from the half spinor minus (Hn′) to be transferred to the positive neighbor of the gauge matrix and the half spinor plus (Hp′) to be transferred to the negative neighbor of the gauge matrix.

A system for calculating lattice quantum chromodynamics (QCD) comprising:

• a loading module for loading an aligned spinor (W) array and gauge matrix (U) into a hardware processor;
• a half-spinor constructor module for constructing a half spinor for plus (Hp) and half-spinor for minus (Hn) from the aligned spinors (W) array;
• a calculation module for calculating a half spinor minus (Hn′) to be transferred to a positive neighbor of the gauge matrix by multiplying the gauge matrix (U) and half-spinor for minus (Hn) before storing in a shared memory;
• a storage module storing the half-spinor for plus (Hp) and minus (Hn′) in the shared memory;
• a calculating module for calculating a half spinor plus (Hp′) to be transferred to a negative neighbor of the gauge matrix by multiplying the gauge matrix (U) and the half-spinor for plus (Hp);
• and
• a calculating module for calculating the sum of half-spinors from the half spinor minus (Hn′) to be transferred to the positive neighbor of the gauge matrix and the half spinor plus (Hp′) to be transferred to the negative neighbor of the gauge matrix.