Propagation of large uncertainty sets in orbital dynamics by automatic domain splitting