add isl_basic_set_max
[platform/upstream/isl.git] / doc / user.pod
index 47c720c..e3277a3 100644 (file)
@@ -73,6 +73,20 @@ of a B<map> as input.  An old call
 C<isl_map_identity(dim)> can be rewritten to
 C<isl_map_identity(isl_dim_map_from_set(dim))>.
 
+=item * The function C<isl_map_power> no longer takes
+a parameter position as input.  Instead, the exponent
+is now expressed as the domain of the resulting relation.
+
+=back
+
+=head3 Changes since isl-0.06
+
+=over
+
+=item * The format of C<isl_printer_print_qpolynomial>'s
+C<ISL_FORMAT_ISL> output has changed.
+Use C<ISL_FORMAT_C> to obtain the old output.
+
 =back
 
 =head1 Installation
@@ -357,7 +371,7 @@ will also release all the objects passed as arguments.
 If the user still wants to use one or more of these arguments
 after the function call, she should pass along a copy of the
 object rather than the object itself.
-The user is then responsible for make sure that the original
+The user is then responsible for making sure that the original
 object gets used somewhere else or is explicitly freed.
 
 The arguments and return values of all documents functions are
@@ -452,6 +466,10 @@ specification of the original object.
        __isl_give isl_dim *isl_union_map_get_dim(
                __isl_keep isl_union_map *umap);
 
+       #include <isl/constraint.h>
+       __isl_give isl_dim *isl_constraint_get_dim(
+               __isl_keep isl_constraint *constraint);
+
        #include <isl/polynomial.h>
        __isl_give isl_dim *isl_qpolynomial_get_dim(
                __isl_keep isl_qpolynomial *qp);
@@ -462,6 +480,10 @@ specification of the original object.
        __isl_give isl_dim *isl_union_pw_qpolynomial_fold_get_dim(
                __isl_keep isl_union_pw_qpolynomial_fold *upwf);
 
+       #include <isl/aff.h>
+       __isl_give isl_dim *isl_aff_get_dim(
+               __isl_keep isl_aff *aff);
+
 The names of the individual dimensions may be set or read off
 using the following functions.
 
@@ -536,10 +558,50 @@ specifications using the following functions.
                enum isl_dim_type type, unsigned first, unsigned n);
        __isl_give isl_dim *isl_dim_map_from_set(
                __isl_take isl_dim *dim);
+       __isl_give isl_dim *isl_dim_zip(__isl_take isl_dim *dim);
 
 Note that if dimensions are added or removed from a space, then
 the name and the internal structure are lost.
 
+=head2 Local Spaces
+
+A local space is essentially a dimension specification with
+zero or more existentially quantified variables.
+The local space of a basic set or relation can be obtained
+using the following functions.
+
+       #include <isl/set.h>
+       __isl_give isl_local_space *isl_basic_set_get_local_space(
+               __isl_keep isl_basic_set *bset);
+
+       #include <isl/map.h>
+       __isl_give isl_local_space *isl_basic_map_get_local_space(
+               __isl_keep isl_basic_map *bmap);
+
+A new local space can be created from a dimension specification using
+
+       #include <isl/local_space.h>
+       __isl_give isl_local_space *isl_local_space_from_dim(
+               __isl_take isl_dim *dim);
+
+They can be inspected, copied and freed using the following functions.
+
+       #include <isl/local_space.h>
+       isl_ctx *isl_local_space_get_ctx(
+               __isl_keep isl_local_space *ls);
+       int isl_local_space_dim(__isl_keep isl_local_space *ls,
+               enum isl_dim_type type);
+       const char *isl_local_space_get_dim_name(
+               __isl_keep isl_local_space *ls,
+               enum isl_dim_type type, unsigned pos);
+       __isl_give isl_dim *isl_local_space_get_dim(
+               __isl_keep isl_local_space *ls);
+       __isl_give isl_div *isl_local_space_get_div(
+               __isl_keep isl_local_space *ls, int pos);
+       __isl_give isl_local_space *isl_local_space_copy(
+               __isl_keep isl_local_space *ls);
+       void *isl_local_space_free(__isl_take isl_local_space *ls);
+
 =head2 Input and Output
 
 C<isl> supports its own input/output format, which is similar
