From bab88bc4f7bd59bbeb4e088e3bd933681e82c2e1 Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Thu, 6 May 2021 23:25:12 +0200
Subject: [PATCH] Fix compile

---
 tests/forces/Test_gp_rect_force.cc   | 11 ++++++-----
 tests/forces/Test_momentum_filter.cc |  5 ++++-
 2 files changed, 10 insertions(+), 6 deletions(-)

diff --git a/tests/forces/Test_gp_rect_force.cc b/tests/forces/Test_gp_rect_force.cc
index 98ebb2fa..627d264e 100644
--- a/tests/forces/Test_gp_rect_force.cc
+++ b/tests/forces/Test_gp_rect_force.cc
@@ -29,7 +29,6 @@ Author: paboyle <paboyle@ph.ed.ac.uk>
 
 using namespace std;
 using namespace Grid;
- ;
 
  
 
@@ -53,12 +52,14 @@ int main (int argc, char ** argv)
   pRNG.SeedFixedIntegers(std::vector<int>({45,12,81,9}));
 
   LatticeGaugeField U(&Grid);
-
   SU<Nc>::HotConfiguration(pRNG,U);
-  
+
   double beta = 1.0;
   double c1   = 0.331;
 
+  std::vector<int> twists(Nd,0);
+  twists[1] = 0;
+  ConjugateGimplD::setDirections(twists);
   ConjugatePlaqPlusRectangleActionR Action(beta,c1);
   //ConjugateWilsonGaugeActionR Action(beta);
   //WilsonGaugeActionR Action(beta);
@@ -67,13 +68,13 @@ int main (int argc, char ** argv)
 
   // get the deriv of phidag MdagM phi with respect to "U"
   LatticeGaugeField UdSdU(&Grid);
-
+  
   Action.deriv(U,UdSdU);
 
   ////////////////////////////////////
   // Modify the gauge field a little 
   ////////////////////////////////////
-  RealD dt = 0.01;
+  RealD dt = 0.001;
 
   LatticeColourMatrix mommu(&Grid); 
   LatticeColourMatrix forcemu(&Grid); 
diff --git a/tests/forces/Test_momentum_filter.cc b/tests/forces/Test_momentum_filter.cc
index 856ea0f2..c7888948 100644
--- a/tests/forces/Test_momentum_filter.cc
+++ b/tests/forces/Test_momentum_filter.cc
@@ -59,9 +59,12 @@ int main (int argc, char ** argv)
   std::cout<<GridLogMessage << "Grid is setup to use "<<threads<<" threads"<<std::endl;
 
   std::vector<int> seeds({1,2,3,4});
+  std::vector<int> serial_seeds({5,6,7,8});
 
+  GridSerialRNG            sRNG; 
   GridParallelRNG          pRNG(&Grid);
   pRNG.SeedFixedIntegers(seeds);
+  sRNG.SeedFixedIntegers(serial_seeds);
 
   typedef PeriodicGimplR Gimpl;
   typedef WilsonGaugeAction<Gimpl> GaugeAction;
@@ -115,7 +118,7 @@ int main (int argc, char ** argv)
   
   integrator.setMomentumFilter(filter);
 
-  integrator.refresh(U, pRNG); //doesn't actually change the gauge field
+  integrator.refresh(U, sRNG, pRNG); //doesn't actually change the gauge field
 
   //Check the momentum is zero on the boundary
   const auto &P = integrator.getMomentum();