25 bool init = Force_G::Factory::Register(
"Force_G_Rectangle", create_object);
35 const string str_vlevel = params.
get_string(
"verbose_level");
40 double beta, c_plaq, c_rect;
86 assert(
m_U->
nin() == Nc * Nc * 2);
88 assert(
m_U->
nex() == Ndim);
90 assert(force.
nin() == Nc * Nc * 2);
91 assert(force.
nvol() == Nvol);
92 assert(force.
nex() == Ndim);
96 Field_G force1(Nvol, 1), force2(Nvol, 1);
98 Field_G Cup1(Nvol, 1), Cup2(Nvol, 1);
99 Field_G Cdn1(Nvol, 1), Cdn2(Nvol, 1);
100 Field_G Umu(Nvol, 1), Unu(Nvol, 1);
101 Field_G v(Nvol, 1), w(Nvol, 1), c(Nvol, 1);
103 for (
int mu = 0; mu < Ndim; ++mu) {
105 for (
int nu = 0; nu < Ndim; ++nu) {
106 if (nu == mu)
continue;
114 force1.addpart_ex(0, Cup1, 0,
m_c_plaq);
115 force1.addpart_ex(0, Cdn1, 0,
m_c_plaq);
121 Umu.setpart_ex(0, *
m_U, mu);
134 force1.addpart_ex(0, c, 0,
m_c_rect);
148 force1.addpart_ex(0, c, 0,
m_c_rect);
160 force1.addpart_ex(0, c, 0,
m_c_rect);
173 force1.addpart_ex(0, c, 0,
m_c_rect);
188 force1.addpart_ex(0, c, 0,
m_c_rect);
201 force1.addpart_ex(0, c, 0,
m_c_rect);
211 double Fave, Fmax, Fdev;
212 force.
stat(Fave, Fmax, Fdev);
static double epsilon_criterion()
void general(const char *format,...)
Container of Field-type object.
int fetch_double(const string &key, double &value) const
void at_Field_G(Field_G &W, const int ex)
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 lower(Field_G &, const Field_G &, const int mu, const int nu)
constructs lower staple in mu-nu plane.
static const std::string class_name
Base class of gauge force calculation.
Bridge::VerboseLevel m_vl
void backward(Field &, const Field &, const int mu)
void upper(Field_G &, const Field_G &, const int mu, const int nu)
constructs upper staple in mu-nu plane.
void set_parameters(const Parameters ¶ms)
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)
HMC force class for rectangular gauge action.
void crucial(const char *format,...)
void stat(double &Fave, double &Fmax, double &Fdev) const
determines the statistics of the field. average, maximum value, and deviation is determined over glob...
void setpart_ex(int ex, const Field &w, int exw)
string get_string(const string &key) const
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)