FLUE_stoutSmearing
Module to do stout-link smearing on SU(3) gaugelinks i.e. Morningstar & Peardon: hep-lat/0311018
Source: src/gauge/stoutSmearing.f90
Dependencies
Contents
Subroutines
StoutSmearLinks
Stout-Smear the gaugelinks in U with a smearing weight of rho Do it nSweeps times output is in USmeared i.e. MorningStar & Peardon
Attributes: pure
fortran
subroutine StoutSmearLinks(U, rho, nSweeps, USmeared)Arguments
| Name | Type | Intent | Attributes | Description |
|---|---|---|---|---|
U | complex(kind=WC) | in | dimension(:, :, :, :, :, :, 😃 | Input gaugefield [3,3,4,nt,nx,ny,nz] |
rho | real(kind=WP) | in | Smearing strength | |
nSweeps | integer | in | Number of sweeps of smearing | |
USmeared | complex(kind=WC) | out | dimension(:, :, :, :, :, :, 😃 | Output smeared gaugefield [3,3,4,nt,nx,ny,nz] |
Call graph
StoutSmearOnce
Does a single sweep of stout-link smearing
Attributes: pure
fortran
subroutine StoutSmearOnce(U_old, rho, U_new)Arguments
| Name | Type | Intent | Attributes | Description |
|---|---|---|---|---|
U_old | complex(kind=WC) | in | dimension(:, :, :, :, :, :, 😃 | Input gaugefield [3,3,4,nt,nx,ny,nz] |
rho | real(kind=WP) | in | Smering strength | |
U_new | complex(kind=WC) | out | dimension(:, :, :, :, :, :, 😃 | Output gaugefield after 1 sweep of smearing [3,3,4,nt,nx,ny,nz] |
Call graph
ComputeSpatialStaples
Computes all of the spatial staples for the gaugefield with factor of smearing strength rho included Computes sum of backward/forward staples at given site/direction
Attributes: pure
fortran
subroutine ComputeSpatialStaples(data, rho, staples)Arguments
| Name | Type | Intent | Attributes | Description |
|---|---|---|---|---|
data | complex(kind=WC) | in | dimension(:, :, :, :, :, :, 😃 | Input gaugefield [3,3,4,nt,nx,ny,nz] |
rho | real(kind=WP) | in | Smearing strength rho | |
staples | complex(kind=WC) | out | dimension(:, :, :, :, :, :, 😃 | Output staples [3,3,3,nt,nx,ny,nz] |
Call graph
Functions
ComputeQMatrix
Computes the Q Matrix required Eqn2 of Morningstar & Peardon
Attributes: pure
Returns: complex(kind=WC)
fortran
function ComputeQMatrix(U, C) result(Q)Arguments
| Name | Type | Intent | Attributes | Description |
|---|---|---|---|---|
U | complex(kind=WC) | in | dimension(3, 3) | Input gauge link and sum of staples |
C | complex(kind=WC) | in | dimension(3, 3) | Input gauge link and sum of staples |
Call graph