diff --git a/lib/FFT.h b/lib/FFT.h index b5b31d82..240f338b 100644 --- a/lib/FFT.h +++ b/lib/FFT.h @@ -244,7 +244,10 @@ namespace Grid { pokeLocalSite(s,pgbuf,cbuf); } } - result = Cshift(result,dim,L); + if (p != processors[dim] - 1) + { + result = Cshift(result,dim,L); + } } // Loop over orthog coords @@ -287,10 +290,10 @@ namespace Grid { cgbuf = clbuf; cgbuf[dim] = clbuf[dim]+L*pc; peekLocalSite(s,pgbuf,cgbuf); - s = s * div; pokeLocalSite(s,result,clbuf); } } + result = result*div; // destroying plan FFTW::fftw_destroy_plan(p);