Bridge++  Ver. 1.3.x
Public Member Functions | Static Public Attributes | Protected Attributes | Private Member Functions | Private Attributes | List of all members
Staples_SF Class Reference

#include <staples_SF.h>

Public Member Functions

 Staples_SF ()
 
void set_parameters (const Parameters &params)
 
void set_parameters (double *phi, double *phipr)
 
void set_parameters (const double *phi, const double *phipr, const double *pomega)
 
void set_parameters (std::vector< double > &phi, std::vector< double > &phipr, std::vector< double > &pomega)
 
void set_parameter_verboselevel (const Bridge::VerboseLevel vl)
 
void upper (Field_G_SF &, const Field_G &, const int, const int)
 
void lower (Field_G_SF &, const Field_G &, const int, const int)
 
double plaq_s (const Field_G &)
 
double plaq_t (const Field_G &)
 
double plaq_t_ct (const Field_G &, double ct)
 
double plaquette (const Field_G &)
 
double plaquette_ct (const Field_G &, double ct)
 
double sf_coupling_plaq (const Field_G &, double ct)
 
double sf_coupling_rect (const Field_G &, double ctr)
 
void staple (Field_G &, const Field_G &, const int)
 
void staple_ct (Field_G &, const Field_G &, const int, double ct)
 
void print_plaquette (const Field_G &)
 

Static Public Attributes

static const std::string class_name = "Staples_SF"
 

Protected Attributes

Bridge::VerboseLevel m_vl
 

Private Member Functions

 Staples_SF (const Staples_SF &)
 
Staples_SF & operator= (const Staples_SF &)
 

Private Attributes

int Nc
 
int Ndim
 
int Nvol
 
int Lvol
 
Field_G_SF Umu
 
Field_G_SF Unu
 
Field_G_SF v
 
Field_G_SF w
 
Index_lex index
 
ShiftField_lex shift
 
int Nx
 
int Ny
 
int Nz
 
int Nt
 
int Lx
 
int Ly
 
int Lz
 
int Lt
 
int NPEt
 
Mat_SU_N wk
 
Mat_SU_N wkpr
 
Mat_SU_N iomega0
 
int initialized
 

Detailed Description

Definition at line 47 of file staples_SF.h.

Constructor & Destructor Documentation

Staples_SF::Staples_SF ( )
inline

Definition at line 71 of file staples_SF.h.

Staples_SF::Staples_SF ( const Staples_SF &  )
private

Member Function Documentation

void Staples_SF::lower ( Field_G_SF &  c,
const Field_G &  U,
const int  mu,
const int  nu 
)

Evaluate lower staple for all the links in mu direction with SF BC. 

      +     +       +
   nu |     |     w=|
     i+-->--+       +
    (1)  mu (2)    (2)
  • lower staple attached to the t=0 boundary spatial link is set to zero.
  • no link exist at t=T boundary in this code. (t=T is identified with t=0 with PBC).

Definition at line 903 of file staples_SF.cpp.

Staples_SF& Staples_SF::operator= ( const Staples_SF &  )
private
double Staples_SF::plaq_s ( const Field_G &  U)

Evaluate summed spatial ReTr plaquette with SF BC.

  • Plaquette is NOT normalized at all.
  • Contribution from the spatial plaquettes at t=0 and t=T are not included since they are always unity.
  • The function upper(U,mu=k,nu) is modified to give zero.

Definition at line 639 of file staples_SF.cpp.

double Staples_SF::plaq_t ( const Field_G &  U)

Evaluate summed temporal ReTr plaquette with SF BC.

  • Plaquette is NOT normalized at all.
  • The temporal plaquettes attached to t=0 and t=T boundary are evaluated with Wk and Wk'.
  • The function lower(U,3,nu) is modified to use Wk and Wk'.

Definition at line 675 of file staples_SF.cpp.

double Staples_SF::plaq_t_ct ( const Field_G &  U,
double  ct 
)

Evaluate summed temporal ReTr plaquette with SF BC.

  • The temporal plaquette attached to the boundary is multiplied with ct. 
         +---+
  ct |   |
 t=0 x---+
  • The temporal plaquettes attached to t=0 and t=T boundary are evaluated with Wk and Wk'.
  • The function lower(U,3,nu) is modified to use Wk and Wk'.
  • Plaquette is NOT normalized at all.

Definition at line 715 of file staples_SF.cpp.

double Staples_SF::plaquette ( const Field_G &  U)

Evaluate summed ReTr plaquette with SF BC.

  • Plaquette is NOT normalized at all.
  • Contribution from the spatial plaquettes at t=0 and t=T are not included since they are always unity.

