View source for 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"

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);
    
    int nx,ny,nz,nt;
    nx=ny=nz=nt=4;
    std::string latname="auto";
       
    latname="c44_b5.3_k0.158_100";
    int bc[4]={1,1,1,-1};
    
    layout_minsurface_eo desc(nx, ny, nz, nt);                                                        
    su3_field links(&desc);                                                                           
                                                                                                    
    read_kentucky_lattice(latname, links);                                                            
                                                                                                    
    apply_boundary(links, bc);
    
    return 0;
  }