Reading a Gauge Configuration

First we need to include the qcd headers for the functions we will need.

    #include "layout_minsurface.h"
    #include "qcd_fields.h"
    #include "boundary_flip.h"

Then for output

    #include <iostream>

The qcd library has its own namespace so to avoid excessive typing of qcd::

    using namespace qcd;

Then we can get started with the code

    int main(int argc, char **argv)
    {
      init_machine(argc, argv);

This function initializes communications. There is more but I don't want to type all this up just to replace it. For now just testing to see if indentation keeps spacing

      int nx,ny,nz,nt;
      std::string latname="auto";

      nx=ny=nz=nt=4;
      latname="c44_b5.3_k0.158_100";
      int bc[4]={1,1,1,-1};

      int this_node = get_node_rank();