cosmo overlap hdr

function [y_in_x, x_in_y] = cosmo_overlap(xs, ys)