portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-04 11:15:55 +01:00
Code Releases Activity

Small problem in compressor fix

Browse Source
This commit is contained in:
Peter Boyle 2016-02-17 19:03:09 -06:00
parent b6f6da923e
commit 22422a84d9
2 changed files with 9 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
lib/qcd/action/fermion/WilsonCompressor.h
View File

@ -212,14 +212,14 @@ namespace QCD {
this->u_comm_offset=0; this->u_comm_offset=0;
int dag = compress.dag; int dag = compress.dag;
static std::vector<int> dirs(8); static std::vector<int> dirs(Nd*2);
for(int mu=0;mu<4;mu++){ for(int mu=0;mu<Nd;mu++){
if ( dag ) { if ( dag ) {
dirs[mu] =mu; dirs[mu] =mu;
dirs[mu+4]=mu+4; dirs[mu+4]=mu+Nd;
} else { } else {
dirs[mu] =mu+4; dirs[mu] =mu+Nd;
dirs[mu+4]=mu; dirs[mu+Nd]=mu;
} }
} }

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

@ -304,8 +304,8 @@ void WilsonFermion5D<Impl>::DhopInternalCommsThenCompute(StencilImpl & st, Lebes
int nwork = U._grid->oSites(); int nwork = U._grid->oSites();
commtime -=usecond(); commtime -=usecond();
auto handle = st.HaloExchangeOptBegin(in,compressor); auto handle = st.HaloExchangeBegin(in,compressor);
st.HaloExchangeOptComplete(handle); st.HaloExchangeComplete(handle);
commtime +=usecond(); commtime +=usecond();
jointime -=usecond(); jointime -=usecond();
@ -440,7 +440,7 @@ void WilsonFermion5D<Impl>::DhopInternalCommsOverlapCompute(StencilImpl & st, Le
int nwork = U._grid->oSites(); int nwork = U._grid->oSites();
commtime -=usecond(); commtime -=usecond();
auto handle = st.HaloExchangeOptBegin(in,compressor); auto handle = st.HaloExchangeBegin(in,compressor);
commtime +=usecond(); commtime +=usecond();
// Dhop takes the 4d grid from U, and makes a 5d index for fermion // Dhop takes the 4d grid from U, and makes a 5d index for fermion
@ -498,7 +498,7 @@ PARALLEL_FOR_LOOP
dslashtime +=usecond(); dslashtime +=usecond();
jointime -=usecond(); jointime -=usecond();
st.HaloExchangeOptComplete(handle); st.HaloExchangeComplete(handle);
jointime +=usecond(); jointime +=usecond();
local = false; local = false;