Updated the comment to be precise on the case

This commit is contained in:
Trung Nguyen
2022-02-13 14:51:00 -06:00
parent 159f107abd
commit 193dea6327

View File

@ -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 this loop apparently iterates
// over neighbors of j that are not intended
// for neigh no within pair hybrid: still need the check below
// 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<cutsq_jk)) continue;
@ -788,8 +788,8 @@ __kernel void k_sw_three_end_vatom(const __global numtyp4 *restrict x_,
numtyp delr2z = kx.z - jx.z;
numtyp rsq2 = delr2x*delr2x + delr2y*delr2y + delr2z*delr2z;
// still need this check because this loop apparently iterates
// over neighbors of j that are not intended
// for neigh no within pair hybrid: still need the check below
// 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<cutsq_jk)) continue;