216 	isl_basic_set *bset;  in rational_universe()  local
218 	bset = isl_basic_set_universe(space);  in rational_universe()
219 	bset = isl_basic_set_set_rational(bset);  in rational_universe()
221 	return bset;  in rational_universe()
237 static __isl_give isl_basic_set *farkas(__isl_take isl_basic_set *bset,  in farkas()  argument
246 	total = isl_basic_set_dim(bset, isl_dim_all);  in farkas()
248 		return isl_basic_set_free(bset);  in farkas()
250 	ctx = isl_basic_set_get_ctx(bset);  in farkas()
252 	if (isl_basic_set_plain_is_empty(bset)) {  in farkas()
253 		isl_basic_set_free(bset);  in farkas()
257 	dual = isl_basic_set_alloc_space(space, bset->n_eq + bset->n_ineq,  in farkas()
258 					total, bset->n_ineq + (shift > 0));  in farkas()
261 	for (i = 0; i < bset->n_eq + bset->n_ineq; ++i) {  in farkas()
274 		for (j = 0; j < bset->n_eq; ++j)  in farkas()
276 				    bset->eq[j][1 + i]);  in farkas()
277 		for (j = 0; j < bset->n_ineq; ++j)  in farkas()
278 			isl_int_set(dual->eq[k][1 + shift + total + bset->n_eq + j],  in farkas()
279 				    bset->ineq[j][1 + i]);  in farkas()
282 	for (i = 0; i < bset->n_ineq; ++i) {  in farkas()
287 			    1 + shift + total + bset->n_eq + bset->n_ineq);  in farkas()
288 		isl_int_set_si(dual->ineq[k][1 + shift + total + bset->n_eq + i], 1);  in farkas()
297 		for (j = 0; j < bset->n_eq; ++j)  in farkas()
299 				    bset->eq[j][0]);  in farkas()
300 		for (j = 0; j < bset->n_ineq; ++j)  in farkas()
301 			isl_int_neg(dual->ineq[k][2 + total + bset->n_eq + j],  in farkas()
302 				    bset->ineq[j][0]);  in farkas()
309 	isl_basic_set_free(bset);  in farkas()
312 	isl_basic_set_free(bset);  in farkas()
322 	__isl_take isl_basic_set *bset)  in isl_basic_set_coefficients_base()  argument
324 	return farkas(bset, 1);  in isl_basic_set_coefficients_base()
414 static int is_ray(__isl_keep isl_basic_set *bset, int ineq)  in is_ray()  argument
416 	return isl_int_is_zero(bset->ineq[ineq][1]);  in is_ray()
426 	__isl_keep isl_basic_set *bset, void *user)  in isl_basic_set_coefficients_factor()  argument
433 	coeff = isl_basic_set_coefficients_base(isl_basic_set_copy(bset));  in isl_basic_set_coefficients_factor()
438 	dim = isl_basic_set_dim(bset, isl_dim_set);  in isl_basic_set_coefficients_factor()
491 static __isl_give isl_basic_set *add_lines(__isl_take isl_basic_set *bset,  in add_lines()  argument
499 		k = isl_basic_set_alloc_equality(bset);  in add_lines()
501 			return isl_basic_set_free(bset);  in add_lines()
502 		clear_entry(bset->eq[k], total);  in add_lines()
503 		set_factor(bset->eq[k], factor->coeff->eq[i], factor);  in add_lines()
506 	return bset;  in add_lines()
513 static __isl_give isl_basic_set *add_rays(__isl_take isl_basic_set *bset,  in add_rays()  argument
525 		k = isl_basic_set_alloc_inequality(bset);  in add_rays()
527 			return isl_basic_set_free(bset);  in add_rays()
528 		clear_entry(bset->ineq[k], total);  in add_rays()
529 		set_factor(bset->ineq[k], data->coeff->ineq[i], data);  in add_rays()
532 	return bset;  in add_rays()
602 static __isl_give isl_basic_set *add_vertex(__isl_take isl_basic_set *bset,  in add_vertex()  argument
609 	k = isl_basic_set_alloc_inequality(bset);  in add_vertex()
611 		return isl_basic_set_free(bset);  in add_vertex()
622 	isl_int_set_si(bset->ineq[k][0], 0);  in add_vertex()
623 	isl_int_set(bset->ineq[k][1], lcm);  in add_vertex()
630 		scale_factor(bset->ineq[k], coeff->ineq[pos], f, factor);  in add_vertex()
636 	return bset;  in add_vertex()
814 	__isl_take isl_basic_set *bset)  in basic_set_coefficients()  argument
819 	nparam = isl_basic_set_dim(bset, isl_dim_param);  in basic_set_coefficients()
821 		return isl_basic_set_free(bset);  in basic_set_coefficients()
822 	bset = isl_basic_set_move_dims(bset, isl_dim_set, 0,  in basic_set_coefficients()
825 	f = isl_basic_set_factorizer(bset);  in basic_set_coefficients()
827 		return isl_basic_set_free(bset);  in basic_set_coefficients()
829 		isl_basic_set_free(bset);  in basic_set_coefficients()
833 	return isl_basic_set_coefficients_base(bset);  in basic_set_coefficients()
840 	__isl_take isl_basic_set *bset)  in isl_basic_set_coefficients()  argument
844 	if (!bset)  in isl_basic_set_coefficients()
846 	if (bset->n_div)  in isl_basic_set_coefficients()
847 		isl_die(bset->ctx, isl_error_invalid,  in isl_basic_set_coefficients()
851 	space = isl_basic_set_get_space(bset);  in isl_basic_set_coefficients()
854 	bset = basic_set_coefficients(bset);  in isl_basic_set_coefficients()
855 	bset = isl_basic_set_reset_space(bset, space);  in isl_basic_set_coefficients()
856 	return bset;  in isl_basic_set_coefficients()
858 	isl_basic_set_free(bset);  in isl_basic_set_coefficients()
867 	__isl_take isl_basic_set *bset)  in isl_basic_set_solutions()  argument
871 	if (!bset)  in isl_basic_set_solutions()
873 	if (bset->n_div)  in isl_basic_set_solutions()
874 		isl_die(bset->ctx, isl_error_invalid,  in isl_basic_set_solutions()
878 	space = isl_basic_set_get_space(bset);  in isl_basic_set_solutions()
881 	bset = farkas(bset, -1);  in isl_basic_set_solutions()
882 	bset = isl_basic_set_reset_space(bset, space);  in isl_basic_set_solutions()
883 	return bset;  in isl_basic_set_solutions()
885 	isl_basic_set_free(bset);  in isl_basic_set_solutions()
909 		isl_basic_set *bset, *coeff_i;  in isl_set_coefficients()  local
910 		bset = isl_basic_set_copy(set->p[i]);  in isl_set_coefficients()
911 		coeff_i = isl_basic_set_coefficients(bset);  in isl_set_coefficients()
923 	__isl_take isl_basic_set *bset, void *user)  in coefficients_wrap()  argument
925 	return isl_basic_set_coefficients(bset);  in coefficients_wrap()
958 		isl_basic_set *bset, *sol_i;  in isl_set_solutions()  local
959 		bset = isl_basic_set_copy(set->p[i]);  in isl_set_solutions()
960 		sol_i = isl_basic_set_solutions(bset);  in isl_set_solutions()