           if((n_stack1 - c_stack0) >= (LMEM_STACK_SIZE << SHIFT))
           {
             //We overwrote our current stack
-#if 0
             apprCount = -1; 
-#endif
             return acc_i;	 
           }
         }

#pragma unroll
          for (int i = 0; i < DIMx; i++)
            acc_i = add_acc(acc_i, pos_i, node_mon0[offs + i], node_mon1[offs+i], eps2);
#if 0
          apprCount += DIMx;
#endif
          __syncthreads();
        }
        __syncthreads();



  //Sum the interaction counters
  float  *sh_ds2 = (float*)&sh_acc[DIM];
  int    *sh_ngb = (int*  )&sh_ds2[DIM];
  sh_ds2[tid] = direCount;
  sh_ngb[tid] = apprCount;

