/*

 * Copyright 2008-2009 Katholieke Universiteit Leuven

 *

 * Use of this software is governed by the MIT license

 *

 * Written by Sven Verdoolaege, K.U.Leuven, Departement

 * Computerwetenschappen, Celestijnenlaan 200A, B-3001 Leuven, Belgium

 */

#ifndef ISL_MAP_PRIVATE_H

#define ISL_MAP_PRIVATE_H

#define isl_basic_set	isl_basic_map

#define isl_set		isl_map

#define isl_basic_set_list	isl_basic_map_list

#define isl_set_list	isl_map_list

#include <isl/list.h>

ISL_DECLARE_LIST(basic_map)

ISL_DECLARE_LIST(map)

#include <isl/set.h>

#include <isl/map.h>

#include <isl_reordering.h>

#include <isl/vec.h>

#include <isl_blk.h>

/* A "basic map" is a relation between two sets of variables,

 * called the "in" and "out" variables.

 * A "basic set" is a basic map with a zero-dimensional

 * domain.

 *

 * It is implemented as a set with two extra fields:

 * n_in is the number of in variables

 * n_out is the number of out variables

 * n_in + n_out should be equal to set.dim

 */

struct isl_basic_map {

	int ref;

#define ISL_BASIC_MAP_FINAL		(1 << 0)

#define ISL_BASIC_MAP_EMPTY		(1 << 1)

#define ISL_BASIC_MAP_NO_IMPLICIT	(1 << 2)

#define ISL_BASIC_MAP_NO_REDUNDANT	(1 << 3)

#define ISL_BASIC_MAP_RATIONAL		(1 << 4)

#define ISL_BASIC_MAP_NORMALIZED	(1 << 5)

#define ISL_BASIC_MAP_NORMALIZED_DIVS	(1 << 6)

#define ISL_BASIC_MAP_ALL_EQUALITIES	(1 << 7)

#define ISL_BASIC_SET_FINAL		(1 << 0)

#define ISL_BASIC_SET_EMPTY		(1 << 1)

#define ISL_BASIC_SET_NO_IMPLICIT	(1 << 2)

#define ISL_BASIC_SET_NO_REDUNDANT	(1 << 3)

#define ISL_BASIC_SET_RATIONAL		(1 << 4)

#define ISL_BASIC_SET_NORMALIZED	(1 << 5)

#define ISL_BASIC_SET_NORMALIZED_DIVS	(1 << 6)

#define ISL_BASIC_SET_ALL_EQUALITIES	(1 << 7)

	unsigned flags;

	struct isl_ctx *ctx;

	isl_space *dim;

	unsigned extra;

	unsigned n_eq;

	unsigned n_ineq;

	size_t c_size;

	isl_int **eq;

	isl_int **ineq;

	unsigned n_div;

	isl_int **div;

	struct isl_vec *sample;

	struct isl_blk block;

	struct isl_blk block2;

};

#undef EL

#define EL isl_basic_set

#include <isl_list_templ.h>

/* A "map" is a (possibly disjoint) union of basic maps.

 * A "set" is a (possibly disjoint) union of basic sets.

 *

 * Currently, the isl_set structure is identical to the isl_map structure

 * and the library depends on this correspondence internally.

 * However, users should not depend on this correspondence.

 */

struct isl_map {

	int ref;

#define ISL_MAP_DISJOINT		(1 << 0)

#define ISL_MAP_NORMALIZED		(1 << 1)

#define ISL_SET_DISJOINT		(1 << 0)

#define ISL_SET_NORMALIZED		(1 << 1)

	unsigned flags;

	struct isl_ctx *ctx;

	isl_space *dim;

	int n;

	size_t size;

	struct isl_basic_map *p[1];

};

#undef EL

#define EL isl_set

#include <isl_list_templ.h>

__isl_give isl_map *isl_map_realign(__isl_take isl_map *map,

	__isl_take isl_reordering *r);

__isl_give isl_set *isl_set_realign(__isl_take isl_set *set,

	__isl_take isl_reordering *r);

__isl_give isl_map *isl_map_reset(__isl_take isl_map *map,

	enum isl_dim_type type);

__isl_give isl_basic_set *isl_basic_set_reset_space(

	__isl_take isl_basic_set *bset, __isl_take isl_space *dim);

