41 const string str_vlevel = params.
get_string(
"verbose_level");
46 std::vector<double> phi, phipr, pomega;
65 std::vector<double>& pomega)
81 assert(phi.size() == 3);
82 assert(phipr.size() == 3);
83 assert(pomega.size() == 3);
93 double c0r, c0i, c1r, c1i, c2r, c2i;
95 c0r = cos(phi[0] /
Lx);
96 c0i = sin(phi[0] /
Lx);
97 c1r = cos(phi[1] /
Lx);
98 c1i = sin(phi[1] /
Lx);
99 c2r = cos(phi[2] /
Lx);
100 c2i = sin(phi[2] /
Lx);
103 wk.
set(0, 0, c0r, c0i);
104 wk.
set(1, 1, c1r, c1i);
105 wk.
set(2, 2, c2r, c2i);
107 c0r = cos(phipr[0] /
Lx);
108 c0i = sin(phipr[0] /
Lx);
109 c1r = cos(phipr[1] /
Lx);
110 c1i = sin(phipr[1] /
Lx);
111 c2r = cos(phipr[2] /
Lx);
112 c2i = sin(phipr[2] /
Lx);
147 double c0r, c0i, c1r, c1i, c2r, c2i;
149 c0r = cos(phi[0] /
Lx);
150 c0i = sin(phi[0] /
Lx);
151 c1r = cos(phi[1] /
Lx);
152 c1i = sin(phi[1] /
Lx);
153 c2r = cos(phi[2] /
Lx);
154 c2i = sin(phi[2] /
Lx);
156 wk.
set(0, 0, c0r, c0i);
157 wk.
set(1, 1, c1r, c1i);
158 wk.
set(2, 2, c2r, c2i);
160 c0r = cos(phipr[0] /
Lx);
161 c0i = sin(phipr[0] /
Lx);
162 c1r = cos(phipr[1] /
Lx);
163 c1i = sin(phipr[1] /
Lx);
164 c2r = cos(phipr[2] /
Lx);
165 c2i = sin(phipr[2] /
Lx);
240 for (
int nu = 0; nu <
Ndim - 1; nu++) {
241 upper(staple, U, 3, nu);
243 for (
int z = 0; z <
Nz; z++) {
244 for (
int y = 0; y <
Ny; y++) {
245 for (
int x = 0; x <
Nx; x++) {
293 vout.
general(
m_vl,
"SF_delSg_plaq, from 0, from T = %.8f %.8f %.8f\n", plaq, plaqt0, plaqtT);
424 for (
int mu = 0; mu <
Ndim - 1; mu++) {
448 upper(Cup2, m_U, nu, mu);
454 for (
int z = 0; z <
Nz; z++) {
455 for (
int y = 0; y <
Ny; y++) {
456 for (
int x = 0; x <
Nx; x++) {
463 wmat = cmat * Unu.
mat_dag(site);
467 wmat = cmat * Unu.
mat_dag(site);
468 rect02 +=
ReTr(
wk * wmat);
475 cmat = Unu.
mat(site) * wmat;
479 cmat = Unu.
mat(site) * wmat;
493 upper(Cup1, m_U, mu, nu);
497 for (
int z = 0; z <
Nz; z++) {
498 for (
int y = 0; y <
Ny; y++) {
499 for (
int x = 0; x <
Nx; x++) {
505 cmat = Unu.
mat(site) * wmat;
506 wmat =
wk * cmat.
dag();
520 lower(Cdn1, m_U, mu, nu);
523 for (
int z = 0; z <
Nz; z++) {
524 for (
int y = 0; y <
Ny; y++) {
525 for (
int x = 0; x <
Nx; x++) {
532 wmat = cmat * Cdn1.mat_dag(site);
533 rectt3 +=
ReTr(Unu.
mat(site) * wmat);
552 rect = -rect01 - rect02 - rect03 + rectt1 + rectt2 + rectt3;
554 vout.
general(
m_vl,
"SF_delSg_rect, at 01, 02, 03, at T1, T2, T3 = %.8f %.8f %.8f %.8f %.8f %.8f %.8f\n",
555 rect, rect01, rect02, rect03, rectt1, rectt2, rectt3);
621 upper(staple, U, 0, 1);
622 for (
int site = 0; site <
Nvol; site++) {
626 upper(staple, U, 1, 2);
627 for (
int site = 0; site <
Nvol; site++) {
631 upper(staple, U, 2, 0);
632 for (
int site = 0; site <
Nvol; site++) {
662 for (
int nu = 0; nu <
Ndim - 1; nu++) {
663 lower(staple, U, 3, nu);
665 for (
int site = 0; site <
Nvol; site++) {
702 for (
int nu = 0; nu <
Ndim - 1; nu++) {
703 lower(staple, U, 3, nu);
711 for (
int site = 0; site <
Nvol; site++) {
751 for (
int nu = 0; nu <
Ndim; nu++) {
753 upper(c_tmp, U, mu, nu);
755 lower(c_tmp, U, mu, nu);
788 for (
int nu = 0; nu <
Ndim; nu++) {
790 upper(staple_upper, U, mu, nu);
791 lower(staple_lower, U, mu, nu);
810 axpy(W, 1.0, staple_upper);
811 axpy(W, 1.0, staple_lower);
895 double plaq2 = plaq + 3 * 3 *
Lx *
Ly *
Lz;
898 plaq / (3 *
Lx *
Ly * Lz * (6 *
Lt - 3)));
900 plaq2 / (3 * 6 *
Lx *
Ly * Lz *
Lt));
SU(N) gauge field class in which a few functions are added for the SF.
static const std::string class_name
void set_boundary_zero()
Set the boundary matrix to 0 for SF bc.
double plaquette_ct(const Field_G &, double ct)
void set_boundary_wkpr(const Mat_SU_N &U)
Set the boundary spatial link at t=Nt-1 for SF bc.
void print_plaquette(const Field_G &)
int fetch_double_vector(const string &key, vector< double > &value) const
void set(const int jin, const int site, const int jex, double v)
int site(const int &x, const int &y, const int &z, const int &t) const
void general(const char *format,...)
Bridge::VerboseLevel m_vl
double sf_coupling_plaq(const Field_G &, double ct)
void staple_ct(Field_G &, const Field_G &, const int, double ct)
void mult_Field_Gdn(Field_G &W, const int ex, const Field_G &U1, const int ex1, const Field_G &U2, const int ex2)
static int ipe(const int dir)
logical coordinate of current proc.
double plaq_s(const Field_G &)
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)
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)
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,...)
double plaq_t(const Field_G &)
double sf_coupling_rect(const Field_G &, double ctr)
void staple(Field_G &, const Field_G &, const int)
Mat_SU_N mat_dag(const int site, const int mn=0) const
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...
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
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 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)
double plaq_t_ct(const Field_G &, double ct)
double plaquette(const Field_G &)