reduce_domain: compute convex union of domain