portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-04 03:05:55 +01:00
Code Releases Activity

Compile fix

Browse Source
This commit is contained in:
Azusa Yamaguchi 2015-06-05 10:29:42 +01:00
parent bb36139fc2
commit 5f33cc3a95
2 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
lib/qcd/utils/WilsonLoops.h
View File

@ -120,6 +120,10 @@ void siteRectangle(GaugeMat &plaq,const std::vector<GaugeMat> &U, const int mu,
*/ */
}; };
typedef WilsonLoops<LatticeColourMatrix,LatticeGaugeField> ColourWilsonLoops;
}} }}

10
tests/Test_GaugeAction.cc
View File

@ -1,12 +1,12 @@
#include <Grid.h> #include <Grid.h>
#include <qcd/utils/CovariantCshift.h>
#include <qcd/utils/WilsonLoops.h>
using namespace std; using namespace std;
using namespace Grid; using namespace Grid;
using namespace Grid::QCD; using namespace Grid::QCD;
int main (int argc, char ** argv) int main (int argc, char ** argv)
{ {
Grid_init(&argc,&argv); Grid_init(&argc,&argv);
@ -96,14 +96,14 @@ int main (int argc, char ** argv)
Complex p = TensorRemove(Tp); Complex p = TensorRemove(Tp);
std::cout << "calculated plaquettes " <<p*PlaqScale<<std::endl; std::cout << "calculated plaquettes " <<p*PlaqScale<<std::endl;
RealD avg_plaq = WilsonLoops::avgPlaquette(Umu); RealD avg_plaq = ColourWilsonLoops::avgPlaquette(Umu);
std::cout << "NEW : calculated real plaquettes " <<avg_plaq<<std::endl; std::cout << "NEW : calculated real plaquettes " <<avg_plaq<<std::endl;
RealD stap_plaq=0.0; RealD stap_plaq=0.0;
LatticeColourMatrix stap(&Fine); LatticeColourMatrix stap(&Fine);
LatticeComplex stap_tr(&Fine); LatticeComplex stap_tr(&Fine);
for(int mu=0;mu<Nd;mu++){ for(int mu=0;mu<Nd;mu++){
WilsonLoops::Staple(stap,Umu,mu); ColourWilsonLoops::Staple(stap,Umu,mu);
stap_tr = trace(stap*adj(U[mu])); stap_tr = trace(stap*adj(U[mu]));
TComplex Ts = sum(stap_tr); TComplex Ts = sum(stap_tr);
Complex s = TensorRemove(Ts); Complex s = TensorRemove(Ts);