#include "hila.h"

Include dependency graph for staples.h:

This graph shows which files directly or indirectly include this file:

Functions
template<typename T >
void	staplesum (const GaugeField< T > &U, Field< T > &staples, Direction d1, const plaqw_t< hila::arithmetic_type< T > > &plaqw, Parity par=ALL)
	Sum the staples of link matrices to direction dir taking into account plaquette weights.

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() [1/2]

template<typename T >

void staplesum	(	const GaugeField< T > &	U,
		Field< T > &	staples,
		Direction	d1,
		const plaqw_t< hila::arithmetic_type< T > > &	plaqw,
		Parity	par = `ALL`
	)

Sum the staples of link matrices to direction dir taking into account plaquette weights.

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]

But the method is computed in a slightly more optimized way

Template Parameters

T

Parameters

U	GaugeField to compute staples for
staples	Filed to compute staplesum into at each lattice point
d1	Direction to compute staplesum for
plaqw	plaquette weights
par	Parity to compute staplesum for

Definition at line 84 of file staples.h.

◆ staplesum() [2/2]