From 10f2c2530b557a601fe50d9e5ac53bec6bca00a2 Mon Sep 17 00:00:00 2001
From: Quadro <Peter Boyle>
Date: Tue, 25 May 2021 21:08:09 -0400
Subject: [PATCH] Temporary debug / display force structure

---
 Grid/qcd/hmc/integrators/Integrator.h | 10 ++++++++++
 1 file changed, 10 insertions(+)

diff --git a/Grid/qcd/hmc/integrators/Integrator.h b/Grid/qcd/hmc/integrators/Integrator.h
index 229f7393..5d875c6a 100644
--- a/Grid/qcd/hmc/integrators/Integrator.h
+++ b/Grid/qcd/hmc/integrators/Integrator.h
@@ -137,6 +137,8 @@ protected:
       force = FieldImplementation::projectForce(force); // Ta for gauge fields
       double end_force = usecond();
 
+      MomFilter->applyFilter(force);
+      
       Real force_abs   = std::sqrt(norm2(force)/U.Grid()->gSites()); //average per-site norm.  nb. norm2(latt) = \sum_x norm2(latt[x]) 
       Real impulse_abs = force_abs * ep * HMC_MOMENTUM_DENOMINATOR;    
 
@@ -153,6 +155,14 @@ protected:
       double time_full  = (end_full - start_full) / 1e3;
       double time_force = (end_force - start_force) / 1e3;
       std::cout << GridLogMessage << "["<<level<<"]["<<a<<"] P update elapsed time: " << time_full << " ms (force: " << time_force << " ms)"  << std::endl;
+
+      DumpSliceNorm("force",force,Nd-1);
+      auto pol = PeekIndex<LorentzIndex>(force,Nd-1);
+      DumpSliceNorm("force_t",pol);
+      pol=Zero();
+      PokeIndex<LorentzIndex>(force,pol,Nd-1);
+      DumpSliceNorm("force_xyz",force);
+      
     }
 
     // Force from the other representations