isl_map_transitive_closure: intersect with domain and range before projection