25 bool init = Force_G::Factory::Register(
"Force_G_Rectangle_SF", create_object);
35 const string str_vlevel = params.
get_string(
"verbose_level");
40 double beta, c_plaq, c_rect;
41 double ct0, ct1, ct2, ctr0, ctr1, ctr2;
42 std::vector<double> phi, phipr;
66 double gg = 6.0 / beta;
68 double ct = ct0 + ct1 * gg + ct2 * gg * gg;
69 double ctr = ctr0 + ctr1 * gg + ctr2 * gg * gg;
91 double *phi,
double *phipr,
double ct,
double ctr)
124 double c0r, c0i, c1r, c1i, c2r, c2i;
125 c0r = cos(phi[0] / Lx);
126 c0i = sin(phi[0] / Lx);
127 c1r = cos(phi[1] / Lx);
128 c1i = sin(phi[1] / Lx);
129 c2r = cos(phi[2] / Lx);
130 c2i = sin(phi[2] / Lx);
136 c0r = cos(phipr[0] / Lx);
137 c0i = sin(phipr[0] / Lx);
138 c1r = cos(phipr[1] / Lx);
139 c1i = sin(phipr[1] / Lx);
140 c2r = cos(phipr[2] / Lx);
141 c2i = sin(phipr[2] / Lx);
215 assert(
m_U->
nex() == Ndim);
218 assert(force.
nvol() == Nvol);
219 assert(force.
nex() == Ndim);
236 for (
int mu = 0; mu < Ndim; ++mu) {
238 for (
int nu = 0; nu < Ndim; ++nu) {
239 if (nu == mu)
continue;
381 double Fave, Fmax, Fdev;
382 force.
stat(Fave, Fmax, Fdev);
SU(N) gauge field class in which a few functions are added for the SF.
void set_boundary_zero()
Set the boundary matrix to 0 for SF bc.
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)
static const std::string class_name
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 set_parameters(const Parameters ¶ms)
static int ipe(const int dir)
logical coordinate of current proc.
void addpart_ex(int ex, const Field &w, int exw)
double m_ct
SF boundary improvement coefficient for the plaquatte action.
void set_parameters(const Parameters ¶ms)
void mult_ct_boundary(int t, double ct)
Multiply the boundary improvement factor ct or ctr to an SU(N) matrix object which belongs to a site ...
void upper(Field_G_SF &, const Field_G &, const int, const int)
Base class of gauge force calculation.
Bridge::VerboseLevel m_vl
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)
Mat_SU_N m_wk
SF boundary condition.
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 stat(double &Fave, double &Fmax, double &Fdev) const
determines the statistics of the field. average, maximum value, and deviation is determined over glob...
HMC force class for rectangular gauge action with the SF BC.
double m_ctr
SF boundary improvement coefficient for the rectangle action.
void set(int c, double re, const double &im)
void lower(Field_G_SF &, const Field_G &, const int, const int)
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)