Definition at line 593 of file staples_SF.cpp.

double Staples_SF::plaquette_ct ( const Field_G &  U,
double  ct 
)

Evaluate summed ReTr plaquette with SF BC.

  • The temporal plaquette attached to the boundary is multiplied with the improvement factor ct. 
         +---+
  ct |   |
 t=0 x---+
  • Plaquette is NOT normalized at all.
  • Contribution from the spatial plaquettes at t=0 and t=T are not included since they are always unity.

Definition at line 619 of file staples_SF.cpp.

void Staples_SF::print_plaquette ( const Field_G &  U)

Print out the plaquette for two definitions.

\[ \frac{1}{3}\frac{1}{6L_xL_yL_zL_t-3L_xL_yL_z}\sum_{p\notin{\rm boundary}}U_p \]

  • Contributions from the boundary spatial plaquettes is not included.
  • Normalized with the number of plaquettes Lx*Ly*Lz*(6Lt-3) and color factor 3.

\[ \frac{1}{3}\frac{1}{6L_xL_yL_zL_t}\sum_{p}U_p \]

  • Contributions from the boundary spatial plaquettes is included.
  • Normalized with the number of plaquettes 6*Lx*Ly*Lz*Lt and color factor 3.

Definition at line 949 of file staples_SF.cpp.

void Staples_SF::set_parameter_verboselevel ( const Bridge::VerboseLevel  vl)
inline

Definition at line 101 of file staples_SF.h.

void Staples_SF::set_parameters ( const Parameters &  params)

Definition at line 62 of file staples_SF.cpp.

void Staples_SF::set_parameters ( double *  phi,
double *  phipr 
)

Definition at line 114 of file staples_SF.cpp.

void Staples_SF::set_parameters ( const double *  phi,
const double *  phipr,
const double *  pomega 
)

Set the SF BC for wk, wkpr and omega0.

\[ U_k(x)|_{t=0}=W_k(\vec{x})=\exp\left(a C_k\right),\quad C_k=\frac{i}{L}\pmatrix{\phi_1\cr&\phi_2\cr&&\phi_3\cr}, \]

\[ U_k(x)|_{t=T}=W_k'(\vec{x})=\exp\left(a C_k'\right),\quad C_k'=\frac{i}{L}\pmatrix{\phi'_1\cr&\phi'_2\cr&&\phi'_3\cr}, \]

omega0 is set to the input value

\[ \verb|omega0|=i\Omega_0 \]

Definition at line 168 of file staples_SF.cpp.

void Staples_SF::set_parameters ( std::vector< double > &  phi,
std::vector< double > &  phipr,
std::vector< double > &  pomega 
)

Definition at line 87 of file staples_SF.cpp.

double Staples_SF::sf_coupling_plaq ( const Field_G &  U,
double  ct 
)

Evaluate boudary plaqutte VEV for the SF coupling. The SF running coupling is defined as

\[ \frac{k}{\overline{g}^2}= \left\langle\frac{\partial S_g}{\partial\eta}\right\rangle +\left\langle\frac{\partial S_f}{\partial\eta}\right\rangle \]

\[ k=12\left(\frac{L}{a}\right)^2 \left(\sin\theta+\sin\left(2\theta\right)\right), \quad \theta=\frac{1}{3}\pi\left(\frac{a^2}{TL}\right) \]

