25 bool init = Action::Factory::Register(
"Action_G_Rectangle_SF", create_object);
36 const string str_vlevel = params.
get_string(
"verbose_level");
41 double beta, c_plaq, c_rect;
42 double ct0, ct1, ct2, ctr0, ctr1, ctr2;
43 std::vector<double> phi, phipr;
67 double gg = 6.0 / beta;
69 double ct = ct0 + ct1 * gg + ct2 * gg * gg;
70 double ctr = ctr0 + ctr1 * gg + ctr2 * gg * gg;
95 double *phi,
double *phipr,
double ct,
double ctr)
128 double c0r, c0i, c1r, c1i, c2r, c2i;
129 c0r = cos(phi[0] / Lx);
130 c0i = sin(phi[0] / Lx);
131 c1r = cos(phi[1] / Lx);
132 c1i = sin(phi[1] / Lx);
133 c2r = cos(phi[2] / Lx);
134 c2i = sin(phi[2] / Lx);
140 c0r = cos(phipr[0] / Lx);
141 c0i = sin(phipr[0] / Lx);
142 c1r = cos(phipr[1] / Lx);
143 c1i = sin(phipr[1] / Lx);
144 c2r = cos(phipr[2] / Lx);
145 c2i = sin(phipr[2] / Lx);
156 double H_U =
calcH();
232 for (
int mu = 0; mu < Ndim; ++mu) {
233 for (
int nu = mu + 1; nu < Ndim; ++nu) {
246 for (
int site = 0; site < Nvol; ++site) {
275 for (
int site = 0; site < Nvol; ++site) {
290 for (
int site = 0; site < Nvol; ++site) {
SU(N) gauge field class in which a few functions are added for the SF.
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,...)
Container of Field-type object.
int fetch_double(const string &key, double &value) const
static int ipe(const int dir)
logical coordinate of current proc.
virtual void force_core(Field &)=0
Base class of HMC action class family.
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 set_parameters(const Parameters ¶ms)
void upper(Field_G_SF &, const Field_G &, const int, const int)
void set_boundary_wk(const Mat_SU_N &U)
Set the boundary spatial link at t=0 for SF bc.
double m_ct
SF boundary improvement coefficient for the plaquatte action.
void backward(Field &, const Field &, const int mu)
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,...)
HMC action class for rectangular gauge action with the SF BC.
Bridge::VerboseLevel m_vl
Base class of random number generators.
double m_ctr
SF boundary improvement coefficient for the rectangle action.
static const std::string class_name
Mat_SU_N mat_dag(const int site, const int mn=0) const
double langevin(RandomNumbers *)
Langevis step.
static int reduce_sum(int count, double *recv_buf, double *send_buf, int pattern=0)
make a global sum of an array of double over the communicator. pattern specifies the dimensions to be...
Mat_SU_N m_wk
SF boundary condition.
virtual void set_parameters(const Parameters &)=0
void set(int c, double re, const double &im)
void setpart_ex(int ex, const Field &w, int exw)
string get_string(const string &key) const
Mat_SU_N mat(const int site, const int mn=0) const
static VerboseLevel set_verbose_level(const std::string &str)
double ReTr(const Mat_SU_N &m)
void mult_Field_Gnd(Field_G &W, const int ex, const Field_G &U1, const int ex1, const Field_G &U2, const int ex2)