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.