void
open_hypercube(int dimension, ppl_Polyhedron_t ph) {
int i;
mpz_t z_one;
mpz_t z_minus_one;
ppl_Coefficient_t coeff_one;
ppl_Coefficient_t coeff_minus_one;
ppl_Linear_Expression_t le;
ppl_Constraint_t c;
ppl_Constraint_System_t cs;
mpz_init_set_si(z_one, 1);
mpz_init_set_si(z_minus_one, -1);
ppl_new_Coefficient(&coeff_one);
ppl_assign_Coefficient_from_mpz_t(coeff_one, z_one);
ppl_new_Coefficient(&coeff_minus_one);
ppl_assign_Coefficient_from_mpz_t(coeff_minus_one, z_minus_one);
ppl_new_Linear_Expression_with_dimension(&le, dimension);
ppl_new_Constraint_System(&cs);
for (i = 0; i < dimension; ++i) {
ppl_Linear_Expression_add_to_coefficient(le, i, coeff_one);
/* Variable(i) > 0 */
ppl_new_Constraint(&c, le, PPL_CONSTRAINT_TYPE_GREATER_THAN);
ppl_Constraint_System_insert_Constraint(cs, c);
ppl_delete_Constraint(c);
/* Variable(i) < 1 */
ppl_Linear_Expression_add_to_inhomogeneous(le, coeff_minus_one);
ppl_new_Constraint(&c, le, PPL_CONSTRAINT_TYPE_LESS_THAN);
ppl_Constraint_System_insert_Constraint(cs, c);
ppl_delete_Constraint(c);
/* Zero `le' */
ppl_Linear_Expression_add_to_coefficient(le, i, coeff_minus_one);
ppl_Linear_Expression_add_to_inhomogeneous(le, coeff_one);
}
ppl_Polyhedron_add_constraints(ph, cs);
ppl_delete_Constraint_System(cs);
ppl_delete_Linear_Expression(le);
ppl_delete_Coefficient(coeff_minus_one);
ppl_delete_Coefficient(coeff_one);
mpz_clear(z_minus_one);
mpz_clear(z_one);
}
static void
new_Constraint_System_from_Cloog_Matrix (ppl_Constraint_System_t *pcs,
CloogMatrix *matrix)
{
int i;
ppl_new_Constraint_System (pcs);
for (i = 0; i < matrix->NbRows; i++)
{
ppl_Constraint_t c = cloog_matrix_to_ppl_constraint (matrix, i);
ppl_Constraint_System_insert_Constraint (*pcs, c);
ppl_delete_Constraint (c);
}
}
static void
add_constraints(ppl_Linear_Expression_t ppl_le,
int type, mpq_t rational_lb, mpq_t rational_ub, mpz_t den_lcm,
ppl_Constraint_System_t ppl_cs) {
ppl_Constraint_t ppl_c;
ppl_Linear_Expression_t ppl_le2;
switch (type) {
case GLP_FR:
break;
case GLP_LO:
mpz_mul(tmp_z, den_lcm, mpq_numref(rational_lb));
mpz_divexact(tmp_z, tmp_z, mpq_denref(rational_lb));
mpz_neg(tmp_z, tmp_z);
ppl_assign_Coefficient_from_mpz_t(ppl_coeff, tmp_z);
ppl_Linear_Expression_add_to_inhomogeneous(ppl_le, ppl_coeff);
ppl_new_Constraint(&ppl_c, ppl_le, PPL_CONSTRAINT_TYPE_GREATER_OR_EQUAL);
if (verbosity >= 4) {
ppl_io_fprint_Constraint(output_file, ppl_c);
fprintf(output_file, "\n");
}
ppl_Constraint_System_insert_Constraint(ppl_cs, ppl_c);
ppl_delete_Constraint(ppl_c);
break;
case GLP_UP:
mpz_mul(tmp_z, den_lcm, mpq_numref(rational_ub));
mpz_divexact(tmp_z, tmp_z, mpq_denref(rational_ub));
mpz_neg(tmp_z, tmp_z);
ppl_assign_Coefficient_from_mpz_t(ppl_coeff, tmp_z);
ppl_Linear_Expression_add_to_inhomogeneous(ppl_le, ppl_coeff);
ppl_new_Constraint(&ppl_c, ppl_le,
PPL_CONSTRAINT_TYPE_LESS_OR_EQUAL);
if (verbosity >= 4) {
ppl_io_fprint_Constraint(output_file, ppl_c);
fprintf(output_file, "\n");
}
ppl_Constraint_System_insert_Constraint(ppl_cs, ppl_c);
ppl_delete_Constraint(ppl_c);
break;
case GLP_DB:
ppl_new_Linear_Expression_from_Linear_Expression(&ppl_le2, ppl_le);
mpz_mul(tmp_z, den_lcm, mpq_numref(rational_lb));
mpz_divexact(tmp_z, tmp_z, mpq_denref(rational_lb));
mpz_neg(tmp_z, tmp_z);
ppl_assign_Coefficient_from_mpz_t(ppl_coeff, tmp_z);
ppl_Linear_Expression_add_to_inhomogeneous(ppl_le, ppl_coeff);
ppl_new_Constraint(&ppl_c, ppl_le, PPL_CONSTRAINT_TYPE_GREATER_OR_EQUAL);
if (verbosity >= 4) {
ppl_io_fprint_Constraint(output_file, ppl_c);
fprintf(output_file, "\n");
}
ppl_Constraint_System_insert_Constraint(ppl_cs, ppl_c);
ppl_delete_Constraint(ppl_c);
mpz_mul(tmp_z, den_lcm, mpq_numref(rational_ub));
mpz_divexact(tmp_z, tmp_z, mpq_denref(rational_ub));
mpz_neg(tmp_z, tmp_z);
ppl_assign_Coefficient_from_mpz_t(ppl_coeff, tmp_z);
ppl_Linear_Expression_add_to_inhomogeneous(ppl_le2, ppl_coeff);
ppl_new_Constraint(&ppl_c, ppl_le2, PPL_CONSTRAINT_TYPE_LESS_OR_EQUAL);
ppl_delete_Linear_Expression(ppl_le2);
if (verbosity >= 4) {
ppl_io_fprint_Constraint(output_file, ppl_c);
fprintf(output_file, "\n");
}
ppl_Constraint_System_insert_Constraint(ppl_cs, ppl_c);
ppl_delete_Constraint(ppl_c);
break;
case GLP_FX:
mpz_mul(tmp_z, den_lcm, mpq_numref(rational_lb));
mpz_divexact(tmp_z, tmp_z, mpq_denref(rational_lb));
mpz_neg(tmp_z, tmp_z);
ppl_assign_Coefficient_from_mpz_t(ppl_coeff, tmp_z);
ppl_Linear_Expression_add_to_inhomogeneous(ppl_le, ppl_coeff);
ppl_new_Constraint(&ppl_c, ppl_le,
PPL_CONSTRAINT_TYPE_EQUAL);
if (verbosity >= 4) {
ppl_io_fprint_Constraint(output_file, ppl_c);
fprintf(output_file, "\n");
}
ppl_Constraint_System_insert_Constraint(ppl_cs, ppl_c);
ppl_delete_Constraint(ppl_c);
break;
default:
fatal("internal error");
break;
}
}
