function [cell_indices, unique_values] = cosmo_index_unique(values)
cosmo independent samples partitioner hdr
cosmo interval neighborhood hdr