function [y_in_x, x_in_y] = cosmo_overlap(xs, ys)
cosmo oddeven partitioner hdr
cosmo parallel get nproc available hdr