Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members | Related Pages

cps::GimprOLSym Class Reference
[Gauge Actions]

A class implementing a lattice with the 1-loop Symanzik improved gauge action. More...

#include <lattice.h>

Inherits cps::Lattice.

Inherited by cps::GimprOLSymFasqtad, cps::GimprOLSymFclover, cps::GimprOLSymFdwf, cps::GimprOLSymFnone, cps::GimprOLSymFp4, cps::GimprOLSymFstag, and cps::GimprOLSymFwilson.

List of all members.

Public Member Functions

 GimprOLSym ()
virtual ~GimprOLSym ()
GclassType Gclass ()
 Returns the type of gauge action.
void GactionGradient (Matrix &grad, int *x, int mu)
 Calculates the partial derivative of the gauge action w.r.t. the link U_mu(x).
void GforceSite (Matrix &force, int *x, int mu)
 Calculates the gauge force at site x and direction mu.
ForceArg EvolveMomGforce (Matrix *mom, Float step_size)
 Molecular dynamics evolution of the conjugate momentum.
Float GhamiltonNode ()
 Computes the pure gauge action on the local sublattice.
void AllStaple (Matrix &stap, const int *x, int mu)
 Computes the sum of all the staples around a link.

Detailed Description

A class implementing a lattice with the 1-loop Symanzik improved gauge action.

This is a sum of the plaquette, rectangle and cube loops:

\[ \sum_x\sum_\mu\sum_{\nu>\mu} \left[\right. U_\mu(x) U_\nu(x+\nu) U^\dagger_\mu(x+\nu) U^\dagger_\nu(x) \]

\[ -\frac{ 1 + 0.4805 \alpha_s}{20 u_0^2} U_\mu(x) U_\mu(x+\mu) U_\nu(x+2\mu) U^\dagger_\mu(x+\mu+\nu) U^\dagger_\mu(x+\nu) U^\dagger_\nu(x) \]

\[ -\frac{ 0.03325 \alpha_s}{u_0^2} \sum_{\rho>\nu} U_\mu(x) U_\nu(x+\mu) U_\rho(x+\mu+\nu) U^\dagger_\mu(x+\mu+\nu+\rho) U^\dagger_\nu(x+\nu+\rho) U^\dagger_\rho(x+\rho) \left.\right] \]

where $u_0$ is the tadpole coefficient and $\alpha_s = -4\log(u_0)/3.06839$.

Definition at line 1507 of file lattice.h.

Constructor & Destructor Documentation

cps::GimprOLSym::GimprOLSym (   ) 
 

Definition at line 78 of file g_impr_OLSym.C.

References cps::GlobalJobParameter::Beta(), cps::ERR, cps::Float, cps::Verbose::Func(), cps::GJP, cps::Error::NotImplemented(), cps::GlobalJobParameter::u0(), cps::VRB, and cps::GlobalJobParameter::XiBare().

cps::GimprOLSym::~GimprOLSym (   )  [virtual]
 

Definition at line 101 of file g_impr_OLSym.C.

References cps::Verbose::Func(), and cps::VRB.

Member Function Documentation

void cps::GimprOLSym::AllStaple (  Matrix &  stap,
const int *  x,
int  mu
)  [virtual]
 

Computes the sum of all the staples around a link.

The staple sum around the link $ U_\mu(x) $ is

\[ \sum_{\pm \nu, |\nu|\neq \mu} \left\{\right. U_\nu(x+\mu) U^\dagger_\mu(x+\nu) U^\dagger_\nu(x) \]

\[ -\frac{ 1 + 0.4805 \alpha_s}{20 u_0^2} \left[\right. U_\mu(x+\mu) U_\nu(x+2\mu) U^\dagger_\mu(x+\mu+\nu) U^\dagger_\mu(x+\nu) U^\dagger_\nu(x) \]

\[ + U_\nu(x+\mu) U^\dagger_\mu(x+\nu) U^\dagger_\mu(x-\mu+\nu) U^\dagger_\nu(x-\mu) U_\mu(x-\mu) \]

\[ + U_\nu(x+\mu) U_\nu(x+\mu+\nu) U^\dagger_\mu(x+2\nu) U^\dagger_\nu(x+\nu) U^\dagger_\nu(x) \left.\right]\]

\[ -\frac{ 0.03325 \alpha_s}{u_0^2} \sum_{\pm \rho, |\rho|\neq \nu, |\rho|\neq \mu} U_\nu(x+\mu) U_\rho(x+\mu+\nu) U_{-\mu}(x+\mu+\nu+\rho) U_{-\nu}(x+\nu+\rho) U_{-\rho}(x+\rho) \left.\right\} \]