__isl_give isl_basic_map *isl_basic_map_reset_space(

	__isl_take isl_basic_map *bmap, __isl_take isl_space *dim);

__isl_give isl_map *isl_map_reset_space(__isl_take isl_map *map,

	__isl_take isl_space *dim);

unsigned isl_basic_map_offset(struct isl_basic_map *bmap,

					enum isl_dim_type type);

unsigned isl_basic_set_offset(struct isl_basic_set *bset,

					enum isl_dim_type type);

int isl_basic_map_may_be_set(__isl_keep isl_basic_map *bmap);

int isl_map_may_be_set(__isl_keep isl_map *map);

int isl_map_compatible_domain(struct isl_map *map, struct isl_set *set);

int isl_basic_map_compatible_domain(struct isl_basic_map *bmap,

		struct isl_basic_set *bset);

int isl_basic_map_compatible_range(struct isl_basic_map *bmap,

		struct isl_basic_set *bset);

struct isl_basic_map *isl_basic_map_extend_space(struct isl_basic_map *base,

		__isl_take isl_space *dim, unsigned extra,

		unsigned n_eq, unsigned n_ineq);

struct isl_basic_set *isl_basic_set_extend_space(struct isl_basic_set *base,

		__isl_take isl_space *dim, unsigned extra,

		unsigned n_eq, unsigned n_ineq);

struct isl_basic_set *isl_basic_set_add_constraints(struct isl_basic_set *bset1,

		struct isl_basic_set *bset2, unsigned pos);

struct isl_map *isl_map_grow(struct isl_map *map, int n);

struct isl_set *isl_set_grow(struct isl_set *set, int n);

int isl_basic_set_contains(struct isl_basic_set *bset, struct isl_vec *vec);

int isl_basic_map_contains(struct isl_basic_map *bmap, struct isl_vec *vec);

__isl_give isl_basic_set *isl_basic_set_alloc_space(__isl_take isl_space *dim,

		unsigned extra, unsigned n_eq, unsigned n_ineq);

__isl_give isl_set *isl_set_alloc_space(__isl_take isl_space *dim, int n,

	unsigned flags);

__isl_give isl_basic_map *isl_basic_map_alloc_space(__isl_take isl_space *dim,

		unsigned extra, unsigned n_eq, unsigned n_ineq);

__isl_give isl_map *isl_map_alloc_space(__isl_take isl_space *dim, int n,

	unsigned flags);

unsigned isl_basic_map_total_dim(const struct isl_basic_map *bmap);

int isl_basic_map_alloc_equality(struct isl_basic_map *bmap);

int isl_basic_set_alloc_equality(struct isl_basic_set *bset);

int isl_basic_set_free_inequality(struct isl_basic_set *bset, unsigned n);

int isl_basic_map_free_equality(struct isl_basic_map *bmap, unsigned n);

int isl_basic_set_free_equality(struct isl_basic_set *bset, unsigned n);

int isl_basic_set_alloc_inequality(struct isl_basic_set *bset);

int isl_basic_map_alloc_inequality(struct isl_basic_map *bmap);

int isl_basic_map_free_inequality(struct isl_basic_map *bmap, unsigned n);

int isl_basic_map_alloc_div(struct isl_basic_map *bmap);

int isl_basic_set_alloc_div(struct isl_basic_set *bset);

int isl_basic_map_free_div(struct isl_basic_map *bmap, unsigned n);

int isl_basic_set_free_div(struct isl_basic_set *bset, unsigned n);

void isl_basic_map_inequality_to_equality(

		struct isl_basic_map *bmap, unsigned pos);

int isl_basic_map_drop_equality(struct isl_basic_map *bmap, unsigned pos);

int isl_basic_set_drop_equality(struct isl_basic_set *bset, unsigned pos);

int isl_basic_set_drop_inequality(struct isl_basic_set *bset, unsigned pos);

int isl_basic_map_drop_inequality(struct isl_basic_map *bmap, unsigned pos);

__isl_give isl_basic_set *isl_basic_set_add_eq(__isl_take isl_basic_set *bset,

	isl_int *eq);