@@ -743,6 +805,10 @@ is only used to specify the parameters.
                __isl_take isl_dim *dim);
        __isl_give isl_map *isl_map_universe(
                __isl_take isl_dim *dim);
+       __isl_give isl_union_set *isl_union_set_universe(
+               __isl_take isl_union_set *uset);
+       __isl_give isl_union_map *isl_union_map_universe(
+               __isl_take isl_union_map *umap);
 
 The sets and relations constructed by the functions above
 contain all integer values, while those constructed by the
@@ -842,7 +908,7 @@ from a universe set or relation, adding equality and/or
 inequality constraints and then projecting out the
 existentially quantified variables, if any.
 Constraints can be constructed, manipulated and
-added to basic sets and relations using the following functions.
+added to (basic) sets and relations using the following functions.
 
        #include <isl/constraint.h>
        __isl_give isl_constraint *isl_equality_alloc(
@@ -860,6 +926,12 @@ added to basic sets and relations using the following functions.
        __isl_give isl_basic_set *isl_basic_set_add_constraint(
                __isl_take isl_basic_set *bset,
                __isl_take isl_constraint *constraint);
+       __isl_give isl_map *isl_map_add_constraint(
+               __isl_take isl_map *map,
+               __isl_take isl_constraint *constraint);
+       __isl_give isl_set *isl_set_add_constraint(
+               __isl_take isl_set *set,
+               __isl_take isl_constraint *constraint);
 
 For example, to create a set containing the even integers
 between 10 and 42, you would use the following code.
@@ -970,6 +1042,8 @@ using the following functions, which compute an overapproximation.
                __isl_take isl_basic_map *bmap);
        __isl_give isl_set *isl_set_remove_divs(
                __isl_take isl_set *set);
+       __isl_give isl_map *isl_map_remove_divs(
+               __isl_take isl_map *map);
 
 To iterate over all the sets or maps in a union set or map, use
 
@@ -1052,6 +1126,9 @@ the following functions.
        void isl_constraint_get_coefficient(
                __isl_keep isl_constraint *constraint,
                enum isl_dim_type type, int pos, isl_int *v);
+       int isl_constraint_involves_dims(
+               __isl_keep isl_constraint *constraint,
+               enum isl_dim_type type, unsigned first, unsigned n);
 
 The explicit representations of the existentially quantified
 variables can be inspected using the following functions.
@@ -1061,6 +1138,7 @@ to C<isl_set_compute_divs> or C<isl_map_compute_divs>.
 
        __isl_give isl_div *isl_constraint_div(
                __isl_keep isl_constraint *constraint, int pos);
