add isl_basic_set_foreach_bound_pair

add isl_qpolynomial_as_polynomial_on_domain

add isl_qpolynomial_substitute

add isl_basic_set_dims_get_sign

isl_ctx: add zero

add isl_qpolynomial_neginfty

add isl_basic_set_drop

isl_qpolynomial_fold_opt_on_domain: fix typo

isl_pw_*_eval: fix handling of empty input

add isl_basic_set_move_dims

export isl_qpolynomial_sub

add isl_qpolynomial_get_dim

isl_map_transitive_closure: coalesce input before checking closedness of input

add map_lex_*_first

isl_map_lex_le: use one basic set less in result

isl_stream: accept lexicographic operators

isl_map_transitive_closure: check if input is closed already

add isl_printer_flush

isl_map_read: parse parametric Omega input

isl_map_read: fix typo

isl_convex_hull.c: is_bound: any inequality is a bound on an empty set

isl_convex_hull.c: add_bounds: propagate errors from add_bound

document isl_map_reverse

isl_stream_next_token: terminate string tokens

isl_map_transitive_closure: treat equality as pair of ineq for mixed constraints

isl_map_transitive_closure, omega-like: only use ?-closure if result is exact

isl_map_transitive_closure: use simple hull in omega-like implementation

isl_map_transitive_closure: try and compute closure incrementally

isl_map_coalesce: allow wrapping in sets that stick out in different directions

Before, we would only wrap in sets that extend beyond exactly
one facet of a given set and such that the wrapped in set itself
has a facet that is adjacent to this facet.
We now extend this process in several ways.
First, the wrapped-in set is allowed to extend beyond several facets.
Second, the wrapped-in set need not have facets adjacent to these facets.
It is sufficient that the adjacent facets are valid for the wrapped-in set.
Finally, the wrapped-in set is allowed to be of a lower dimension.

isl_map_coalesce: allowing cutting equalities in adjacent facet

From the start, we have been very careful not to allow any cutting
equalities anywhere.  However, it makes sense to wrap in a piece
of a hyperplane adjacent to a facet, even if this piece is
lower-dimensional.

add some more coalescing test cases

isl_basic_map_gist: be more aggressive in removing constraints

In particular, detect equalities in intersection of input and context
and check for each constraint if negating it would lead to a contradiction.

isl_stream: accept string tokens

isl_tab_rollback: restore rows that are no longer marked redundant

When a row is marked redundant, it may have a negative (but strictly
larger than -1) sample value.  When undoing this marking, we need
to make sure the sample value is restored to a non-negative value.

doc: document simple hull

doc: emphasize that we are dealing with integer sets

isl_transitive_closure.c: path_along_delta: only avoid splitting off identity

Commit 39345fc (isl_transitive_closure.c: path_along_delta: try and avoid
splitting off identity) attempted to avoid splitting off the identity mapping.
However, it would only check if relaxing the length to zero would
not introduce any extra elements.  The path may have constraints
on the parameters, though, and so we may actually be removing elements.
This means that we need to check for equality with the identity relation.
Perhaps this makes the whole optimization pointless.

isl_map_transitive_closure: compute divs of input

isl_map_compute_divs: make sure result is gaussed and finalized

In particular, when the variables in the input can be compressed,
basic_set_append_equalities should gauss and finalize the result
of adding the equalities back in.

isl_printer_print_pw_qpolynomial: properly print constraints in C format

In particular, use an explicit multiplication in affine expressions.

isl_pip and isl_cat: terminate isl_arg array with ISL_ARG_END

On most systems, leaving it out apparently just happens to work,
but not so on MacOS X.

Problem reported and analyzed by Tomofumi Yuki <yuki@cs.colostate.edu>

doc: document setting and getting of dimension names

doc: document *_get_dim functions

isl_map_convex_hull: handle non full-dimensional pairs of basic sets

isl_convex_hull.c: modulo_affine_hull: drop redundant argument

isl_convex_hull.c: convex_hull_pair: skip direction change if result is bounded

isl_printer_print_basic_set: fix up printing to string printer

In particular, only print prefix and suffix (and call strlen on them)
if they are not NULL.

Problem reported by Tomofumi Yuki <yuki@cs.colostate.edu>

add isl_set_neg and isl_set_sum

isl_transitive_closure.c: fix documentation of construct_component

isl_map_coalesce: make sure result still contains set after wrapping in a facet

isl_map_transitive_closure: Floyd-Warshall: properly keep track of current group

isl_map_transitive_closure: use Floyd-Warshall on disjoint domains and ranges

isl_transitive_closure.c: extract out isl_set_overlaps

isl_map_transitive_closure: project out length on each component separately

There is no need to keep track of the length over multiple components.

isl_transitive_closure.c: construct_component: don't include identity in result

isl_map_transitive_closure: don't add length parameter

