19 namespace InlineQQbarEnv
52 XMLReader paramtop(xml, path);
55 read(paramtop,
"version", version);
64 QDPIO::cerr <<
"Input parameter version " << version <<
" unsupported." << std::endl;
77 write(xml,
"version", version);
86 XMLReader inputtop(xml, path);
113 XMLReader paramtop(xml_in, path);
115 if (paramtop.count(
"Frequency") == 1)
128 QDPIO::cerr << __func__ <<
": Caught Exception reading XML: " << e << std::endl;
160 XMLBufferWriter gauge_xml;
166 catch( std::bad_cast )
178 const multi1d<LatticeColorMatrix>&
u =
181 push(xml_out,
"qqbar");
182 write(xml_out,
"update_no", update_no);
195 write(xml_out,
"Config_info", gauge_xml);
197 push(xml_out,
"Output_version");
198 write(xml_out,
"out_version", 1);
202 MesPlq(xml_out,
"Observables",
u);
205 const int Nprops = 2;
208 QDPIO::cerr <<
InlineQQbarEnv::name <<
": Error on input params - expecting 2 filenames" << std::endl;
218 multi1d<LatticePropagator> quark_propagator(Nprops);
222 for(
int i=0;
i < Nprops; ++
i)
227 quark_propagator[
i] =
231 XMLReader prop_file_xml, prop_record_xml;
243 catch( std::bad_cast )
259 write(xml_out,
"Propagator_input", qqbar);
269 for(
int i=0;
i < Nprops; ++
i)
271 multi1d<Double> prop_corr = sumMulti(localNorm2(quark_propagator[
i]),
274 push(xml_out,
"SinkSmearedProp_correlator");
275 write(xml_out,
"correlator_num",
i);
276 write(xml_out,
"sink_smeared_prop_corr", prop_corr);
283 multiNd<Complex> mesprop;
293 LatticePropagator q_tmp;
294 for(
int i=0;
i < Nprops; ++
i)
296 q_tmp = adj(
U) * quark_propagator[
i] *
U;
297 quark_propagator[
i] = q_tmp;
308 multi1d<Complex> mesprop_1d;
314 XMLBufferWriter file_xml;
315 push(file_xml,
"qqbar");
319 XMLBufferWriter record_xml;
320 push(record_xml,
"QQbar");
321 write(record_xml,
".", qqbar);
322 write(record_xml,
"Config_info", gauge_xml);
327 QDPIO_SINGLEFILE, QDPIO_SERIAL, QDPIO_OPEN);
328 write(to,record_xml,mesprop_1d);
337 << snoop.getTimeInSeconds()
338 <<
" secs" << std::endl;
Inline measurement factory.
Inline measurement of quark-antiquark 2-pt correlators.
void operator()(const unsigned long update_no, XMLWriter &xml_out)
Do the measurement.
Fourier transform phase factor support.
const Set & getSet() const
The set to be used in sumMulti.
Basis rotation matrix from Dirac to Degrand-Rossi (and reverse)
SpinMatrixD DiracToDRMat()
The Dirac to Degrand-Rossi spin transformation matrix (and reverse)
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 mescomp(multiNd< Complex > &mesprop, const LatticePropagator &quark_propagator_1, const LatticePropagator &quark_propagator_2, const SftMom &phases, int t0)
Construct all components of a meson propagator.
std::string uniqueId()
Return a unique id.
void proginfo(XMLWriter &xml)
Print out basic information about this program.
Inline construction of qqbar.
Construct all components of a meson propagator.
Named object function std::map.
static bool registered
Local registration flag.
bool registerAll()
Register all the factories.
Asqtad Staggered-Dirac operator.
static multi1d< LatticeColorMatrix > u
push(xml_out,"Condensates")
void convertMescomp(multi1d< Complex > &mesprop_1d, const multiNd< Complex > &mesprop, const int j_decay)
Convert generalized correlator object.
void MesPlq(const multi1d< LatticeColorMatrixF3 > &u, multi2d< Double > &plane_plaq, Double &link)
Print out basic info about this program.
Fourier transform phase factor support.
multi1d< std::string > prop_ids
void write(XMLWriter &xml_out, const std::string &path)
struct Chroma::InlineQQbarParams::Param_t param
struct Chroma::InlineQQbarParams::NamedObject_t named_obj
Mega structure holding QQbar props.
multi1d< ForwardProp_t > forward_props
multi1d< LatticeColorMatrix > U