portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Dirichlet fix

Browse Source
This commit is contained in:
Peter Boyle 2022-06-15 00:05:20 -04:00
parent 501bb117bf
commit fdef7a1a8c
1 changed files with 9 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
Grid/qcd/action/fermion/implementation/WilsonFermion5DImplementation.h
View File

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