Introduce halfspace_matrix_pair(ZZ, ...)
#4709
Draft
halfspace_matrix_pair(ZZ, ...)
#4709