From 137f1902580cede3251b590df62668456432f4cb Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Thu, 30 Jun 2022 13:45:07 -0400
Subject: [PATCH] Dirichlet implementation

---
 .../implementation/WilsonFermion5DImplementation.h     | 10 +++++++---
 1 file changed, 7 insertions(+), 3 deletions(-)

diff --git a/Grid/qcd/action/fermion/implementation/WilsonFermion5DImplementation.h b/Grid/qcd/action/fermion/implementation/WilsonFermion5DImplementation.h
index 681a6914..51c7df57 100644
--- a/Grid/qcd/action/fermion/implementation/WilsonFermion5DImplementation.h
+++ b/Grid/qcd/action/fermion/implementation/WilsonFermion5DImplementation.h
@@ -92,15 +92,19 @@ WilsonFermion5D<Impl>::WilsonFermion5D(GaugeField &_Umu,
     assert(FourDimRedBlackGrid._simd_layout[d]  ==FourDimGrid._simd_layout[d]);
   }
 
-  if ( p.dirichlet.size() ) {
+  if ( p.dirichlet.size() == Nd+1) {
     Coordinate block = p.dirichlet;
-    assert(block.size()==Nd+1);
     if ( block[0] || block[1] || block[2] || block[3] || block[4] ){
       Dirichlet = 1;
       Block = block;
     }
+  } else {
+    Coordinate block(Nd+1,0);
+    Block = block;
   }
-  
+
+  ZeroCounters();
+
   if (Impl::LsVectorised) { 
 
     int nsimd = Simd::Nsimd();