* 91893 Orsay, France
*/
+#include <isl_ctx_private.h>
#include <isl_map_private.h>
#include <isl_bound.h>
#include <isl_bernstein.h>
#include <isl_range.h>
#include <isl_polynomial_private.h>
+#include <isl_options_private.h>
/* Compute a bound on the polynomial defined over the parametric polytope
* using either range propagation or bernstein expansion and
struct isl_bound *bound = (struct isl_bound *)user;
isl_pw_qpolynomial_fold *top_pwf;
isl_pw_qpolynomial_fold *top_pwf_tight;
- isl_dim *dim;
+ isl_space *dim;
isl_morph *morph;
- unsigned orig_nvar, final_nvar;
int r;
bset = isl_basic_set_detect_equalities(bset);
if (bset->n_eq == 0)
return compressed_guarded_poly_bound(bset, poly, user);
- orig_nvar = isl_basic_set_dim(bset, isl_dim_set);
-
morph = isl_basic_set_full_compression(bset);
bset = isl_morph_basic_set(isl_morph_copy(morph), bset);
- poly = isl_qpolynomial_morph(poly, isl_morph_copy(morph));
-
- final_nvar = isl_basic_set_dim(bset, isl_dim_set);
+ poly = isl_qpolynomial_morph_domain(poly, isl_morph_copy(morph));
- dim = isl_morph_get_ran_dim(morph);
- dim = isl_dim_drop(dim, isl_dim_set, 0, isl_dim_size(dim, isl_dim_set));
+ dim = isl_morph_get_ran_space(morph);
+ dim = isl_space_params(dim);
top_pwf = bound->pwf;
top_pwf_tight = bound->pwf_tight;
- bound->pwf = isl_pw_qpolynomial_fold_zero(isl_dim_copy(dim),
+ dim = isl_space_from_domain(dim);
+ dim = isl_space_add_dims(dim, isl_dim_out, 1);
+ bound->pwf = isl_pw_qpolynomial_fold_zero(isl_space_copy(dim),
bound->type);
bound->pwf_tight = isl_pw_qpolynomial_fold_zero(dim, bound->type);
r = compressed_guarded_poly_bound(bset, poly, user);
- morph = isl_morph_remove_dom_dims(morph, isl_dim_set, 0, orig_nvar);
- morph = isl_morph_remove_ran_dims(morph, isl_dim_set, 0, final_nvar);
+ morph = isl_morph_dom_params(morph);
+ morph = isl_morph_ran_params(morph);
morph = isl_morph_inverse(morph);
- bound->pwf = isl_pw_qpolynomial_fold_morph(bound->pwf,
+ bound->pwf = isl_pw_qpolynomial_fold_morph_domain(bound->pwf,
isl_morph_copy(morph));
- bound->pwf_tight = isl_pw_qpolynomial_fold_morph(bound->pwf_tight, morph);
+ bound->pwf_tight = isl_pw_qpolynomial_fold_morph_domain(
+ bound->pwf_tight, morph);
bound->pwf = isl_pw_qpolynomial_fold_fold(top_pwf, bound->pwf);
bound->pwf_tight = isl_pw_qpolynomial_fold_fold(top_pwf_tight,
__isl_take isl_qpolynomial *poly, void *user)
{
struct isl_bound *bound = (struct isl_bound *)user;
- isl_dim *dim;
+ isl_space *dim;
isl_pw_qpolynomial_fold *top_pwf;
isl_pw_qpolynomial_fold *top_pwf_tight;
int nparam;
if (!bound->wrapping)
return unwrapped_guarded_poly_bound(bset, poly, user);
- nparam = isl_dim_size(bound->dim, isl_dim_param);
- n_in = isl_dim_size(bound->dim, isl_dim_set);
+ nparam = isl_space_dim(bound->dim, isl_dim_param);
+ n_in = isl_space_dim(bound->dim, isl_dim_in);
bset = isl_basic_set_move_dims(bset, isl_dim_param, nparam,
isl_dim_set, 0, n_in);
poly = isl_qpolynomial_move_dims(poly, isl_dim_param, nparam,
- isl_dim_set, 0, n_in);
+ isl_dim_in, 0, n_in);
- dim = isl_basic_set_get_dim(bset);
- dim = isl_dim_drop(dim, isl_dim_set, 0, isl_dim_size(dim, isl_dim_set));
+ dim = isl_basic_set_get_space(bset);
+ dim = isl_space_params(dim);
top_pwf = bound->pwf;
top_pwf_tight = bound->pwf_tight;
- bound->pwf = isl_pw_qpolynomial_fold_zero(isl_dim_copy(dim),
+ dim = isl_space_from_domain(dim);
+ dim = isl_space_add_dims(dim, isl_dim_out, 1);
+ bound->pwf = isl_pw_qpolynomial_fold_zero(isl_space_copy(dim),
bound->type);
bound->pwf_tight = isl_pw_qpolynomial_fold_zero(dim, bound->type);
r = unwrapped_guarded_poly_bound(bset, poly, user);
- bound->pwf = isl_pw_qpolynomial_fold_reset_dim(bound->pwf,
- isl_dim_copy(bound->dim));
- bound->pwf_tight = isl_pw_qpolynomial_fold_reset_dim(bound->pwf_tight,
- isl_dim_copy(bound->dim));
+ bound->pwf = isl_pw_qpolynomial_fold_reset_space(bound->pwf,
+ isl_space_copy(bound->dim));
+ bound->pwf_tight = isl_pw_qpolynomial_fold_reset_space(bound->pwf_tight,
+ isl_space_copy(bound->dim));
bound->pwf = isl_pw_qpolynomial_fold_fold(top_pwf, bound->pwf);
bound->pwf_tight = isl_pw_qpolynomial_fold_fold(top_pwf_tight,
if (!pwf)
return NULL;
- bound.dim = isl_pw_qpolynomial_fold_get_dim(pwf);
- nvar = isl_dim_size(bound.dim, isl_dim_set);
+ bound.dim = isl_pw_qpolynomial_fold_get_domain_space(pwf);
- bound.wrapping = isl_dim_is_wrapping(bound.dim);
- if (bound.wrapping) {
- bound.dim = isl_dim_unwrap(bound.dim);
- nvar = isl_dim_size(bound.dim, isl_dim_out);
- bound.dim = isl_dim_domain(bound.dim);
- } else
- bound.dim = isl_dim_drop(bound.dim, isl_dim_set, 0, nvar);
+ bound.wrapping = isl_space_is_wrapping(bound.dim);
+ if (bound.wrapping)
+ bound.dim = isl_space_unwrap(bound.dim);
+ nvar = isl_space_dim(bound.dim, isl_dim_out);
+ bound.dim = isl_space_domain(bound.dim);
+ bound.dim = isl_space_from_domain(bound.dim);
+ bound.dim = isl_space_add_dims(bound.dim, isl_dim_out, 1);
if (nvar == 0) {
if (tight)
*tight = 1;
- return isl_pw_qpolynomial_fold_reset_dim(pwf, bound.dim);
+ return isl_pw_qpolynomial_fold_reset_space(pwf, bound.dim);
}
if (isl_pw_qpolynomial_fold_is_zero(pwf)) {
return isl_pw_qpolynomial_fold_zero(bound.dim, type);
}
- bound.pwf = isl_pw_qpolynomial_fold_zero(isl_dim_copy(bound.dim),
+ bound.pwf = isl_pw_qpolynomial_fold_zero(isl_space_copy(bound.dim),
pwf->type);
- bound.pwf_tight = isl_pw_qpolynomial_fold_zero(isl_dim_copy(bound.dim),
+ bound.pwf_tight = isl_pw_qpolynomial_fold_zero(isl_space_copy(bound.dim),
pwf->type);
bound.check_tight = !!tight;
if (tight)
*tight = covers;
- isl_dim_free(bound.dim);
+ isl_space_free(bound.dim);
isl_pw_qpolynomial_fold_free(pwf);
if (covers) {
isl_pw_qpolynomial_fold_free(bound.pwf_tight);
isl_pw_qpolynomial_fold_free(bound.pwf);
isl_pw_qpolynomial_fold_free(pwf);
- isl_dim_free(bound.dim);
+ isl_space_free(bound.dim);
return NULL;
}
__isl_take isl_union_pw_qpolynomial *upwqp,
enum isl_fold type, int *tight)
{
- isl_dim *dim;
+ isl_space *dim;
struct isl_union_bound_data data = { type, 1, NULL };
if (!upwqp)
if (!tight)
data.tight = 0;
- dim = isl_union_pw_qpolynomial_get_dim(upwqp);
+ dim = isl_union_pw_qpolynomial_get_space(upwqp);
data.res = isl_union_pw_qpolynomial_fold_zero(dim, type);
if (isl_union_pw_qpolynomial_foreach_pw_qpolynomial(upwqp,
&bound_pw, &data) < 0)