Reading a Gauge Configuration

From Gw-qcd-wiki
Revision as of 15:02, 6 May 2019 by Chris (talk | contribs)
Jump to: navigation, search

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);

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;
   nx=ny=nz=nt=4;
   std::string latname="auto";
      
   latname="c44_b5.3_k0.158_100";
   int bc[4]={1,1,1,-1};
   int this_node = get_node_rank();

```