From 491a7082258fb61309e54bde0906d2f30d20fab8 Mon Sep 17 00:00:00 2001
From: paboyle <paboyle@ph.ed.ac.uk>
Date: Sat, 2 Jan 2016 14:02:41 +0000
Subject: [PATCH] Twist params set up correctly in gparity even odd 2f DWF +
 Wilson HMC test

---
 tests/Test_hmc_EODWFRatio_Gparity.cc | 10 ++++++++--
 1 file changed, 8 insertions(+), 2 deletions(-)

diff --git a/tests/Test_hmc_EODWFRatio_Gparity.cc b/tests/Test_hmc_EODWFRatio_Gparity.cc
index 521c496c..b75b3bb5 100644
--- a/tests/Test_hmc_EODWFRatio_Gparity.cc
+++ b/tests/Test_hmc_EODWFRatio_Gparity.cc
@@ -33,11 +33,17 @@ public:
     // Gauge action
     ConjugateWilsonGaugeActionR Waction(5.6);
 
+    // Fermion action
+    const int nu = 3;
+    std::vector<int> twists(Nd,0);
+    twists[nu] = 1;
+    FermionAction::ImplParams params;
+    params.twists = twists;
     Real mass=0.04;
     Real pv  =1.0;
     RealD M5=1.5;
-    FermionAction DenOp(U,*FGrid,*FrbGrid,*UGrid,*UrbGrid,mass,M5);
-    FermionAction NumOp(U,*FGrid,*FrbGrid,*UGrid,*UrbGrid,pv,M5);
+    FermionAction DenOp(U,*FGrid,*FrbGrid,*UGrid,*UrbGrid,mass,M5,params);
+    FermionAction NumOp(U,*FGrid,*FrbGrid,*UGrid,*UrbGrid,pv,M5,params);
   
     ConjugateGradient<FermionField>  CG(1.0e-8,10000);
     TwoFlavourEvenOddRatioPseudoFermionAction<ImplPolicy> Nf2(NumOp, DenOp,CG,CG);