From b92060f511942f0ad1331b0fd2d3ed1b766df854 Mon Sep 17 00:00:00 2001 From: Peter Boyle Date: Tue, 9 Jun 2015 22:41:59 +0100 Subject: [PATCH] 5d OpDir direction interface refers to the 5d dims, not 4d to present a sensible and consistent external interface. --- lib/qcd/action/fermion/WilsonFermion5D.cc | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/lib/qcd/action/fermion/WilsonFermion5D.cc b/lib/qcd/action/fermion/WilsonFermion5D.cc index 7d9ffe22..7bbeeef9 100644 --- a/lib/qcd/action/fermion/WilsonFermion5D.cc +++ b/lib/qcd/action/fermion/WilsonFermion5D.cc @@ -72,6 +72,7 @@ namespace QCD { } void WilsonFermion5D::DoubleStore(LatticeDoubledGaugeField &Uds,const LatticeGaugeField &Umu) { + assert(GaugeGrid()->_ndimension==4); conformable(Uds._grid,GaugeGrid()); conformable(Umu._grid,GaugeGrid()); LatticeColourMatrix U(GaugeGrid()); @@ -82,8 +83,10 @@ void WilsonFermion5D::DoubleStore(LatticeDoubledGaugeField &Uds,const LatticeGau pokeIndex(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( (dir>=0)&&(dir<4) ); //must do x,y,z or t;