From 178376f24bb99a99f268573a0239f88cd897e3b1 Mon Sep 17 00:00:00 2001
From: Alessandro Lupo <Alessandro.Lupo@ed.ac.uk>
Date: Thu, 6 Apr 2023 12:08:17 +0100
Subject: [PATCH] minor stylistic changes

---
 Grid/qcd/action/gauge/GaugeImplTypes.h | 21 ++++-----------------
 Grid/qcd/utils/GaugeGroup.h            |  2 +-
 tests/core/Test_reunitarise.cc         | 12 ++++--------
 tests/sp2n/Test_project_on_Sp.cc       | 14 ++++++--------
 4 files changed, 15 insertions(+), 34 deletions(-)

diff --git a/Grid/qcd/action/gauge/GaugeImplTypes.h b/Grid/qcd/action/gauge/GaugeImplTypes.h
index 17bee86a..e2d98b85 100644
--- a/Grid/qcd/action/gauge/GaugeImplTypes.h
+++ b/Grid/qcd/action/gauge/GaugeImplTypes.h
@@ -133,19 +133,6 @@ public:
     
   static inline Field projectForce(Field &P)
     {
-        /*  // this works
-        iScalar<iScalar<iMatrix<Complex, Nc> > > gen;
-        auto Psum = P;
-        Psum = Zero();
-        for (int a = 0; a < Group::AlgebraDimension; a++)
-            {
-                Group::generator(a, gen);
-                auto coeff = 2. * trace(P * gen);
-                Psum += coeff * gen;
-            
-                   
-        return Psum;}*/
-        //this doesnt
         Field ret(P.Grid());
         Group::taProj(P, ret);
         return ret;
@@ -159,7 +146,7 @@ public:
     autoView(P_v,P,AcceleratorRead);
     accelerator_for(ss, P.Grid()->oSites(),1,{
       for (int mu = 0; mu < Nd; mu++) {
-          U_v[ss](mu) = Exponentiate(P_v[ss](mu), ep, Nexp) * U_v[ss](mu);  //
+          U_v[ss](mu) = Exponentiate(P_v[ss](mu), ep, Nexp) * U_v[ss](mu);
           Group::ProjectOnGaugeGroup(U_v[ss](mu));
       }
     });
@@ -206,9 +193,9 @@ typedef GaugeImplTypes<vComplex, Nc> GimplTypesR;
 typedef GaugeImplTypes<vComplexF, Nc> GimplTypesF;
 typedef GaugeImplTypes<vComplexD, Nc> GimplTypesD;
 
-typedef GaugeImplTypes<vComplex, Nc, 12, Sp<Nc>> SymplGimplTypesR;
-typedef GaugeImplTypes<vComplexF, Nc, 12, Sp<Nc>> SymplGimplTypesF;
-typedef GaugeImplTypes<vComplexD, Nc, 12, Sp<Nc>> SymplGimplTypesD;
+typedef GaugeImplTypes<vComplex, Nc, 12, Sp<Nc> > SymplGimplTypesR;
+typedef GaugeImplTypes<vComplexF, Nc, 12, Sp<Nc> > SymplGimplTypesF;
+typedef GaugeImplTypes<vComplexD, Nc, 12, Sp<Nc> > SymplGimplTypesD;
 
 typedef GaugeImplTypes<vComplex, SU<Nc>::AdjointDimension> GimplAdjointTypesR;
 typedef GaugeImplTypes<vComplexF, SU<Nc>::AdjointDimension> GimplAdjointTypesF;
diff --git a/Grid/qcd/utils/GaugeGroup.h b/Grid/qcd/utils/GaugeGroup.h
index 8e6c5ee4..0e8fbabe 100644
--- a/Grid/qcd/utils/GaugeGroup.h
+++ b/Grid/qcd/utils/GaugeGroup.h
@@ -170,7 +170,7 @@ class GaugeGroup {
     }
   }
 
