10 int main(
int argc,
char *argv[])
17 const int foo[] = {16,16,16,32};
20 multi1d<int> nrow(
Nd);
22 Layout::setLattSize(nrow);
25 XMLFileWriter xml(
"t_fuzwilp.xml");
26 push(xml,
"t_fuzwilp");
31 write(xml,
"nrow",nrow);
36 multi1d<LatticeColorMatrix>
u(
Nd);
49 gaugeFilename=
"/raidd/aci/configs/getlat.1200";
53 QDPIO::cout <<
"Reading config: " << gaugeFilename << std::endl;
54 readArchiv(gauge_xml,
u, gaugeFilename);
57 for(
int m=0;
m <
u.size(); ++
m)
65 QDPIO::cout <<
"Measure the plaquette and link " << std::endl;
73 double time_in_sec = swatch.getTimeInSeconds();
75 QDPIO::cout <<
"w_plaq = " <<
w_plaq << std::endl;
76 QDPIO::cout <<
"s_plaq = " <<
s_plaq << std::endl;
77 QDPIO::cout <<
"t_plaq = " <<
t_plaq << std::endl;
78 QDPIO::cout <<
"link = " <<
link << std::endl;
80 QDPIO::cout <<
"Measurement took " << time_in_sec
81 <<
" secs" << std::endl;
84 push(xml,
"observables");
91 QDPIO::cout <<
"Measure ape-smeared Wilson loops" << std::endl;
98 Real BlkAccu = 1.0e-5;
101 QDPIO::cout <<
"sm_fact: " << sm_fact << std::endl;
102 QDPIO::cout <<
"n_smear: " << n_smear << std::endl;
108 "Fuzzed_Wilson_Loops");
111 time_in_sec = swatch.getTimeInSeconds();
112 QDPIO::cout <<
"Loop measurements complete" << std::endl;
113 QDPIO::cout <<
"Measurement took " << time_in_sec
114 <<
" secs" << std::endl;
Primary include file for CHROMA in application codes.
void write(XMLWriter &xml, const std::string &path, const AsqtadFermActParams ¶m)
Writer parameters.
void fuzwilp(const multi1d< LatticeColorMatrix > &u, int j_decay, int tmax, int n_smear, const Real &sm_fact, const Real &BlkAccu, int BlkMax, XMLWriter &xml, const std::string &xml_group)
Calculate ape-fuzzed Wilson loops.
Asqtad Staggered-Dirac operator.
static multi1d< LatticeColorMatrix > u
push(xml_out,"Condensates")
void initialize(int *argc, char ***argv)
Chroma initialisation routine.
void finalize(void)
Chroma finalization routine.
void reunit(LatticeColorMatrixF3 &xa)
void MesPlq(const multi1d< LatticeColorMatrixF3 > &u, multi2d< Double > &plane_plaq, Double &link)
FloatingPoint< double > Double
int main(int argc, char *argv[])