mark some functions as requiring use of return value

isl_tab: improved error handling

isl_tab_pip.c: incrementally build recession cone of gbr context

When checking the feasibility of a constraint in the gbr context,
we compute a sample.  The first step of this computation builds
a recession cone since we need to know which directions are unbounded.
By incrementally building up the recession cone, we can avoid
the complete recomputation.  This is especially useful when
the context is/becomes bounded, since we can then avoid the
recession cone computation altogether.

isl_tab_pip.c: propagate some equalities from gbr context to main tableau

If a new equality in the context can be used to eliminate
a parameter in the main tableau, we perform the elimination
as this allows us to more easily determine the signs of rows
and it may avoid the introduction of some redundant integer
divisions.

isl_tab_pip.c: detect equalities in gbr context on each parametric cut

A new div that was introduced to express a parametric cut may,
within the current context, be an affine combination of other divs.
If so, we want to find these equalities as soon as possible
as this may allow us to more easily determine feasibility of
a constraint.

add isl_tab_detect_equalities

add isl_tab_set_initial_basis_with_cone

isl_tab_pip: use generalized basis reduction based context by default

isl_tab_pip: add generalized basis reduction based context handling

isl_tab_pip.c: extract out context handling

isl_tab_pip.c: add_div: use more reliable way to test whether div is nonneg

Before, we used the fact that if all parameters are non-negative,
then the big parameter is not used.
However, it is perfectly valid to have negative parameters and
still not use any big parameter.

isl_tab_basic_map_partial_lexopt: remove samples that are no longer useful

When we test an inequality in row_sign, it is useful to store any
samples found, so we can use them to speed up the computation
of future row signs.
However, when we return from a branch, then all the samples found
since we entered the branch all satisfy the constraints on that
branch and so they will not satisfy the constraints on any other
branch.  We therefore save the current number of samples when
we enter a branch and then discard all subsequently added samples
when we return from the branch.

isl_tab.c: extract out samples handling from isl_tab_pip.c

isl_affine_hull.c: uset_affine_hull_bounded: use tableaus to perform computation

isl_basic_map_detect_equalities: keep track of sample

isl_basic_map_detect_equalities: only compute recession cone once

Before, the recession cone was computed both during the computation
of the initial sample and then later again to project out this
recession cone.

isl_tab_sample: handle unbounded directions in initial basis

isl_tab_compute_reduced_basis: handle unbounded directions in initial basis

exploit equalities in isl_tab_sample

isl_tab: keep (in)equalities of bset (if any) in sync

isl_tab_sample: be more verbose about unbounded directions

separate out isl_tab_sample from sample_bounded

sample_bounded: reimplement to work directly on a tableau

In other words, adapt the implementation of isl_basic_set_samples
in polytope_scan.c
This avoids the reconstruction of tableaus and allows for a refactoring
that takes a tableau as input.

isl_tab_compute_reduced_basis: work with affine basis instead of linear basis

isl_tab_compute_reduced_basis: allow incremental computation

isl_polytope_scan: use isl_tab_from_basic_set

separate out isl_tab_reduced_basis from isl_basic_set_reduced_basis

This allows us to compute a reduced basis on a set given in tableau form.

basis_reduction_tab.c: use isl_tab_product to construct product tableau

Before, we would manually construct a tableau for the Cartesian
product of a set with itself.  Now, we first create a tableau
for the set and then construct the product tableau.
This should be more efficient and allows us to refactor
the code to pass in a tableau instead of a set.

add isl_tab_product

isl_basic_set_reduced_basis: fix value in directions with only one integer value

If two or more directions have a width that is smaller than one,
i.e., such that the only one integer value can be attained,
then the standard algorithm would still try to order these
directions.  However, if the first directions only admit one
value, then it doesn't matter in which order they are put.
We therefore keep track of such directions explicitly and
also set the value in the tableau to the fixed value, rather
than keeping the rational interval with only one integer value.

isl_basic_set_reduced_basis: fix up documentation

basis_reduction_tab.c: keep track of con_offset instead of n_ineq

This is a minor cosmetic change.

add isl_tab_add_eq

isl_tab_add_valid_eq: add special treatment for manifestly zero rows

isl_tab_extend_cons: check tab argument

rename isl_tab_detect_equalities to isl_tab_detect_implicit_equalities

isl_tab_from_recession_cone: take basic set instead of basic map as argument

isl_ilp.c: separate out solve_ilp_search

isl_sample.c: basic_set_sample: remember boundedness of basic set

add isl_polyhedron_detect_equalities test application

add isl_basic_set_detect_equalities

isl_tab_allocate_con: add extra assertion

isl_sample.c: basic_set_reduced: fix typo preventing use of gbr_only_first

isl_tab_basic_map_partial_lexopt: simplify result

privately export isl_basic_set_add_{in,}eq

privately export isl_basic_set_sample_with_cone

add isl_basic_set_underlying_set

isl_tab.c: close_row: push undo record for setting row to zero

add isl_mat_vec_inverse_product

add isl_mat_concat

add isl_int_divexact_ui

add isl_set_dim_residue_class

fix long standing bug in isl_mat_inverse_product

add isl_basic_map_foreach_lexmin

add isl_set_project_out

add isl_set_detect_equalities

add isl_set_foreach_basic_set

add isl_set_lifting

isl_dim_size: check argument

export isl_mat header

add isl_basic_map_lexmin

add isl_map_foreach_basic_map

configure.ac: fix cut-and-paste error in original commit

add isl_int_get_si

export isl_set_sample

add isl_basic_set_universe_like_set

add isl_set_is_strict_subset

add isl_set_universe_like

AX_SUBMODULE: set PKG_CONFIG_PATH

create pkg-config file

AX_SUBMODULE: drop options that are meaningless given possible choices

isl_basic_map_from_constraint: return copy of bmap in constraint if possible

Commit fa1cd4b (isl_basic_map_add_constraint: handle constraints obtained from
basic maps) fixed isl_basic_map_add_constraint to also work with constraints
obtained from iterating over the constraints of a basic map by
calling isl_basic_map_from_constraint.
However, this broke the case where a new constraint is created with
existentially quantified variables.  In this case, the divs may be
unknown and we can't compute them because that may result in a union.
So, try to detect this special case in isl_basic_map_from_constraint
and return the basic map stored inside the constraint.

doc: fix description of lexicograhpic order relations

add isl_map_lex_le and isl_map_lex_ge

export isl_basic_map_gist

add isl_basic_map_sample and isl_map_sample

isl_sample.c: move isl_basic_set_from_vec from isl_affine_hull.c

rename isl_basic_set_sample to isl_basic_set_sample_vec

make some internal functions static

isl_basic_map_add_constraint: handle constraints obtained from basic maps

The bmap field of a constraint may contain a basic map that is bounded
by more than one constraint.  The constraint can therefore not be
added by simply taking the intersection with this bmap.
Instead, first create a basic map that exactly corresponds to the constraint
and then take the intersection with that basic map.

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