25 bool init = Smear::Factory::Register(
"APE_spatial", create_object);
36 const string str_vlevel = params.
get_string(
"verbose_level");
75 assert(U.
nvol() == Nvol);
78 assert(Usmear.
nvol() == Nvol);
83 Field_G c_tmp(Nvol, 1), u_tmp(Nvol, 1), u_tmp2(Nvol, 1);
87 double plaq = staple.
plaq_s(U);
94 for (
int mu = 0; mu < Ndim_spc; ++mu) {
96 u_tmp.setpart_ex(0, U, mu);
98 for (
int nu = 0; nu < Ndim_spc; ++nu) {
100 staple.
upper(u_tmp2, U, mu, nu);
101 c_tmp.addpart_ex(0, u_tmp2, 0,
m_rho);
103 staple.
lower(u_tmp2, U, mu, nu);
104 c_tmp.addpart_ex(0, u_tmp2, 0,
m_rho);
115 plaq = staple.
plaq_s(Usmear);
117 plaq = staple.
plaq_t(Usmear);
Projection * m_proj
projector to group element.
APE type smearing of spatial link variables.
void set(const int jin, const int site, const int jex, double v)
void general(const char *format,...)
int fetch_double(const string &key, double &value) const
double plaq_t(const Field_G &)
calculates temporal plaquette value.
void smear(Field_G &Usmear, const Field_G &U)
Smearing of a given gauge field.
void lower(Field_G &, const Field_G &, const int mu, const int nu)
constructs lower staple in mu-nu plane.
static const std::string class_name
void set_parameters(const Parameters ¶ms)
Setting parameters with Parameters object.
void staple(Field_G &, const Field_G &, const Field_G &, int mu, int nu)
Staple construction.
void upper(Field_G &, const Field_G &, const int mu, const int nu)
constructs upper staple in mu-nu plane.
void crucial(const char *format,...)
base class for projection operator into gauge group.
double m_rho
smearing parameter
double plaq_s(const Field_G &)
calculates spatial plaquette value.
base class for smearing of link variables.
virtual void project(Field_G &v, double alpha, const Field_G &C, const Field_G &U)=0
projection V = P[alpha, C, U]
void setpart_ex(int ex, const Field &w, int exw)
string get_string(const string &key) const
Bridge::VerboseLevel m_vl
int m_Ndim
spacetime dimension
static VerboseLevel set_verbose_level(const std::string &str)