16 #ifdef USE_PARAMETERS_FACTORY
30 bool init = Smear::Factory::Register(
"HYP_SF", create_object);
48 #ifdef USE_PARAMETERS_FACTORY
63 const string str_vlevel = params.
get_string(
"verbose_level");
68 double alpha1, alpha2, alpha3;
69 valarray<double> phi, phipr;
90 double *phi,
double *phipr)
140 for (
int mu = 0; mu <
m_Ndim; ++mu) {
141 m_U[mu].setpart_ex(0, U, mu);
159 for (
int mu = 0; mu <
m_Ndim; ++mu) {
160 for (
int nu = 0; nu <
m_Ndim; ++nu) {
161 if (nu == mu)
continue;
163 for (
int rho = nu + 1; rho <
m_Ndim; ++rho) {
164 if (rho == mu)
continue;
165 int sig = 6 - mu - nu - rho;
182 for (
int mu = 0; mu <
m_Ndim; ++mu) {
183 for (
int nu = 0; nu <
m_Ndim; ++nu) {
184 if (nu == mu)
continue;
187 for (
int rho = 0; rho <
m_Ndim; ++rho) {
188 if ((rho != mu) && (rho != nu)) {
191 c2.addpart_ex(0, u_tmp, 0);
209 for (
int mu = 0; mu <
m_Ndim; ++mu) {
212 for (
int nu = 0; nu <
m_Ndim; ++nu) {
216 c3.addpart_ex(0, u_tmp, 0);
void staple(Field_G &, const Field_G &, const Field_G &, int mu, int nu)
void set_boundary_wkpr(const Mat_SU_N &U)
Set the boundary spatial link at t=Nt-1 for SF bc.
int index_v2(int mu, int nu)
void Register_string(const string &, const string &)
void mult_Field_Gdn(Field_G &w, const int ex, const Field_G &u1, const int ex1, const Field_G &u2, const int ex2)
int index_v1(int mu, int nu, int rho)
void general(const char *format,...)
void smear(Field_G &Usmear, const Field_G &U)
void addpart_ex(int ex, const Field &w, int exw)
void set_parameters(const Parameters ¶ms)
int fetch_double_vector(const string &key, std::valarray< double > &val) const
void mult_Field_Gnd(Field_G &w, const int ex, const Field_G &u1, const int ex1, const Field_G &u2, const int ex2)
std::valarray< Field_G > m_v1
Parameters_Smear_HYP_SF()
void set_boundary_wk(const Mat_SU_N &U)
Set the boundary spatial link at t=0 for SF bc.
void mult_Field_Gnn(Field_G &w, const int ex, const Field_G &u1, const int ex1, const Field_G &u2, const int ex2)
void backward(Field &, const Field &, const int mu)
void crucial(const char *format,...)
base class for projection operator into gauge group.
void set_parameters(const Parameters_Field_G_SF ¶ms)
Set the parameter with Parameters_Field_G_SF class.
static const std::string class_name
static bool Register(const std::string &realm, const creator_callback &cb)
std::valarray< Field_G > m_v2
void Register_double_vector(const string &, const std::valarray< double > &)
std::valarray< Field_G > m_U
void Register_double(const string &, const double)
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)
int fetch_double(const string &key, double &val) const
string get_string(const string &key) const
Bridge::VerboseLevel m_vl
static VerboseLevel set_verbose_level(const std::string &str)
void forward(Field &, const Field &, const int mu)