where $u_0$ is the tadpole coefficient and $\alpha_s = -4\log(u_0)/3.06839$.

Parameters:
stap The computed staple sum.
x The coordinates of the lattice site
mu The link direction

Implements cps::Lattice.

Definition at line 287 of file g_impr_OLSym.C.

References cps::Lattice::BufferedCubeStaple(), cps::Lattice::BufferedRectStaple(), cps::Lattice::BufferedStaple(), cps::Verbose::Func(), cps::IFloat, MATRIX_SIZE, cps::vecAddEquVec(), cps::vecTimesEquFloat(), and cps::VRB.

Referenced by GactionGradient(), and GforceSite().

ForceArg cps::GimprOLSym::EvolveMomGforce (  Matrix *  mom,
Float  step_size
)  [virtual]
 

Molecular dynamics evolution of the conjugate momentum.

The momentum is evolved for a single molecular dynamics timestep using the force from the pure gauge action.

Parameters:
mom The momentum matrices on all links.
step_size The molecular dynamics timestep used in the numerical integration.
Postcondition:
mom is assigned the value of the momentum after the molecular dynamics evolution.

Implements cps::Lattice.

Definition at line 19 of file g_impr_OLSym_force.C.

References cps::Matrix::Dagger(), cps::dclock(), cps::ffree(), cps::Float, cps::fmalloc(), cps::fTimesV1PlusV2(), cps::Verbose::Func(), cps::GJP, cps::glb_max(), cps::glb_sum(), cps::Lattice::GsiteOffset(), cps::IFloat, cps::N, cps::Matrix::norm(), cps::print_flops(), cps::ParTransGauge::run(), cps::GlobalJobParameter::TnodeSites(), cps::Matrix::TrLessAntiHermMatrix(), cps::vaxpy3_m(), cps::GlobalJobParameter::VolNodeSites(), cps::GlobalJobParameter::VolSites(), cps::VRB, cps::GlobalJobParameter::XnodeSites(), cps::GlobalJobParameter::YnodeSites(), and cps::GlobalJobParameter::ZnodeSites().

void cps::GimprOLSym::GactionGradient (  Matrix &  grad,
int *  x,
int  mu
)  [virtual]
 

Calculates the partial derivative of the gauge action w.r.t. the link U_mu(x).

Parameters:
grad The computed gradient.
x The coordinates of the lattice site.
mu The direction of the link.

Implements cps::Lattice.

Definition at line 233 of file g_impr_OLSym.C.

References AllStaple(), cps::Verbose::Func(), cps::IFloat, MATRIX_SIZE, cps::vecTimesEquFloat(), and cps::VRB.

GclassType cps::GimprOLSym::Gclass (   )  [virtual]
 

Returns the type of gauge action.

Implements cps::Lattice.

Definition at line 112 of file g_impr_OLSym.C.

References cps::GclassType.

void cps::GimprOLSym::GforceSite (  Matrix &  force,
int *  x,
int  mu
) 
 

Calculates the gauge force at site x and direction mu.

Definition at line 125 of file g_impr_OLSym.C.

References AllStaple(), cps::CBUF_MODE4, cps::Matrix::Dagger(), cps::Verbose::Func(), cps::Lattice::GaugeField(), cps::Lattice::GsiteOffset(), cps::IFloat, MATRIX_SIZE, cps::mDotMEqual(), cps::moveMem(), cps::setCbufCntrlReg(), cps::Matrix::TrLessAntiHermMatrix(), cps::vecTimesEquFloat(), and cps::VRB.

Float cps::GimprOLSym::GhamiltonNode (   )  [virtual]
 

Computes the pure gauge action on the local sublattice.

Returns:
The value of the pure gauge action on this node.

Implements cps::Lattice.

Definition at line 213 of file g_impr_OLSym.C.

References cps::Float, cps::Verbose::Func(), cps::Lattice::SumReTrCubeNode(), cps::Lattice::SumReTrPlaqNode(), cps::Lattice::SumReTrRectNode(), and cps::VRB.

The documentation for this class was generated from the following files:
Generated on Sat Oct 10 14:12:12 2009 for Columbia Physics System by  doxygen 1.3.9.1