26 bool init = ForceSmear::Factory::Register(
"HYP", create_object);
37 const string str_vlevel = params.
get_string(
"verbose_level");
42 double alpha1, alpha2, alpha3;
106 for (
int mu = 0; mu <
m_Ndim; ++mu) {
107 m_U[mu].setpart_ex(0, U, mu);
136 for (
int mu = 0; mu <
m_Ndim; ++mu) {
139 for (
int nu = 0; nu <
m_Ndim; ++nu) {
140 if (nu == mu)
continue;
142 C.addpart_ex(0, c_tmp, 0);
146 Sigmap_tmp.setpart_ex(0, Sigmap, mu);
152 for (
int mu = 0; mu <
m_Ndim; ++mu) {
153 for (
int nu = 0; nu <
m_Ndim; ++nu) {
154 if (nu == mu)
continue;
171 for (
int mu = 0; mu <
m_Ndim; ++mu) {
172 for (
int nu = 0; nu <
m_Ndim; ++nu) {
173 if (nu == mu)
continue;
177 for (
int rho = 0; rho <
m_Ndim; ++rho) {
178 if ((rho == mu) || (rho == nu))
continue;
181 C.addpart_ex(0, c_tmp, 0);
191 for (
int mu = 0; mu <
m_Ndim; ++mu) {
192 for (
int nu = 0; nu <
m_Ndim; ++nu) {
193 if (nu == mu)
continue;
194 for (
int rho = 0; rho <
m_Ndim; ++rho) {
195 if ((rho == mu) || (rho == nu))
continue;
211 for (
int mu = 0; mu <
m_Ndim; ++mu) {
212 for (
int nu = 0; nu <
m_Ndim; ++nu) {
213 if (nu == mu)
continue;
214 for (
int rho = 0; rho <
m_Ndim; ++rho) {
215 if ((rho == mu) || (rho == nu))
continue;
217 int sig = 6 - mu - nu - rho;
229 for (
int mu = 0; mu <
m_Ndim; ++mu) {
230 for (
int nu = 0; nu <
m_Ndim; ++nu) {
231 if (nu == mu)
continue;
232 for (
int rho = 0; rho <
m_Ndim; ++rho) {
233 if ((rho == mu) || (rho == nu))
continue;
234 int sig = 6 - mu - nu - rho;
265 axpy(Sigma_mu, 1.0, vt3);
270 axpy(Sigma_mu, 1.0, vt3);
280 axpy(Sigma_mu, 1.0, vt2);
294 for (
int mu = 0; mu <
m_Ndim; ++mu) {
295 for (
int nu = 0; nu <
m_Ndim; ++nu) {
296 if (nu == mu)
continue;
297 for (
int rho = nu + 1; rho <
m_Ndim; ++rho) {
298 if (rho == mu)
continue;
299 int sig = 6 - mu - nu - rho;
314 for (
int mu = 0; mu <
m_Ndim; ++mu) {
315 for (
int nu = 0; nu <
m_Ndim; ++nu) {
316 if (nu == mu)
continue;
320 for (
int rho = 0; rho <
m_Ndim; ++rho) {
321 if ((rho != mu) && (rho != nu)) {
324 c2.addpart_ex(0, u_tmp, 0);
void scal(Field &x, const double a)
scal(x, a): x = a * x
Bridge::VerboseLevel m_vl
void set_parameters(const Parameters ¶ms)
void set(const int jin, const int site, const int jex, double v)
void general(const char *format,...)
std::vector< Field_G > m_U
void force_step3(Field_G &, const Field_G &)
void multadd_Field_Gnd(Field_G &W, const int ex, const Field_G &U1, const int ex1, const Field_G &U2, const int ex2, const double ff)
int fetch_double(const string &key, double &value) const
void mult_Field_Gdn(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_iTheta2
void force_each(Field_G &, const Field_G &, const Field_G &, const Field_G &, const Field_G &, int mu, int nu)
void addpart_ex(int ex, const Field &w, int exw)
void staple(Field_G &, const Field_G &, const Field_G &, int mu, int nu)
std::vector< Field_G > m_v2
Base class for force calculation of smeared operators.
virtual void force_recursive(Field_G &Xi, Field_G &iTheta, double alpha, const Field_G &Sigmap, const Field_G &C, const Field_G &U)=0
determination of fields for force calculation
void mult_Field_Gdd(Field_G &W, const int ex, const Field_G &U1, const int ex1, const Field_G &U2, const int ex2)
static const std::string class_name
void backward(Field &, const Field &, const int mu)
Recursive calculation of HYP smeared fermion force.
int idx1(int mu, int nu, int rho)
void axpy(Field &y, const double a, const Field &x)
axpy(y, a, x): y := a * x + y
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 crucial(const char *format,...)
int idx1b(int mu, int nu, int rho)
void force_udiv(Field_G &Sigma, const Field_G &Sigma_p, const Field_G &U)
base class for projection operator into gauge group.
std::vector< Field_G > m_iTheta1
std::vector< Field_G > m_v1
virtual void project(Field_G &v, double alpha, const Field_G &C, const Field_G &U)=0
projection V = P[alpha, C, U]
std::vector< Field_G > m_Sigma2
void force_step1(Field_G &)
void force_step2(Field_G &)
std::vector< Field_G > m_iTheta3
string get_string(const string &key) const
std::vector< Field_G > m_Sigma3
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)