9 int main(
int argc,
char *argv[])
15 const int foo[] = {4,4,4,8};
16 multi1d<int> nrow(
Nd);
18 Layout::setLattSize(nrow);
21 XMLFileWriter xml(
"t_hypsmear.xml");
22 push(xml,
"t_hypsmear");
27 write(xml,
"nrow", nrow);
32 multi1d<LatticeColorMatrix>
u(
Nd);
38 QDPIO::cout <<
"Start gaussian\n";
39 for(
int m=0;
m <
u.size(); ++
m)
43 for(
int m=0;
m <
u.size(); ++
m)
48 QDPIO::cout <<
"w_plaq = " <<
w_plaq << std::endl;
49 QDPIO::cout <<
"link = " <<
link << std::endl;
52 push(xml,
"observables");
58 multi1d<LatticeColorMatrix> u_hyp(
Nd);
59 Real BlkAccu = 1.0e-5;
64 Hyp_Smear(
u, u_hyp, alpha1, alpha2, alpha3, BlkAccu, BlkMax);
68 QDPIO::cout <<
"w_plaq = " <<
w_plaq << std::endl;
69 QDPIO::cout <<
"link = " <<
link << std::endl;
72 push(xml,
"HYP_observables");
Primary include file for CHROMA in application codes.
void write(XMLWriter &xml, const std::string &path, const AsqtadFermActParams ¶m)
Writer parameters.
void Hyp_Smear(const multi1d< LatticeColorMatrix > &u, multi1d< LatticeColorMatrix > &u_hyp, const Real &alpha1, const Real &alpha2, const Real &alpha3, const Real &BlkAccu, int BlkMax)
Construct the "hyp-smeared" links of Anna Hasenfratz.
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[])