2 * Copyright 2008-2009 Katholieke Universiteit Leuven
4 * Use of this software is governed by the GNU LGPLv2.1 license
6 * Written by Sven Verdoolaege, K.U.Leuven, Departement
7 * Computerwetenschappen, Celestijnenlaan 200A, B-3001 Leuven, Belgium
13 struct isl_dim
*isl_dim_alloc(struct isl_ctx
*ctx
,
14 unsigned nparam
, unsigned n_in
, unsigned n_out
)
18 dim
= isl_alloc_type(ctx
, struct isl_dim
);
35 struct isl_dim
*isl_dim_set_alloc(struct isl_ctx
*ctx
,
36 unsigned nparam
, unsigned dim
)
38 return isl_dim_alloc(ctx
, nparam
, 0, dim
);
41 static unsigned global_pos(struct isl_dim
*dim
,
42 enum isl_dim_type type
, unsigned pos
)
44 struct isl_ctx
*ctx
= dim
->ctx
;
48 isl_assert(ctx
, pos
< dim
->nparam
, return isl_dim_total(dim
));
51 isl_assert(ctx
, pos
< dim
->n_in
, return isl_dim_total(dim
));
52 return pos
+ dim
->nparam
;
54 isl_assert(ctx
, pos
< dim
->n_out
, return isl_dim_total(dim
));
55 return pos
+ dim
->nparam
+ dim
->n_in
;
57 isl_assert(ctx
, 0, return isl_dim_total(dim
));
59 return isl_dim_total(dim
);
62 /* Extend length of names array to the total number of dimensions.
64 static __isl_give isl_dim
*extend_names(__isl_take isl_dim
*dim
)
66 struct isl_name
**names
;
69 if (isl_dim_total(dim
) <= dim
->n_name
)
73 dim
->names
= isl_calloc_array(dim
->ctx
,
74 struct isl_name
*, isl_dim_total(dim
));
78 names
= isl_realloc_array(dim
->ctx
, dim
->names
,
79 struct isl_name
*, isl_dim_total(dim
));
83 for (i
= dim
->n_name
; i
< isl_dim_total(dim
); ++i
)
87 dim
->n_name
= isl_dim_total(dim
);
95 static struct isl_dim
*set_name(struct isl_dim
*dim
,
96 enum isl_dim_type type
, unsigned pos
,
97 struct isl_name
*name
)
99 struct isl_ctx
*ctx
= dim
->ctx
;
100 dim
= isl_dim_cow(dim
);
105 pos
= global_pos(dim
, type
, pos
);
106 isl_assert(ctx
, pos
!= isl_dim_total(dim
), goto error
);
108 if (pos
>= dim
->n_name
) {
111 dim
= extend_names(dim
);
116 dim
->names
[pos
] = name
;
120 isl_name_free(ctx
, name
);
125 static struct isl_name
*get_name(struct isl_dim
*dim
,
126 enum isl_dim_type type
, unsigned pos
)
131 pos
= global_pos(dim
, type
, pos
);
132 if (pos
== isl_dim_total(dim
))
134 if (pos
>= dim
->n_name
)
136 return dim
->names
[pos
];
139 static unsigned offset(struct isl_dim
*dim
, enum isl_dim_type type
)
142 case isl_dim_param
: return 0;
143 case isl_dim_in
: return dim
->nparam
;
144 case isl_dim_out
: return dim
->nparam
+ dim
->n_in
;
148 static unsigned n(struct isl_dim
*dim
, enum isl_dim_type type
)
151 case isl_dim_param
: return dim
->nparam
;
152 case isl_dim_in
: return dim
->n_in
;
153 case isl_dim_out
: return dim
->n_out
;
157 unsigned isl_dim_size(struct isl_dim
*dim
, enum isl_dim_type type
)
164 static struct isl_dim
*copy_names(struct isl_dim
*dst
,
165 enum isl_dim_type dst_type
, unsigned offset
, struct isl_dim
*src
,
166 enum isl_dim_type src_type
)
169 struct isl_name
*name
;
171 for (i
= 0; i
< n(src
, src_type
); ++i
) {
172 name
= get_name(src
, src_type
, i
);
175 dst
= set_name(dst
, dst_type
, offset
+ i
,
176 isl_name_copy(dst
->ctx
, name
));
183 struct isl_dim
*isl_dim_dup(struct isl_dim
*dim
)
186 dup
= isl_dim_alloc(dim
->ctx
, dim
->nparam
, dim
->n_in
, dim
->n_out
);
189 dup
= copy_names(dup
, isl_dim_param
, 0, dim
, isl_dim_param
);
190 dup
= copy_names(dup
, isl_dim_in
, 0, dim
, isl_dim_in
);
191 dup
= copy_names(dup
, isl_dim_out
, 0, dim
, isl_dim_out
);
195 struct isl_dim
*isl_dim_cow(struct isl_dim
*dim
)
203 return isl_dim_dup(dim
);
206 struct isl_dim
*isl_dim_copy(struct isl_dim
*dim
)
215 void isl_dim_free(struct isl_dim
*dim
)
225 for (i
= 0; i
< dim
->n_name
; ++i
)
226 isl_name_free(dim
->ctx
, dim
->names
[i
]);
228 isl_ctx_deref(dim
->ctx
);
233 struct isl_dim
*isl_dim_set_name(struct isl_dim
*dim
,
234 enum isl_dim_type type
, unsigned pos
,
237 struct isl_name
*name
;
240 name
= isl_name_get(dim
->ctx
, s
);
243 return set_name(dim
, type
, pos
, name
);
249 const char *isl_dim_get_name(struct isl_dim
*dim
,
250 enum isl_dim_type type
, unsigned pos
)
252 struct isl_name
*name
= get_name(dim
, type
, pos
);
253 return name
? name
->name
: NULL
;
256 static int match(struct isl_dim
*dim1
, enum isl_dim_type dim1_type
,
257 struct isl_dim
*dim2
, enum isl_dim_type dim2_type
)
261 if (n(dim1
, dim1_type
) != n(dim2
, dim2_type
))
264 if (!dim1
->names
&& !dim2
->names
)
267 for (i
= 0; i
< n(dim1
, dim1_type
); ++i
) {
268 if (get_name(dim1
, dim1_type
, i
) !=
269 get_name(dim2
, dim2_type
, i
))
275 int isl_dim_match(struct isl_dim
*dim1
, enum isl_dim_type dim1_type
,
276 struct isl_dim
*dim2
, enum isl_dim_type dim2_type
)
278 return match(dim1
, dim1_type
, dim2
, dim2_type
);
281 static void get_names(struct isl_dim
*dim
, enum isl_dim_type type
,
282 unsigned first
, unsigned n
, struct isl_name
**names
)
286 for (i
= 0; i
< n
; ++i
)
287 names
[i
] = get_name(dim
, type
, first
+i
);
290 struct isl_dim
*isl_dim_extend(struct isl_dim
*dim
,
291 unsigned nparam
, unsigned n_in
, unsigned n_out
)
293 struct isl_name
**names
= NULL
;
297 if (dim
->nparam
== nparam
&& dim
->n_in
== n_in
&& dim
->n_out
== n_out
)
300 isl_assert(dim
->ctx
, dim
->nparam
<= nparam
, goto error
);
301 isl_assert(dim
->ctx
, dim
->n_in
<= n_in
, goto error
);
302 isl_assert(dim
->ctx
, dim
->n_out
<= n_out
, goto error
);
304 dim
= isl_dim_cow(dim
);
307 names
= isl_calloc_array(dim
->ctx
, struct isl_name
*,
308 nparam
+ n_in
+ n_out
);
311 get_names(dim
, isl_dim_param
, 0, dim
->nparam
, names
);
312 get_names(dim
, isl_dim_in
, 0, dim
->n_in
, names
+ nparam
);
313 get_names(dim
, isl_dim_out
, 0, dim
->n_out
,
314 names
+ nparam
+ n_in
);
317 dim
->n_name
= nparam
+ n_in
+ n_out
;
319 dim
->nparam
= nparam
;
330 struct isl_dim
*isl_dim_add(struct isl_dim
*dim
, enum isl_dim_type type
,
335 return isl_dim_extend(dim
,
336 dim
->nparam
+ n
, dim
->n_in
, dim
->n_out
);
338 return isl_dim_extend(dim
,
339 dim
->nparam
, dim
->n_in
+ n
, dim
->n_out
);
341 return isl_dim_extend(dim
,
342 dim
->nparam
, dim
->n_in
, dim
->n_out
+ n
);
347 __isl_give isl_dim
*isl_dim_insert(__isl_take isl_dim
*dim
,
348 enum isl_dim_type type
, unsigned pos
, unsigned n
)
350 struct isl_name
**names
= NULL
;
357 isl_assert(dim
->ctx
, pos
<= isl_dim_size(dim
, type
), goto error
);
359 dim
= isl_dim_cow(dim
);
367 names
= isl_calloc_array(dim
->ctx
, struct isl_name
*,
368 dim
->nparam
+ dim
->n_in
+ dim
->n_out
+ n
);
372 size
[isl_dim_param
] = dim
->nparam
;
373 size
[isl_dim_in
] = dim
->n_in
;
374 size
[isl_dim_out
] = dim
->n_out
;
375 for (t
= isl_dim_param
; t
<= isl_dim_out
; ++t
) {
377 get_names(dim
, t
, 0, size
[t
], names
+ off
);
380 get_names(dim
, t
, 0, pos
, names
+ off
);
382 get_names(dim
, t
, pos
, size
[t
]-pos
, names
+off
);
383 off
+= size
[t
] - pos
;
388 dim
->n_name
= dim
->nparam
+ dim
->n_in
+ dim
->n_out
+ n
;
391 case isl_dim_param
: dim
->nparam
+= n
; break;
392 case isl_dim_in
: dim
->n_in
+= n
; break;
393 case isl_dim_out
: dim
->n_out
+= n
; break;
402 __isl_give isl_dim
*isl_dim_move(__isl_take isl_dim
*dim
,
403 enum isl_dim_type dst_type
, unsigned dst_pos
,
404 enum isl_dim_type src_type
, unsigned src_pos
, unsigned n
)
411 isl_assert(dim
->ctx
, src_pos
+ n
<= isl_dim_size(dim
, src_type
),
414 if (dst_type
== src_type
&& dst_pos
== src_pos
)
417 isl_assert(dim
->ctx
, dst_type
!= src_type
, goto error
);
419 dim
= isl_dim_cow(dim
);
424 struct isl_name
**names
;
428 names
= isl_calloc_array(dim
->ctx
, struct isl_name
*,
429 dim
->nparam
+ dim
->n_in
+ dim
->n_out
);
433 size
[isl_dim_param
] = dim
->nparam
;
434 size
[isl_dim_in
] = dim
->n_in
;
435 size
[isl_dim_out
] = dim
->n_out
;
436 for (t
= isl_dim_param
; t
<= isl_dim_out
; ++t
) {
438 get_names(dim
, t
, 0, dst_pos
, names
+ off
);
440 get_names(dim
, src_type
, src_pos
, n
, names
+off
);
442 get_names(dim
, t
, dst_pos
, size
[t
] - dst_pos
,
444 off
+= size
[t
] - dst_pos
;
445 } else if (t
== src_type
) {
446 get_names(dim
, t
, 0, src_pos
, names
+ off
);
448 get_names(dim
, t
, src_pos
+ n
,
449 size
[t
] - src_pos
- n
, names
+ off
);
450 off
+= size
[t
] - src_pos
- n
;
452 get_names(dim
, t
, 0, size
[t
], names
+ off
);
458 dim
->n_name
= dim
->nparam
+ dim
->n_in
+ dim
->n_out
;
462 case isl_dim_param
: dim
->nparam
+= n
; break;
463 case isl_dim_in
: dim
->n_in
+= n
; break;
464 case isl_dim_out
: dim
->n_out
+= n
; break;
468 case isl_dim_param
: dim
->nparam
-= n
; break;
469 case isl_dim_in
: dim
->n_in
-= n
; break;
470 case isl_dim_out
: dim
->n_out
-= n
; break;
479 struct isl_dim
*isl_dim_join(struct isl_dim
*left
, struct isl_dim
*right
)
486 isl_assert(left
->ctx
, match(left
, isl_dim_param
, right
, isl_dim_param
),
488 isl_assert(left
->ctx
, n(left
, isl_dim_out
) == n(right
, isl_dim_in
),
491 dim
= isl_dim_alloc(left
->ctx
, left
->nparam
, left
->n_in
, right
->n_out
);
495 dim
= copy_names(dim
, isl_dim_param
, 0, left
, isl_dim_param
);
496 dim
= copy_names(dim
, isl_dim_in
, 0, left
, isl_dim_in
);
497 dim
= copy_names(dim
, isl_dim_out
, 0, right
, isl_dim_out
);
509 struct isl_dim
*isl_dim_product(struct isl_dim
*left
, struct isl_dim
*right
)
516 isl_assert(left
->ctx
, match(left
, isl_dim_param
, right
, isl_dim_param
),
519 dim
= isl_dim_alloc(left
->ctx
, left
->nparam
,
520 left
->n_in
+ right
->n_in
, left
->n_out
+ right
->n_out
);
524 dim
= copy_names(dim
, isl_dim_param
, 0, left
, isl_dim_param
);
525 dim
= copy_names(dim
, isl_dim_in
, 0, left
, isl_dim_in
);
526 dim
= copy_names(dim
, isl_dim_in
, left
->n_in
, right
, isl_dim_in
);
527 dim
= copy_names(dim
, isl_dim_out
, 0, left
, isl_dim_out
);
528 dim
= copy_names(dim
, isl_dim_out
, left
->n_out
, right
, isl_dim_out
);
540 struct isl_dim
*isl_dim_map(struct isl_dim
*dim
)
542 struct isl_name
**names
= NULL
;
546 isl_assert(dim
->ctx
, dim
->n_in
== 0, goto error
);
549 dim
= isl_dim_cow(dim
);
553 names
= isl_calloc_array(dim
->ctx
, struct isl_name
*,
554 dim
->nparam
+ dim
->n_out
+ dim
->n_out
);
557 get_names(dim
, isl_dim_param
, 0, dim
->nparam
, names
);
558 get_names(dim
, isl_dim_out
, 0, dim
->n_out
, names
+ dim
->nparam
);
560 dim
->n_in
= dim
->n_out
;
564 dim
->n_name
= dim
->nparam
+ dim
->n_out
+ dim
->n_out
;
565 dim
= copy_names(dim
, isl_dim_out
, 0, dim
, isl_dim_in
);
573 static struct isl_dim
*set_names(struct isl_dim
*dim
, enum isl_dim_type type
,
574 unsigned first
, unsigned n
, struct isl_name
**names
)
578 for (i
= 0; i
< n
; ++i
)
579 dim
= set_name(dim
, type
, first
+i
, names
[i
]);
584 struct isl_dim
*isl_dim_reverse(struct isl_dim
*dim
)
587 struct isl_name
**names
= NULL
;
591 if (match(dim
, isl_dim_in
, dim
, isl_dim_out
))
594 dim
= isl_dim_cow(dim
);
599 names
= isl_alloc_array(dim
->ctx
, struct isl_name
*,
600 dim
->n_in
+ dim
->n_out
);
603 get_names(dim
, isl_dim_in
, 0, dim
->n_in
, names
);
604 get_names(dim
, isl_dim_out
, 0, dim
->n_out
, names
+ dim
->n_in
);
608 dim
->n_in
= dim
->n_out
;
612 dim
= set_names(dim
, isl_dim_out
, 0, dim
->n_out
, names
);
613 dim
= set_names(dim
, isl_dim_in
, 0, dim
->n_in
, names
+ dim
->n_out
);
624 struct isl_dim
*isl_dim_drop(struct isl_dim
*dim
, enum isl_dim_type type
,
625 unsigned first
, unsigned num
)
635 isl_assert(dim
->ctx
, first
+ num
<= n(dim
, type
), goto error
);
636 dim
= isl_dim_cow(dim
);
640 dim
= extend_names(dim
);
643 for (i
= 0; i
< num
; ++i
)
644 isl_name_free(dim
->ctx
, get_name(dim
, type
, first
+i
));
645 for (i
= first
+num
; i
< n(dim
, type
); ++i
)
646 set_name(dim
, type
, i
- num
, get_name(dim
, type
, i
));
649 get_names(dim
, isl_dim_in
, 0, dim
->n_in
,
650 dim
->names
+ offset(dim
, isl_dim_in
) - num
);
652 get_names(dim
, isl_dim_out
, 0, dim
->n_out
,
653 dim
->names
+ offset(dim
, isl_dim_out
) - num
);
660 case isl_dim_param
: dim
->nparam
-= num
; break;
661 case isl_dim_in
: dim
->n_in
-= num
; break;
662 case isl_dim_out
: dim
->n_out
-= num
; break;
670 struct isl_dim
*isl_dim_drop_inputs(struct isl_dim
*dim
,
671 unsigned first
, unsigned n
)
673 return isl_dim_drop(dim
, isl_dim_in
, first
, n
);
676 struct isl_dim
*isl_dim_drop_outputs(struct isl_dim
*dim
,
677 unsigned first
, unsigned n
)
679 return isl_dim_drop(dim
, isl_dim_out
, first
, n
);
682 struct isl_dim
*isl_dim_domain(struct isl_dim
*dim
)
686 dim
= isl_dim_drop_outputs(dim
, 0, dim
->n_out
);
687 return isl_dim_reverse(dim
);
690 struct isl_dim
*isl_dim_range(struct isl_dim
*dim
)
694 return isl_dim_drop_inputs(dim
, 0, dim
->n_in
);
697 struct isl_dim
*isl_dim_underlying(struct isl_dim
*dim
, unsigned n_div
)
704 dim
->nparam
== 0 && dim
->n_in
== 0 && dim
->n_name
== 0)
706 dim
= isl_dim_cow(dim
);
709 dim
->n_out
+= dim
->nparam
+ dim
->n_in
+ n_div
;
713 for (i
= 0; i
< dim
->n_name
; ++i
)
714 isl_name_free(dim
->ctx
, get_name(dim
, isl_dim_out
, i
));
720 unsigned isl_dim_total(struct isl_dim
*dim
)
722 return dim
->nparam
+ dim
->n_in
+ dim
->n_out
;
725 int isl_dim_equal(struct isl_dim
*dim1
, struct isl_dim
*dim2
)
727 return match(dim1
, isl_dim_param
, dim2
, isl_dim_param
) &&
728 n(dim1
, isl_dim_in
) == n(dim2
, isl_dim_in
) &&
729 n(dim1
, isl_dim_out
) == n(dim2
, isl_dim_out
);
732 int isl_dim_compatible(struct isl_dim
*dim1
, struct isl_dim
*dim2
)
734 return dim1
->nparam
== dim2
->nparam
&&
735 dim1
->n_in
+ dim1
->n_out
== dim2
->n_in
+ dim2
->n_out
;