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

isl_map_transitive_closure: coalesce after each step

Without coalescing, the number of disjuncts in the result
is exponential in the number of disjuncts in the input.

export isl_pw_qpolynomial_coalesce

isl_tab.c: unrelax: restore row if variable is non-negative

Unrelaxing a constraint means that the constraint is tightened.
The sample value may therefore become negative and then have to
restore the row again.  This is guaranteed to succeed because
the tableau was supposed to be in a valid state before the constraint
was relaxed.

isl_tab.c: cut_to_hyperplane: add extra sanity check

add isl_pw_qpolynomial_fold_size

add isl_pw_qpolynomial_split_dims

isl_qpolynomial_max: fix detection of inactive variables

Without this patch, some active set dimensions would be treated
as inactive and effectively be removed, resulting in an incorrect maximum.

doc: fix transitive closure example

isl_convex_hull.c: valid_direction: properly normalize valid direction

The original code would only normalize all but one of the coordinates
of the valid direction, possibly resulting in a scaled direction, which
may no longer be valid.  Normalize the whole vector, resulting in only
a change of representation and not a change of direction.

isl_convex_hull.c: simplify computation of initial facet constraint

The old implementation was needlessly complicated and apparently
incorrect (witness the extra test case).
The simplified version passes the new test case.

The new test case is based on a bug report for CLooG/isl
by Muthu Manikandan <muthumanikandan@gmail.com>.

export isl_pw_qpolynomial_is_zero

add isl_pw_qpolynomial_fold_dim

add isl_int_get_d

add isl_pw_qpolynomial_min

add isl_pw_qpolynomial_fix_dim

add isl_pw_qpolynomial_add_dims

export isl_dim_move

isl_pw_qpolynomial_fold_add: remove redundant quasipolynomials from folds

isl_qpolynomial_alloc: take initial struct isl_upoly as argument

add isl_set_solve_lp

isl_polynomial.c: separate out fold functionality to isl_fold.c

isl_pw_qpolynomial_{fold_,}print: set output format

improve isl_pw_qpolynomial_move

add isl_mat_move_cols

add isl_pw_qpolynomial_drop_dims

add isl_set_drop