Cut routines
Here is an overview of the cutting routines supported by the code.
Gluon fields can be computed from the logarithm or approximation to the logarithm of the SU(NC)-valued parallel transport matrices. By default the code uses the normal lattice approximation but can also use the logarithmic definition although for NC > 3, this becomes slow and unstable as series expansions are used to compute this thing.
Gluon fields are not uniquely defined unless the gauge is fixed, once it is though we can define our gauge fields from the parallel transport matrices,
The momentum-space gluon fields are matrices as are their Fourier transforms,