Just project out the final coordinates, intead of equating the difference
to a parameter and then projecting out the difference.

doc: fix typos

isl_set_gist: make context argument an isl_set for consistency

isl_map_transitive_closure: coalesce input map

isl_transitive_closure.c: construct_component: always include identity

Commit e06f258 (isl_map_transitive_closure: break early if input map
doesn't compose with itself) added a special case to construct_component,
but neglected to include the identity map in the result.

doc: cite Omega library documentation on computation of underapproximations

add isl_stream_read_map

isl_arg.h: add missing include

isl_basic_map_move_dims: normalize equalities

rename isl_map_move to isl_map_move_dims

isl_options_parse: only parse known options by default

That is, only barf on unknown options when the user explicitly
asks for this behavior.

isl_map_closure: optionally use Omega-like algorithm

isl_closure: take options

isl_map_affine_hull: handle non-obviously empty input

In particular, detect equalities before checking whether the input is empty.
Otherwise, we could be accessing the first element of an empty list.

isl_map_compute_divs: handle non-obviously empty input

In particular, properly handle the case where variable compression
detects that the input map is actually empty.

isl_basic_set_preimage: add extra sanity check

isl_basic_map_set_to_empty: add missing assignment

isl_equalities.c: isl_mat_variable_compression: clarify result

isl_sample.c: initial_basis: set n_unbounded and n_zero

n_unbounded and n_zero are initialized to zero, but they can
be set to non-zero values in isl_tab_set_initial_basis_with_cone.
If isl_tab_sample is later called on an updated tab, but with
the old basis removed, then these values, in particular n_unbounded,
may not be valid anymore.
In particular, this happens in isl_tab_pip.c's gbr_get_sample.
Arguably, gbr_get_sample should reset the values of n_zero and n_unbounded
when it drops the basis, but it doens't hurt for initial_basis
to reset those values too.

isl_map_coalesce: only use non-redundant constraints during wrapping

Constraints that have been marked redundant are ignored and may
no longer be valid, so they shouldn't be taken into account during wrapping.

isl_map_print: allow empty constraints

A basic map that is passed to the user should never contain a constraint
"0 >= 0".  However, internally, such constraints can exist for a while
and for debugging purposes, it's convenient to be able to print
basic maps with such constraints.

isl_options.h: allow inclusion from C++

isl_arg.h: allow inclusion from C++

isl_map_transitive_closure: handle existentials

isl_map_deltas: preserve dimension names

isl_map_print: support LaTeX output format

isl_map_simplify: don't remove any div definitions if all divs are known

isl_map_simplify is called from within compute_divs, so we need
to make sure it doesn't remove any divs just computed.
Unfortunately, being more strict about not removing any div definition
means that we break some test cases, so they are disabled now.
Perhaps we can think of a better way of making the tests pass later.

isl_map_intersect: add special cases for empty input maps

In particular, if both input maps consist of a single basic map
and that of the first map turns out to be empty, then we need
to take into account that after cowing the map, it may not
have any basic maps anymore.

isl_transitive_closure.c: path_along_delta: share code for handling {in,}equalities

isl_transitive_closure.c: path_along_delta: try and avoid splitting off identity

Usually, path_along_delta returns a union of the identity mapping and
of paths of length at least 1.  In some cases, the paths of length
0 correspond exactly to the identity mapping and then it is more
efficient to return a single mapping with paths of length at least 0.

export isl_qpolynomial_move_dims

add isl_qpolynomial_get_ctx

add isl_pw_qpolynomial_fold_get_dim

add some qpolynomial test functions

add isl_pw_qpolynomial_fold_{max,min}

isl_tab_min_at_most_neg_one: restore sample value of non-redundant constraints

isl_map_coalesce: avoid reconsidering pairs considered before on change

Before, whenever anything changed (two elements get fused or one
element gets dropped), we would start all over from scratch.
Pairs that were already considered before and that are not related
to the change do not need to be reconsidered.
We may still reconsider some pairs, even after this commit.
In particular, if one element gets dropped, then we will still reconsider
all combinations with the element that did not get dropped.

isl_tab.c: update debugging aid check_table to presence of big parameter

isl_tab_relax: make sure no non-negative rows get a negative sample value

isl_tab_pip.c: always set row sign of found row to isl_tab_row_neg

first_neg returns row that is obviously negative and usually
sets the sign of the row to isl_tab_row_neg.  However, if
the row happens to be obviously negative because the big parameter
has a negative coefficient, then it would fail to set the row sign,
while update_row_sign check that the row sign is effectively set
to isl_tab_row_neg.

isl_printer_print_basic_{set,map}: add closing "}"

They were removed by accident in commit 5aef299 (add isl_printer).

export isl_set_split_dims

add isl_set_recession_cone

add isl_set_dim_is_bounded

isl_map_transitive_closure: break early if input map doesn't compose with itself