function [nh_split, masks] = cosmo_neighborhood_split(nh, varargin)
cosmo nchoosek partitioner hdr
cosmo nfold partitioner hdr