Fixed ordering of offsets when constructing vector of nearest neighbours.
Setting nearest neighbour construction to work when there is more than one neighbour.
Changed way in which nearest neighbour dimension is stored and used.
Fixed wrong index in construction of columnLink array in Dirac matrix construction.
Changed to cold-start initial configuration of gauge field to allow for free-field tests of propagator.
Added flexibility in spatial and time extents to todo list.
Reversed order of gamma matrix addition and subtration in dirac matrix comutation.
Added if statement so that only non-zero elements are added to Dirac matrix
Replaced 1/2 with 0.5 to make code succinct.
Added choice between CG and BiCGSTAB solvers (BiCGSTAB fails for free fermions)
Removed momentum projection as this must be done after the correlator has beenc computed.
Updated all lattice.cpp-dependent files to reflect changes from previous commits.
No point in writing gamma0 explicitly, as it's just gamma4.
Added a couple more things to the todo list.
Updated interfaces to offer choice over solver used in propagator.
Added propagator calculation and storage to run script
Updated readme to reflect changes in the source files