HILA
Loading...
Searching...
No Matches
staples.h File Reference
#include "hila.h"
Include dependency graph for staples.h:
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions

template<typename T >
void staplesum (const GaugeField< T > &U, Field< T > &staples, Direction d1, Parity par=ALL)
 Sum the staples of link matrices to direction dir.
 

Function Documentation

◆ staplesum()

template<typename T >
void staplesum ( const GaugeField< T > &  U,
Field< T > &  staples,
Direction  d1,
Parity  par = ALL 
)

Sum the staples of link matrices to direction dir.

Naive method is to compute:

foralldir(d2) if (d2 != d1)
stapes[par] += U[d2][X]*U[d1][X+d2]*U[d2][X+d1].dagger() +
U[d2][X-d2].dagger()*U[d1][X-d2]*U[d2][X-d2+d1]
#define foralldir(d)
Macro to loop over (all) Direction(s)
Definition coordinates.h:78

But the method is computed in a slightly more optimized way

Template Parameters
T
Parameters
UGaugeField to compute staples for
staplesFiled to compute staplesum into at each lattice point
d1Direction to compute staplesum for
parParity to compute staplesum for

Definition at line 28 of file staples.h.