\[ \frac{\partial S_g}{\partial\eta}= -\frac{2}{g_0^2}\sum_{\vec{x}} c_t\sum_{k=1}^{3}\frac{a}{L}{\rm Re\ }{\rm tr}\Bigl( i\Omega_0W_k(\vec{x})U_0(\vec{x}+\hat{k},0)U_k^\dagger(\vec{x},1) U_0^\dagger(\vec{x},0) -i\Omega_0W_k'U_0^\dagger(\vec{x}+\hat{k},T-1)U_k^\dagger(\vec{x},T-1) U_0(\vec{x},T-1) \Bigr) \]

\[ \frac{\partial S_g}{\partial\eta}= -\frac{2}{g_0^2}\sum_{\vec{x}} c_t\sum_{k=1}^{3}\frac{a}{L}{\rm Re\ }{\rm tr}\Bigl( i\Omega_0\verb|upper|(3,k)(\vec{x},0)U_0^\dagger(\vec{x},0) -i\Omega_0\verb|upper|(3,k)(\vec{x},T-1)^\dagger U_0(\vec{x},T-1) \Bigr) \]

  • Evaluete the following quantity and print the result

    \[ -c_t\sum_{\vec{x}}\sum_{k=1}^{3}\frac{a}{L}{\rm Re\ }{\rm tr}\Bigl( i\Omega_0\verb|upper|(3,k)(\vec{x},0)U_0^\dagger(\vec{x},0) -i\Omega_0\verb|upper|(3,k)(\vec{x},T-1)^\dagger U_0(\vec{x},T-1) \Bigr) \]

  • Clover term contribution $\partial S_f/\partial\eta$ is not evaluated here.

Definition at line 248 of file staples_SF.cpp.

double Staples_SF::sf_coupling_rect ( const Field_G &  m_U,
double  ctr 
)

Evaluate the temporal rectangle VEV at the boundary for the SF running coupling.

The SF running coupling is given by

\[ \frac{k}{\overline{g}^2}= \left\langle\frac{\partial S_g^{\rm plaq}}{\partial\eta}\right\rangle +\left\langle\frac{\partial S_g^{\rm rect}}{\partial\eta}\right\rangle +\left\langle\frac{\partial S_f}{\partial\eta}\right\rangle \]

\[ \frac{\partial S_g}{\partial\eta}= -\frac{2}{g_0^2}c_1\sum_{\vec{x}}\sum_{k=1}^{3} \frac{ia}{L}{\rm Re}{\rm Tr}\Bigl( c_t^R\left( \Omega S_{kk0}(\vec{x},0) +\Omega S_{0kk}(\vec{x},0)\right) +\Omega S_{00k}(\vec{x},0) \Bigr) \]

\[ +\frac{2}{g_0^2}c_1\sum_{\vec{x}}\sum_{k=1}^{3} \frac{ia}{L}{\rm Re}{\rm Tr}\Bigl( c_t^R\left( \Omega'S_{kk0}(\vec{x},T) +\Omega'S_{0kk}(\vec{x},T)\right) +\Omega'S_{00k}(\vec{x},T) \Bigr) \]

\[ S_{kk0}(\vec{x},0)=W_kW_kU_0(\vec{x}+2\hat{k},0) U_k^\dagger(\vec{x}+\hat{k},1)U_k^\dagger(\vec{x},1)U_0^\dagger(\vec{x},0) \]

\[ S_{0kk}(\vec{x},0)=W_kU_0(\vec{x}+\hat{k},0)U_k^\dagger(\vec{x},1) U_k^\dagger(\vec{x}-\hat{k},1)U_0^\dagger(\vec{x}-\hat{k},0)W_k \]

\[ S_{00k}(\vec{x},0)=W_kU_0(\vec{x}+\hat{k},0)U_0(\vec{x}+\hat{k},1) U_k^\dagger(\vec{x},2)U_0^\dagger(\vec{x},1)U_0^\dagger(\vec{x},0) \]

\[ S_{kk0}(\vec{x},T)=W_k'W_k'U_0^\dagger(\vec{x}+2\hat{k},T) U_k^\dagger(\vec{x}+\hat{k},T-1)U_k^\dagger(\vec{x},T-1)U_0(\vec{x},T-1) \]

\[ S_{0kk}(\vec{x},T)=W_k'U_0^\dagger(\vec{x}+\hat{k},T-1) U_k^\dagger(\vec{x},T-1)U_k^\dagger(\vec{x}-\hat{k},T-1) U_0(\vec{x}-\hat{k},T-1)W_k' \]

\[ S_{00k}(\vec{x},T)=W_k'U_0^\dagger(\vec{x}+\hat{k},T-1) U_0^\dagger(\vec{x}+\hat{k},T-2)U_k^\dagger(\vec{x},T-2)U_0(\vec{x},T-2) U_0(\vec{x},T-1) \]

In this function we evaluate

\[ -\sum_{\vec{x}}\sum_{k=1}^{3} \frac{ia}{L}{\rm Re}{\rm Tr}\Bigl( c_t^R\left( \Omega S_{kk0}(\vec{x},0) +\Omega S_{0kk}(\vec{x},0)\right) +\Omega S_{00k}(\vec{x},0) \Bigr) \]

\[ +\sum_{\vec{x}}\sum_{k=1}^{3} \frac{ia}{L}{\rm Re}{\rm Tr}\Bigl( c_t^R\left( \Omega'S_{kk0}(\vec{x},T) +\Omega'S_{0kk}(\vec{x},T)\right) +\Omega'S_{00k}(\vec{x},T) \Bigr) \]

The following quantities are also printed out. 

  rect01         rect02
      <—<—+      <—<—+
      |       |      |       |
  t=0 x—>—+  t=0 +—x—+
  omega0               omega0
  rectt1         rectt2
  omega0              omega0
 t=Nt x—>—+ t=Nt +—x—+
      |       |      |       |
      +—<—+      <—<—+
  rect03      rectt3
              omega0
      +—+  t=Nt x—>
      |   |       |   |
      v   ^       ^   v
      |   |       |   |
  t=0 x—>       +—+
   omega0

Definition at line 422 of file staples_SF.cpp.

void Staples_SF::staple ( Field_G &  W,
const Field_G &  U,
const int  mu 
)

Evaluate staple for all the links in mu direction with SF BC. 

   (1)  mu (2)
      +-->--+
   nu |     |
     i+     +
      +     +
   nu |     |
     i+-->--+
    (1)  mu (2)
  • We should notice that all the staple attached to the boundary spatial link is set to zero.
  • This staple at the boundary shall not be used since the corresponding force for the boundary spatial link is always set to zero.

Definition at line 765 of file staples_SF.cpp.

void Staples_SF::staple_ct ( Field_G &  W,
const Field_G &  U,
const int  mu,
double  ct 
)

Evaluate staple for all the links in mu direction with SF BC and boundary improvement factor ct. 

   (1)  mu (2)
      +-->--+
   nu |     |
     i+     +
      +     +
   nu |     |
     i+-->--+
    (1)  mu (2)
  • staples attached to spatial boundary link is set to zero.
  • contributions from the temporal plaquette attached to the boundary is multiplied with ct.

Definition at line 805 of file staples_SF.cpp.

void Staples_SF::upper ( Field_G_SF &  c,
const Field_G &  U,
const int  mu,
const int  nu 
)

Evaluate upper staple for all the links in mu direction with SF BC. 

   (1)  mu (2)                 (2)
      +-->--+    c=(+-->--+)    +
   nu |     |                 v=|
     i+     +                   +
  • all upper staple attached to the t=0 boundary spatial link is set to zero.
    • if(mu!=3) c.set_boundary_zero()
  • no link exist at t=T boundary in this code. (t=T is identified with t=0 with PBC).
  • Implementation of the boundary link is given by set_boundary_wk(wk) and set_boundary_wkpr(wkpr).

Definition at line 860 of file staples_SF.cpp.

Member Data Documentation

const std::string Staples_SF::class_name = "Staples_SF"
static

Definition at line 50 of file staples_SF.h.

Index_lex Staples_SF::index
private

Definition at line 60 of file staples_SF.h.

int Staples_SF::initialized
private

Definition at line 67 of file staples_SF.h.

Mat_SU_N Staples_SF::iomega0
private

Definition at line 66 of file staples_SF.h.

int Staples_SF::Lt
private

Definition at line 64 of file staples_SF.h.

int Staples_SF::Lvol
private

Definition at line 58 of file staples_SF.h.

int Staples_SF::Lx
private

Definition at line 64 of file staples_SF.h.

int Staples_SF::Ly
private

Definition at line 64 of file staples_SF.h.

int Staples_SF::Lz
private

Definition at line 64 of file staples_SF.h.

Bridge::VerboseLevel Staples_SF::m_vl
protected

Definition at line 53 of file staples_SF.h.

int Staples_SF::Nc
private

Definition at line 56 of file staples_SF.h.

int Staples_SF::Ndim
private

Definition at line 57 of file staples_SF.h.

int Staples_SF::NPEt
private

Definition at line 65 of file staples_SF.h.

int Staples_SF::Nt
private

Definition at line 63 of file staples_SF.h.

int Staples_SF::Nvol
private

Definition at line 58 of file staples_SF.h.

int Staples_SF::Nx
private

Definition at line 63 of file staples_SF.h.

int Staples_SF::Ny
private

Definition at line 63 of file staples_SF.h.

int Staples_SF::Nz
private

Definition at line 63 of file staples_SF.h.

ShiftField_lex Staples_SF::shift
private

Definition at line 61 of file staples_SF.h.

Field_G_SF Staples_SF::Umu
private

Definition at line 59 of file staples_SF.h.

Field_G_SF Staples_SF::Unu
private

Definition at line 59 of file staples_SF.h.

Field_G_SF Staples_SF::v
private

Definition at line 59 of file staples_SF.h.

Field_G_SF Staples_SF::w
private

Definition at line 59 of file staples_SF.h.

Mat_SU_N Staples_SF::wk
private

Definition at line 66 of file staples_SF.h.

Mat_SU_N Staples_SF::wkpr
private

Definition at line 66 of file staples_SF.h.

The documentation for this class was generated from the following files: