From 8b070ae54ce83f39cda2ac89aa1bd059abe8d1ce Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Wed, 19 Aug 2015 11:26:01 +0100
Subject: [PATCH] Gparity now accepting twists through constructor

---
 tests/Test_gparity.cc | 9 +++++++--
 1 file changed, 7 insertions(+), 2 deletions(-)

diff --git a/tests/Test_gparity.cc b/tests/Test_gparity.cc
index 3a5a5dd8..a301861e 100644
--- a/tests/Test_gparity.cc
+++ b/tests/Test_gparity.cc
@@ -13,6 +13,9 @@ int main (int argc, char ** argv)
 
   Grid_init(&argc,&argv);
 
+  std::vector<int> twists(Nd,0);
+  twists[nu] = 1;
+
   const int Ls=4;
   const int L =4;
   std::vector<int> latt_2f(Nd,L);
@@ -95,8 +98,10 @@ int main (int argc, char ** argv)
   SchurDiagMooeeOperator<DomainWallFermionR,LatticeFermion> HermOpEO(Ddwf);
   ConjugateGradient<LatticeFermion> CG(1.0e-8,10000);
   CG(HermOpEO,src_o_1f,result_o_1f);
-
-  GparityDomainWallFermionR GPDdwf(Umu_2f,*FGrid_2f,*FrbGrid_2f,*UGrid_2f,*UrbGrid_2f,mass,M5);
+  
+  GparityDomainWallFermionR::ImplParams params;
+  params.twists = twists;
+  GparityDomainWallFermionR GPDdwf(Umu_2f,*FGrid_2f,*FrbGrid_2f,*UGrid_2f,*UrbGrid_2f,mass,M5,params);
 
   for(int disp=-1;disp<=1;disp+=2)
   for(int mu=0;mu<5;mu++)