unsigned int excl_bit;
real dx, dy, dz;
- get_nbl_exclusions_1(nbl, cj4_ind, j>>2, &excl);
+ const int jHalf = j/(c_nbnxnGpuClusterSize/c_nbnxnGpuClusterpairSplit);
+ get_nbl_exclusions_1(nbl, cj4_ind, jHalf, &excl);
excl_pair = a_mod_wj(j)*nbl->na_ci + i;
excl_bit = (1U << (gcj*c_gpuNumClusterPerCell + c));