26 bool init = Smear::Factory::Register(
"APE_SF", create_object);
37 const string str_vlevel = params.
get_string(
"verbose_level");
43 std::vector<double> phi, phipr;
61 double *phi,
double *phipr)
80 for (
int mu = 0; mu <
m_Ndim; ++mu) {
81 for (
int nu = 0; nu <
m_Ndim; ++nu) {
86 for (
int i = 0; i < 3; ++i) {
98 for (
int mu = 0; mu <
m_Ndim; ++mu) {
112 assert(rho.size() == m_Ndim *
m_Ndim);
116 for (
int mu = 0; mu <
m_Ndim; ++mu) {
117 for (
int nu = 0; nu <
m_Ndim; ++nu) {
122 for (
int i = 0; i < 3; ++i) {
134 assert(U.
nvol() == Nvol);
136 assert(Usmear.
nvol() == Nvol);
139 Field_G c_tmp(Nvol, 1), u_tmp(Nvol, 1);
149 for (
int mu = 0; mu <
m_Ndim; ++mu) {
151 u_tmp.setpart_ex(0, U, mu);
153 for (
int nu = 0; nu <
m_Ndim; ++nu) {
155 double rho =
m_rho[mu + m_Ndim * nu];
156 staple.
upper(u_tmp2, U, mu, nu);
157 c_tmp.addpart_ex(0, u_tmp2, 0, rho);
159 staple.
lower(u_tmp2, U, mu, nu);
160 c_tmp.addpart_ex(0, u_tmp2, 0, rho);
164 double rho0 =
m_rho[mu + m_Ndim * mu];
SU(N) gauge field class in which a few functions are added for the SF.
int fetch_double_vector(const string &key, vector< double > &value) const
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
APE type smearing of link variables.
double m_phi[3]
SF boundary condition at t=0.
static const std::string class_name
void set_parameters(const Parameters ¶ms)
void upper(Field_G_SF &, const Field_G &, const int, const int)
void set_boundary_wk(const Mat_SU_N &U)
Set the boundary spatial link at t=0 for SF bc.
void crucial(const char *format,...)
base class for projection operator into gauge group.
void set_parameters(const Parameters ¶ms)
double m_phipr[3]
SF boundary condition at t=Nt.
std::valarray< double > m_rho
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 lower(Field_G_SF &, const Field_G &, const int, const int)
void setpart_ex(int ex, const Field &w, int exw)
string get_string(const string &key) const
Bridge::VerboseLevel m_vl
void smear(Field_G &Usmear, const Field_G &U)
static VerboseLevel set_verbose_level(const std::string &str)