__isl_give isl_basic_map *isl_basic_map_add_eq(__isl_take isl_basic_map *bmap,

	isl_int *eq);

__isl_give isl_basic_set *isl_basic_set_add_ineq(__isl_take isl_basic_set *bset,

	isl_int *ineq);

__isl_give isl_basic_map *isl_basic_map_add_ineq(__isl_take isl_basic_map *bmap,

	isl_int *ineq);

int isl_inequality_negate(struct isl_basic_map *bmap, unsigned pos);

struct isl_basic_set *isl_basic_set_cow(struct isl_basic_set *bset);

struct isl_basic_map *isl_basic_map_cow(struct isl_basic_map *bmap);

struct isl_set *isl_set_cow(struct isl_set *set);

struct isl_map *isl_map_cow(struct isl_map *map);

struct isl_basic_map *isl_basic_map_set_to_empty(struct isl_basic_map *bmap);

struct isl_basic_set *isl_basic_set_set_to_empty(struct isl_basic_set *bset);

struct isl_basic_set *isl_basic_set_order_divs(struct isl_basic_set *bset);

void isl_basic_map_swap_div(struct isl_basic_map *bmap, int a, int b);

struct isl_basic_map *isl_basic_map_order_divs(struct isl_basic_map *bmap);

__isl_give isl_map *isl_map_order_divs(__isl_take isl_map *map);

struct isl_basic_map *isl_basic_map_align_divs(

		struct isl_basic_map *dst, struct isl_basic_map *src);

struct isl_basic_set *isl_basic_set_align_divs(

		struct isl_basic_set *dst, struct isl_basic_set *src);

__isl_give isl_set *isl_set_align_divs_to_basic_set_list(

	__isl_take isl_set *set, __isl_keep isl_basic_set_list *list);

__isl_give isl_basic_set_list *isl_basic_set_list_align_divs_to_basic_set(

	__isl_take isl_basic_set_list *list, __isl_keep isl_basic_set *bset);

__isl_give isl_basic_map *isl_basic_map_sort_divs(

	__isl_take isl_basic_map *bmap);

__isl_give isl_map *isl_map_sort_divs(__isl_take isl_map *map);

struct isl_basic_map *isl_basic_map_gauss(

	struct isl_basic_map *bmap, int *progress);

struct isl_basic_set *isl_basic_set_gauss(

	struct isl_basic_set *bset, int *progress);

__isl_give isl_basic_set *isl_basic_set_sort_constraints(

	__isl_take isl_basic_set *bset);

int isl_basic_map_plain_cmp(const __isl_keep isl_basic_map *bmap1,

	const __isl_keep isl_basic_map *bmap2);

int isl_basic_map_plain_is_equal(__isl_keep isl_basic_map *bmap1,

	__isl_keep isl_basic_map *bmap2);

struct isl_basic_map *isl_basic_map_normalize_constraints(

	struct isl_basic_map *bmap);

struct isl_basic_set *isl_basic_set_normalize_constraints(

	struct isl_basic_set *bset);

struct isl_basic_map *isl_basic_map_implicit_equalities(

						struct isl_basic_map *bmap);

struct isl_basic_set *isl_basic_map_underlying_set(struct isl_basic_map *bmap);

__isl_give isl_basic_set *isl_basic_set_underlying_set(

		__isl_take isl_basic_set *bset);

__isl_give isl_basic_set_list *isl_basic_set_list_underlying_set(

	__isl_take isl_basic_set_list *list);

struct isl_set *isl_map_underlying_set(struct isl_map *map);

struct isl_basic_map *isl_basic_map_overlying_set(struct isl_basic_set *bset,

	struct isl_basic_map *like);

__isl_give isl_basic_set *isl_basic_set_drop_constraints_involving(

	__isl_take isl_basic_set *bset, unsigned first, unsigned n);

__isl_give isl_basic_set *isl_basic_set_drop(__isl_take isl_basic_set *bset,

	enum isl_dim_type type, unsigned first, unsigned n);

struct isl_basic_map *isl_basic_map_drop(struct isl_basic_map *bmap,

	enum isl_dim_type type, unsigned first, unsigned n);

struct isl_set *isl_set_drop(struct isl_set *set,

	enum isl_dim_type type, unsigned first, unsigned n);

