add isl_basic_map_from_constraint

isl_map_simplify.c: break if set turns out to be empty during div coalescing

add isl_basic_map_universe_like

add isl_map_lex_lt and isl_map_lex_gt

isl_tab.c: cut_to_hyperplane: do nothing if selected constraint is an equality

reverse meaning of isl_basic_map_{less,more}_at

Before, these functions were interpreted as returning maps that map elements
to elements that are less/more.
Now, they are interpreted as returning relations that express
the concept of being less/more, i.e., the elements in the domain
are now less/more than the elements in the range.

isl_tab_dup: avoid out-of-bounds array access

tab->n_col is between 0 and mat->n_col - 2 - M
and may not be equal to n_var after some columns have been dropped.

isl 0.01

add a rudimentary manual

add some memory management annotations

export isl_basic_set_project_out

isl_basic_set_normalize_constraints: add missing return

isl_map_simplify.c: remove unused variables

isl_map.c: remove unused variables

isl_basic_set_get_hash: initialize hash value

isl_affine_hull.c: remove unused variable

isl_hash.c: remove unused variable

isl_lp.c: add missing include

isl_equalities.c: drop redundant error handling

isl_equalities.c: remove unused variable

isl_convex_hull.c: remove unused variables

isl_constraint.c: remove unused variables

isl_tab.c: sign_of_min: only pivot back if we performed any pivot

isl_tab.c: isl_tab_mark_redundant: fix up error return

isl_tab.c: remove unused variables

isl_tab_pip.c: remove unused variables

isl_vec.c: remove unused variable

isl_test.c: remove unused variable

isl_assert: validate all arguments and fix up fallout

polyhedron_sample.c: add missing include

isl_pip.c: check return value of fscanf

polyhedron_minimize.c: handle all enumeration values in switch

configure.ac: set maximal optimization compiler flags

isl_set_copy_basic_set: add missing return

isl_set_drop_basic_set: add missing return

isl_set_coalesce: add missing return

isl_coalesce.c: add missing include

isl_test: check srcdir has been set

privately export isl_tab_kill_col

isl_tab.c: fix up to_col

isl_tab_pip.c: add missing include

isl_vec.c: add missing include

export isl_seq_cmp

fix return type of isl_constraint_free

export isl_token_free and isl_stream_error

isl_list.c: add missing include

isl_div.c: add missing include

isl_constraint.c: add missing include

basis_reduction_tab.c: add missing include

isl_map_no_piplib.c: add missing include

ax_create_stdint_h.m4: protect some limits from redefinition

polytope_scan.c: add missing include

fix return type of isl_div_free

isl_map_simplify.c: drop return from void function

isl_map_simplify.c: add missing include

isl_convex_hull.c: use isl_seq_get_hash instead of isl_seq_hash

declare flexarrays of size 1 to silence sun compiler

change isl_basic_map_empty interface for consistency

AX_SUBMODULE: mention default type of library to use

configure.ac: use AX_SUBMODULE for gmp for consistency

separate out config header from isl_ctx.h

pilp solver: don't ignore feasibility test on context

In the initial implementation, the sample point of the context was
kept integer.  This was later changed in such a way that now the
sample point is just lexico-smallest rational value.
This means that after checking for an integer point, we rollback
until the state before we performed this check.
However, the results of the feasibility test were only checked
after the rollback, completely ignoring the actual results of
the feasibility test.

isl_basic_map_gauss: try not to remove any div definitions

After computing divs, the result may get simplified and during
this simplification, isl_basic_map_gauss may actually remove some
divs because it is trying to play safe wrt circular div definitions.
By keeping the divs ordered, divs can safely be eliminated without
introducing any circular definitions.

add isl_map_sum

add isl_map_floordiv

add isl_map_neg

add isl_map_is_strict_subset

add isl_map_fast_is_fixed

add isl_map_identity_like

rename isl_map_identity_like to isl_map_identity_like_basic_map

add isl_polytope_scan application

