diff options
Diffstat (limited to 'cloog-0.16.3/isl/isl_point.c')
| -rw-r--r-- | cloog-0.16.3/isl/isl_point.c | 488 |
1 files changed, 0 insertions, 488 deletions
diff --git a/cloog-0.16.3/isl/isl_point.c b/cloog-0.16.3/isl/isl_point.c deleted file mode 100644 index e840162..0000000 --- a/cloog-0.16.3/isl/isl_point.c +++ /dev/null @@ -1,488 +0,0 @@ -#include <isl_map_private.h> -#include <isl_point_private.h> -#include <isl/set.h> -#include <isl_sample.h> -#include <isl_scan.h> -#include <isl/seq.h> -#include <isl_dim_private.h> - -isl_ctx *isl_point_get_ctx(__isl_keep isl_point *pnt) -{ - return pnt ? isl_dim_get_ctx(pnt->dim) : NULL; -} - -__isl_give isl_dim *isl_point_get_dim(__isl_keep isl_point *pnt) -{ - return pnt ? isl_dim_copy(pnt->dim) : NULL; -} - -__isl_give isl_point *isl_point_alloc(__isl_take isl_dim *dim, - __isl_take isl_vec *vec) -{ - struct isl_point *pnt; - - if (!dim || !vec) - goto error; - - if (vec->size > 1 + isl_dim_total(dim)) { - vec = isl_vec_cow(vec); - if (!vec) - goto error; - vec->size = 1 + isl_dim_total(dim); - } - - pnt = isl_alloc_type(dim->ctx, struct isl_point); - if (!pnt) - goto error; - - pnt->ref = 1; - pnt->dim = dim; - pnt->vec = vec; - - return pnt; -error: - isl_dim_free(dim); - isl_vec_free(vec); - return NULL; -} - -__isl_give isl_point *isl_point_zero(__isl_take isl_dim *dim) -{ - isl_vec *vec; - - if (!dim) - return NULL; - vec = isl_vec_alloc(dim->ctx, 1 + isl_dim_total(dim)); - if (!vec) - goto error; - isl_int_set_si(vec->el[0], 1); - isl_seq_clr(vec->el + 1, vec->size - 1); - return isl_point_alloc(dim, vec); -error: - isl_dim_free(dim); - return NULL; -} - -__isl_give isl_point *isl_point_dup(__isl_keep isl_point *pnt) -{ - struct isl_point *pnt2; - - if (!pnt) - return NULL; - pnt2 = isl_point_alloc(isl_dim_copy(pnt->dim), isl_vec_copy(pnt->vec)); - return pnt2; -} - -__isl_give isl_point *isl_point_cow(__isl_take isl_point *pnt) -{ - struct isl_point *pnt2; - if (!pnt) - return NULL; - - if (pnt->ref == 1) - return pnt; - - pnt2 = isl_point_dup(pnt); - isl_point_free(pnt); - return pnt2; -} - -__isl_give isl_point *isl_point_copy(__isl_keep isl_point *pnt) -{ - if (!pnt) - return NULL; - - pnt->ref++; - return pnt; -} - -void isl_point_free(__isl_take isl_point *pnt) -{ - if (!pnt) - return; - - if (--pnt->ref > 0) - return; - - isl_dim_free(pnt->dim); - isl_vec_free(pnt->vec); - free(pnt); -} - -__isl_give isl_point *isl_point_void(__isl_take isl_dim *dim) -{ - if (!dim) - return NULL; - - return isl_point_alloc(dim, isl_vec_alloc(dim->ctx, 0)); -} - -int isl_point_is_void(__isl_keep isl_point *pnt) -{ - if (!pnt) - return -1; - - return pnt->vec->size == 0; -} - -void isl_point_get_coordinate(__isl_keep isl_point *pnt, - enum isl_dim_type type, int pos, isl_int *v) -{ - if (!pnt || isl_point_is_void(pnt)) - return; - if (type == isl_dim_set) - pos += isl_dim_size(pnt->dim, isl_dim_param); - isl_int_set(*v, pnt->vec->el[1 + pos]); -} - -__isl_give isl_point *isl_point_set_coordinate(__isl_take isl_point *pnt, - enum isl_dim_type type, int pos, isl_int v) -{ - if (!pnt || isl_point_is_void(pnt)) - return pnt; - - pnt = isl_point_cow(pnt); - if (!pnt) - return NULL; - pnt->vec = isl_vec_cow(pnt->vec); - if (!pnt->vec) - goto error; - - if (type == isl_dim_set) - pos += isl_dim_size(pnt->dim, isl_dim_param); - - isl_int_set(pnt->vec->el[1 + pos], v); - - return pnt; -error: - isl_point_free(pnt); - return NULL; -} - -__isl_give isl_point *isl_point_add_ui(__isl_take isl_point *pnt, - enum isl_dim_type type, int pos, unsigned val) -{ - if (!pnt || isl_point_is_void(pnt)) - return pnt; - - pnt = isl_point_cow(pnt); - if (!pnt) - return NULL; - pnt->vec = isl_vec_cow(pnt->vec); - if (!pnt->vec) - goto error; - - if (type == isl_dim_set) - pos += isl_dim_size(pnt->dim, isl_dim_param); - - isl_int_add_ui(pnt->vec->el[1 + pos], pnt->vec->el[1 + pos], val); - - return pnt; -error: - isl_point_free(pnt); - return NULL; -} - -__isl_give isl_point *isl_point_sub_ui(__isl_take isl_point *pnt, - enum isl_dim_type type, int pos, unsigned val) -{ - if (!pnt || isl_point_is_void(pnt)) - return pnt; - - pnt = isl_point_cow(pnt); - if (!pnt) - return NULL; - pnt->vec = isl_vec_cow(pnt->vec); - if (!pnt->vec) - goto error; - - if (type == isl_dim_set) - pos += isl_dim_size(pnt->dim, isl_dim_param); - - isl_int_sub_ui(pnt->vec->el[1 + pos], pnt->vec->el[1 + pos], val); - - return pnt; -error: - isl_point_free(pnt); - return NULL; -} - -struct isl_foreach_point { - struct isl_scan_callback callback; - int (*fn)(__isl_take isl_point *pnt, void *user); - void *user; - isl_dim *dim; -}; - -static int foreach_point(struct isl_scan_callback *cb, __isl_take isl_vec *sample) -{ - struct isl_foreach_point *fp = (struct isl_foreach_point *)cb; - isl_point *pnt; - - pnt = isl_point_alloc(isl_dim_copy(fp->dim), sample); - - return fp->fn(pnt, fp->user); -} - -int isl_set_foreach_point(__isl_keep isl_set *set, - int (*fn)(__isl_take isl_point *pnt, void *user), void *user) -{ - struct isl_foreach_point fp = { { &foreach_point }, fn, user }; - int i; - - if (!set) - return -1; - - fp.dim = isl_set_get_dim(set); - if (!fp.dim) - return -1; - - set = isl_set_copy(set); - set = isl_set_cow(set); - set = isl_set_make_disjoint(set); - set = isl_set_compute_divs(set); - if (!set) - goto error; - - for (i = 0; i < set->n; ++i) - if (isl_basic_set_scan(isl_basic_set_copy(set->p[i]), - &fp.callback) < 0) - goto error; - - isl_set_free(set); - isl_dim_free(fp.dim); - - return 0; -error: - isl_set_free(set); - isl_dim_free(fp.dim); - return -1; -} - -/* Return 1 if "bmap" contains the point "point". - * "bmap" is assumed to have known divs. - * The point is first extended with the divs and then passed - * to basic_map_contains. - */ -int isl_basic_map_contains_point(__isl_keep isl_basic_map *bmap, - __isl_keep isl_point *point) -{ - int i; - struct isl_vec *vec; - unsigned dim; - int contains; - - if (!bmap || !point) - return -1; - isl_assert(bmap->ctx, isl_dim_equal(bmap->dim, point->dim), return -1); - if (bmap->n_div == 0) - return isl_basic_map_contains(bmap, point->vec); - - dim = isl_basic_map_total_dim(bmap) - bmap->n_div; - vec = isl_vec_alloc(bmap->ctx, 1 + dim + bmap->n_div); - if (!vec) - return -1; - - isl_seq_cpy(vec->el, point->vec->el, point->vec->size); - for (i = 0; i < bmap->n_div; ++i) { - isl_seq_inner_product(bmap->div[i] + 1, vec->el, - 1 + dim + i, &vec->el[1+dim+i]); - isl_int_fdiv_q(vec->el[1+dim+i], vec->el[1+dim+i], - bmap->div[i][0]); - } - - contains = isl_basic_map_contains(bmap, vec); - - isl_vec_free(vec); - return contains; -} - -int isl_map_contains_point(__isl_keep isl_map *map, __isl_keep isl_point *point) -{ - int i; - int found = 0; - - if (!map || !point) - return -1; - - map = isl_map_copy(map); - map = isl_map_compute_divs(map); - if (!map) - return -1; - - for (i = 0; i < map->n; ++i) { - found = isl_basic_map_contains_point(map->p[i], point); - if (found < 0) - goto error; - if (found) - break; - } - isl_map_free(map); - - return found; -error: - isl_map_free(map); - return -1; -} - -int isl_set_contains_point(__isl_keep isl_set *set, __isl_keep isl_point *point) -{ - return isl_map_contains_point((isl_map *)set, point); -} - -__isl_give isl_basic_set *isl_basic_set_from_point(__isl_take isl_point *pnt) -{ - isl_basic_set *bset; - isl_basic_set *model; - - model = isl_basic_set_empty(isl_dim_copy(pnt->dim)); - bset = isl_basic_set_from_vec(isl_vec_copy(pnt->vec)); - bset = isl_basic_set_from_underlying_set(bset, model); - isl_point_free(pnt); - - return bset; -} - -__isl_give isl_set *isl_set_from_point(__isl_take isl_point *pnt) -{ - isl_basic_set *bset; - bset = isl_basic_set_from_point(pnt); - return isl_set_from_basic_set(bset); -} - -__isl_give isl_basic_set *isl_basic_set_box_from_points( - __isl_take isl_point *pnt1, __isl_take isl_point *pnt2) -{ - isl_basic_set *bset; - unsigned total; - int i; - int k; - isl_int t; - - isl_int_init(t); - - if (!pnt1 || !pnt2) - goto error; - - isl_assert(pnt1->dim->ctx, - isl_dim_equal(pnt1->dim, pnt2->dim), goto error); - - if (isl_point_is_void(pnt1) && isl_point_is_void(pnt2)) { - isl_dim *dim = isl_dim_copy(pnt1->dim); - isl_point_free(pnt1); - isl_point_free(pnt2); - isl_int_clear(t); - return isl_basic_set_empty(dim); - } - if (isl_point_is_void(pnt1)) { - isl_point_free(pnt1); - isl_int_clear(t); - return isl_basic_set_from_point(pnt2); - } - if (isl_point_is_void(pnt2)) { - isl_point_free(pnt2); - isl_int_clear(t); - return isl_basic_set_from_point(pnt1); - } - - total = isl_dim_total(pnt1->dim); - bset = isl_basic_set_alloc_dim(isl_dim_copy(pnt1->dim), 0, 0, 2 * total); - - for (i = 0; i < total; ++i) { - isl_int_mul(t, pnt1->vec->el[1 + i], pnt2->vec->el[0]); - isl_int_submul(t, pnt2->vec->el[1 + i], pnt1->vec->el[0]); - - k = isl_basic_set_alloc_inequality(bset); - if (k < 0) - goto error; - isl_seq_clr(bset->ineq[k] + 1, total); - if (isl_int_is_pos(t)) { - isl_int_set_si(bset->ineq[k][1 + i], -1); - isl_int_set(bset->ineq[k][0], pnt1->vec->el[1 + i]); - } else { - isl_int_set_si(bset->ineq[k][1 + i], 1); - isl_int_neg(bset->ineq[k][0], pnt1->vec->el[1 + i]); - } - isl_int_fdiv_q(bset->ineq[k][0], bset->ineq[k][0], pnt1->vec->el[0]); - - k = isl_basic_set_alloc_inequality(bset); - if (k < 0) - goto error; - isl_seq_clr(bset->ineq[k] + 1, total); - if (isl_int_is_pos(t)) { - isl_int_set_si(bset->ineq[k][1 + i], 1); - isl_int_neg(bset->ineq[k][0], pnt2->vec->el[1 + i]); - } else { - isl_int_set_si(bset->ineq[k][1 + i], -1); - isl_int_set(bset->ineq[k][0], pnt2->vec->el[1 + i]); - } - isl_int_fdiv_q(bset->ineq[k][0], bset->ineq[k][0], pnt2->vec->el[0]); - } - - bset = isl_basic_set_finalize(bset); - - isl_point_free(pnt1); - isl_point_free(pnt2); - - isl_int_clear(t); - - return bset; -error: - isl_point_free(pnt1); - isl_point_free(pnt2); - isl_int_clear(t); - return NULL; -} - -__isl_give isl_set *isl_set_box_from_points(__isl_take isl_point *pnt1, - __isl_take isl_point *pnt2) -{ - isl_basic_set *bset; - bset = isl_basic_set_box_from_points(pnt1, pnt2); - return isl_set_from_basic_set(bset); -} - -void isl_point_print(__isl_keep isl_point *pnt, FILE *out) -{ - int i; - unsigned nparam; - unsigned dim; - - if (!pnt) - return; - if (isl_point_is_void(pnt)) { - fprintf(out, "void\n"); - return; - } - - nparam = isl_dim_size(pnt->dim, isl_dim_param); - dim = isl_dim_size(pnt->dim, isl_dim_set); - if (nparam > 0) { - fprintf(out, "["); - for (i = 0; i < nparam; ++i) { - const char *name; - if (i) - fprintf(out, ", "); - name = isl_dim_get_name(pnt->dim, isl_dim_param, i); - if (name) - fprintf(out, "%s = ", name); - isl_int_print(out, pnt->vec->el[1 + i], 0); - if (!isl_int_is_one(pnt->vec->el[0])) { - fprintf(out, "/"); - isl_int_print(out, pnt->vec->el[0], 0); - } - } - fprintf(out, "] -> "); - } - fprintf(out, "["); - for (i = 0; i < dim; ++i) { - if (i) - fprintf(out, ", "); - isl_int_print(out, pnt->vec->el[1 + nparam + i], 0); - if (!isl_int_is_one(pnt->vec->el[0])) { - fprintf(out, "/"); - isl_int_print(out, pnt->vec->el[0], 0); - } - } - fprintf(out, "]\n"); -} |