struct isl_basic_set *isl_basic_set_drop_dims(

		struct isl_basic_set *bset, unsigned first, unsigned n);

struct isl_set *isl_set_drop_dims(

		struct isl_set *set, unsigned first, unsigned n);

struct isl_map *isl_map_drop_inputs(

		struct isl_map *map, unsigned first, unsigned n);

struct isl_map *isl_map_drop(struct isl_map *map,

	enum isl_dim_type type, unsigned first, unsigned n);

__isl_give isl_basic_map *isl_basic_map_remove_duplicate_constraints(

	__isl_take isl_basic_map *bmap, int *progress, int detect_divs);

struct isl_map *isl_map_remove_empty_parts(struct isl_map *map);

struct isl_set *isl_set_remove_empty_parts(struct isl_set *set);

__isl_give isl_map *isl_map_remove_obvious_duplicates(__isl_take isl_map *map);

struct isl_set *isl_set_normalize(struct isl_set *set);

struct isl_set *isl_set_drop_vars(

		struct isl_set *set, unsigned first, unsigned n);

struct isl_basic_map *isl_basic_map_eliminate_vars(

	struct isl_basic_map *bmap, unsigned pos, unsigned n);

struct isl_basic_set *isl_basic_set_eliminate_vars(

	struct isl_basic_set *bset, unsigned pos, unsigned n);

__isl_give isl_map *isl_map_eliminate(__isl_take isl_map *map,

	enum isl_dim_type type, unsigned first, unsigned n);

__isl_give isl_set *isl_set_eliminate(__isl_take isl_set *set,

	enum isl_dim_type type, unsigned first, unsigned n);

int isl_basic_set_constraint_is_redundant(struct isl_basic_set **bset,

	isl_int *c, isl_int *opt_n, isl_int *opt_d);

int isl_basic_map_add_div_constraint(__isl_keep isl_basic_map *bmap,

	unsigned div, int sign);

int isl_basic_map_add_div_constraints(struct isl_basic_map *bmap, unsigned div);

__isl_give isl_basic_map *isl_basic_map_add_known_div_constraints(

	__isl_take isl_basic_map *bmap);

struct isl_basic_map *isl_basic_map_drop_redundant_divs(

	struct isl_basic_map *bmap);

struct isl_basic_set *isl_basic_set_drop_redundant_divs(

	struct isl_basic_set *bset);

struct isl_basic_set *isl_basic_set_recession_cone(struct isl_basic_set *bset);

struct isl_basic_set *isl_basic_set_lineality_space(struct isl_basic_set *bset);

struct isl_basic_set *isl_basic_set_set_rational(struct isl_basic_set *bset);

__isl_give isl_set *isl_set_set_rational(__isl_take isl_set *set);

__isl_give isl_basic_map *isl_basic_map_set_rational(

	__isl_take isl_basic_map *bmap);

__isl_give isl_map *isl_map_set_rational(__isl_take isl_map *map);

int isl_map_has_rational(__isl_keep isl_map *map);

int isl_set_has_rational(__isl_keep isl_set *set);

struct isl_mat;

struct isl_basic_set *isl_basic_set_preimage(struct isl_basic_set *bset,

	struct isl_mat *mat);

struct isl_set *isl_set_preimage(struct isl_set *set, struct isl_mat *mat);

__isl_give isl_basic_set *isl_basic_set_transform_dims(

	__isl_take isl_basic_set *bset, enum isl_dim_type type, unsigned first,

	__isl_take isl_mat *trans);

isl_int *isl_set_wrap_facet(__isl_keep isl_set *set,

	isl_int *facet, isl_int *ridge);

int isl_basic_map_contains_point(__isl_keep isl_basic_map *bmap,

	__isl_keep isl_point *point);

int isl_set_contains_point(__isl_keep isl_set *set, __isl_keep isl_point *point);

int isl_basic_set_vars_get_sign(__isl_keep isl_basic_set *bset,

	unsigned first, unsigned n, int *signs);

int isl_set_foreach_orthant(__isl_keep isl_set *set,

	int (*fn)(__isl_take isl_set *orthant, int *signs, void *user),

	void *user);

int isl_basic_map_add_div_constraints_var(__isl_keep isl_basic_map *bmap,

	unsigned pos, isl_int *div);