+       isl_ctx *isl_div_get_ctx(__isl_keep isl_div *div);
        void isl_div_get_constant(__isl_keep isl_div *div,
                isl_int *v);
        void isl_div_get_denominator(__isl_keep isl_div *div,
@@ -1117,6 +1195,9 @@ the following functions.
        const char *isl_constraint_get_dim_name(
                __isl_keep isl_constraint *constraint,
                enum isl_dim_type type, unsigned pos);
+       const char *isl_basic_set_get_dim_name(
+               __isl_keep isl_basic_set *bset,
+               enum isl_dim_type type, unsigned pos);
        const char *isl_set_get_dim_name(
                __isl_keep isl_set *set,
                enum isl_dim_type type, unsigned pos);
@@ -1139,17 +1220,18 @@ of the parameters.
 =item * Emptiness
 
 The following functions test whether the given set or relation
-contains any integer points.  The ``fast'' variants do not perform
+contains any integer points.  The ``plain'' variants do not perform
 any computations, but simply check if the given set or relation
 is already known to be empty.
 
-       int isl_basic_set_fast_is_empty(__isl_keep isl_basic_set *bset);
+       int isl_basic_set_plain_is_empty(__isl_keep isl_basic_set *bset);
        int isl_basic_set_is_empty(__isl_keep isl_basic_set *bset);
+       int isl_set_plain_is_empty(__isl_keep isl_set *set);
        int isl_set_is_empty(__isl_keep isl_set *set);
        int isl_union_set_is_empty(__isl_keep isl_union_set *uset);
-       int isl_basic_map_fast_is_empty(__isl_keep isl_basic_map *bmap);
+       int isl_basic_map_plain_is_empty(__isl_keep isl_basic_map *bmap);
        int isl_basic_map_is_empty(__isl_keep isl_basic_map *bmap);
-       int isl_map_fast_is_empty(__isl_keep isl_map *map);
+       int isl_map_plain_is_empty(__isl_keep isl_map *map);
        int isl_map_is_empty(__isl_keep isl_map *map);
        int isl_union_map_is_empty(__isl_keep isl_union_map *umap);
 
@@ -1157,25 +1239,46 @@ is already known to be empty.
 
        int isl_basic_set_is_universe(__isl_keep isl_basic_set *bset);
        int isl_basic_map_is_universe(__isl_keep isl_basic_map *bmap);
-       int isl_set_fast_is_universe(__isl_keep isl_set *set);
+       int isl_set_plain_is_universe(__isl_keep isl_set *set);
 
 =item * Single-valuedness
 
        int isl_map_is_single_valued(__isl_keep isl_map *map);
+       int isl_union_map_is_single_valued(__isl_keep isl_union_map *umap);
+
+=item * Injectivity
+
+       int isl_map_plain_is_injective(__isl_keep isl_map *map);
+       int isl_map_is_injective(__isl_keep isl_map *map);
+       int isl_union_map_plain_is_injective(
+               __isl_keep isl_union_map *umap);
+       int isl_union_map_is_injective(
+               __isl_keep isl_union_map *umap);
 
 =item * Bijectivity
 
        int isl_map_is_bijective(__isl_keep isl_map *map);
+       int isl_union_map_is_bijective(__isl_keep isl_union_map *umap);
 
 =item * Wrapping
 
-The followning functions check whether the domain of the given
+The following functions check whether the domain of the given
 (basic) set is a wrapped relation.
 
        int isl_basic_set_is_wrapping(
                __isl_keep isl_basic_set *bset);
        int isl_set_is_wrapping(__isl_keep isl_set *set);
 
+=item * Internal Product
+
+       int isl_basic_map_can_zip(
+               __isl_keep isl_basic_map *bmap);
+       int isl_map_can_zip(__isl_keep isl_map *map);
+
+Check whether the product of domain and range of the given relation
+can be computed,
+i.e., whether both domain and range are nested relations.
+
 =back
 
 =head3 Binary Properties
@@ -1184,7 +1287,7 @@ The followning functions check whether the domain of the given
 
 =item * Equality
 
-       int isl_set_fast_is_equal(__isl_keep isl_set *set1,
+       int isl_set_plain_is_equal(__isl_keep isl_set *set1,
                __isl_keep isl_set *set2);
        int isl_set_is_equal(__isl_keep isl_set *set1,
                __isl_keep isl_set *set2);
@@ -1196,7 +1299,7 @@ The followning functions check whether the domain of the given
                __isl_keep isl_basic_map *bmap2);
        int isl_map_is_equal(__isl_keep isl_map *map1,
                __isl_keep isl_map *map2);
-       int isl_map_fast_is_equal(__isl_keep isl_map *map1,
+       int isl_map_plain_is_equal(__isl_keep isl_map *map1,
                __isl_keep isl_map *map2);
        int isl_union_map_is_equal(
                __isl_keep isl_union_map *umap1,
@@ -1204,7 +1307,7 @@ The followning functions check whether the domain of the given
 
 =item * Disjointness
 
-       int isl_set_fast_is_disjoint(__isl_keep isl_set *set1,
+       int isl_set_plain_is_disjoint(__isl_keep isl_set *set1,
                __isl_keep isl_set *set2);
 
 =item * Subset
@@ -1298,6 +1401,15 @@ The followning functions check whether the domain of the given
 The functions above construct a (basic, regular or union) relation
 that maps (a wrapped version of) the input relation to its domain or range.
 
+=item * Elimination
+
+       __isl_give isl_set *isl_set_eliminate(
+               __isl_take isl_set *set, enum isl_dim_type type,
+               unsigned first, unsigned n);
+
+Eliminate the coefficients for the given dimensions from the constraints,
+without removing the dimensions.
+
 =item * Identity
 
        __isl_give isl_map *isl_set_identity(
@@ -1318,6 +1430,16 @@ Construct an identity relation on the given (union) set.
 These functions return a (basic) set containing the differences
 between image elements and corresponding domain elements in the input.
 
+       __isl_give isl_basic_map *isl_basic_map_deltas_map(
+               __isl_take isl_basic_map *bmap);
+       __isl_give isl_map *isl_map_deltas_map(
+               __isl_take isl_map *map);
+       __isl_give isl_union_map *isl_union_map_deltas_map(
+               __isl_take isl_union_map *umap);
+
+The functions above construct a (basic, regular or union) relation
+that maps (a wrapped version of) the input relation to its delta set.
+
 =item * Coalescing
 
 Simplify the representation of a set or relation by trying
@@ -1349,6 +1471,13 @@ basic set or relation.
 Simplify the representation of a set or relation by detecting implicit
 equalities.
 
+=item * Removing redundant constraints
+
+       __isl_give isl_basic_set *isl_basic_set_remove_redundancies(
+               __isl_take isl_basic_set *bset);
+       __isl_give isl_basic_map *isl_basic_map_remove_redundancies(
+               __isl_take isl_basic_map *bmap);
+
 =item * Convex hull
 
        __isl_give isl_basic_set *isl_set_convex_hull(
@@ -1414,17 +1543,58 @@ that contains the whole input set or relation.
 In case of union sets and relations, the polyhedral hull is computed
 per space.
 
+=item * Optimization
+
+       #include <isl/ilp.h>
+       enum isl_lp_result isl_basic_set_max(
+               __isl_keep isl_basic_set *bset,
+               __isl_keep isl_aff *obj, isl_int *opt)
+       enum isl_lp_result isl_set_max(__isl_keep isl_set *set,
+               __isl_keep isl_aff *obj, isl_int *opt);
+
+Compute the maximum of the integer affine expression C<obj>
+over the points in C<set>, returning the result in C<opt>.
+The return value may be one of C<isl_lp_error>,
+C<isl_lp_ok>, C<isl_lp_unbounded> or C<isl_lp_empty>.
+
+=item * Dual
+
+The following functions compute either the set of (rational) coefficient
+values of valid constraints for the given set or the set of (rational)
+values satisfying the constraints with coefficients from the given set.
+Internally, these two sets of functions perform essentially the
+same operations, except that the set of coefficients is assumed to
+be a cone, while the set of values may be any polyhedron.
+The current implementation is based on the Farkas lemma and
+Fourier-Motzkin elimination, but this may change or be made optional
+in future.  In particular, future implementations may use different
+dualization algorithms or skip the elimination step.
+
+       __isl_give isl_basic_set *isl_basic_set_coefficients(
+               __isl_take isl_basic_set *bset);
+       __isl_give isl_basic_set *isl_set_coefficients(
+               __isl_take isl_set *set);
+       __isl_give isl_union_set *isl_union_set_coefficients(
+               __isl_take isl_union_set *bset);
+       __isl_give isl_basic_set *isl_basic_set_solutions(
+               __isl_take isl_basic_set *bset);
+       __isl_give isl_basic_set *isl_set_solutions(
+               __isl_take isl_set *set);
+       __isl_give isl_union_set *isl_union_set_solutions(
+               __isl_take isl_union_set *bset);
+
 =item * Power
 
        __isl_give isl_map *isl_map_power(__isl_take isl_map *map,
-               unsigned param, int *exact);
+               int *exact);
+       __isl_give isl_union_map *isl_union_map_power(
+               __isl_take isl_union_map *umap, int *exact);
 
 Compute a parametric representation for all positive powers I<k> of C<map>.
-The power I<k> is equated to the parameter at position C<param>.
-The result may be an overapproximation.  If the result is exact,
+The result maps I<k> to a nested relation corresponding to the
+I<k>th power of C<map>.
+The result may be an overapproximation.  If the result is known to be exact,
 then C<*exact> is set to C<1>.
-The current implementation only produces exact results for particular
-cases of piecewise translations (i.e., piecewise uniform dependences).
 
 =item * Transitive closure
 
@@ -1436,8 +1606,6 @@ cases of piecewise translations (i.e., piecewise uniform dependences).
 Compute the transitive closure of C<map>.
 The result may be an overapproximation.  If the result is known to be exact,
 then C<*exact> is set to C<1>.
-The current implementation only produces exact results for particular
-cases of piecewise translations (i.e., piecewise uniform dependences).
 
 =item * Reaching path lengths
 
@@ -1494,6 +1662,47 @@ then the name of the space is also removed.
 The function above constructs a relation
 that maps the input set to a flattened version of the set.
 
+=item * Lifting
+
+Lift the input set to a space with extra dimensions corresponding
+to the existentially quantified variables in the input.
+In particular, the result lives in a wrapped map where the domain
+is the original space and the range corresponds to the original
+existentially quantified variables.
+
+       __isl_give isl_basic_set *isl_basic_set_lift(
+               __isl_take isl_basic_set *bset);
+       __isl_give isl_set *isl_set_lift(
+               __isl_take isl_set *set);
+       __isl_give isl_union_set *isl_union_set_lift(
+               __isl_take isl_union_set *uset);
+
+=item * Internal Product
+
+       __isl_give isl_basic_map *isl_basic_map_zip(
+               __isl_take isl_basic_map *bmap);
+       __isl_give isl_map *isl_map_zip(
+               __isl_take isl_map *map);
+       __isl_give isl_union_map *isl_union_map_zip(
+               __isl_take isl_union_map *umap);
+
+Given a relation with nested relations for domain and range,
+interchange the range of the domain with the domain of the range.
+
+=item * Aligning parameters
+
+       __isl_give isl_set *isl_set_align_params(
+               __isl_take isl_set *set,
+               __isl_take isl_dim *model);
+       __isl_give isl_map *isl_map_align_params(
+               __isl_take isl_map *map,
+               __isl_take isl_dim *model);
+
+Change the order of the parameters of the given set or relation
+such that the first parameters match those of C<model>.
+This may involve the introduction of extra parameters.
+All parameters need to be named.
+
 =item * Dimension manipulation
 
        __isl_give isl_set *isl_set_add_dims(
@@ -1506,7 +1715,9 @@ that maps the input set to a flattened version of the set.
 It is usually not advisable to directly change the (input or output)
 space of a set or a relation as this removes the name and the internal
 structure of the space.  However, the above functions can be useful
-to add new parameters.
+to add new parameters, assuming
+C<isl_set_align_params> and C<isl_map_align_params>
+are not sufficient.
 
 =back
 
@@ -1803,12 +2014,15 @@ Matrices can be created, copied and freed using the following functions.
 Note that the elements of a newly created matrix may have arbitrary values.
 The elements can be changed and inspected using the following functions.
 
+       isl_ctx *isl_mat_get_ctx(__isl_keep isl_mat *mat);
        int isl_mat_rows(__isl_keep isl_mat *mat);
        int isl_mat_cols(__isl_keep isl_mat *mat);
        int isl_mat_get_element(__isl_keep isl_mat *mat,
                int row, int col, isl_int *v);
        __isl_give isl_mat *isl_mat_set_element(__isl_take isl_mat *mat,
                int row, int col, isl_int v);
+       __isl_give isl_mat *isl_mat_set_element_si(__isl_take isl_mat *mat,
+               int row, int col, int v);
 
 C<isl_mat_get_element> will return a negative value if anything went wrong.
 In that case, the value of C<*v> is undefined.
@@ -1826,6 +2040,92 @@ the original and the kernel (in that order) is the zero matrix.
 
        __isl_give isl_mat *isl_mat_right_kernel(__isl_take isl_mat *mat);
 
+=head2 Quasi Affine Expressions
+
+The zero quasi affine expression can be created using
+
+       __isl_give isl_aff *isl_aff_zero(
+               __isl_take isl_local_space *ls);
+
+Quasi affine expressions can be copied and free using
+
+       #include <isl/aff.h>
+       __isl_give isl_aff *isl_aff_copy(__isl_keep isl_aff *aff);
+       void *isl_aff_free(__isl_take isl_aff *aff);
+
+A (rational) bound on a dimension can be extracted from an C<isl_constraint>
+using the following function.  The constraint is required to have
+a non-zero coefficient for the specified dimension.
+
+       #include <isl/constraint.h>
+       __isl_give isl_aff *isl_constraint_get_bound(
+               __isl_keep isl_constraint *constraint,
+               enum isl_dim_type type, int pos);
+
+Conversely, an equality constraint can be constructed, equating
+the affine expression to zero, using
+
+       __isl_give isl_constraint *isl_equality_from_aff(
+               __isl_take isl_aff *aff);
+
+The expression can be inspected using
+
+       #include <isl/aff.h>
+       isl_ctx *isl_aff_get_ctx(__isl_keep isl_aff *aff);
+       int isl_aff_dim(__isl_keep isl_aff *aff,
+               enum isl_dim_type type);
+       __isl_give isl_local_space *isl_aff_get_local_space(
+               __isl_keep isl_aff *aff);
+       const char *isl_aff_get_dim_name(__isl_keep isl_aff *aff,
+               enum isl_dim_type type, unsigned pos);
+       int isl_aff_get_constant(__isl_keep isl_aff *aff,
+               isl_int *v);
+       int isl_aff_get_coefficient(__isl_keep isl_aff *aff,
+               enum isl_dim_type type, int pos, isl_int *v);
+       int isl_aff_get_denominator(__isl_keep isl_aff *aff,
+               isl_int *v);
+       __isl_give isl_div *isl_aff_get_div(
+               __isl_keep isl_aff *aff, int pos);
+
+It can be modified using
+
+       #include <isl/aff.h>
+       __isl_give isl_aff *isl_aff_set_constant(
+               __isl_take isl_aff *aff, isl_int v);
+       __isl_give isl_aff *isl_aff_set_constant_si(
+               __isl_take isl_aff *aff, int v);
+       __isl_give isl_aff *isl_aff_set_coefficient(
+               __isl_take isl_aff *aff,
+               enum isl_dim_type type, int pos, isl_int v);
+       __isl_give isl_aff *isl_aff_set_coefficient_si(
+               __isl_take isl_aff *aff,
+               enum isl_dim_type type, int pos, int v);
+       __isl_give isl_aff *isl_aff_set_denominator(
+               __isl_take isl_aff *aff, isl_int v);
+
+       __isl_give isl_aff *isl_aff_add_constant(
+               __isl_take isl_aff *aff, isl_int v);
+       __isl_give isl_aff *isl_aff_add_coefficient_si(
+               __isl_take isl_aff *aff,
+               enum isl_dim_type type, int pos, int v);
+
+Note that the C<set_constant> and C<set_coefficient> functions
+set the I<numerator> of the constant or coefficient, while
+C<add_constant> and C<add_coefficient> add an integer value to
+the possibly rational constant or coefficient.
+
+Operations include
+
+       #include <isl/aff.h>
+       __isl_give isl_aff *isl_aff_neg(__isl_take isl_aff *aff);
+       __isl_give isl_aff *isl_aff_ceil(__isl_take isl_aff *aff);
+
+An expression can be printed using
+
+       #include <isl/aff.h>
+       __isl_give isl_printer *isl_printer_print_aff(
+               __isl_take isl_printer *p, __isl_keep isl_aff *aff);
+
 =head2 Points
 
 Points are elements of a set.  They can be used to construct
@@ -1991,6 +2291,8 @@ on the resulting quasipolynomials
        __isl_give isl_qpolynomial *isl_qpolynomial_var(
                __isl_take isl_dim *dim,
                enum isl_dim_type type, unsigned pos);
+       __isl_give isl_qpolynomial *isl_qpolynomial_from_aff(
+               __isl_take isl_aff *aff);
 
 The zero piecewise quasipolynomial or a piecewise quasipolynomial
 with a single cell can be created using the following functions.
@@ -2176,9 +2478,17 @@ are returned in C<*n> and C<*d>, respectively.
                __isl_take isl_union_pw_qpolynomial *upwpq,
                __isl_take isl_union_set *uset);
 
+       __isl_give isl_qpolynomial *isl_qpolynomial_align_params(
+               __isl_take isl_qpolynomial *qp,
+               __isl_take isl_dim *model);
+
        __isl_give isl_union_pw_qpolynomial *isl_union_pw_qpolynomial_coalesce(
                __isl_take isl_union_pw_qpolynomial *upwqp);
 
+       __isl_give isl_qpolynomial *isl_qpolynomial_gist(
+               __isl_take isl_qpolynomial *qp,
+               __isl_take isl_set *context);
+
        __isl_give isl_pw_qpolynomial *isl_pw_qpolynomial_gist(
                __isl_take isl_pw_qpolynomial *pwqp,
                __isl_take isl_set *context);
@@ -2515,6 +2825,51 @@ Any of C<must_dep>, C<may_dep>, C<must_no_source>
 or C<may_no_source> may be C<NULL>, but a C<NULL> value for
 any of the other arguments is treated as an error.
 
+=head2 Scheduling
+
+B<The functionality described in this section is fairly new
+and may be subject to change.>
+
+The following function can be used to compute a schedule
+for a union of domains.  The generated schedule respects
+all C<validity> dependences.  That is, all dependence distances
+over these dependences in the scheduled space are lexicographically
+positive.  The generated schedule schedule also tries to minimize
+the dependence distances over C<proximity> dependences.
+Moreover, it tries to obtain sequences (bands) of schedule dimensions
+for groups of domains where the dependence distances have only
+non-negative values.
+The algorithm used to construct the schedule is similar to that
+of C<Pluto>.
+
+       #include <isl/schedule.h>
+       __isl_give isl_schedule *isl_union_set_compute_schedule(
+               __isl_take isl_union_set *domain,
+               __isl_take isl_union_map *validity,
+               __isl_take isl_union_map *proximity);
+       void *isl_schedule_free(__isl_take isl_schedule *sched);
+
+A mapping from the domains to the scheduled space can be obtained
+from an C<isl_schedule> using the following function.
+
+       __isl_give isl_union_map *isl_schedule_get_map(
+               __isl_keep isl_schedule *sched);
+
+This mapping can also be obtained in pieces using the following functions.
+
+       int isl_schedule_n_band(__isl_keep isl_schedule *sched);
+       __isl_give isl_union_map *isl_schedule_get_band(
+               __isl_keep isl_schedule *sched, unsigned band);
+
+C<isl_schedule_n_band> returns the maximal number of bands.
+C<isl_schedule_get_band> returns a union of mappings from a domain to
+the band of consecutive schedule dimensions with the given sequence
+number for that domain.  Bands with the same sequence number but for
+different domains may be completely unrelated.
+Within a band, the corresponding coordinates of the distance vectors
+are all non-negative, assuming that the coordinates for all previous
+bands are all zero.
+
 =head2 Parametric Vertex Enumeration
 
 The parametric vertex enumeration described in this section
@@ -2617,27 +2972,3 @@ is given, then the constant should appear in the last column.
 
 Given a polytope, C<isl_polytope_scan> prints
 all integer points in the polytope.
-
-=head1 C<isl-polylib>
-
-The C<isl-polylib> library provides the following functions for converting
-between C<isl> objects and C<PolyLib> objects.
-The library is distributed separately for licensing reasons.
-
-       #include <isl_set_polylib.h>
-       __isl_give isl_basic_set *isl_basic_set_new_from_polylib(
-               Polyhedron *P, __isl_take isl_dim *dim);
-       Polyhedron *isl_basic_set_to_polylib(
-               __isl_keep isl_basic_set *bset);
-       __isl_give isl_set *isl_set_new_from_polylib(Polyhedron *D,
-               __isl_take isl_dim *dim);
-       Polyhedron *isl_set_to_polylib(__isl_keep isl_set *set);
-
-       #include <isl_map_polylib.h>
-       __isl_give isl_basic_map *isl_basic_map_new_from_polylib(
-               Polyhedron *P, __isl_take isl_dim *dim);
-       __isl_give isl_map *isl_map_new_from_polylib(Polyhedron *D,
-               __isl_take isl_dim *dim);
-       Polyhedron *isl_basic_map_to_polylib(
-               __isl_keep isl_basic_map *bmap);
-       Polyhedron *isl_map_to_polylib(__isl_keep isl_map *map);