diff --git a/lib/gpu/lal_sw.cu b/lib/gpu/lal_sw.cu index 5c13caaa43..0e4c236f1a 100644 --- a/lib/gpu/lal_sw.cu +++ b/lib/gpu/lal_sw.cu @@ -666,8 +666,8 @@ __kernel void k_sw_three_end(const __global numtyp4 *restrict x_, numtyp delr2z = kx.z - jx.z; numtyp rsq2 = delr2x*delr2x + delr2y*delr2y + delr2z*delr2z; - // still need this check because j and k can be ghost atoms - // and short neighbor lists are built for local atoms only + // still need this check because this loop apparently iterates + // over neighbors of j that are not intended const numtyp cutsq_jk=sw_cut_ik*sw_cut_ik; if (!(rsq2>(numtyp)0 && rsq2(numtyp)0 && rsq2