int isl_basic_set_add_div_constraints_var(__isl_keep isl_basic_set *bset,

	unsigned pos, isl_int *div);

int isl_basic_map_is_div_constraint(__isl_keep isl_basic_map *bmap,

	isl_int *constraint, unsigned div);

int isl_basic_set_is_div_constraint(__isl_keep isl_basic_set *bset,

	isl_int *constraint, unsigned div);

__isl_give isl_basic_set *isl_basic_set_from_local_space(

	__isl_take isl_local_space *ls);

__isl_give isl_basic_map *isl_basic_map_from_local_space(

	__isl_take isl_local_space *ls);

__isl_give isl_basic_set *isl_basic_set_expand_divs(

	__isl_take isl_basic_set *bset, __isl_take isl_mat *div, int *exp);

int isl_basic_map_divs_known(__isl_keep isl_basic_map *bmap);

__isl_give isl_mat *isl_basic_set_get_divs(__isl_keep isl_basic_set *bset);

__isl_give isl_mat *isl_basic_map_get_divs(__isl_keep isl_basic_map *bmap);

__isl_give isl_map *isl_map_inline_foreach_basic_map(__isl_take isl_map *map,

	__isl_give isl_basic_map *(*fn)(__isl_take isl_basic_map *bmap));

__isl_give isl_map *isl_map_align_params_map_map_and(

	__isl_take isl_map *map1, __isl_take isl_map *map2,

	__isl_give isl_map *(*fn)(__isl_take isl_map *map1,

				    __isl_take isl_map *map2));

int isl_map_align_params_map_map_and_test(__isl_keep isl_map *map1,

	__isl_keep isl_map *map2,

	int (*fn)(__isl_keep isl_map *map1, __isl_keep isl_map *map2));

int isl_basic_map_foreach_lexopt(__isl_keep isl_basic_map *bmap, int max,

	int (*fn)(__isl_take isl_basic_set *dom, __isl_take isl_aff_list *list,

		  void *user),

	void *user);

int isl_basic_set_foreach_lexopt(__isl_keep isl_basic_set *bset, int max,

	int (*fn)(__isl_take isl_basic_set *dom, __isl_take isl_aff_list *list,

		  void *user),

	void *user);

__isl_give isl_set *isl_set_substitute(__isl_take isl_set *set,

	enum isl_dim_type type, unsigned pos, __isl_keep isl_aff *subs);

__isl_give isl_set *isl_set_gist_params_basic_set(__isl_take isl_set *set,

	__isl_take isl_basic_set *context);

int isl_map_compatible_range(__isl_keep isl_map *map, __isl_keep isl_set *set);

int isl_basic_map_plain_is_single_valued(__isl_keep isl_basic_map *bmap);

int isl_map_is_set(__isl_keep isl_map *map);

int isl_basic_set_plain_dim_is_fixed(__isl_keep isl_basic_set *bset,

	unsigned dim, isl_int *val);

int isl_basic_set_dim_residue_class(struct isl_basic_set *bset,

	int pos, isl_int *modulo, isl_int *residue);

int isl_set_dim_residue_class(struct isl_set *set,

	int pos, isl_int *modulo, isl_int *residue);

__isl_give isl_basic_set *isl_basic_set_fix(__isl_take isl_basic_set *bset,

	enum isl_dim_type type, unsigned pos, isl_int value);

__isl_give isl_set *isl_set_fix(__isl_take isl_set *set,

	enum isl_dim_type type, unsigned pos, isl_int value);

int isl_map_plain_is_fixed(__isl_keep isl_map *map,

	enum isl_dim_type type, unsigned pos, isl_int *val);

int isl_basic_map_output_defining_equality(__isl_keep isl_basic_map *bmap,

	int pos);

__isl_give isl_basic_set_list *isl_set_get_basic_set_list(

	__isl_keep isl_set *set);

__isl_give isl_map *isl_map_fixed_power(__isl_take isl_map *map, isl_int exp);

int isl_basic_set_count_upto(__isl_keep isl_basic_set *bset,

	isl_int max, isl_int *count);

int isl_set_count_upto(__isl_keep isl_set *set, isl_int max, isl_int *count);

#endif