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::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();