From 80176b1b394e343ade8baf4e879caf6ebe906739 Mon Sep 17 00:00:00 2001
From: Christopher Kelly <ckelly@bnl.gov>
Date: Tue, 4 May 2021 13:12:23 -0400
Subject: [PATCH] RHMC now outputs some initial norms to the logs Fixed DWF+I
 Gparity binaries not correctly assigning twist directions (thanks Peter!)

---
 Grid/algorithms/iterative/ConjugateGradientMultiShift.h | 3 +++
 HMC/DWF2p1fIwasakiGparity.cc                            | 2 +-
 HMC/DWF2p1fIwasakiGparityRHMCdouble.cc                  | 2 +-
 3 files changed, 5 insertions(+), 2 deletions(-)

diff --git a/Grid/algorithms/iterative/ConjugateGradientMultiShift.h b/Grid/algorithms/iterative/ConjugateGradientMultiShift.h
index 252e4506..3b079e99 100644
--- a/Grid/algorithms/iterative/ConjugateGradientMultiShift.h
+++ b/Grid/algorithms/iterative/ConjugateGradientMultiShift.h
@@ -182,6 +182,9 @@ public:
     for(int s=0;s<nshift;s++) {
       axpby(psi[s],0.,-bs[s]*alpha[s],src,src);
     }
+
+    std::cout << GridLogIterative << "ConjugateGradientMultiShift: initial rn (|src|^2) =" << rn << " qq (|MdagM src|^2) =" << qq << " d ( dot(src, [MdagM + m_0]src) ) =" << d << " c=" << c << std::endl;
+    
   
   ///////////////////////////////////////
   // Timers
diff --git a/HMC/DWF2p1fIwasakiGparity.cc b/HMC/DWF2p1fIwasakiGparity.cc
index e9865a54..c02ada51 100644
--- a/HMC/DWF2p1fIwasakiGparity.cc
+++ b/HMC/DWF2p1fIwasakiGparity.cc
@@ -354,7 +354,7 @@ int main(int argc, char **argv) {
  
   //Setup the BCs
   FermionActionD::ImplParams Params;
-  for(int i=0;i<Nd-1;i++) Params.twists = user_params.GparityDirs[i]; //G-parity directions
+  for(int i=0;i<Nd-1;i++) Params.twists[i] = user_params.GparityDirs[i]; //G-parity directions
   Params.twists[Nd-1] = 1; //APBC in time direction
 
   std::vector<int> dirs4(Nd);
diff --git a/HMC/DWF2p1fIwasakiGparityRHMCdouble.cc b/HMC/DWF2p1fIwasakiGparityRHMCdouble.cc
index 533ad8a2..58295d5e 100644
--- a/HMC/DWF2p1fIwasakiGparityRHMCdouble.cc
+++ b/HMC/DWF2p1fIwasakiGparityRHMCdouble.cc
@@ -354,7 +354,7 @@ int main(int argc, char **argv) {
  
   //Setup the BCs
   FermionActionD::ImplParams Params;
-  for(int i=0;i<Nd-1;i++) Params.twists = user_params.GparityDirs[i]; //G-parity directions
+  for(int i=0;i<Nd-1;i++) Params.twists[i] = user_params.GparityDirs[i]; //G-parity directions
   Params.twists[Nd-1] = 1; //APBC in time direction
 
   std::vector<int> dirs4(Nd);