// $Id: t_mesplq.cc,v 1.3 2003/02/20 19:37:32 dgr Exp $ #include #include #include "qdp.h" #include "proto.h" // This is needed for readArchiv #include "tutorial.h" using namespace QDP; int main(int argc, char *argv[]) { // Put the machine into a known state QDP_initialize(&argc, &argv); // Setup the layout const int foo[] = {4,4,4,8}; multi1d nrow(Nd); nrow = foo; // Use only Nd elements Layout::setLattSize(nrow); Layout::create(); //! Example of calling a plaquette routine multi1d u(Nd); Double w_plaq, s_plaq, t_plaq, link; readArchiv(u, argv[1]); // Now the business cerr << "Start mesplq\n"; MesPlq(u, w_plaq, s_plaq, t_plaq, link); cerr << "w_plaq = " << w_plaq << endl; cerr << "link = " << link << endl; // Time to bolt QDP_finalize(); return 0; }