diff --git a/lib/parallelIO/BinaryIO.h b/lib/parallelIO/BinaryIO.h
index 3265ba54..9dab4d24 100644
--- a/lib/parallelIO/BinaryIO.h
+++ b/lib/parallelIO/BinaryIO.h
@@ -372,6 +372,8 @@ class BinaryIO {
     if (grid->IsBoss()) 
       fout.close();
 
+    std::cout << GridLogMessage << "RNG file checksum " << csum << std::endl;
+
     return csum;
   }
   static inline uint32_t readRNGSerial(GridSerialRNG &serial,GridParallelRNG &parallel,std::string file,int offset)
@@ -421,6 +423,8 @@ class BinaryIO {
       Uint32Checksum((uint32_t *)&saved[0],bytes,csum);
     }
 
+    std::cout << GridLogMessage << "RNG file checksum " << csum << std::endl;
+
     grid->Broadcast(0,(void *)&csum,sizeof(csum));
 
     return csum;
diff --git a/lib/qcd/action/gauge/GaugeImpl.h b/lib/qcd/action/gauge/GaugeImpl.h
index 08fbfa68..3448062b 100644
--- a/lib/qcd/action/gauge/GaugeImpl.h
+++ b/lib/qcd/action/gauge/GaugeImpl.h
@@ -85,6 +85,7 @@ public:
     LinkField Pmu(P._grid);
     Pmu = zero;
     for (int mu = 0; mu < Nd; mu++) {
+      std::cout << "generating momenta " << mu << "\n"; 
       SU<Nrepresentation>::GaussianFundamentalLieAlgebraMatrix(pRNG, Pmu);
       PokeIndex<LorentzIndex>(P, Pmu, mu);
     }
diff --git a/tests/hmc/Test_hmc_EODWFRatio_Binary.cc b/tests/hmc/Test_hmc_EODWFRatio_Binary.cc
index a7fb4db1..0b912c66 100644
--- a/tests/hmc/Test_hmc_EODWFRatio_Binary.cc
+++ b/tests/hmc/Test_hmc_EODWFRatio_Binary.cc
@@ -161,4 +161,6 @@ int main(int argc, char **argv) {
   TheHMC.MDparameters.set(20, 1.0);// MDsteps, traj length
 
   TheHMC.BuildTheAction(argc, argv);
+
+  Grid_finalize();
 }
diff --git a/tests/hmc/Test_hmc_WilsonGauge_Binary.cc b/tests/hmc/Test_hmc_WilsonGauge_Binary.cc
index d6041679..e1f1e8df 100644
--- a/tests/hmc/Test_hmc_WilsonGauge_Binary.cc
+++ b/tests/hmc/Test_hmc_WilsonGauge_Binary.cc
@@ -133,4 +133,6 @@ int main(int argc, char **argv) {
   TheHMC.MDparameters.set(20, 1.0);// MDsteps, traj length
 
   TheHMC.BuildTheAction(argc, argv);
+
+  Grid_finalize();
 }