4 struct isl_dim *isl_dim_alloc(struct isl_ctx *ctx,
5 unsigned nparam, unsigned n_in, unsigned n_out)
9 dim = isl_alloc_type(ctx, struct isl_dim);
26 struct isl_dim *isl_dim_set_alloc(struct isl_ctx *ctx,
27 unsigned nparam, unsigned dim)
29 return isl_dim_alloc(ctx, nparam, 0, dim);
32 static unsigned global_pos(struct isl_dim *dim,
33 enum isl_dim_type type, unsigned pos)
35 struct isl_ctx *ctx = dim->ctx;
39 isl_assert(ctx, pos < dim->nparam, return isl_dim_total(dim));
42 isl_assert(ctx, pos < dim->n_in, return isl_dim_total(dim));
43 return pos + dim->nparam;
45 isl_assert(ctx, pos < dim->n_out, return isl_dim_total(dim));
46 return pos + dim->nparam + dim->n_in;
48 isl_assert(ctx, 0, goto error);
50 return isl_dim_total(dim);
53 static struct isl_dim *set_name(struct isl_dim *dim,
54 enum isl_dim_type type, unsigned pos,
55 struct isl_name *name)
57 struct isl_ctx *ctx = dim->ctx;
58 dim = isl_dim_cow(dim);
63 pos = global_pos(dim, type, pos);
64 isl_assert(ctx, pos != isl_dim_total(dim), goto error);
66 if (pos >= dim->n_name) {
70 dim->names = isl_calloc_array(dim->ctx,
71 struct isl_name *, isl_dim_total(dim));
76 dim->names = isl_realloc_array(dim->ctx, dim->names,
77 struct isl_name *, isl_dim_total(dim));
80 for (i = dim->n_name; i < isl_dim_total(dim); ++i)
83 dim->n_name = isl_dim_total(dim);
86 dim->names[pos] = name;
90 isl_name_free(ctx, name);
95 static struct isl_name *get_name(struct isl_dim *dim,
96 enum isl_dim_type type, unsigned pos)
101 pos = global_pos(dim, type, pos);
102 if (pos == isl_dim_total(dim))
104 if (pos >= dim->n_name)
106 return dim->names[pos];
109 static unsigned n(struct isl_dim *dim, enum isl_dim_type type)
112 case isl_dim_param: return dim->nparam;
113 case isl_dim_in: return dim->n_in;
114 case isl_dim_out: return dim->n_out;
118 unsigned isl_dim_size(struct isl_dim *dim, enum isl_dim_type type)
123 static struct isl_dim *copy_names(struct isl_dim *dst,
124 enum isl_dim_type dst_type, struct isl_dim *src,
125 enum isl_dim_type src_type)
128 struct isl_name *name;
130 for (i = 0; i < n(dst, dst_type); ++i) {
131 name = get_name(src, src_type, i);
134 dst = set_name(dst, dst_type, i, isl_name_copy(dst->ctx, name));
141 struct isl_dim *isl_dim_dup(struct isl_dim *dim)
144 dup = isl_dim_alloc(dim->ctx, dim->nparam, dim->n_in, dim->n_out);
147 dup = copy_names(dup, isl_dim_param, dim, isl_dim_param);
148 dup = copy_names(dup, isl_dim_in, dim, isl_dim_in);
149 dup = copy_names(dup, isl_dim_out, dim, isl_dim_out);
153 struct isl_dim *isl_dim_cow(struct isl_dim *dim)
161 return isl_dim_dup(dim);
164 struct isl_dim *isl_dim_copy(struct isl_dim *dim)
173 void isl_dim_free(struct isl_dim *dim)
183 for (i = 0; i < dim->n_name; ++i)
184 isl_name_free(dim->ctx, dim->names[i]);
186 isl_ctx_deref(dim->ctx);
191 struct isl_dim *isl_dim_set_name(struct isl_dim *dim,
192 enum isl_dim_type type, unsigned pos,
195 struct isl_name *name;
198 name = isl_name_get(dim->ctx, s);
201 return set_name(dim, type, pos, name);
207 const char *isl_dim_get_name(struct isl_dim *dim,
208 enum isl_dim_type type, unsigned pos)
210 struct isl_name *name = get_name(dim, type, pos);
211 return name ? name->name : NULL;
214 static int match(struct isl_dim *dim1, enum isl_dim_type dim1_type,
215 struct isl_dim *dim2, enum isl_dim_type dim2_type)
219 if (n(dim1, dim1_type) != n(dim2, dim2_type))
222 if (!dim1->names && !dim2->names)
225 for (i = 0; i < n(dim1, dim1_type); ++i) {
226 if (get_name(dim1, dim1_type, i) !=
227 get_name(dim2, dim2_type, i))
233 int isl_dim_match(struct isl_dim *dim1, enum isl_dim_type dim1_type,
234 struct isl_dim *dim2, enum isl_dim_type dim2_type)
236 return match(dim1, dim1_type, dim2, dim2_type);
239 static void get_names(struct isl_dim *dim, enum isl_dim_type type,
240 unsigned first, unsigned n, struct isl_name **names)
244 for (i = 0; i < n ; ++i)
245 names[i] = get_name(dim, type, first+i);
248 struct isl_dim *isl_dim_extend(struct isl_dim *dim,
249 unsigned nparam, unsigned n_in, unsigned n_out)
251 struct isl_name **names = NULL;
255 if (dim->nparam == nparam && dim->n_in == n_in && dim->n_out == n_out)
258 isl_assert(dim->ctx, dim->nparam <= nparam, goto error);
259 isl_assert(dim->ctx, dim->n_in <= n_in, goto error);
260 isl_assert(dim->ctx, dim->n_out <= n_out, goto error);
262 dim = isl_dim_cow(dim);
265 names = isl_calloc_array(dim->ctx, struct isl_name *,
266 nparam + n_in + n_out);
269 get_names(dim, isl_dim_param, 0, dim->nparam, names);
270 get_names(dim, isl_dim_in, 0, dim->n_in, names + nparam);
271 get_names(dim, isl_dim_out, 0, dim->n_out,
272 names + nparam + n_in);
275 dim->n_name = nparam + n_in + n_out;
277 dim->nparam = nparam;
288 struct isl_dim *isl_dim_join(struct isl_dim *left, struct isl_dim *right)
295 isl_assert(left->ctx, match(left, isl_dim_param, right, isl_dim_param),
297 isl_assert(left->ctx, match(left, isl_dim_out, right, isl_dim_in),
300 dim = isl_dim_alloc(left->ctx, left->nparam, left->n_in, right->n_out);
304 dim = copy_names(dim, isl_dim_param, left, isl_dim_param);
305 dim = copy_names(dim, isl_dim_in, left, isl_dim_in);
306 dim = copy_names(dim, isl_dim_out, right, isl_dim_out);
318 struct isl_dim *isl_dim_map(struct isl_dim *dim)
320 struct isl_name **names = NULL;
324 isl_assert(dim->ctx, dim->n_in == 0, goto error);
327 dim = isl_dim_cow(dim);
331 names = isl_calloc_array(dim->ctx, struct isl_name *,
332 dim->nparam + dim->n_out + dim->n_out);
335 get_names(dim, isl_dim_param, 0, dim->nparam, names);
336 get_names(dim, isl_dim_out, 0, dim->n_out, names + dim->nparam);
338 dim->n_in = dim->n_out;
340 copy_names(dim, isl_dim_out, dim, isl_dim_in);
343 dim->n_name = dim->nparam + dim->n_out + dim->n_out;
351 static struct isl_dim *set_names(struct isl_dim *dim, enum isl_dim_type type,
352 unsigned first, unsigned n, struct isl_name **names)
356 for (i = 0; i < n ; ++i)
357 dim = set_name(dim, type, first+i, names[i]);
362 struct isl_dim *isl_dim_reverse(struct isl_dim *dim)
365 struct isl_name **names = NULL;
369 if (match(dim, isl_dim_in, dim, isl_dim_out))
372 dim = isl_dim_cow(dim);
377 names = isl_alloc_array(dim->ctx, struct isl_name *,
378 dim->n_in + dim->n_out);
381 get_names(dim, isl_dim_in, 0, dim->n_in, names);
382 get_names(dim, isl_dim_out, 0, dim->n_out, names + dim->n_in);
386 dim->n_in = dim->n_out;
390 dim = set_names(dim, isl_dim_out, 0, dim->n_out, names);
391 dim = set_names(dim, isl_dim_in, 0, dim->n_in, names + dim->n_out);
402 struct isl_dim *isl_dim_drop_inputs(struct isl_dim *dim,
403 unsigned first, unsigned n)
413 isl_assert(dim->ctx, first + n <= dim->n_in, goto error);
414 dim = isl_dim_cow(dim);
418 for (i = 0; i < n; ++i) {
419 isl_name_free(dim->ctx,
420 get_name(dim, isl_dim_in, first+i));
422 for (i = first+n; i < dim->n_in; ++i)
423 set_name(dim, isl_dim_in, i - n,
424 get_name(dim, isl_dim_in, i));
425 get_names(dim, isl_dim_out, 0, dim->n_out,
426 dim->names + dim->nparam + dim->n_in - n);
435 struct isl_dim *isl_dim_drop_outputs(struct isl_dim *dim,
436 unsigned first, unsigned n)
446 isl_assert(dim->ctx, first + n <= dim->n_out, goto error);
447 dim = isl_dim_cow(dim);
451 for (i = 0; i < n; ++i) {
452 isl_name_free(dim->ctx,
453 get_name(dim, isl_dim_out, first+i));
455 for (i = first+n; i < dim->n_out; ++i)
456 set_name(dim, isl_dim_out, i - n,
457 get_name(dim, isl_dim_out, i));
466 struct isl_dim *isl_dim_domain(struct isl_dim *dim)
470 dim = isl_dim_drop_outputs(dim, 0, dim->n_out);
471 return isl_dim_reverse(dim);
474 struct isl_dim *isl_dim_underlying(struct isl_dim *dim, unsigned n_div)
481 dim->nparam == 0 && dim->n_in == 0 && dim->n_name == 0)
483 dim = isl_dim_cow(dim);
486 dim->n_out += dim->nparam + dim->n_in + n_div;
490 for (i = 0; i < dim->n_name; ++i)
491 isl_name_free(dim->ctx, get_name(dim, isl_dim_out, i));
497 unsigned isl_dim_total(struct isl_dim *dim)
499 return dim->nparam + dim->n_in + dim->n_out;
502 int isl_dim_equal(struct isl_dim *dim1, struct isl_dim *dim2)
504 return match(dim1, isl_dim_param, dim2, isl_dim_param) &&
505 match(dim1, isl_dim_in, dim2, isl_dim_in) &&
506 match(dim1, isl_dim_out, dim2, isl_dim_out);
509 int isl_dim_compatible(struct isl_dim *dim1, struct isl_dim *dim2)
511 return dim1->nparam == dim2->nparam &&
512 dim1->n_in + dim1->n_out == dim2->n_in + dim2->n_out;