- template <typename LatticeMatrixType>
+  template <typename LatticeMatrixType>
   static void LieRandomize(GridParallelRNG &pRNG, LatticeMatrixType &out,
                            double scale = 1.0) {
     GridBase *grid = out.Grid();
diff --git a/tests/core/Test_reunitarise.cc b/tests/core/Test_reunitarise.cc
index e421a2a8..110ebabf 100644
--- a/tests/core/Test_reunitarise.cc
+++ b/tests/core/Test_reunitarise.cc
@@ -122,15 +122,16 @@ int main (int argc, char ** argv)
   detU=detU-1.0;
   std::cout << "Determinant defect before projection " <<norm2(detU)<<std::endl;
   tmp = U*adj(U) - ident;
-  std::cout << "Unitarity check before projection    " << norm2(tmp)<<std::endl; 
-  #if Nc==3
+  std::cout << "Unitarity check before projection    " << norm2(tmp)<<std::endl;
+    
+#if Nc==3
   ProjectSU3(U);
   detU= Determinant(U) ;
   detU= detU -1.0;
   std::cout << "Determinant ProjectSU3 defect " <<norm2(detU)<<std::endl;
   tmp = U*adj(U) - ident;
   std::cout << "Unitarity check after projection    " << norm2(tmp)<<std::endl; 
-  #endif
+#endif
   ProjectSUn(UU);
   detUU= Determinant(UU);
   detUU= detUU -1.0;
@@ -140,8 +141,3 @@ int main (int argc, char ** argv)
   
   Grid_finalize();
 }
-
-
-
-
-
diff --git a/tests/sp2n/Test_project_on_Sp.cc b/tests/sp2n/Test_project_on_Sp.cc
index 516a7f93..13c2bc17 100644
--- a/tests/sp2n/Test_project_on_Sp.cc
+++ b/tests/sp2n/Test_project_on_Sp.cc
@@ -31,6 +31,7 @@ int main (int argc, char **argv)
     
     identity = 1.0;
     RealD epsilon = 0.01;
+    RealD Delta = 666.;
     Complex i(0., 1.);
     RealD u = 0.;
     double vol = Umu.Grid()->gSites();
@@ -118,7 +119,7 @@ int main (int argc, char **argv)
     
     
     std::cout << GridLogMessage << "Testing ProjectOnGaugeGroup" << std::endl;
-    U = U + 2932.111*identity;
+    U = U + Delta*identity;
     std::cout << GridLogMessage << "Apply ProjectOnGaugeGroup to deformed matrix" << std::endl;
     Sp<Nc>::ProjectOnGaugeGroup(U);
     aux = U*adj(U) - identity;
@@ -162,7 +163,7 @@ int main (int argc, char **argv)
     }
     
     std::cout << GridLogMessage << "Testing ProjectGn" << std::endl;
-    U = U + 2932.111*identity;
+    U = U + Delta*identity;
     std::cout << GridLogMessage << "Apply ProjectGn to deformed matrix" << std::endl;
     Sp<Nc>::ProjectGn(U);
     aux = U*adj(U) - identity;
@@ -214,7 +215,7 @@ int main (int argc, char **argv)
     std::cout << GridLogMessage << "Testing SpTa" << std::endl;
     
     U = PeekIndex<LorentzIndex>(Umu,1);
-    U = U + 666.*identity;
+    U = U + Delta*identity;
     std::cout << GridLogMessage << "Matrix deformed " << std::endl;
     std::cout << GridLogMessage << "Apply SpTa to deformed matrix" << std::endl;
     U = SpTa(U);
@@ -279,7 +280,7 @@ int main (int argc, char **argv)
     
     // test taProj
     std::cout << GridLogMessage << "Testing taProj" << std::endl;
-    U = U + 666.*identity;
+    U = U + Delta*identity;
     std::cout << GridLogMessage << "Matrix deformed " << std::endl;
     std::cout << GridLogMessage << "Apply taProj to deformed matrix" << std::endl;
     Sp<Nc>::taProj(U, Up);
@@ -333,10 +334,7 @@ int main (int argc, char **argv)
             assert(  amizeroo.real() < 10e-6 );
         }
     }
-    
-    
+        
     Grid_finalize();
 
-
 }
-