isl_seq_normalize: use pre-allocated temporary variable in isl_ctx

This requires that we pass along an isl_ctx, but it removes
the need for allocating a temporary variable every time
isl_seq_normalize is called.
In some cases, up to 20% of the execution time of isl_seq_normalize
was wasted on initializing and clearing this temporary variable.

add isl_vec_normalize

isl_seq_normalize: no need to scale down by one

add isl_polyhedron_minimize application

add generalized basis reduction based ILP solver

add isl_vec_mat_product

export isl_vec_ceil

isl_basic_set_sample: only perform basis reduction once

On larger problem, much time was spent on repeated basis reductions
with little or no reduction on the number of points visited.
The initial basis reduction seems to be sufficient in most cases.

rename isl_solve_lp to isl_basic_{map,set}_solve_lp

isl_solve_lp: optionally return solution point

The caller may now not be interested in the optimal value anymore,
so make returning this optimal value optional.

export isl_vec header

add isl_int_fdiv_q_ui

add isl_vec_read_from_file

isl_basic_map_drop_redundant_divs: also investigate divs that have a definition

They may have become redundant after having been defined.

isl_basic_map_apply_{domain,range}: drop redundant divs in result

isl_map_simplify.c: fix typo in comment

isl_basic_set_project_out: drop redundant divs in result

isl_basic_map_remove: only drop divs if basic map did not turn out to be empty

If isl_basic_map_eliminate_vars finds that its argument is empty,
it resets its constraints to the canonical representation of an empty
set and drops all divs.

isl_basic_map_detect_equalities: explicitly keep track of any equalities found

The way extend_affine_hull works is that it starts off from an
under-approximation of the affine hull and then extends it with
each extra point found outside of the current approximation.
The affine hull is then determined by the equalities that remain.

However, during the computation, we may find an equality in
the current approximation such that there is no point in the set
outside of this equality.  If the approximation is later extended
with a point found outside some other equality, then this first
equality would get tested again.  If this is a difficult computation,
then we needlessly go through it again.

By explicitly storing the equality that we have detected in the set,
this computation becomes trivial, to there is no real harm in doing
again.  An alternative would be to keep track of all equalities found
in some additional data structure, but adding it to the original set
seems preferable since it may also simplify tests for points outside
other equalities.

Perhaps the calling convention should be changed such that
extend_affine_hull would return the original set intersected
with its affine hull, rather than the affine hull itself.

isl_basic_map_extend_dim: keep hold of sample if dimension doesn't change

isl_tab_min: read off all information from tableau before rollback

As it happens, the rollback has little effect in this case and
the information is still available, but we shouldn't depend on this.

isl_tab_add_valid_eq: keep track of whether equality is negated

This is important during generalized basis reduction, because
this procedure requires the dual variables of some of the
equalities added to the tableau.  If the original equality get
negated then the dual of the original equality will be the
opposite of the dual of the equality that was actually added.

Without this fix, gbr could in some cases make completely
inappropriate decisions.

configure.ac: no longer use piplib by default

add an internal parametric integer linear program solver

This removes the dependence on piplib.

The solver is very similar to the one in piplib.
In particular, as in piplib, the lexicographic solver is
used recursively on the context.  This means that the current
solver will run into the same problems the piplib solver runs into
while looking for an integer point in a context that no longer has any.

In the future, we will replace this recursive call by a call
to an incremental solver based on the generalized basis reduction solver.

isl_mat_extend: make sure the number of rows never decreases

isl_affine_hull.c: only construct affine hull in bounded directions

The affine hull of a set necessarily contains the directions
of the recession cone of the set.
Before we exploited this fact by adding these directions to
our initial approximation of the affine hull.
The set for which the affine hull was computed was not changed, however.
Now, we project out the recession cone and only construct an
affine hull in the remaining part.

isl_tab: add support for keeping track of samples

isl_tab: optionally keep track of row signs

isl_tab: add isl_basic_set field for optionally keeping track of inequalities