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

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

    }
  }
  __syncthreads();

  return acc_i;
}
