diff --git a/lib/qcd/utils/WilsonLoops.h b/lib/qcd/utils/WilsonLoops.h
index f2d6d4de..82230aa8 100644
--- a/lib/qcd/utils/WilsonLoops.h
+++ b/lib/qcd/utils/WilsonLoops.h
@@ -120,6 +120,10 @@ void siteRectangle(GaugeMat &plaq,const std::vector<GaugeMat> &U, const int mu,
   */
 
 };
+
+
+ typedef WilsonLoops<LatticeColourMatrix,LatticeGaugeField> ColourWilsonLoops;
+
 }}
 
 
diff --git a/tests/Test_GaugeAction.cc b/tests/Test_GaugeAction.cc
index 07e7a4eb..0dde6d8f 100644
--- a/tests/Test_GaugeAction.cc
+++ b/tests/Test_GaugeAction.cc
@@ -1,12 +1,12 @@
 #include <Grid.h>
 
+#include <qcd/utils/CovariantCshift.h>
+#include <qcd/utils/WilsonLoops.h>
+
 using namespace std;
 using namespace Grid;
 using namespace Grid::QCD;
 
-
-
-
 int main (int argc, char ** argv)
 {
   Grid_init(&argc,&argv);
@@ -96,14 +96,14 @@ int main (int argc, char ** argv)
   Complex p  = TensorRemove(Tp);
   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;
 
   RealD stap_plaq=0.0;
   LatticeColourMatrix stap(&Fine);
   LatticeComplex stap_tr(&Fine);
   for(int mu=0;mu<Nd;mu++){
-    WilsonLoops::Staple(stap,Umu,mu);
+    ColourWilsonLoops::Staple(stap,Umu,mu);
     stap_tr = trace(stap*adj(U[mu]));
     TComplex Ts = sum(stap_tr);
     Complex s  = TensorRemove(Ts);