16 #ifdef USE_FACTORY_AUTOREGISTER
18 bool init = Smear_HYP_SF::register_factory();
27 const string str_vlevel = params.
get_string(
"verbose_level");
32 double alpha1, alpha2, alpha3;
33 std::vector<double> phi, phipr;
53 double *phi,
double *phipr)
103 for (
int mu = 0; mu <
m_Ndim; ++mu) {
104 m_U[mu].setpart_ex(0, U, mu);
120 for (
int mu = 0; mu <
m_Ndim; ++mu) {
121 for (
int nu = 0; nu <
m_Ndim; ++nu) {
122 if (nu == mu)
continue;
124 for (
int rho = nu + 1; rho <
m_Ndim; ++rho) {
125 if (rho == mu)
continue;
127 int sig = 6 - mu - nu - rho;
145 for (
int mu = 0; mu <
m_Ndim; ++mu) {
146 for (
int nu = 0; nu <
m_Ndim; ++nu) {
147 if (nu == mu)
continue;
152 for (
int rho = 0; rho <
m_Ndim; ++rho) {
153 if ((rho != mu) && (rho != nu)) {
174 for (
int mu = 0; mu <
m_Ndim; ++mu) {
180 for (
int nu = 0; nu <
m_Ndim; ++nu) {
202 const int mu,
const int nu)
void scal(Field &x, const double a)
scal(x, a): x = a * x
void set_boundary_wkpr(const Mat_SU_N &U)
Set the boundary spatial link at t=Nt-1 for SF bc.
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
void smear(Field_G &Usmear, const Field_G &U)
void mult_Field_Gdn(Field_G &W, const int ex, const Field_G &U1, const int ex1, const Field_G &U2, const int ex2)
void addpart_ex(int ex, const Field &w, int exw)
void set_parameters(const Parameters ¶ms)
virtual void project(Field_G &v, const double alpha, const Field_G &C, const Field_G &U)=0
projection V = P[alpha, C, U]
std::vector< Field_G > m_v1
void set_parameters(const std::vector< double > &phi, const std::vector< double > &phipr)
Set the parameter by giving vector objects.
int index_v2(const int mu, int nu)
void set_boundary_wk(const Mat_SU_N &U)
Set the boundary spatial link at t=0 for SF bc.
void backward(Field &, const Field &, const int mu)
int index_v1(const int mu, const int nu, const int rho)
void mult_Field_Gnn(Field_G &W, const int ex, const Field_G &U1, const int ex1, const Field_G &U2, const int ex2)
std::vector< Field_G > m_U
void crucial(const char *format,...)
static const std::string class_name
void staple(Field_G &, const Field_G &, const Field_G &, const int mu, const int nu)
void setpart_ex(int ex, const Field &w, int exw)
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)
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::vector< Field_G > m_v2