From fdef7a1a8cedad899a073d4ad46c58c7514d069d Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Wed, 15 Jun 2022 00:05:20 -0400
Subject: [PATCH] Dirichlet fix

---
 .../implementation/WilsonFermion5DImplementation.h       | 9 +++++++++
 1 file changed, 9 insertions(+)

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