 #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

      exit(EXIT_FAILURE);
    }
    unsigned int block_rng = min_multiplicator(i_to-i_from, block_size);
    for (unsigned int i=0; i*block_size < n_rows; ++i) {
      population_count <<< block_rng
                         , block_size

                                                       , i_to);
    }
    cudaDeviceSynchronize();
    check_error("after kernel loop");
    // initialize neighborhoods
    Neighborhood nh;
    Neighborhood nhhd;
