From a9a3248cb571f5c7306043ea4bdcff2d70146cfe Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Sun, 24 Feb 2019 15:21:19 +0000
Subject: [PATCH] More precision in rect force test

---
 tests/forces/Test_rect_force.cc | 14 +++++++++++---
 1 file changed, 11 insertions(+), 3 deletions(-)

diff --git a/tests/forces/Test_rect_force.cc b/tests/forces/Test_rect_force.cc
index a846af30..e0ffd28c 100644
--- a/tests/forces/Test_rect_force.cc
+++ b/tests/forces/Test_rect_force.cc
@@ -57,9 +57,10 @@ int main (int argc, char ** argv)
   SU3::HotConfiguration(pRNG,U);
   
   double beta = 1.0;
-  double c1   = 0.331;
+  double c1   = -0.331;
 
-  PlaqPlusRectangleActionR Action(beta,c1);
+  IwasakiGaugeActionR Action(beta);
+  //  PlaqPlusRectangleActionR Action(beta,c1);
   //  WilsonGaugeActionR Action(beta);
 
   ComplexD S    = Action.S(U);
@@ -87,7 +88,13 @@ int main (int argc, char ** argv)
 
     // fourth order exponential approx
     parallel_for(auto i=mom.begin();i<mom.end();i++){ // exp(pmu dt) * Umu
-      Uprime[i](mu) = U[i](mu) + mom[i](mu)*U[i](mu)*dt ;
+      Uprime[i](mu) = U[i](mu) + mom[i](mu)*U[i](mu)*dt 
+	+ mom[i](mu) *mom[i](mu) *U[i](mu)*(dt*dt/2.0)
+	+ mom[i](mu) *mom[i](mu) *mom[i](mu) *U[i](mu)*(dt*dt*dt/6.0)
+	+ mom[i](mu) *mom[i](mu) *mom[i](mu) *mom[i](mu) *U[i](mu)*(dt*dt*dt*dt/24.0)
+	+ mom[i](mu) *mom[i](mu) *mom[i](mu) *mom[i](mu) *mom[i](mu) *U[i](mu)*(dt*dt*dt*dt*dt/120.0)
+	+ mom[i](mu) *mom[i](mu) *mom[i](mu) *mom[i](mu) *mom[i](mu) *mom[i](mu) *U[i](mu)*(dt*dt*dt*dt*dt*dt/720.0);
+
     }
   }
 
@@ -114,6 +121,7 @@ int main (int argc, char ** argv)
   }
   ComplexD dSpred    = sum(dS);
 
+  std::cout << std::setprecision(15)<<std::endl;
   std::cout << GridLogMessage << " S      "<<S<<std::endl;
   std::cout << GridLogMessage << " Sprime "<<Sprime<<std::endl;
   std::cout << GridLogMessage << "dS      "<<Sprime-S<<std::endl;