26 bool init = ForceSmear::Factory::Register(
"APE", create_object);
37 const string str_vlevel = params.
get_string(
"verbose_level");
69 for (
int mu = 0; mu <
m_Ndim; ++mu) {
70 for (
int nu = 0; nu <
m_Ndim; ++nu) {
82 for (
int mu = 0; mu <
m_Ndim; ++mu) {
88 assert(rho.size() == m_Ndim *
m_Ndim);
92 for (
int mu = 0; mu <
m_Ndim; ++mu) {
93 for (
int nu = 0; nu <
m_Ndim; ++nu) {
110 for (
int mu = 0; mu <
m_Ndim; ++mu) {
111 for (
int nu = 0; nu <
m_Ndim; ++nu) {
122 int NinG = 2 * Nc * Nc;
124 assert(Sigmap.
nin() == NinG);
131 for (
int mu = 0; mu <
m_Ndim; ++mu) {
132 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;
141 double rho =
m_rho[mu + m_Ndim * nu];
143 C.addpart_ex(0, c_tmp, 0, rho);
146 sigmap_tmp.setpart_ex(0, Sigmap, mu);
148 double alpha =
m_rho[mu + m_Ndim * mu];
150 alpha, sigmap_tmp, C,
m_U[mu]);
155 for (
int mu = 0; mu <
m_Ndim; ++mu) {
156 for (
int nu = 0; nu <
m_Ndim; ++nu) {
157 if (nu == mu)
continue;
158 double rho =
m_rho[mu + m_Ndim * nu];
186 axpy(Sigma_mu, 1.0, vt3);
191 axpy(Sigma_mu, 1.0, vt3);
201 axpy(Sigma_mu, 1.0, vt2);
Bridge::VerboseLevel m_vl
void set(const int jin, const int site, const int jex, double v)
void force_each(Field_G &, const Field_G &, const Field_G &, const Field_G &, const Field_G &, int mu, int nu)
Recursive calculation for APE smeared fermion force.
void general(const char *format,...)
std::vector< Field_G > m_iTheta
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)
void addpart_ex(int ex, const Field &w, int exw)
Base class for force calculation of smeared operators.
void staple(Field_G &, const Field_G &, const Field_G &, int mu, int nu)
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)
void set_parameters(const Parameters ¶ms)
void backward(Field &, const Field &, const int mu)
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,...)
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_U
double rho(int mu, int nu)
void setpart_ex(int ex, const Field &w, int exw)
string get_string(const string &key) const
std::vector< double > m_rho
static const std::string class_name
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)