Heatbath

We have a NC and ND generic heatbath-overrelxation routine that works for even and odd lattice lengths when using the normal Wilson plaquette gauge action.

Like everything in the library the NC and ND have to be compiled in and the heatbath and overrelaxation routines scale poorly with NC (actually like NC^2).

We also have the opportunity for so-called improved gauge actions such as the Iwasaki, Symanzik, DBW2 or generic (C_1) actions, but these will only work with even dimensional lattices that are divisible by 4. Be careful about the beta to get the regime you want! These have to be set at configure time again with the

     --with-IMPROVED_STAPLE={IWASAKI,SYMANZIK,DBW2}

Be careful with this if you wish to measure Wilson-flow or Topology or do anything as it will use these rectangles too! I suggest to have separate compiles with these things on or off.