 #include "lts_cuda_kernels.cuh"
 
 // for pops
-#define BSIZE_POPS 64
+//#define BSIZE_POPS 128
+#define BSIZE_POPS 1024
 
 // for neighborhood search
 #define BSIZE_NH 128

                                         , i_to);
    }
    cudaDeviceSynchronize();
    check_error();
    // get partial results from GPU
    std::vector<unsigned int> partial_pops(n_rows*n_radii);
    cudaMemcpy(partial_pops.data()


  ////

  void check_error(std::string msg="") {
    cudaError_t err = cudaGetLastError();
    if (err != cudaSuccess) {
      std::cerr << "CUDA error: "
                << msg << "\n"
                << cudaGetErrorString(err) << std::endl;
      exit(EXIT_FAILURE);
    }
  }
