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

5d OpDir direction interface refers to the 5d dims, not 4d to present a

Browse Source 
sensible and consistent external interface.
This commit is contained in:
Peter Boyle 2015-06-09 22:41:59 +01:00
parent 3a52bc4ce1
commit 4963f7356a
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
lib/qcd/action/fermion/WilsonFermion5D.cc
View File

@ -72,6 +72,7 @@ namespace QCD {
} }
void WilsonFermion5D::DoubleStore(LatticeDoubledGaugeField &Uds,const LatticeGaugeField &Umu) void WilsonFermion5D::DoubleStore(LatticeDoubledGaugeField &Uds,const LatticeGaugeField &Umu)
{ {
assert(GaugeGrid()->_ndimension==4);
conformable(Uds._grid,GaugeGrid()); conformable(Uds._grid,GaugeGrid());
conformable(Umu._grid,GaugeGrid()); conformable(Umu._grid,GaugeGrid());
LatticeColourMatrix U(GaugeGrid()); LatticeColourMatrix U(GaugeGrid());
@ -82,8 +83,10 @@ void WilsonFermion5D::DoubleStore(LatticeDoubledGaugeField &Uds,const LatticeGau
pokeIndex<LorentzIndex>(Uds,U,mu+4); pokeIndex<LorentzIndex>(Uds,U,mu+4);
} }
} }
void WilsonFermion5D::DhopDir(const LatticeFermion &in, LatticeFermion &out,int dir,int disp) void WilsonFermion5D::DhopDir(const LatticeFermion &in, LatticeFermion &out,int dir5,int disp)
{ {
int dir = dir5-1; // Maps to the ordering above in "directions" that is passed to stencil
// we drop off the innermost fifth dimension
assert( (disp==1)||(disp==-1) ); assert( (disp==1)||(disp==-1) );
assert( (dir>=0)&&(dir<4) ); //must do x,y,z or t; assert( (dir>=0)&&(dir<4) ); //must do x,y,z or t;