2 * Copyright 2008-2009 Katholieke Universiteit Leuven
3 * Copyright 2011 INRIA Saclay
4 * Copyright 2012-2013 Ecole Normale Superieure
5 * Copyright 2017 Sven Verdoolaege
7 * Use of this software is governed by the MIT license
9 * Written by Sven Verdoolaege, K.U.Leuven, Departement
10 * Computerwetenschappen, Celestijnenlaan 200A, B-3001 Leuven, Belgium
11 * and INRIA Saclay - Ile-de-France, Parc Club Orsay Universite,
12 * ZAC des vignes, 4 rue Jacques Monod, 91893 Orsay, France
13 * and Ecole Normale Superieure, 45 rue d’Ulm, 75230 Paris, France
17 #include <isl_tarjan.h>
18 #include <isl/printer.h>
20 #include <isl_list_macro.h>
22 #define xS(TYPE,NAME) struct TYPE ## _ ## NAME
23 #define S(TYPE,NAME) xS(TYPE,NAME)
25 isl_ctx
*FN(LIST(EL
),get_ctx
)(__isl_keep
LIST(EL
) *list
)
27 return list
? list
->ctx
: NULL
;
30 __isl_give
LIST(EL
) *FN(LIST(EL
),alloc
)(isl_ctx
*ctx
, int n
)
35 isl_die(ctx
, isl_error_invalid
,
36 "cannot create list of negative length",
38 list
= isl_alloc(ctx
, LIST(EL
),
39 sizeof(LIST(EL
)) + (n
- 1) * sizeof(struct EL
*));
51 __isl_give
LIST(EL
) *FN(LIST(EL
),copy
)(__isl_keep
LIST(EL
) *list
)
60 __isl_give
LIST(EL
) *FN(LIST(EL
),dup
)(__isl_keep
LIST(EL
) *list
)
68 dup
= FN(LIST(EL
),alloc
)(FN(LIST(EL
),get_ctx
)(list
), list
->n
);
71 for (i
= 0; i
< list
->n
; ++i
)
72 dup
= FN(LIST(EL
),add
)(dup
, FN(EL
,copy
)(list
->p
[i
]));
76 __isl_give
LIST(EL
) *FN(LIST(EL
),cow
)(__isl_take
LIST(EL
) *list
)
84 return FN(LIST(EL
),dup
)(list
);
87 /* Make sure "list" has room for at least "n" more pieces.
88 * Always return a list with a single reference.
90 * If there is only one reference to list, we extend it in place.
91 * Otherwise, we create a new LIST(EL) and copy the elements.
93 static __isl_give
LIST(EL
) *FN(LIST(EL
),grow
)(__isl_take
LIST(EL
) *list
, int n
)
101 if (list
->ref
== 1 && list
->n
+ n
<= list
->size
)
104 ctx
= FN(LIST(EL
),get_ctx
)(list
);
105 new_size
= ((list
->n
+ n
+ 1) * 3) / 2;
106 if (list
->ref
== 1) {
107 res
= isl_realloc(ctx
, list
, LIST(EL
),
108 sizeof(LIST(EL
)) + (new_size
- 1) * sizeof(EL
*));
110 return FN(LIST(EL
),free
)(list
);
111 res
->size
= new_size
;
115 if (list
->n
+ n
<= list
->size
&& list
->size
< new_size
)
116 new_size
= list
->size
;
118 res
= FN(LIST(EL
),alloc
)(ctx
, new_size
);
120 return FN(LIST(EL
),free
)(list
);
122 for (i
= 0; i
< list
->n
; ++i
)
123 res
= FN(LIST(EL
),add
)(res
, FN(EL
,copy
)(list
->p
[i
]));
125 FN(LIST(EL
),free
)(list
);
129 /* Check that "index" is a valid position in "list".
131 static isl_stat
FN(LIST(EL
),check_index
)(__isl_keep
LIST(EL
) *list
, int index
)
134 return isl_stat_error
;
135 if (index
< 0 || index
>= list
->n
)
136 isl_die(FN(LIST(EL
),get_ctx
)(list
), isl_error_invalid
,
137 "index out of bounds", return isl_stat_error
);
141 __isl_give
LIST(EL
) *FN(LIST(EL
),add
)(__isl_take
LIST(EL
) *list
,
142 __isl_take
struct EL
*el
)
144 list
= FN(LIST(EL
),grow
)(list
, 1);
147 list
->p
[list
->n
] = el
;
152 FN(LIST(EL
),free
)(list
);
156 /* Remove the "n" elements starting at "first" from "list".
158 __isl_give
LIST(EL
) *FN(LIST(EL
),drop
)(__isl_take
LIST(EL
) *list
,
159 unsigned first
, unsigned n
)
165 if (first
+ n
> list
->n
|| first
+ n
< first
)
166 isl_die(list
->ctx
, isl_error_invalid
,
167 "index out of bounds", return FN(LIST(EL
),free
)(list
));
170 list
= FN(LIST(EL
),cow
)(list
);
173 for (i
= 0; i
< n
; ++i
)
174 FN(EL
,free
)(list
->p
[first
+ i
]);
175 for (i
= first
; i
+ n
< list
->n
; ++i
)
176 list
->p
[i
] = list
->p
[i
+ n
];
181 /* Remove all elements from "list".
183 __isl_give
LIST(EL
) *FN(LIST(EL
),clear
)(__isl_take
LIST(EL
) *list
)
187 return FN(LIST(EL
),drop
)(list
, 0, list
->n
);
190 /* Insert "el" at position "pos" in "list".
192 * If there is only one reference to "list" and if it already has space
193 * for one extra element, we insert it directly into "list".
194 * Otherwise, we create a new list consisting of "el" and copied
195 * elements from "list".
197 __isl_give
LIST(EL
) *FN(LIST(EL
),insert
)(__isl_take
LIST(EL
) *list
,
198 unsigned pos
, __isl_take
struct EL
*el
)
206 ctx
= FN(LIST(EL
),get_ctx
)(list
);
208 isl_die(ctx
, isl_error_invalid
,
209 "index out of bounds", goto error
);
211 if (list
->ref
== 1 && list
->size
> list
->n
) {
212 for (i
= list
->n
; i
> pos
; --i
)
213 list
->p
[i
] = list
->p
[i
- 1];
219 res
= FN(LIST(EL
),alloc
)(ctx
, list
->n
+ 1);
220 for (i
= 0; i
< pos
; ++i
)
221 res
= FN(LIST(EL
),add
)(res
, FN(EL
,copy
)(list
->p
[i
]));
222 res
= FN(LIST(EL
),add
)(res
, el
);
223 for (i
= pos
; i
< list
->n
; ++i
)
224 res
= FN(LIST(EL
),add
)(res
, FN(EL
,copy
)(list
->p
[i
]));
225 FN(LIST(EL
),free
)(list
);
230 FN(LIST(EL
),free
)(list
);
234 __isl_null
LIST(EL
) *FN(LIST(EL
),free
)(__isl_take
LIST(EL
) *list
)
244 isl_ctx_deref(list
->ctx
);
245 for (i
= 0; i
< list
->n
; ++i
)
246 FN(EL
,free
)(list
->p
[i
]);
252 /* Return the number of elements in "list".
254 isl_size
FN(LIST(EL
),size
)(__isl_keep
LIST(EL
) *list
)
256 return list
? list
->n
: isl_size_error
;
259 /* This is an alternative name for the function above.
261 isl_size
FN(FN(LIST(EL
),n
),EL_BASE
)(__isl_keep
LIST(EL
) *list
)
263 return FN(LIST(EL
),size
)(list
);
266 /* Return the element at position "index" in "list".
268 __isl_keep EL
*FN(LIST(EL
),peek
)(__isl_keep
LIST(EL
) *list
, int index
)
270 if (FN(LIST(EL
),check_index
)(list
, index
) < 0)
272 return list
->p
[index
];
275 /* Return a copy of the element at position "index" in "list".
277 __isl_give EL
*FN(LIST(EL
),get_at
)(__isl_keep
LIST(EL
) *list
, int index
)
279 return FN(EL
,copy
)(FN(LIST(EL
),peek
)(list
, index
));
282 /* This is an alternative name for the function above.
284 __isl_give EL
*FN(FN(LIST(EL
),get
),EL_BASE
)(__isl_keep
LIST(EL
) *list
,
287 return FN(LIST(EL
),get_at
)(list
, index
);
290 /* Replace the element at position "index" in "list" by "el".
292 __isl_give
LIST(EL
) *FN(LIST(EL
),set_at
)(__isl_take
LIST(EL
) *list
,
293 int index
, __isl_take EL
*el
)
297 if (FN(LIST(EL
),check_index
)(list
, index
) < 0)
299 if (list
->p
[index
] == el
) {
303 list
= FN(LIST(EL
),cow
)(list
);
306 FN(EL
,free
)(list
->p
[index
]);
311 FN(LIST(EL
),free
)(list
);
315 /* This is an alternative name for the function above.
317 __isl_give
LIST(EL
) *FN(FN(LIST(EL
),set
),EL_BASE
)(__isl_take
LIST(EL
) *list
,
318 int index
, __isl_take EL
*el
)
320 return FN(LIST(EL
),set_at
)(list
, index
, el
);
323 /* Return the element at position "index" of "list".
324 * This may be either a copy or the element itself
325 * if there is only one reference to "list".
326 * This allows the element to be modified inplace
327 * if both the list and the element have only a single reference.
328 * The caller is not allowed to modify "list" between
329 * this call to isl_list_*_take_* and a subsequent call
330 * to isl_list_*_restore_*.
331 * The only exception is that isl_list_*_free can be called instead.
333 static __isl_give EL
*FN(FN(LIST(EL
),take
),EL_BASE
)(__isl_keep
LIST(EL
) *list
,
338 if (FN(LIST(EL
),check_index
)(list
, index
) < 0)
341 return FN(FN(LIST(EL
),get
),EL_BASE
)(list
, index
);
343 list
->p
[index
] = NULL
;
347 /* Set the element at position "index" of "list" to "el",
348 * where the position may be empty due to a previous call
349 * to isl_list_*_take_*.
351 static __isl_give
LIST(EL
) *FN(FN(LIST(EL
),restore
),EL_BASE
)(
352 __isl_take
LIST(EL
) *list
, int index
, __isl_take EL
*el
)
354 return FN(FN(LIST(EL
),set
),EL_BASE
)(list
, index
, el
);
357 /* Swap the elements of "list" in positions "pos1" and "pos2".
359 __isl_give
LIST(EL
) *FN(LIST(EL
),swap
)(__isl_take
LIST(EL
) *list
,
360 unsigned pos1
, unsigned pos2
)
366 el1
= FN(FN(LIST(EL
),take
),EL_BASE
)(list
, pos1
);
367 el2
= FN(FN(LIST(EL
),take
),EL_BASE
)(list
, pos2
);
368 list
= FN(FN(LIST(EL
),restore
),EL_BASE
)(list
, pos1
, el2
);
369 list
= FN(FN(LIST(EL
),restore
),EL_BASE
)(list
, pos2
, el1
);
373 /* Reverse the elements of "list".
375 __isl_give
LIST(EL
) *FN(LIST(EL
),reverse
)(__isl_take
LIST(EL
) *list
)
379 n
= FN(LIST(EL
),size
)(list
);
380 for (i
= 0; i
< n
- 1 - i
; ++i
)
381 list
= FN(LIST(EL
),swap
)(list
, i
, n
- 1 - i
);
385 isl_stat
FN(LIST(EL
),foreach
)(__isl_keep
LIST(EL
) *list
,
386 isl_stat (*fn
)(__isl_take EL
*el
, void *user
), void *user
)
391 return isl_stat_error
;
393 for (i
= 0; i
< list
->n
; ++i
) {
394 EL
*el
= FN(EL
,copy
)(list
->p
[i
]);
396 return isl_stat_error
;
397 if (fn(el
, user
) < 0)
398 return isl_stat_error
;
404 /* Does "test" succeed on every element of "list"?
406 isl_bool
FN(LIST(EL
),every
)(__isl_keep
LIST(EL
) *list
,
407 isl_bool (*test
)(__isl_keep EL
*el
, void *user
), void *user
)
412 return isl_bool_error
;
414 for (i
= 0; i
< list
->n
; ++i
) {
417 r
= test(list
->p
[i
], user
);
422 return isl_bool_true
;
425 /* Replace each element in "list" by the result of calling "fn"
428 __isl_give
LIST(EL
) *FN(LIST(EL
),map
)(__isl_keep
LIST(EL
) *list
,
429 __isl_give EL
*(*fn
)(__isl_take EL
*el
, void *user
), void *user
)
437 for (i
= 0; i
< n
; ++i
) {
438 EL
*el
= FN(FN(LIST(EL
),take
),EL_BASE
)(list
, i
);
440 return FN(LIST(EL
),free
)(list
);
442 list
= FN(FN(LIST(EL
),restore
),EL_BASE
)(list
, i
, el
);
448 /* Internal data structure for isl_*_list_sort.
450 * "cmp" is the original comparison function.
451 * "user" is a user provided pointer that should be passed to "cmp".
453 S(LIST(EL
),sort_data
) {
454 int (*cmp
)(__isl_keep EL
*a
, __isl_keep EL
*b
, void *user
);
458 /* Compare two entries of an isl_*_list based on the user provided
459 * comparison function on pairs of isl_* objects.
461 static int FN(LIST(EL
),cmp
)(const void *a
, const void *b
, void *user
)
463 S(LIST(EL
),sort_data
) *data
= user
;
467 return data
->cmp(*el1
, *el2
, data
->user
);
470 /* Sort the elements of "list" in ascending order according to
471 * comparison function "cmp".
473 __isl_give
LIST(EL
) *FN(LIST(EL
),sort
)(__isl_take
LIST(EL
) *list
,
474 int (*cmp
)(__isl_keep EL
*a
, __isl_keep EL
*b
, void *user
), void *user
)
476 S(LIST(EL
),sort_data
) data
= { cmp
, user
};
482 list
= FN(LIST(EL
),cow
)(list
);
486 if (isl_sort(list
->p
, list
->n
, sizeof(list
->p
[0]),
487 &FN(LIST(EL
),cmp
), &data
) < 0)
488 return FN(LIST(EL
),free
)(list
);
493 /* Internal data structure for isl_*_list_foreach_scc.
495 * "list" is the original list.
496 * "follows" is the user provided callback that defines the edges of the graph.
498 S(LIST(EL
),foreach_scc_data
) {
500 isl_bool (*follows
)(__isl_keep EL
*a
, __isl_keep EL
*b
, void *user
);
504 /* Does element i of data->list follow element j?
506 * Use the user provided callback to find out.
508 static isl_bool
FN(LIST(EL
),follows
)(int i
, int j
, void *user
)
510 S(LIST(EL
),foreach_scc_data
) *data
= user
;
512 return data
->follows(data
->list
->p
[i
], data
->list
->p
[j
],
516 /* Call "fn" on the sublist of "list" that consists of the elements
517 * with indices specified by the "n" elements of "pos".
519 static isl_stat
FN(LIST(EL
),call_on_scc
)(__isl_keep
LIST(EL
) *list
, int *pos
,
520 int n
, isl_stat (*fn
)(__isl_take
LIST(EL
) *scc
, void *user
), void *user
)
526 ctx
= FN(LIST(EL
),get_ctx
)(list
);
527 slice
= FN(LIST(EL
),alloc
)(ctx
, n
);
528 for (i
= 0; i
< n
; ++i
) {
531 el
= FN(EL
,copy
)(list
->p
[pos
[i
]]);
532 slice
= FN(LIST(EL
),add
)(slice
, el
);
535 return fn(slice
, user
);
538 /* Call "fn" on each of the strongly connected components (SCCs) of
539 * the graph with as vertices the elements of "list" and
540 * a directed edge from node b to node a iff follows(a, b)
541 * returns 1. follows should return -1 on error.
543 * If SCC a contains a node i that follows a node j in another SCC b
544 * (i.e., follows(i, j, user) returns 1), then fn will be called on SCC a
545 * after being called on SCC b.
547 * We simply call isl_tarjan_graph_init, extract the SCCs from the result and
548 * call fn on each of them.
550 isl_stat
FN(LIST(EL
),foreach_scc
)(__isl_keep
LIST(EL
) *list
,
551 isl_bool (*follows
)(__isl_keep EL
*a
, __isl_keep EL
*b
, void *user
),
553 isl_stat (*fn
)(__isl_take
LIST(EL
) *scc
, void *user
), void *fn_user
)
555 S(LIST(EL
),foreach_scc_data
) data
= { list
, follows
, follows_user
};
558 struct isl_tarjan_graph
*g
;
561 return isl_stat_error
;
565 return fn(FN(LIST(EL
),copy
)(list
), fn_user
);
567 ctx
= FN(LIST(EL
),get_ctx
)(list
);
569 g
= isl_tarjan_graph_init(ctx
, n
, &FN(LIST(EL
),follows
), &data
);
571 return isl_stat_error
;
577 if (g
->order
[i
] == -1)
578 isl_die(ctx
, isl_error_internal
, "cannot happen",
581 while (g
->order
[i
] != -1) {
584 if (first
== 0 && n
== 0) {
585 isl_tarjan_graph_free(g
);
586 return fn(FN(LIST(EL
),copy
)(list
), fn_user
);
588 if (FN(LIST(EL
),call_on_scc
)(list
, g
->order
+ first
, i
- first
,
594 isl_tarjan_graph_free(g
);
596 return n
> 0 ? isl_stat_error
: isl_stat_ok
;
599 __isl_give
LIST(EL
) *FN(FN(LIST(EL
),from
),EL_BASE
)(__isl_take EL
*el
)
606 ctx
= FN(EL
,get_ctx
)(el
);
607 list
= FN(LIST(EL
),alloc
)(ctx
, 1);
610 list
= FN(LIST(EL
),add
)(list
, el
);
617 /* This function performs the same operation as isl_*_list_from_*,
618 * but is considered as a function on the element when exported.
620 __isl_give
LIST(EL
) *FN(EL
,to_list
)(__isl_take EL
*el
)
622 return FN(FN(LIST(EL
),from
),EL_BASE
)(el
);
625 /* Append the elements of "list2" to "list1", where "list1" is known
626 * to have only a single reference and enough room to hold
627 * the extra elements.
629 static __isl_give
LIST(EL
) *FN(LIST(EL
),concat_inplace
)(
630 __isl_take
LIST(EL
) *list1
, __isl_take
LIST(EL
) *list2
)
634 for (i
= 0; i
< list2
->n
; ++i
)
635 list1
= FN(LIST(EL
),add
)(list1
, FN(EL
,copy
)(list2
->p
[i
]));
636 FN(LIST(EL
),free
)(list2
);
640 /* Concatenate "list1" and "list2".
641 * If "list1" has only one reference and has enough room
642 * for the elements of "list2", the add the elements to "list1" itself.
643 * Otherwise, create a new list to store the result.
645 __isl_give
LIST(EL
) *FN(LIST(EL
),concat
)(__isl_take
LIST(EL
) *list1
,
646 __isl_take
LIST(EL
) *list2
)
652 if (!list1
|| !list2
)
655 if (list1
->ref
== 1 && list1
->n
+ list2
->n
<= list1
->size
)
656 return FN(LIST(EL
),concat_inplace
)(list1
, list2
);
658 ctx
= FN(LIST(EL
),get_ctx
)(list1
);
659 res
= FN(LIST(EL
),alloc
)(ctx
, list1
->n
+ list2
->n
);
660 for (i
= 0; i
< list1
->n
; ++i
)
661 res
= FN(LIST(EL
),add
)(res
, FN(EL
,copy
)(list1
->p
[i
]));
662 for (i
= 0; i
< list2
->n
; ++i
)
663 res
= FN(LIST(EL
),add
)(res
, FN(EL
,copy
)(list2
->p
[i
]));
665 FN(LIST(EL
),free
)(list1
);
666 FN(LIST(EL
),free
)(list2
);
669 FN(LIST(EL
),free
)(list1
);
670 FN(LIST(EL
),free
)(list2
);
674 __isl_give isl_printer
*CAT(isl_printer_print_
,LIST(EL_BASE
))(
675 __isl_take isl_printer
*p
, __isl_keep
LIST(EL
) *list
)
681 p
= isl_printer_print_str(p
, "(");
682 for (i
= 0; i
< list
->n
; ++i
) {
684 p
= isl_printer_print_str(p
, ",");
685 p
= CAT(isl_printer_print_
,EL_BASE
)(p
, list
->p
[i
]);
687 p
= isl_printer_print_str(p
, ")");
695 #define BASE LIST(EL_BASE)
697 #define PRINT_DUMP_DEFAULT 0
698 #include "print_templ.c"
699 #undef PRINT_DUMP_DEFAULT