?  search:
summary | log | graphiclog1 | graphiclog2 | commit | commitdiff | tree | refs | edit | fork
blame | history | raw | HEAD
drop isl_map_drop_redundant_divs
[isl.git] / isl_multi_templ.c
blob166a3390a5e8341751103d3682f2782f6784acb8
1 /*
2 * Copyright 2011 Sven Verdoolaege
3 * Copyright 2012-2013 Ecole Normale Superieure
4 *
5 * Use of this software is governed by the MIT license
6 *
7 * Written by Sven Verdoolaege,
8 * Ecole Normale Superieure, 45 rue d’Ulm, 75230 Paris, France
9 */
10
11 #include <isl_space_private.h>
12 #include <isl/set.h>
13 #include <isl_reordering.h>
14
15 #include <isl_multi_macro.h>
16
17 #define MULTI_NAME(BASE) "isl_multi_" #BASE
18 #define xLIST(EL) EL ## _list
19 #define LIST(EL) xLIST(EL)
20
21 isl_ctx *FN(MULTI(BASE),get_ctx)(__isl_keep MULTI(BASE) *multi)
22 {
23 return multi ? isl_space_get_ctx(multi->space) : NULL;
24 }
25
26 __isl_give isl_space *FN(MULTI(BASE),get_space)(__isl_keep MULTI(BASE) *multi)
27 {
28 return multi ? isl_space_copy(multi->space) : NULL;
29 }
30
31 /* Return the position of the dimension of the given type and name
32 * in "multi".
33 * Return -1 if no such dimension can be found.
34 */
35 int FN(MULTI(BASE),find_dim_by_name)(__isl_keep MULTI(BASE) *multi,
36 enum isl_dim_type type, const char *name)
37 {
38 if (!multi)
39 return -1;
40 return isl_space_find_dim_by_name(multi->space, type, name);
41 }
42
43 __isl_give isl_space *FN(MULTI(BASE),get_domain_space)(
44 __isl_keep MULTI(BASE) *multi)
45 {
46 return multi ? isl_space_domain(isl_space_copy(multi->space)) : NULL;
47 }
48
49 __isl_give MULTI(BASE) *FN(MULTI(BASE),alloc)(__isl_take isl_space *space)
50 {
51 isl_ctx *ctx;
52 int n;
53 MULTI(BASE) *multi;
54
55 if (!space)
56 return NULL;
57
58 ctx = isl_space_get_ctx(space);
59 n = isl_space_dim(space, isl_dim_out);
60 multi = isl_calloc(ctx, MULTI(BASE),
61 sizeof(MULTI(BASE)) + (n - 1) * sizeof(struct EL *));
62 if (!multi)
63 goto error;
64
65 multi->space = space;
66 multi->n = n;
67 multi->ref = 1;
68 return multi;
69 error:
70 isl_space_free(space);
71 return NULL;
72 }
73
74 __isl_give MULTI(BASE) *FN(MULTI(BASE),dup)(__isl_keep MULTI(BASE) *multi)
75 {
76 int i;
77 MULTI(BASE) *dup;
78
79 if (!multi)
80 return NULL;
81
82 dup = FN(MULTI(BASE),alloc)(isl_space_copy(multi->space));
83 if (!dup)
84 return NULL;
85
86 for (i = 0; i < multi->n; ++i)
87 dup = FN(FN(MULTI(BASE),set),BASE)(dup, i,
88 FN(EL,copy)(multi->p[i]));
89
90 return dup;
91 }
92
93 __isl_give MULTI(BASE) *FN(MULTI(BASE),cow)(__isl_take MULTI(BASE) *multi)
94 {
95 if (!multi)
96 return NULL;
97
98 if (multi->ref == 1)
99 return multi;
100
101 multi->ref--;
102 return FN(MULTI(BASE),dup)(multi);
103 }
104
105 __isl_give MULTI(BASE) *FN(MULTI(BASE),copy)(__isl_keep MULTI(BASE) *multi)
106 {
107 if (!multi)
108 return NULL;
109
110 multi->ref++;
111 return multi;
112 }
113
114 __isl_null MULTI(BASE) *FN(MULTI(BASE),free)(__isl_take MULTI(BASE) *multi)
115 {
116 int i;
117
118 if (!multi)
119 return NULL;
120
121 if (--multi->ref > 0)
122 return NULL;
123
124 isl_space_free(multi->space);
125 for (i = 0; i < multi->n; ++i)
126 FN(EL,free)(multi->p[i]);
127 free(multi);
128
129 return NULL;
130 }
131
132 #ifndef NO_DIMS
133 /* Check whether "multi" has non-zero coefficients for any dimension
134 * in the given range or if any of these dimensions appear
135 * with non-zero coefficients in any of the integer divisions involved.
136 */
137 isl_bool FN(MULTI(BASE),involves_dims)(__isl_keep MULTI(BASE) *multi,
138 enum isl_dim_type type, unsigned first, unsigned n)
139 {
140 int i;
141
142 if (!multi)
143 return isl_bool_error;
144 if (multi->n == 0 || n == 0)
145 return isl_bool_false;
146
147 for (i = 0; i < multi->n; ++i) {
148 isl_bool involves;
149
150 involves = FN(EL,involves_dims)(multi->p[i], type, first, n);
151 if (involves < 0 || involves)
152 return involves;
153 }
154
155 return isl_bool_false;
156 }
157
158 __isl_give MULTI(BASE) *FN(MULTI(BASE),insert_dims)(
159 __isl_take MULTI(BASE) *multi,
160 enum isl_dim_type type, unsigned first, unsigned n)
161 {
162 int i;
163
164 if (!multi)
165 return NULL;
166 if (type == isl_dim_out)
167 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
168 "cannot insert output/set dimensions",
169 return FN(MULTI(BASE),free)(multi));
170 if (n == 0 && !isl_space_is_named_or_nested(multi->space, type))
171 return multi;
172
173 multi = FN(MULTI(BASE),cow)(multi);
174 if (!multi)
175 return NULL;
176
177 multi->space = isl_space_insert_dims(multi->space, type, first, n);
178 if (!multi->space)
179 return FN(MULTI(BASE),free)(multi);
180
181 for (i = 0; i < multi->n; ++i) {
182 multi->p[i] = FN(EL,insert_dims)(multi->p[i], type, first, n);
183 if (!multi->p[i])
184 return FN(MULTI(BASE),free)(multi);
185 }
186
187 return multi;
188 }
189
190 __isl_give MULTI(BASE) *FN(MULTI(BASE),add_dims)(__isl_take MULTI(BASE) *multi,
191 enum isl_dim_type type, unsigned n)
192 {
193 unsigned pos;
194
195 pos = FN(MULTI(BASE),dim)(multi, type);
196
197 return FN(MULTI(BASE),insert_dims)(multi, type, pos, n);
198 }
199 #endif
200
201 unsigned FN(MULTI(BASE),dim)(__isl_keep MULTI(BASE) *multi,
202 enum isl_dim_type type)
203 {
204 return multi ? isl_space_dim(multi->space, type) : 0;
205 }
206
207 /* Return the position of the first dimension of "type" with id "id".
208 * Return -1 if there is no such dimension.
209 */
210 int FN(MULTI(BASE),find_dim_by_id)(__isl_keep MULTI(BASE) *multi,
211 enum isl_dim_type type, __isl_keep isl_id *id)
212 {
213 if (!multi)
214 return -1;
215 return isl_space_find_dim_by_id(multi->space, type, id);
216 }
217
218 /* Return the id of the given dimension.
219 */
220 __isl_give isl_id *FN(MULTI(BASE),get_dim_id)(__isl_keep MULTI(BASE) *multi,
221 enum isl_dim_type type, unsigned pos)
222 {
223 return multi ? isl_space_get_dim_id(multi->space, type, pos) : NULL;
224 }
225
226 __isl_give MULTI(BASE) *FN(MULTI(BASE),set_dim_name)(
227 __isl_take MULTI(BASE) *multi,
228 enum isl_dim_type type, unsigned pos, const char *s)
229 {
230 int i;
231
232 multi = FN(MULTI(BASE),cow)(multi);
233 if (!multi)
234 return NULL;
235
236 multi->space = isl_space_set_dim_name(multi->space, type, pos, s);
237 if (!multi->space)
238 return FN(MULTI(BASE),free)(multi);
239
240 if (type == isl_dim_out)
241 return multi;
242 for (i = 0; i < multi->n; ++i) {
243 multi->p[i] = FN(EL,set_dim_name)(multi->p[i], type, pos, s);
244 if (!multi->p[i])
245 return FN(MULTI(BASE),free)(multi);
246 }
247
248 return multi;
249 }
250
251 const char *FN(MULTI(BASE),get_tuple_name)(__isl_keep MULTI(BASE) *multi,
252 enum isl_dim_type type)
253 {
254 return multi ? isl_space_get_tuple_name(multi->space, type) : NULL;
255 }
256
257 /* Does the specified tuple have an id?
258 */
259 isl_bool FN(MULTI(BASE),has_tuple_id)(__isl_keep MULTI(BASE) *multi,
260 enum isl_dim_type type)
261 {
262 if (!multi)
263 return isl_bool_error;
264 return isl_space_has_tuple_id(multi->space, type);
265 }
266
267 /* Return the id of the specified tuple.
268 */
269 __isl_give isl_id *FN(MULTI(BASE),get_tuple_id)(__isl_keep MULTI(BASE) *multi,
270 enum isl_dim_type type)
271 {
272 return multi ? isl_space_get_tuple_id(multi->space, type) : NULL;
273 }
274
275 __isl_give EL *FN(FN(MULTI(BASE),get),BASE)(__isl_keep MULTI(BASE) *multi,
276 int pos)
277 {
278 isl_ctx *ctx;
279
280 if (!multi)
281 return NULL;
282 ctx = FN(MULTI(BASE),get_ctx)(multi);
283 if (pos < 0 || pos >= multi->n)
284 isl_die(ctx, isl_error_invalid,
285 "index out of bounds", return NULL);
286 return FN(EL,copy)(multi->p[pos]);
287 }
288
289 __isl_give MULTI(BASE) *FN(FN(MULTI(BASE),set),BASE)(
290 __isl_take MULTI(BASE) *multi, int pos, __isl_take EL *el)
291 {
292 isl_space *multi_space = NULL;
293 isl_space *el_space = NULL;
294 isl_bool match;
295
296 multi = FN(MULTI(BASE),cow)(multi);
297 if (!multi || !el)
298 goto error;
299
300 multi_space = FN(MULTI(BASE),get_space)(multi);
301 match = FN(EL,matching_params)(el, multi_space);
302 if (match < 0)
303 goto error;
304 if (!match) {
305 multi = FN(MULTI(BASE),align_params)(multi,
306 FN(EL,get_space)(el));
307 isl_space_free(multi_space);
308 multi_space = FN(MULTI(BASE),get_space)(multi);
309 el = FN(EL,align_params)(el, isl_space_copy(multi_space));
310 }
311 if (FN(EL,check_match_domain_space)(el, multi_space) < 0)
312 goto error;
313
314 if (pos < 0 || pos >= multi->n)
315 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
316 "index out of bounds", goto error);
317
318 FN(EL,free)(multi->p[pos]);
319 multi->p[pos] = el;
320
321 isl_space_free(multi_space);
322 isl_space_free(el_space);
323
324 return multi;
325 error:
326 FN(MULTI(BASE),free)(multi);
327 FN(EL,free)(el);
328 isl_space_free(multi_space);
329 isl_space_free(el_space);
330 return NULL;
331 }
332
333 /* Reset the space of "multi". This function is called from isl_pw_templ.c
334 * and doesn't know if the space of an element object is represented
335 * directly or through its domain. It therefore passes along both,
336 * which we pass along to the element function since we don't how
337 * that is represented either.
338 */
339 __isl_give MULTI(BASE) *FN(MULTI(BASE),reset_space_and_domain)(
340 __isl_take MULTI(BASE) *multi, __isl_take isl_space *space,
341 __isl_take isl_space *domain)
342 {
343 int i;
344
345 multi = FN(MULTI(BASE),cow)(multi);
346 if (!multi || !space || !domain)
347 goto error;
348
349 for (i = 0; i < multi->n; ++i) {
350 multi->p[i] = FN(EL,reset_domain_space)(multi->p[i],
351 isl_space_copy(domain));
352 if (!multi->p[i])
353 goto error;
354 }
355 isl_space_free(domain);
356 isl_space_free(multi->space);
357 multi->space = space;
358
359 return multi;
360 error:
361 isl_space_free(domain);
362 isl_space_free(space);
363 FN(MULTI(BASE),free)(multi);
364 return NULL;
365 }
366
367 __isl_give MULTI(BASE) *FN(MULTI(BASE),reset_domain_space)(
368 __isl_take MULTI(BASE) *multi, __isl_take isl_space *domain)
369 {
370 isl_space *space;
371
372 space = isl_space_extend_domain_with_range(isl_space_copy(domain),
373 isl_space_copy(multi->space));
374 return FN(MULTI(BASE),reset_space_and_domain)(multi, space, domain);
375 }
376
377 __isl_give MULTI(BASE) *FN(MULTI(BASE),reset_space)(
378 __isl_take MULTI(BASE) *multi, __isl_take isl_space *space)
379 {
380 isl_space *domain;
381
382 domain = isl_space_domain(isl_space_copy(space));
383 return FN(MULTI(BASE),reset_space_and_domain)(multi, space, domain);
384 }
385
386 /* Set the id of the given dimension of "multi" to "id".
387 */
388 __isl_give MULTI(BASE) *FN(MULTI(BASE),set_dim_id)(
389 __isl_take MULTI(BASE) *multi,
390 enum isl_dim_type type, unsigned pos, __isl_take isl_id *id)
391 {
392 isl_space *space;
393
394 multi = FN(MULTI(BASE),cow)(multi);
395 if (!multi || !id)
396 goto error;
397
398 space = FN(MULTI(BASE),get_space)(multi);
399 space = isl_space_set_dim_id(space, type, pos, id);
400
401 return FN(MULTI(BASE),reset_space)(multi, space);
402 error:
403 isl_id_free(id);
404 FN(MULTI(BASE),free)(multi);
405 return NULL;
406 }
407
408 __isl_give MULTI(BASE) *FN(MULTI(BASE),set_tuple_name)(
409 __isl_keep MULTI(BASE) *multi, enum isl_dim_type type,
410 const char *s)
411 {
412 isl_space *space;
413
414 multi = FN(MULTI(BASE),cow)(multi);
415 if (!multi)
416 return NULL;
417
418 space = FN(MULTI(BASE),get_space)(multi);
419 space = isl_space_set_tuple_name(space, type, s);
420
421 return FN(MULTI(BASE),reset_space)(multi, space);
422 }
423
424 __isl_give MULTI(BASE) *FN(MULTI(BASE),set_tuple_id)(
425 __isl_take MULTI(BASE) *multi, enum isl_dim_type type,
426 __isl_take isl_id *id)
427 {
428 isl_space *space;
429
430 multi = FN(MULTI(BASE),cow)(multi);
431 if (!multi)
432 goto error;
433
434 space = FN(MULTI(BASE),get_space)(multi);
435 space = isl_space_set_tuple_id(space, type, id);
436
437 return FN(MULTI(BASE),reset_space)(multi, space);
438 error:
439 isl_id_free(id);
440 return NULL;
441 }
442
443 /* Drop the id on the specified tuple.
444 */
445 __isl_give MULTI(BASE) *FN(MULTI(BASE),reset_tuple_id)(
446 __isl_take MULTI(BASE) *multi, enum isl_dim_type type)
447 {
448 isl_space *space;
449
450 if (!multi)
451 return NULL;
452 if (!FN(MULTI(BASE),has_tuple_id)(multi, type))
453 return multi;
454
455 multi = FN(MULTI(BASE),cow)(multi);
456 if (!multi)
457 return NULL;
458
459 space = FN(MULTI(BASE),get_space)(multi);
460 space = isl_space_reset_tuple_id(space, type);
461
462 return FN(MULTI(BASE),reset_space)(multi, space);
463 }
464
465 /* Reset the user pointer on all identifiers of parameters and tuples
466 * of the space of "multi".
467 */
468 __isl_give MULTI(BASE) *FN(MULTI(BASE),reset_user)(
469 __isl_take MULTI(BASE) *multi)
470 {
471 isl_space *space;
472
473 space = FN(MULTI(BASE),get_space)(multi);
474 space = isl_space_reset_user(space);
475
476 return FN(MULTI(BASE),reset_space)(multi, space);
477 }
478
479 __isl_give MULTI(BASE) *FN(MULTI(BASE),realign_domain)(
480 __isl_take MULTI(BASE) *multi, __isl_take isl_reordering *exp)
481 {
482 int i;
483
484 multi = FN(MULTI(BASE),cow)(multi);
485 if (!multi || !exp)
486 goto error;
487
488 for (i = 0; i < multi->n; ++i) {
489 multi->p[i] = FN(EL,realign_domain)(multi->p[i],
490 isl_reordering_copy(exp));
491 if (!multi->p[i])
492 goto error;
493 }
494
495 multi = FN(MULTI(BASE),reset_domain_space)(multi,
496 isl_space_copy(exp->dim));
497
498 isl_reordering_free(exp);
499 return multi;
500 error:
501 isl_reordering_free(exp);
502 FN(MULTI(BASE),free)(multi);
503 return NULL;
504 }
505
506 /* Align the parameters of "multi" to those of "model".
507 */
508 __isl_give MULTI(BASE) *FN(MULTI(BASE),align_params)(
509 __isl_take MULTI(BASE) *multi, __isl_take isl_space *model)
510 {
511 isl_ctx *ctx;
512 isl_reordering *exp;
513
514 if (!multi || !model)
515 goto error;
516
517 if (isl_space_match(multi->space, isl_dim_param,
518 model, isl_dim_param)) {
519 isl_space_free(model);
520 return multi;
521 }
522
523 ctx = isl_space_get_ctx(model);
524 if (!isl_space_has_named_params(model))
525 isl_die(ctx, isl_error_invalid,
526 "model has unnamed parameters", goto error);
527 if (!isl_space_has_named_params(multi->space))
528 isl_die(ctx, isl_error_invalid,
529 "input has unnamed parameters", goto error);
530
531 model = isl_space_params(model);
532 exp = isl_parameter_alignment_reordering(multi->space, model);
533 exp = isl_reordering_extend_space(exp,
534 FN(MULTI(BASE),get_domain_space)(multi));
535 multi = FN(MULTI(BASE),realign_domain)(multi, exp);
536
537 isl_space_free(model);
538 return multi;
539 error:
540 isl_space_free(model);
541 FN(MULTI(BASE),free)(multi);
542 return NULL;
543 }
544
545 __isl_give MULTI(BASE) *FN(FN(MULTI(BASE),from),LIST(BASE))(
546 __isl_take isl_space *space, __isl_take LIST(EL) *list)
547 {
548 int i;
549 int n;
550 isl_ctx *ctx;
551 MULTI(BASE) *multi;
552
553 if (!space || !list)
554 goto error;
555
556 ctx = isl_space_get_ctx(space);
557 n = FN(FN(LIST(EL),n),BASE)(list);
558 if (n != isl_space_dim(space, isl_dim_out))
559 isl_die(ctx, isl_error_invalid,
560 "invalid number of elements in list", goto error);
561
562 multi = FN(MULTI(BASE),alloc)(isl_space_copy(space));
563 for (i = 0; i < n; ++i) {
564 multi = FN(FN(MULTI(BASE),set),BASE)(multi, i,
565 FN(FN(LIST(EL),get),BASE)(list, i));
566 }
567
568 isl_space_free(space);
569 FN(LIST(EL),free)(list);
570 return multi;
571 error:
572 isl_space_free(space);
573 FN(LIST(EL),free)(list);
574 return NULL;
575 }
576
577 #ifndef NO_IDENTITY
578 /* Create a multi expression in the given space that maps each
579 * input dimension to the corresponding output dimension.
580 */
581 __isl_give MULTI(BASE) *FN(MULTI(BASE),identity)(__isl_take isl_space *space)
582 {
583 int i, n;
584 isl_local_space *ls;
585 MULTI(BASE) *multi;
586
587 if (!space)
588 return NULL;
589
590 if (isl_space_is_set(space))
591 isl_die(isl_space_get_ctx(space), isl_error_invalid,
592 "expecting map space", goto error);
593
594 n = isl_space_dim(space, isl_dim_out);
595 if (n != isl_space_dim(space, isl_dim_in))
596 isl_die(isl_space_get_ctx(space), isl_error_invalid,
597 "number of input and output dimensions needs to be "
598 "the same", goto error);
599
600 multi = FN(MULTI(BASE),alloc)(isl_space_copy(space));
601
602 if (!n) {
603 isl_space_free(space);
604 return multi;
605 }
606
607 space = isl_space_domain(space);
608 ls = isl_local_space_from_space(space);
609
610 for (i = 0; i < n; ++i) {
611 EL *el;
612 el = FN(EL,var_on_domain)(isl_local_space_copy(ls),
613 isl_dim_set, i);
614 multi = FN(FN(MULTI(BASE),set),BASE)(multi, i, el);
615 }
616
617 isl_local_space_free(ls);
618
619 return multi;
620 error:
621 isl_space_free(space);
622 return NULL;
623 }
624 #endif
625
626 #ifndef NO_ZERO
627 /* Construct a multi expression in the given space with value zero in
628 * each of the output dimensions.
629 */
630 __isl_give MULTI(BASE) *FN(MULTI(BASE),zero)(__isl_take isl_space *space)
631 {
632 int n;
633 MULTI(BASE) *multi;
634
635 if (!space)
636 return NULL;
637
638 n = isl_space_dim(space , isl_dim_out);
639 multi = FN(MULTI(BASE),alloc)(isl_space_copy(space));
640
641 if (!n)
642 isl_space_free(space);
643 else {
644 int i;
645 isl_local_space *ls;
646 EL *el;
647
648 space = isl_space_domain(space);
649 ls = isl_local_space_from_space(space);
650 el = FN(EL,zero_on_domain)(ls);
651
652 for (i = 0; i < n; ++i)
653 multi = FN(FN(MULTI(BASE),set),BASE)(multi, i,
654 FN(EL,copy)(el));
655
656 FN(EL,free)(el);
657 }
658
659 return multi;
660 }
661 #endif
662
663 #ifndef NO_FROM_BASE
664 /* Create a multiple expression with a single output/set dimension
665 * equal to "el".
666 * For most multiple expression types, the base type has a single
667 * output/set dimension and the space of the result is therefore
668 * the same as the space of the input.
669 * In the case of isl_multi_union_pw_aff, however, the base type
670 * lives in a parameter space and we therefore need to add
671 * a single set dimension.
672 */
673 __isl_give MULTI(BASE) *FN(FN(MULTI(BASE),from),BASE)(__isl_take EL *el)
674 {
675 isl_space *space;
676 MULTI(BASE) *multi;
677
678 space = FN(EL,get_space(el));
679 if (isl_space_is_params(space)) {
680 space = isl_space_set_from_params(space);
681 space = isl_space_add_dims(space, isl_dim_set, 1);
682 }
683 multi = FN(MULTI(BASE),alloc)(space);
684 multi = FN(FN(MULTI(BASE),set),BASE)(multi, 0, el);
685
686 return multi;
687 }
688 #endif
689
690 __isl_give MULTI(BASE) *FN(MULTI(BASE),drop_dims)(
691 __isl_take MULTI(BASE) *multi,
692 enum isl_dim_type type, unsigned first, unsigned n)
693 {
694 int i;
695 unsigned dim;
696
697 multi = FN(MULTI(BASE),cow)(multi);
698 if (!multi)
699 return NULL;
700
701 dim = FN(MULTI(BASE),dim)(multi, type);
702 if (first + n > dim || first + n < first)
703 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
704 "index out of bounds",
705 return FN(MULTI(BASE),cow)(multi));
706
707 multi->space = isl_space_drop_dims(multi->space, type, first, n);
708 if (!multi->space)
709 return FN(MULTI(BASE),cow)(multi);
710
711 if (type == isl_dim_out) {
712 for (i = 0; i < n; ++i)
713 FN(EL,free)(multi->p[first + i]);
714 for (i = first; i + n < multi->n; ++i)
715 multi->p[i] = multi->p[i + n];
716 multi->n -= n;
717
718 return multi;
719 }
720
721 for (i = 0; i < multi->n; ++i) {
722 multi->p[i] = FN(EL,drop_dims)(multi->p[i], type, first, n);
723 if (!multi->p[i])
724 return FN(MULTI(BASE),cow)(multi);
725 }
726
727 return multi;
728 }
729
730 /* Align the parameters of "multi1" and "multi2" (if needed) and call "fn".
731 */
732 static __isl_give MULTI(BASE) *FN(MULTI(BASE),align_params_multi_multi_and)(
733 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2,
734 __isl_give MULTI(BASE) *(*fn)(__isl_take MULTI(BASE) *multi1,
735 __isl_take MULTI(BASE) *multi2))
736 {
737 isl_ctx *ctx;
738
739 if (!multi1 || !multi2)
740 goto error;
741 if (isl_space_match(multi1->space, isl_dim_param,
742 multi2->space, isl_dim_param))
743 return fn(multi1, multi2);
744 ctx = FN(MULTI(BASE),get_ctx)(multi1);
745 if (!isl_space_has_named_params(multi1->space) ||
746 !isl_space_has_named_params(multi2->space))
747 isl_die(ctx, isl_error_invalid,
748 "unaligned unnamed parameters", goto error);
749 multi1 = FN(MULTI(BASE),align_params)(multi1,
750 FN(MULTI(BASE),get_space)(multi2));
751 multi2 = FN(MULTI(BASE),align_params)(multi2,
752 FN(MULTI(BASE),get_space)(multi1));
753 return fn(multi1, multi2);
754 error:
755 FN(MULTI(BASE),free)(multi1);
756 FN(MULTI(BASE),free)(multi2);
757 return NULL;
758 }
759
760 /* Given two MULTI(BASE)s A -> B and C -> D,
761 * construct a MULTI(BASE) (A * C) -> [B -> D].
762 *
763 * The parameters are assumed to have been aligned.
764 */
765 static __isl_give MULTI(BASE) *FN(MULTI(BASE),range_product_aligned)(
766 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2)
767 {
768 int i, n1, n2;
769 EL *el;
770 isl_space *space;
771 MULTI(BASE) *res;
772
773 if (!multi1 || !multi2)
774 goto error;
775
776 space = isl_space_range_product(FN(MULTI(BASE),get_space)(multi1),
777 FN(MULTI(BASE),get_space)(multi2));
778 res = FN(MULTI(BASE),alloc)(space);
779
780 n1 = FN(MULTI(BASE),dim)(multi1, isl_dim_out);
781 n2 = FN(MULTI(BASE),dim)(multi2, isl_dim_out);
782
783 for (i = 0; i < n1; ++i) {
784 el = FN(FN(MULTI(BASE),get),BASE)(multi1, i);
785 res = FN(FN(MULTI(BASE),set),BASE)(res, i, el);
786 }
787
788 for (i = 0; i < n2; ++i) {
789 el = FN(FN(MULTI(BASE),get),BASE)(multi2, i);
790 res = FN(FN(MULTI(BASE),set),BASE)(res, n1 + i, el);
791 }
792
793 FN(MULTI(BASE),free)(multi1);
794 FN(MULTI(BASE),free)(multi2);
795 return res;
796 error:
797 FN(MULTI(BASE),free)(multi1);
798 FN(MULTI(BASE),free)(multi2);
799 return NULL;
800 }
801
802 /* Given two MULTI(BASE)s A -> B and C -> D,
803 * construct a MULTI(BASE) (A * C) -> [B -> D].
804 */
805 __isl_give MULTI(BASE) *FN(MULTI(BASE),range_product)(
806 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2)
807 {
808 return FN(MULTI(BASE),align_params_multi_multi_and)(multi1, multi2,
809 &FN(MULTI(BASE),range_product_aligned));
810 }
811
812 /* Is the range of "multi" a wrapped relation?
813 */
814 isl_bool FN(MULTI(BASE),range_is_wrapping)(__isl_keep MULTI(BASE) *multi)
815 {
816 if (!multi)
817 return isl_bool_error;
818 return isl_space_range_is_wrapping(multi->space);
819 }
820
821 /* Given a function A -> [B -> C], extract the function A -> B.
822 */
823 __isl_give MULTI(BASE) *FN(MULTI(BASE),range_factor_domain)(
824 __isl_take MULTI(BASE) *multi)
825 {
826 isl_space *space;
827 int total, keep;
828
829 if (!multi)
830 return NULL;
831 if (!isl_space_range_is_wrapping(multi->space))
832 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
833 "range is not a product",
834 return FN(MULTI(BASE),free)(multi));
835
836 space = FN(MULTI(BASE),get_space)(multi);
837 total = isl_space_dim(space, isl_dim_out);
838 space = isl_space_range_factor_domain(space);
839 keep = isl_space_dim(space, isl_dim_out);
840 multi = FN(MULTI(BASE),drop_dims)(multi,
841 isl_dim_out, keep, total - keep);
842 multi = FN(MULTI(BASE),reset_space)(multi, space);
843
844 return multi;
845 }
846
847 /* Given a function A -> [B -> C], extract the function A -> C.
848 */
849 __isl_give MULTI(BASE) *FN(MULTI(BASE),range_factor_range)(
850 __isl_take MULTI(BASE) *multi)
851 {
852 isl_space *space;
853 int total, keep;
854
855 if (!multi)
856 return NULL;
857 if (!isl_space_range_is_wrapping(multi->space))
858 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
859 "range is not a product",
860 return FN(MULTI(BASE),free)(multi));
861
862 space = FN(MULTI(BASE),get_space)(multi);
863 total = isl_space_dim(space, isl_dim_out);
864 space = isl_space_range_factor_range(space);
865 keep = isl_space_dim(space, isl_dim_out);
866 multi = FN(MULTI(BASE),drop_dims)(multi, isl_dim_out, 0, total - keep);
867 multi = FN(MULTI(BASE),reset_space)(multi, space);
868
869 return multi;
870 }
871
872 /* Given a function [B -> C], extract the function C.
873 */
874 __isl_give MULTI(BASE) *FN(MULTI(BASE),factor_range)(
875 __isl_take MULTI(BASE) *multi)
876 {
877 isl_space *space;
878 int total, keep;
879
880 if (!multi)
881 return NULL;
882 if (!isl_space_is_wrapping(multi->space))
883 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
884 "not a product", return FN(MULTI(BASE),free)(multi));
885
886 space = FN(MULTI(BASE),get_space)(multi);
887 total = isl_space_dim(space, isl_dim_out);
888 space = isl_space_factor_range(space);
889 keep = isl_space_dim(space, isl_dim_out);
890 multi = FN(MULTI(BASE),drop_dims)(multi, isl_dim_out, 0, total - keep);
891 multi = FN(MULTI(BASE),reset_space)(multi, space);
892
893 return multi;
894 }
895
896 #ifndef NO_PRODUCT
897 /* Given two MULTI(BASE)s A -> B and C -> D,
898 * construct a MULTI(BASE) [A -> C] -> [B -> D].
899 *
900 * The parameters are assumed to have been aligned.
901 */
902 __isl_give MULTI(BASE) *FN(MULTI(BASE),product_aligned)(
903 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2)
904 {
905 int i;
906 EL *el;
907 isl_space *space;
908 MULTI(BASE) *res;
909 int in1, in2, out1, out2;
910
911 in1 = FN(MULTI(BASE),dim)(multi1, isl_dim_in);
912 in2 = FN(MULTI(BASE),dim)(multi2, isl_dim_in);
913 out1 = FN(MULTI(BASE),dim)(multi1, isl_dim_out);
914 out2 = FN(MULTI(BASE),dim)(multi2, isl_dim_out);
915 space = isl_space_product(FN(MULTI(BASE),get_space)(multi1),
916 FN(MULTI(BASE),get_space)(multi2));
917 res = FN(MULTI(BASE),alloc)(isl_space_copy(space));
918 space = isl_space_domain(space);
919
920 for (i = 0; i < out1; ++i) {
921 el = FN(FN(MULTI(BASE),get),BASE)(multi1, i);
922 el = FN(EL,insert_dims)(el, isl_dim_in, in1, in2);
923 el = FN(EL,reset_domain_space)(el, isl_space_copy(space));
924 res = FN(FN(MULTI(BASE),set),BASE)(res, i, el);
925 }
926
927 for (i = 0; i < out2; ++i) {
928 el = FN(FN(MULTI(BASE),get),BASE)(multi2, i);
929 el = FN(EL,insert_dims)(el, isl_dim_in, 0, in1);
930 el = FN(EL,reset_domain_space)(el, isl_space_copy(space));
931 res = FN(FN(MULTI(BASE),set),BASE)(res, out1 + i, el);
932 }
933
934 isl_space_free(space);
935 FN(MULTI(BASE),free)(multi1);
936 FN(MULTI(BASE),free)(multi2);
937 return res;
938 }
939
940 /* Given two MULTI(BASE)s A -> B and C -> D,
941 * construct a MULTI(BASE) [A -> C] -> [B -> D].
942 */
943 __isl_give MULTI(BASE) *FN(MULTI(BASE),product)(
944 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2)
945 {
946 return FN(MULTI(BASE),align_params_multi_multi_and)(multi1, multi2,
947 &FN(MULTI(BASE),product_aligned));
948 }
949 #endif
950
951 __isl_give MULTI(BASE) *FN(MULTI(BASE),flatten_range)(
952 __isl_take MULTI(BASE) *multi)
953 {
954 if (!multi)
955 return NULL;
956
957 if (!multi->space->nested[1])
958 return multi;
959
960 multi = FN(MULTI(BASE),cow)(multi);
961 if (!multi)
962 return NULL;
963
964 multi->space = isl_space_flatten_range(multi->space);
965 if (!multi->space)
966 return FN(MULTI(BASE),free)(multi);
967
968 return multi;
969 }
970
971 /* Given two MULTI(BASE)s A -> B and C -> D,
972 * construct a MULTI(BASE) (A * C) -> (B, D).
973 */
974 __isl_give MULTI(BASE) *FN(MULTI(BASE),flat_range_product)(
975 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2)
976 {
977 MULTI(BASE) *multi;
978
979 multi = FN(MULTI(BASE),range_product)(multi1, multi2);
980 multi = FN(MULTI(BASE),flatten_range)(multi);
981 return multi;
982 }
983
984 /* Given two multi expressions, "multi1"
985 *
986 * [A] -> [B1 B2]
987 *
988 * where B2 starts at position "pos", and "multi2"
989 *
990 * [A] -> [D]
991 *
992 * return the multi expression
993 *
994 * [A] -> [B1 D B2]
995 */
996 __isl_give MULTI(BASE) *FN(MULTI(BASE),range_splice)(
997 __isl_take MULTI(BASE) *multi1, unsigned pos,
998 __isl_take MULTI(BASE) *multi2)
999 {
1000 MULTI(BASE) *res;
1001 unsigned dim;
1002
1003 if (!multi1 || !multi2)
1004 goto error;
1005
1006 dim = FN(MULTI(BASE),dim)(multi1, isl_dim_out);
1007 if (pos > dim)
1008 isl_die(FN(MULTI(BASE),get_ctx)(multi1), isl_error_invalid,
1009 "index out of bounds", goto error);
1010
1011 res = FN(MULTI(BASE),copy)(multi1);
1012 res = FN(MULTI(BASE),drop_dims)(res, isl_dim_out, pos, dim - pos);
1013 multi1 = FN(MULTI(BASE),drop_dims)(multi1, isl_dim_out, 0, pos);
1014
1015 res = FN(MULTI(BASE),flat_range_product)(res, multi2);
1016 res = FN(MULTI(BASE),flat_range_product)(res, multi1);
1017
1018 return res;
1019 error:
1020 FN(MULTI(BASE),free)(multi1);
1021 FN(MULTI(BASE),free)(multi2);
1022 return NULL;
1023 }
1024
1025 #ifndef NO_SPLICE
1026 /* Given two multi expressions, "multi1"
1027 *
1028 * [A1 A2] -> [B1 B2]
1029 *
1030 * where A2 starts at position "in_pos" and B2 starts at position "out_pos",
1031 * and "multi2"
1032 *
1033 * [C] -> [D]
1034 *
1035 * return the multi expression
1036 *
1037 * [A1 C A2] -> [B1 D B2]
1038 *
1039 * We first insert input dimensions to obtain
1040 *
1041 * [A1 C A2] -> [B1 B2]
1042 *
1043 * and
1044 *
1045 * [A1 C A2] -> [D]
1046 *
1047 * and then apply range_splice.
1048 */
1049 __isl_give MULTI(BASE) *FN(MULTI(BASE),splice)(
1050 __isl_take MULTI(BASE) *multi1, unsigned in_pos, unsigned out_pos,
1051 __isl_take MULTI(BASE) *multi2)
1052 {
1053 unsigned n_in1;
1054 unsigned n_in2;
1055
1056 if (!multi1 || !multi2)
1057 goto error;
1058
1059 n_in1 = FN(MULTI(BASE),dim)(multi1, isl_dim_in);
1060 if (in_pos > n_in1)
1061 isl_die(FN(MULTI(BASE),get_ctx)(multi1), isl_error_invalid,
1062 "index out of bounds", goto error);
1063
1064 n_in2 = FN(MULTI(BASE),dim)(multi2, isl_dim_in);
1065
1066 multi1 = FN(MULTI(BASE),insert_dims)(multi1, isl_dim_in, in_pos, n_in2);
1067 multi2 = FN(MULTI(BASE),insert_dims)(multi2, isl_dim_in, n_in2,
1068 n_in1 - in_pos);
1069 multi2 = FN(MULTI(BASE),insert_dims)(multi2, isl_dim_in, 0, in_pos);
1070
1071 return FN(MULTI(BASE),range_splice)(multi1, out_pos, multi2);
1072 error:
1073 FN(MULTI(BASE),free)(multi1);
1074 FN(MULTI(BASE),free)(multi2);
1075 return NULL;
1076 }
1077 #endif
1078
1079 /* This function is currently only used from isl_aff.c
1080 */
1081 static __isl_give MULTI(BASE) *FN(MULTI(BASE),bin_op)(
1082 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2,
1083 __isl_give EL *(*fn)(__isl_take EL *, __isl_take EL *))
1084 __attribute__ ((unused));
1085
1086 /* Pairwise perform "fn" to the elements of "multi1" and "multi2" and
1087 * return the result.
1088 */
1089 static __isl_give MULTI(BASE) *FN(MULTI(BASE),bin_op)(
1090 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2,
1091 __isl_give EL *(*fn)(__isl_take EL *, __isl_take EL *))
1092 {
1093 int i;
1094 isl_ctx *ctx;
1095
1096 multi1 = FN(MULTI(BASE),cow)(multi1);
1097 if (!multi1 || !multi2)
1098 goto error;
1099
1100 ctx = FN(MULTI(BASE),get_ctx)(multi1);
1101 if (!isl_space_is_equal(multi1->space, multi2->space))
1102 isl_die(ctx, isl_error_invalid,
1103 "spaces don't match", goto error);
1104
1105 for (i = 0; i < multi1->n; ++i) {
1106 multi1->p[i] = fn(multi1->p[i], FN(EL,copy)(multi2->p[i]));
1107 if (!multi1->p[i])
1108 goto error;
1109 }
1110
1111 FN(MULTI(BASE),free)(multi2);
1112 return multi1;
1113 error:
1114 FN(MULTI(BASE),free)(multi1);
1115 FN(MULTI(BASE),free)(multi2);
1116 return NULL;
1117 }
1118
1119 /* Add "multi2" from "multi1" and return the result.
1120 *
1121 * The parameters of "multi1" and "multi2" are assumed to have been aligned.
1122 */
1123 static __isl_give MULTI(BASE) *FN(MULTI(BASE),add_aligned)(
1124 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2)
1125 {
1126 return FN(MULTI(BASE),bin_op)(multi1, multi2, &FN(EL,add));
1127 }
1128
1129 /* Add "multi2" from "multi1" and return the result.
1130 */
1131 __isl_give MULTI(BASE) *FN(MULTI(BASE),add)(__isl_take MULTI(BASE) *multi1,
1132 __isl_take MULTI(BASE) *multi2)
1133 {
1134 return FN(MULTI(BASE),align_params_multi_multi_and)(multi1, multi2,
1135 &FN(MULTI(BASE),add_aligned));
1136 }
1137
1138 /* Subtract "multi2" from "multi1" and return the result.
1139 *
1140 * The parameters of "multi1" and "multi2" are assumed to have been aligned.
1141 */
1142 static __isl_give MULTI(BASE) *FN(MULTI(BASE),sub_aligned)(
1143 __isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2)
1144 {
1145 return FN(MULTI(BASE),bin_op)(multi1, multi2, &FN(EL,sub));
1146 }
1147
1148 /* Subtract "multi2" from "multi1" and return the result.
1149 */
1150 __isl_give MULTI(BASE) *FN(MULTI(BASE),sub)(__isl_take MULTI(BASE) *multi1,
1151 __isl_take MULTI(BASE) *multi2)
1152 {
1153 return FN(MULTI(BASE),align_params_multi_multi_and)(multi1, multi2,
1154 &FN(MULTI(BASE),sub_aligned));
1155 }
1156
1157 /* Multiply the elements of "multi" by "v" and return the result.
1158 */
1159 __isl_give MULTI(BASE) *FN(MULTI(BASE),scale_val)(__isl_take MULTI(BASE) *multi,
1160 __isl_take isl_val *v)
1161 {
1162 int i;
1163
1164 if (!multi || !v)
1165 goto error;
1166
1167 if (isl_val_is_one(v)) {
1168 isl_val_free(v);
1169 return multi;
1170 }
1171
1172 if (!isl_val_is_rat(v))
1173 isl_die(isl_val_get_ctx(v), isl_error_invalid,
1174 "expecting rational factor", goto error);
1175
1176 multi = FN(MULTI(BASE),cow)(multi);
1177 if (!multi)
1178 return NULL;
1179
1180 for (i = 0; i < multi->n; ++i) {
1181 multi->p[i] = FN(EL,scale_val)(multi->p[i], isl_val_copy(v));
1182 if (!multi->p[i])
1183 goto error;
1184 }
1185
1186 isl_val_free(v);
1187 return multi;
1188 error:
1189 isl_val_free(v);
1190 return FN(MULTI(BASE),free)(multi);
1191 }
1192
1193 /* Divide the elements of "multi" by "v" and return the result.
1194 */
1195 __isl_give MULTI(BASE) *FN(MULTI(BASE),scale_down_val)(
1196 __isl_take MULTI(BASE) *multi, __isl_take isl_val *v)
1197 {
1198 int i;
1199
1200 if (!multi || !v)
1201 goto error;
1202
1203 if (isl_val_is_one(v)) {
1204 isl_val_free(v);
1205 return multi;
1206 }
1207
1208 if (!isl_val_is_rat(v))
1209 isl_die(isl_val_get_ctx(v), isl_error_invalid,
1210 "expecting rational factor", goto error);
1211 if (isl_val_is_zero(v))
1212 isl_die(isl_val_get_ctx(v), isl_error_invalid,
1213 "cannot scale down by zero", goto error);
1214
1215 multi = FN(MULTI(BASE),cow)(multi);
1216 if (!multi)
1217 return NULL;
1218
1219 for (i = 0; i < multi->n; ++i) {
1220 multi->p[i] = FN(EL,scale_down_val)(multi->p[i],
1221 isl_val_copy(v));
1222 if (!multi->p[i])
1223 goto error;
1224 }
1225
1226 isl_val_free(v);
1227 return multi;
1228 error:
1229 isl_val_free(v);
1230 return FN(MULTI(BASE),free)(multi);
1231 }
1232
1233 /* Multiply the elements of "multi" by the corresponding element of "mv"
1234 * and return the result.
1235 */
1236 __isl_give MULTI(BASE) *FN(MULTI(BASE),scale_multi_val)(
1237 __isl_take MULTI(BASE) *multi, __isl_take isl_multi_val *mv)
1238 {
1239 int i;
1240
1241 if (!multi || !mv)
1242 goto error;
1243
1244 if (!isl_space_tuple_is_equal(multi->space, isl_dim_out,
1245 mv->space, isl_dim_set))
1246 isl_die(isl_multi_val_get_ctx(mv), isl_error_invalid,
1247 "spaces don't match", goto error);
1248
1249 multi = FN(MULTI(BASE),cow)(multi);
1250 if (!multi)
1251 goto error;
1252
1253 for (i = 0; i < multi->n; ++i) {
1254 isl_val *v;
1255
1256 v = isl_multi_val_get_val(mv, i);
1257 multi->p[i] = FN(EL,scale_val)(multi->p[i], v);
1258 if (!multi->p[i])
1259 goto error;
1260 }
1261
1262 isl_multi_val_free(mv);
1263 return multi;
1264 error:
1265 isl_multi_val_free(mv);
1266 return FN(MULTI(BASE),free)(multi);
1267 }
1268
1269 /* Divide the elements of "multi" by the corresponding element of "mv"
1270 * and return the result.
1271 */
1272 __isl_give MULTI(BASE) *FN(MULTI(BASE),scale_down_multi_val)(
1273 __isl_take MULTI(BASE) *multi, __isl_take isl_multi_val *mv)
1274 {
1275 int i;
1276
1277 if (!multi || !mv)
1278 goto error;
1279
1280 if (!isl_space_tuple_is_equal(multi->space, isl_dim_out,
1281 mv->space, isl_dim_set))
1282 isl_die(isl_multi_val_get_ctx(mv), isl_error_invalid,
1283 "spaces don't match", goto error);
1284
1285 multi = FN(MULTI(BASE),cow)(multi);
1286 if (!multi)
1287 return NULL;
1288
1289 for (i = 0; i < multi->n; ++i) {
1290 isl_val *v;
1291
1292 v = isl_multi_val_get_val(mv, i);
1293 multi->p[i] = FN(EL,scale_down_val)(multi->p[i], v);
1294 if (!multi->p[i])
1295 goto error;
1296 }
1297
1298 isl_multi_val_free(mv);
1299 return multi;
1300 error:
1301 isl_multi_val_free(mv);
1302 return FN(MULTI(BASE),free)(multi);
1303 }
1304
1305 /* Compute the residues of the elements of "multi" modulo
1306 * the corresponding element of "mv" and return the result.
1307 */
1308 __isl_give MULTI(BASE) *FN(MULTI(BASE),mod_multi_val)(
1309 __isl_take MULTI(BASE) *multi, __isl_take isl_multi_val *mv)
1310 {
1311 int i;
1312
1313 if (!multi || !mv)
1314 goto error;
1315
1316 if (!isl_space_tuple_is_equal(multi->space, isl_dim_out,
1317 mv->space, isl_dim_set))
1318 isl_die(isl_multi_val_get_ctx(mv), isl_error_invalid,
1319 "spaces don't match", goto error);
1320
1321 multi = FN(MULTI(BASE),cow)(multi);
1322 if (!multi)
1323 return NULL;
1324
1325 for (i = 0; i < multi->n; ++i) {
1326 isl_val *v;
1327
1328 v = isl_multi_val_get_val(mv, i);
1329 multi->p[i] = FN(EL,mod_val)(multi->p[i], v);
1330 if (!multi->p[i])
1331 goto error;
1332 }
1333
1334 isl_multi_val_free(mv);
1335 return multi;
1336 error:
1337 isl_multi_val_free(mv);
1338 return FN(MULTI(BASE),free)(multi);
1339 }
1340
1341 #ifndef NO_MOVE_DIMS
1342 /* Move the "n" dimensions of "src_type" starting at "src_pos" of "multi"
1343 * to dimensions of "dst_type" at "dst_pos".
1344 *
1345 * We only support moving input dimensions to parameters and vice versa.
1346 */
1347 __isl_give MULTI(BASE) *FN(MULTI(BASE),move_dims)(__isl_take MULTI(BASE) *multi,
1348 enum isl_dim_type dst_type, unsigned dst_pos,
1349 enum isl_dim_type src_type, unsigned src_pos, unsigned n)
1350 {
1351 int i;
1352
1353 if (!multi)
1354 return NULL;
1355
1356 if (n == 0 &&
1357 !isl_space_is_named_or_nested(multi->space, src_type) &&
1358 !isl_space_is_named_or_nested(multi->space, dst_type))
1359 return multi;
1360
1361 if (dst_type == isl_dim_out || src_type == isl_dim_out)
1362 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
1363 "cannot move output/set dimension",
1364 return FN(MULTI(BASE),free)(multi));
1365 if (dst_type == isl_dim_div || src_type == isl_dim_div)
1366 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
1367 "cannot move divs",
1368 return FN(MULTI(BASE),free)(multi));
1369 if (src_pos + n > isl_space_dim(multi->space, src_type))
1370 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
1371 "range out of bounds",
1372 return FN(MULTI(BASE),free)(multi));
1373 if (dst_type == src_type)
1374 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_unsupported,
1375 "moving dims within the same type not supported",
1376 return FN(MULTI(BASE),free)(multi));
1377
1378 multi = FN(MULTI(BASE),cow)(multi);
1379 if (!multi)
1380 return NULL;
1381
1382 multi->space = isl_space_move_dims(multi->space, dst_type, dst_pos,
1383 src_type, src_pos, n);
1384 if (!multi->space)
1385 return FN(MULTI(BASE),free)(multi);
1386
1387 for (i = 0; i < multi->n; ++i) {
1388 multi->p[i] = FN(EL,move_dims)(multi->p[i], dst_type, dst_pos,
1389 src_type, src_pos, n);
1390 if (!multi->p[i])
1391 return FN(MULTI(BASE),free)(multi);
1392 }
1393
1394 return multi;
1395 }
1396 #endif
1397
1398 /* Convert a multiple expression defined over a parameter domain
1399 * into one that is defined over a zero-dimensional set.
1400 */
1401 __isl_give MULTI(BASE) *FN(MULTI(BASE),from_range)(
1402 __isl_take MULTI(BASE) *multi)
1403 {
1404 isl_space *space;
1405
1406 if (!multi)
1407 return NULL;
1408 if (!isl_space_is_set(multi->space))
1409 isl_die(FN(MULTI(BASE),get_ctx)(multi), isl_error_invalid,
1410 "not living in a set space",
1411 return FN(MULTI(BASE),free)(multi));
1412
1413 space = FN(MULTI(BASE),get_space)(multi);
1414 space = isl_space_from_range(space);
1415 multi = FN(MULTI(BASE),reset_space)(multi, space);
1416
1417 return multi;
1418 }
1419
1420 /* Are "multi1" and "multi2" obviously equal?
1421 */
1422 isl_bool FN(MULTI(BASE),plain_is_equal)(__isl_keep MULTI(BASE) *multi1,
1423 __isl_keep MULTI(BASE) *multi2)
1424 {
1425 int i;
1426 isl_bool equal;
1427
1428 if (!multi1 || !multi2)
1429 return isl_bool_error;
1430 if (multi1->n != multi2->n)
1431 return isl_bool_false;
1432 equal = isl_space_is_equal(multi1->space, multi2->space);
1433 if (equal < 0 || !equal)
1434 return equal;
1435
1436 for (i = 0; i < multi1->n; ++i) {
1437 equal = FN(EL,plain_is_equal)(multi1->p[i], multi2->p[i]);
1438 if (equal < 0 || !equal)
1439 return equal;
1440 }
1441
1442 return isl_bool_true;
1443 }
1444
1445 /* Does "multi" involve any NaNs?
1446 */
1447 isl_bool FN(MULTI(BASE),involves_nan)(__isl_keep MULTI(BASE) *multi)
1448 {
1449 int i;
1450
1451 if (!multi)
1452 return isl_bool_error;
1453 if (multi->n == 0)
1454 return isl_bool_false;
1455
1456 for (i = 0; i < multi->n; ++i) {
1457 isl_bool has_nan = FN(EL,involves_nan)(multi->p[i]);
1458 if (has_nan < 0 || has_nan)
1459 return has_nan;
1460 }
1461
1462 return isl_bool_false;
1463 }
1464
1465 #ifndef NO_DOMAIN
1466 /* Return the shared domain of the elements of "multi".
1467 */
1468 __isl_give isl_set *FN(MULTI(BASE),domain)(__isl_take MULTI(BASE) *multi)
1469 {
1470 int i;
1471 isl_set *dom;
1472
1473 if (!multi)
1474 return NULL;
1475
1476 dom = isl_set_universe(FN(MULTI(BASE),get_domain_space)(multi));
1477 for (i = 0; i < multi->n; ++i) {
1478 isl_set *dom_i;
1479
1480 dom_i = FN(EL,domain)(FN(FN(MULTI(BASE),get),BASE)(multi, i));
1481 dom = isl_set_intersect(dom, dom_i);
1482 }
1483
1484 FN(MULTI(BASE),free)(multi);
1485 return dom;
1486 }
1487 #endif
1488
1489 #ifndef NO_NEG
1490 /* Return the opposite of "multi".
1491 */
1492 __isl_give MULTI(BASE) *FN(MULTI(BASE),neg)(__isl_take MULTI(BASE) *multi)
1493 {
1494 int i;
1495
1496 multi = FN(MULTI(BASE),cow)(multi);
1497 if (!multi)
1498 return NULL;
1499
1500 for (i = 0; i < multi->n; ++i) {
1501 multi->p[i] = FN(EL,neg)(multi->p[i]);
1502 if (!multi->p[i])
1503 return FN(MULTI(BASE),free)(multi);
1504 }
1505
1506 return multi;
1507 }
1508 #endif
Integer set library