37 multi1d<int> boundary;
39 bool do_gauge_transform ;
61 XMLReader inputtop(xml, path);
75 QDPIO::cerr <<
"Error reading IO version: " << e << std::endl;
79 XMLReader paramtop(inputtop,
"param");
92 QDPIO::cerr <<
"Error parameters: " << e << std::endl;
105 read(inputtop,
"Cfg", input.
cfg);
109 QDPIO::cerr <<
"Error gauge name data: " << e << std::endl;
124 int main(
int argc,
char **argv)
129 QDPIO::cout <<
"Measure the Wilson flow " << std::endl;
130 QDPIO::cout <<
"Calculation for SU(" << Nc <<
")" << std::endl;
133 TextFileReader text_in ;
137 text_in.open(queuefile_name);
141 QDPIO::cerr <<
"Error reading input file " << queuefile_name << std::endl;
142 QDPIO::cerr <<
"The input file name can be passed via the -i flag " << std::endl;
143 QDPIO::cerr <<
"The default name is ./DATA" << std::endl;
150 bool successful_read=
true;
152 while(successful_read){
156 successful_read=!text_in.fail( );
157 if (successful_read){
162 QDPIO::cout <<
"Starting run with in_name:\t"<< in_name << std::endl;
163 QDPIO::cout <<
"Output name = out_name:\t"<< out_name << std::endl ;
173 xml_in.open(in_name);
177 QDPIO::cerr <<
"Error reading input file " << in_name << std::endl;
178 QDPIO::cerr <<
"The input file name can be passed via the -i flag " << std::endl;
184 read(xml_in,
"/WilsonFlow", input);
192 multi1d<LatticeColorMatrix>
u(
Nd);
194 XMLReader gauge_file_xml, gauge_xml;
203 XMLFileWriter xml_out(out_name);
204 push(xml_out,
"wilsonFlow_output");
207 write(xml_out,
"Input", xml_in);
210 push(xml_out,
"Output_version");
211 write(xml_out,
"out_version", 1);
218 MesPlq(xml_out,
"Observables",
u);
224 QDPIO::cout <<
"FLOW for " << input.
cfg.
cfg_file <<
"\n" ;
238 QDPIO::cout <<
"End of measurements " << std::endl;
Primary include file for CHROMA in application codes.
void read(XMLReader &xml, const std::string &path, AsqtadFermActParams ¶m)
Read parameters.
void write(XMLWriter &xml, const std::string &path, const AsqtadFermActParams ¶m)
Writer parameters.
void unitarityCheck(const multi1d< LatticeColorMatrixF3 > &u)
Check the unitarity of color matrix in SU(N)
void gaugeStartup(XMLReader &gauge_file_xml, XMLReader &gauge_xml, multi1d< LatticeColorMatrix > &u, Cfg_t &cfg)
Initialize the gauge fields.
void wilson_flow(XMLWriter &xml, multi1d< LatticeColorMatrix > &u, int nstep, Real wflow_eps, int t_dir, const multi1d< bool > &smear_dirs)
Compute the Wilson flow.
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.
std::string getXMLInputFileName()
Get input file name.
void MesPlq(const multi1d< LatticeColorMatrixF3 > &u, multi2d< Double > &plane_plaq, Double &link)
void setXMLInputFileName(const std::string &name)
Set input file name.
Gauge configuration structure.
Parameters for running program.
int main(int argc, char **argv)