     }
     return clustering;
   }
-
+*/
+
+
+
+/* TODO: remove
   std::vector<unsigned int>
   lumped_clusters__old(std::vector<unsigned int> clustering
                 , std::vector<unsigned int> ref_states

                         , 0);
    check_error("getting max shared mem size");
    unsigned int shared_mem = 2 * BSIZE_SCR * n_cols * sizeof(float);
    unsigned int block_rng, i_from, i_to, i, i_gpu;
    // lump microstates until nothing changes
    bool microstates_lumped = true;

    }
    return clustering;
  }
*/

  std::vector<std::size_t>
  initial_density_clustering(const std::vector<float>& free_energy
