123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684 |
- /*
- * Copyright 2011 INRIA Saclay
- * Copyright 2012-2014 Ecole Normale Superieure
- * Copyright 2016 Sven Verdoolaege
- *
- * Use of this software is governed by the MIT license
- *
- * Written by Sven Verdoolaege, INRIA Saclay - Ile-de-France,
- * Parc Club Orsay Universite, ZAC des vignes, 4 rue Jacques Monod,
- * 91893 Orsay, France
- * and Ecole Normale Superieure, 45 rue d'Ulm, 75230 Paris, France
- */
- #include <isl/ctx.h>
- #include <isl/val.h>
- #include <isl_aff_private.h>
- #include <isl/map.h>
- #include <isl/set.h>
- #include <isl/schedule.h>
- #include <isl/schedule_node.h>
- #include <isl_sort.h>
- #include <isl/printer.h>
- #include <isl_schedule_private.h>
- #include <isl_schedule_tree.h>
- #include <isl_schedule_node_private.h>
- /* Return a schedule encapsulating the given schedule tree.
- *
- * We currently only allow schedule trees with a domain or extension as root.
- *
- * The leaf field is initialized as a leaf node so that it can be
- * used to represent leaves in the constructed schedule.
- * The reference count is set to -1 since the isl_schedule_tree
- * should never be freed. It is up to the (internal) users of
- * these leaves to ensure that they are only used while the schedule
- * is still alive.
- */
- __isl_give isl_schedule *isl_schedule_from_schedule_tree(isl_ctx *ctx,
- __isl_take isl_schedule_tree *tree)
- {
- enum isl_schedule_node_type type;
- isl_schedule *schedule;
- if (!tree)
- return NULL;
- type = isl_schedule_tree_get_type(tree);
- if (type != isl_schedule_node_domain &&
- type != isl_schedule_node_extension)
- isl_die(isl_schedule_tree_get_ctx(tree), isl_error_unsupported,
- "root of schedule tree should be a domain or extension",
- goto error);
- schedule = isl_calloc_type(ctx, isl_schedule);
- if (!schedule)
- goto error;
- schedule->ref = 1;
- schedule->root = tree;
- schedule->leaf = isl_schedule_tree_leaf(ctx);
- if (!schedule->leaf)
- return isl_schedule_free(schedule);
- return schedule;
- error:
- isl_schedule_tree_free(tree);
- return NULL;
- }
- /* Return a pointer to a schedule with as single node
- * a domain node with the given domain.
- */
- __isl_give isl_schedule *isl_schedule_from_domain(
- __isl_take isl_union_set *domain)
- {
- isl_ctx *ctx;
- isl_schedule_tree *tree;
- ctx = isl_union_set_get_ctx(domain);
- tree = isl_schedule_tree_from_domain(domain);
- return isl_schedule_from_schedule_tree(ctx, tree);
- }
- /* Return a pointer to a schedule with as single node
- * a domain node with an empty domain.
- */
- __isl_give isl_schedule *isl_schedule_empty(__isl_take isl_space *space)
- {
- return isl_schedule_from_domain(isl_union_set_empty(space));
- }
- /* Return a new reference to "sched".
- */
- __isl_give isl_schedule *isl_schedule_copy(__isl_keep isl_schedule *sched)
- {
- if (!sched)
- return NULL;
- sched->ref++;
- return sched;
- }
- /* Return an isl_schedule that is equal to "schedule" and that has only
- * a single reference.
- */
- __isl_give isl_schedule *isl_schedule_cow(__isl_take isl_schedule *schedule)
- {
- isl_ctx *ctx;
- isl_schedule_tree *tree;
- if (!schedule)
- return NULL;
- if (schedule->ref == 1)
- return schedule;
- ctx = isl_schedule_get_ctx(schedule);
- schedule->ref--;
- tree = isl_schedule_tree_copy(schedule->root);
- return isl_schedule_from_schedule_tree(ctx, tree);
- }
- __isl_null isl_schedule *isl_schedule_free(__isl_take isl_schedule *sched)
- {
- if (!sched)
- return NULL;
- if (--sched->ref > 0)
- return NULL;
- isl_schedule_tree_free(sched->root);
- isl_schedule_tree_free(sched->leaf);
- free(sched);
- return NULL;
- }
- /* Replace the root of "schedule" by "tree".
- */
- __isl_give isl_schedule *isl_schedule_set_root(
- __isl_take isl_schedule *schedule, __isl_take isl_schedule_tree *tree)
- {
- if (!schedule || !tree)
- goto error;
- if (schedule->root == tree) {
- isl_schedule_tree_free(tree);
- return schedule;
- }
- schedule = isl_schedule_cow(schedule);
- if (!schedule)
- goto error;
- isl_schedule_tree_free(schedule->root);
- schedule->root = tree;
- return schedule;
- error:
- isl_schedule_free(schedule);
- isl_schedule_tree_free(tree);
- return NULL;
- }
- isl_ctx *isl_schedule_get_ctx(__isl_keep isl_schedule *schedule)
- {
- return schedule ? isl_schedule_tree_get_ctx(schedule->leaf) : NULL;
- }
- /* Return a pointer to the leaf of "schedule".
- */
- __isl_keep isl_schedule_tree *isl_schedule_peek_leaf(
- __isl_keep isl_schedule *schedule)
- {
- return schedule ? schedule->leaf : NULL;
- }
- /* Are "schedule1" and "schedule2" obviously equal to each other?
- */
- isl_bool isl_schedule_plain_is_equal(__isl_keep isl_schedule *schedule1,
- __isl_keep isl_schedule *schedule2)
- {
- if (!schedule1 || !schedule2)
- return isl_bool_error;
- if (schedule1 == schedule2)
- return isl_bool_true;
- return isl_schedule_tree_plain_is_equal(schedule1->root,
- schedule2->root);
- }
- /* Return the (parameter) space of the schedule, i.e., the space
- * of the root domain.
- */
- __isl_give isl_space *isl_schedule_get_space(
- __isl_keep isl_schedule *schedule)
- {
- enum isl_schedule_node_type type;
- isl_space *space;
- isl_union_set *domain;
- if (!schedule)
- return NULL;
- type = isl_schedule_tree_get_type(schedule->root);
- if (type != isl_schedule_node_domain)
- isl_die(isl_schedule_get_ctx(schedule), isl_error_internal,
- "root node not a domain node", return NULL);
- domain = isl_schedule_tree_domain_get_domain(schedule->root);
- space = isl_union_set_get_space(domain);
- isl_union_set_free(domain);
- return space;
- }
- /* Return a pointer to the root of "schedule".
- */
- __isl_give isl_schedule_node *isl_schedule_get_root(
- __isl_keep isl_schedule *schedule)
- {
- isl_ctx *ctx;
- isl_schedule_tree *tree;
- isl_schedule_tree_list *ancestors;
- if (!schedule)
- return NULL;
- ctx = isl_schedule_get_ctx(schedule);
- tree = isl_schedule_tree_copy(schedule->root);
- schedule = isl_schedule_copy(schedule);
- ancestors = isl_schedule_tree_list_alloc(ctx, 0);
- return isl_schedule_node_alloc(schedule, tree, ancestors, NULL);
- }
- /* Return the domain of the root domain node of "schedule".
- */
- __isl_give isl_union_set *isl_schedule_get_domain(
- __isl_keep isl_schedule *schedule)
- {
- if (!schedule)
- return NULL;
- return isl_schedule_tree_domain_get_domain(schedule->root);
- }
- /* Traverse all nodes of "sched" in depth first preorder.
- *
- * If "fn" returns -1 on any of the nodes, then the traversal is aborted.
- * If "fn" returns 0 on any of the nodes, then the subtree rooted
- * at that node is skipped.
- *
- * Return 0 on success and -1 on failure.
- */
- isl_stat isl_schedule_foreach_schedule_node_top_down(
- __isl_keep isl_schedule *sched,
- isl_bool (*fn)(__isl_keep isl_schedule_node *node, void *user),
- void *user)
- {
- isl_schedule_node *node;
- isl_stat r;
- if (!sched)
- return isl_stat_error;
- node = isl_schedule_get_root(sched);
- r = isl_schedule_node_foreach_descendant_top_down(node, fn, user);
- isl_schedule_node_free(node);
- return r;
- }
- /* Traverse the node of "sched" in depth first postorder,
- * allowing the user to modify the visited node.
- * The traversal continues from the node returned by the callback function.
- * It is the responsibility of the user to ensure that this does not
- * lead to an infinite loop. It is safest to always return a pointer
- * to the same position (same ancestors and child positions) as the input node.
- */
- __isl_give isl_schedule *isl_schedule_map_schedule_node_bottom_up(
- __isl_take isl_schedule *schedule,
- __isl_give isl_schedule_node *(*fn)(
- __isl_take isl_schedule_node *node, void *user), void *user)
- {
- isl_schedule_node *node;
- node = isl_schedule_get_root(schedule);
- isl_schedule_free(schedule);
- node = isl_schedule_node_map_descendant_bottom_up(node, fn, user);
- schedule = isl_schedule_node_get_schedule(node);
- isl_schedule_node_free(node);
- return schedule;
- }
- /* Wrapper around isl_schedule_node_reset_user for use as
- * an isl_schedule_map_schedule_node_bottom_up callback.
- */
- static __isl_give isl_schedule_node *reset_user(
- __isl_take isl_schedule_node *node, void *user)
- {
- return isl_schedule_node_reset_user(node);
- }
- /* Reset the user pointer on all identifiers of parameters and tuples
- * in the schedule "schedule".
- */
- __isl_give isl_schedule *isl_schedule_reset_user(
- __isl_take isl_schedule *schedule)
- {
- return isl_schedule_map_schedule_node_bottom_up(schedule, &reset_user,
- NULL);
- }
- /* Wrapper around isl_schedule_node_align_params for use as
- * an isl_schedule_map_schedule_node_bottom_up callback.
- */
- static __isl_give isl_schedule_node *align_params(
- __isl_take isl_schedule_node *node, void *user)
- {
- isl_space *space = user;
- return isl_schedule_node_align_params(node, isl_space_copy(space));
- }
- /* Align the parameters of all nodes in schedule "schedule"
- * to those of "space".
- */
- __isl_give isl_schedule *isl_schedule_align_params(
- __isl_take isl_schedule *schedule, __isl_take isl_space *space)
- {
- schedule = isl_schedule_map_schedule_node_bottom_up(schedule,
- &align_params, space);
- isl_space_free(space);
- return schedule;
- }
- /* Wrapper around isl_schedule_node_pullback_union_pw_multi_aff for use as
- * an isl_schedule_map_schedule_node_bottom_up callback.
- */
- static __isl_give isl_schedule_node *pullback_upma(
- __isl_take isl_schedule_node *node, void *user)
- {
- isl_union_pw_multi_aff *upma = user;
- return isl_schedule_node_pullback_union_pw_multi_aff(node,
- isl_union_pw_multi_aff_copy(upma));
- }
- /* Compute the pullback of "schedule" by the function represented by "upma".
- * In other words, plug in "upma" in the iteration domains of "schedule".
- *
- * The schedule tree is not allowed to contain any expansion nodes.
- */
- __isl_give isl_schedule *isl_schedule_pullback_union_pw_multi_aff(
- __isl_take isl_schedule *schedule,
- __isl_take isl_union_pw_multi_aff *upma)
- {
- schedule = isl_schedule_map_schedule_node_bottom_up(schedule,
- &pullback_upma, upma);
- isl_union_pw_multi_aff_free(upma);
- return schedule;
- }
- /* Expand the schedule "schedule" by extending all leaves
- * with an expansion node with as subtree the tree of "expansion".
- * The expansion of the expansion node is determined by "contraction"
- * and the domain of "expansion". That is, the domain of "expansion"
- * is contracted according to "contraction".
- *
- * Call isl_schedule_node_expand after extracting the required
- * information from "expansion".
- */
- __isl_give isl_schedule *isl_schedule_expand(__isl_take isl_schedule *schedule,
- __isl_take isl_union_pw_multi_aff *contraction,
- __isl_take isl_schedule *expansion)
- {
- isl_union_set *domain;
- isl_schedule_node *node;
- isl_schedule_tree *tree;
- domain = isl_schedule_get_domain(expansion);
- node = isl_schedule_get_root(expansion);
- node = isl_schedule_node_child(node, 0);
- tree = isl_schedule_node_get_tree(node);
- isl_schedule_node_free(node);
- isl_schedule_free(expansion);
- node = isl_schedule_get_root(schedule);
- isl_schedule_free(schedule);
- node = isl_schedule_node_expand(node, contraction, domain, tree);
- schedule = isl_schedule_node_get_schedule(node);
- isl_schedule_node_free(node);
- return schedule;
- }
- /* Intersect the domain of the schedule "schedule" with "domain".
- * The root of "schedule" is required to be a domain node.
- */
- __isl_give isl_schedule *isl_schedule_intersect_domain(
- __isl_take isl_schedule *schedule, __isl_take isl_union_set *domain)
- {
- enum isl_schedule_node_type root_type;
- isl_schedule_node *node;
- if (!schedule || !domain)
- goto error;
- root_type = isl_schedule_tree_get_type(schedule->root);
- if (root_type != isl_schedule_node_domain)
- isl_die(isl_schedule_get_ctx(schedule), isl_error_invalid,
- "root node must be a domain node", goto error);
- node = isl_schedule_get_root(schedule);
- isl_schedule_free(schedule);
- node = isl_schedule_node_domain_intersect_domain(node, domain);
- schedule = isl_schedule_node_get_schedule(node);
- isl_schedule_node_free(node);
- return schedule;
- error:
- isl_schedule_free(schedule);
- isl_union_set_free(domain);
- return NULL;
- }
- /* Replace the domain of the schedule "schedule" with the gist
- * of the original domain with respect to the parameter domain "context".
- */
- __isl_give isl_schedule *isl_schedule_gist_domain_params(
- __isl_take isl_schedule *schedule, __isl_take isl_set *context)
- {
- enum isl_schedule_node_type root_type;
- isl_schedule_node *node;
- if (!schedule || !context)
- goto error;
- root_type = isl_schedule_tree_get_type(schedule->root);
- if (root_type != isl_schedule_node_domain)
- isl_die(isl_schedule_get_ctx(schedule), isl_error_invalid,
- "root node must be a domain node", goto error);
- node = isl_schedule_get_root(schedule);
- isl_schedule_free(schedule);
- node = isl_schedule_node_domain_gist_params(node, context);
- schedule = isl_schedule_node_get_schedule(node);
- isl_schedule_node_free(node);
- return schedule;
- error:
- isl_schedule_free(schedule);
- isl_set_free(context);
- return NULL;
- }
- /* Return an isl_union_map representation of the schedule. In particular,
- * return an isl_union_map corresponding to the subtree schedule of the child
- * of the root domain node. That is, we do not intersect the domain
- * of the returned isl_union_map with the domain constraints.
- */
- __isl_give isl_union_map *isl_schedule_get_map(__isl_keep isl_schedule *sched)
- {
- enum isl_schedule_node_type type;
- isl_schedule_node *node;
- isl_union_map *umap;
- if (!sched)
- return NULL;
- type = isl_schedule_tree_get_type(sched->root);
- if (type != isl_schedule_node_domain)
- isl_die(isl_schedule_get_ctx(sched), isl_error_internal,
- "root node not a domain node", return NULL);
- node = isl_schedule_get_root(sched);
- node = isl_schedule_node_child(node, 0);
- umap = isl_schedule_node_get_subtree_schedule_union_map(node);
- isl_schedule_node_free(node);
- return umap;
- }
- /* Insert a band node with partial schedule "partial" between the domain
- * root node of "schedule" and its single child.
- * Return a pointer to the updated schedule.
- *
- * If any of the nodes in the tree depend on the set of outer band nodes
- * then we refuse to insert the band node.
- */
- __isl_give isl_schedule *isl_schedule_insert_partial_schedule(
- __isl_take isl_schedule *schedule,
- __isl_take isl_multi_union_pw_aff *partial)
- {
- isl_schedule_node *node;
- int anchored;
- node = isl_schedule_get_root(schedule);
- isl_schedule_free(schedule);
- if (!node)
- goto error;
- if (isl_schedule_node_get_type(node) != isl_schedule_node_domain)
- isl_die(isl_schedule_node_get_ctx(node), isl_error_internal,
- "root node not a domain node", goto error);
- node = isl_schedule_node_child(node, 0);
- anchored = isl_schedule_node_is_subtree_anchored(node);
- if (anchored < 0)
- goto error;
- if (anchored)
- isl_die(isl_schedule_node_get_ctx(node), isl_error_invalid,
- "cannot insert band node in anchored subtree",
- goto error);
- node = isl_schedule_node_insert_partial_schedule(node, partial);
- schedule = isl_schedule_node_get_schedule(node);
- isl_schedule_node_free(node);
- return schedule;
- error:
- isl_schedule_node_free(node);
- isl_multi_union_pw_aff_free(partial);
- return NULL;
- }
- /* Insert a context node with constraints "context" between the domain
- * root node of "schedule" and its single child.
- * Return a pointer to the updated schedule.
- */
- __isl_give isl_schedule *isl_schedule_insert_context(
- __isl_take isl_schedule *schedule, __isl_take isl_set *context)
- {
- isl_schedule_node *node;
- node = isl_schedule_get_root(schedule);
- isl_schedule_free(schedule);
- node = isl_schedule_node_child(node, 0);
- node = isl_schedule_node_insert_context(node, context);
- schedule = isl_schedule_node_get_schedule(node);
- isl_schedule_node_free(node);
- return schedule;
- }
- /* Insert a guard node with constraints "guard" between the domain
- * root node of "schedule" and its single child.
- * Return a pointer to the updated schedule.
- */
- __isl_give isl_schedule *isl_schedule_insert_guard(
- __isl_take isl_schedule *schedule, __isl_take isl_set *guard)
- {
- isl_schedule_node *node;
- node = isl_schedule_get_root(schedule);
- isl_schedule_free(schedule);
- node = isl_schedule_node_child(node, 0);
- node = isl_schedule_node_insert_guard(node, guard);
- schedule = isl_schedule_node_get_schedule(node);
- isl_schedule_node_free(node);
- return schedule;
- }
- /* Return a tree with as top-level node a filter corresponding to "filter" and
- * as child, the (single) child of "tree".
- * However, if this single child is of type "type", then the filter is inserted
- * in the children of this single child instead.
- */
- static __isl_give isl_schedule_tree *insert_filter_in_child_of_type(
- __isl_take isl_schedule_tree *tree, __isl_take isl_union_set *filter,
- enum isl_schedule_node_type type)
- {
- if (!isl_schedule_tree_has_children(tree)) {
- isl_schedule_tree_free(tree);
- return isl_schedule_tree_from_filter(filter);
- } else {
- tree = isl_schedule_tree_child(tree, 0);
- }
- if (isl_schedule_tree_get_type(tree) == type)
- tree = isl_schedule_tree_children_insert_filter(tree, filter);
- else
- tree = isl_schedule_tree_insert_filter(tree, filter);
- return tree;
- }
- /* Construct a schedule that combines the schedules "schedule1" and "schedule2"
- * with a top-level node (underneath the domain node) of type "type",
- * either isl_schedule_node_sequence or isl_schedule_node_set.
- * The domains of the two schedules are assumed to be disjoint.
- *
- * The new schedule has as domain the union of the domains of the two
- * schedules. The child of the domain node is a node of type "type"
- * with two filters corresponding to the domains of the input schedules.
- * If one (or both) of the top-level nodes of the two schedules is itself
- * of type "type", then the filter is pushed into the children of that
- * node and the sequence or set is flattened.
- */
- __isl_give isl_schedule *isl_schedule_pair(enum isl_schedule_node_type type,
- __isl_take isl_schedule *schedule1, __isl_take isl_schedule *schedule2)
- {
- int disjoint;
- isl_ctx *ctx;
- enum isl_schedule_node_type root_type;
- isl_schedule_tree *tree1, *tree2;
- isl_union_set *filter1, *filter2, *domain;
- if (!schedule1 || !schedule2)
- goto error;
- root_type = isl_schedule_tree_get_type(schedule1->root);
- if (root_type != isl_schedule_node_domain)
- isl_die(isl_schedule_get_ctx(schedule1), isl_error_internal,
- "root node not a domain node", goto error);
- root_type = isl_schedule_tree_get_type(schedule2->root);
- if (root_type != isl_schedule_node_domain)
- isl_die(isl_schedule_get_ctx(schedule1), isl_error_internal,
- "root node not a domain node", goto error);
- ctx = isl_schedule_get_ctx(schedule1);
- tree1 = isl_schedule_tree_copy(schedule1->root);
- filter1 = isl_schedule_tree_domain_get_domain(tree1);
- tree2 = isl_schedule_tree_copy(schedule2->root);
- filter2 = isl_schedule_tree_domain_get_domain(tree2);
- isl_schedule_free(schedule1);
- isl_schedule_free(schedule2);
- disjoint = isl_union_set_is_disjoint(filter1, filter2);
- if (disjoint < 0)
- filter1 = isl_union_set_free(filter1);
- if (!disjoint)
- isl_die(ctx, isl_error_invalid,
- "schedule domains not disjoint",
- filter1 = isl_union_set_free(filter1));
- domain = isl_union_set_union(isl_union_set_copy(filter1),
- isl_union_set_copy(filter2));
- filter1 = isl_union_set_gist(filter1, isl_union_set_copy(domain));
- filter2 = isl_union_set_gist(filter2, isl_union_set_copy(domain));
- tree1 = insert_filter_in_child_of_type(tree1, filter1, type);
- tree2 = insert_filter_in_child_of_type(tree2, filter2, type);
- tree1 = isl_schedule_tree_from_pair(type, tree1, tree2);
- tree1 = isl_schedule_tree_insert_domain(tree1, domain);
- return isl_schedule_from_schedule_tree(ctx, tree1);
- error:
- isl_schedule_free(schedule1);
- isl_schedule_free(schedule2);
- return NULL;
- }
- /* Construct a schedule that combines the schedules "schedule1" and "schedule2"
- * through a sequence node.
- * The domains of the input schedules are assumed to be disjoint.
- */
- __isl_give isl_schedule *isl_schedule_sequence(
- __isl_take isl_schedule *schedule1, __isl_take isl_schedule *schedule2)
- {
- return isl_schedule_pair(isl_schedule_node_sequence,
- schedule1, schedule2);
- }
- /* Construct a schedule that combines the schedules "schedule1" and "schedule2"
- * through a set node.
- * The domains of the input schedules are assumed to be disjoint.
- */
- __isl_give isl_schedule *isl_schedule_set(
- __isl_take isl_schedule *schedule1, __isl_take isl_schedule *schedule2)
- {
- return isl_schedule_pair(isl_schedule_node_set, schedule1, schedule2);
- }
- /* Print "schedule" to "p".
- */
- __isl_give isl_printer *isl_printer_print_schedule(__isl_take isl_printer *p,
- __isl_keep isl_schedule *schedule)
- {
- if (!schedule)
- return isl_printer_free(p);
- return isl_printer_print_schedule_tree(p, schedule->root);
- }
- #undef BASE
- #define BASE schedule
- #include <print_templ_yaml.c>
|