More bug fixes to modified PPPM charge spread kernel.

2011-03-09 08:58:08 -05:00
parent 8e03c58917
commit 111de61716
2 changed files with 2 additions and 2 deletions
--- a/lib/gpu/pppm_gpu_memory.cpp
+++ b/lib/gpu/pppm_gpu_memory.cpp
@ -258,7 +258,7 @@ int PPPMGPUMemoryT::compute(const int ago, const int nlocal, const int nall,
  f_brick_z=boxlo[2]+(_nzlo_out-_nlower+shift)/delzinv;

  BX=block_size();
-  GX=_npts_y*_npts_z*8;
+  GX=static_cast<int>(ceil(static_cast<double>(_npts_y*_npts_z)/8));
  k_make_rho.set_size(GX,BX);
  k_make_rho.run(&atom->dev_x.begin(), &atom->dev_q.begin(),
                 &d_brick_counts.begin(), &d_brick_atoms.begin(),