     
     // normally, data in HOOMD is not perfectly sorted:
     //for (unsigned int i = 0; i < 100; i++)
-    //  tweak_data();
+      //tweak_data();
     
     // prepare the binned data
     rebin_particles_host(gh_idxlist_coord, gh_idxlist_coord_trans, gh_bin_size, gh_pos, g_N, g_Lx, g_Ly, g_Lz, g_Mx, g_My, g_Mz, g_Nmax);

    g_Mx = int((g_Lx) / (g_rcut));
    g_My = int((g_Ly) / (g_rcut));
    g_Mz = int((g_Lz) / (g_rcut));
    g_Mx = std::min(g_Mx, (unsigned)30);
    g_My = std::min(g_My, (unsigned)30);
    g_Mz = std::min(g_Mz, (unsigned)30);

    // allocate bins
    unsigned int Nbins = g_Mx * g_My * g_Mz;

float g_rcut;
const unsigned int g_Nmax = 64;    // Maximum number of particles each cell can hold
const float tweak_dist = 0.1f;
const unsigned int g_neigh_max = 192;

//*************** data structures
float4 *gh_pos, *gd_pos;            // particle positions
