4 #include <barvinok/polylib.h>
6 #if defined(__cplusplus)
10 #define polynomial new_polynomial
11 #define periodic new_periodic
12 #define evector new_evector
13 #define enode_type new_enode_type
14 #define enode_node new_enode_node
15 #define _enumeration _new_enumeration
16 #define Enumeration _new_Enumeration
17 #define enode _new_enode
18 #define _enode _new__enode
19 #define evalue _new_evalue
20 #define _evalue _new__evalue
21 #define eadd _new_eadd
22 #define ecopy _new_ecopy
23 #define new_enode _new_new_enode
24 #define free_evalue_refs _new_free_evalue_refs
25 #define print_evalue _new_print_evalue
26 #define print_enode _new_print_enode
27 #define reduce_evalue _new_reduce
28 #define compute_evalue _new_compute_evalue
29 #define compute_poly _new_compute_poly
30 #define in_domain _new_in_domain
31 #define Enumeration_Free _new_Enumeration_Free
32 #define evalue_div _new_evalue_div
34 typedef enum { polynomial
, periodic
, evector
, fractional
, relation
,
35 partition
, flooring
} enode_type
;
37 typedef struct _evalue
{
38 Value d
; /* denominator */
40 Value n
; /* numerator (if denominator > 0) */
41 struct _enode
*p
; /* pointer (if denominator == 0) */
42 Polyhedron
*D
; /* domain (if denominator == -1) */
46 #define EVALUE_DOMAIN(ev) ((ev).x.D)
47 #define EVALUE_SET_DOMAIN(ev, D) \
49 value_set_si((ev).d, -1); \
50 EVALUE_DOMAIN(ev) = D; \
52 #define EVALUE_IS_DOMAIN(ev) (value_mone_p((ev).d))
54 #define EVALUE_IS_ZERO(ev) (value_pos_p((ev).d) && value_zero_p((ev).x.n))
55 #define EVALUE_IS_ONE(ev) (value_one_p((ev).d) && value_one_p((ev).x.n))
57 #define EVALUE_IS_NAN(ev) (value_cmp_si((ev).d, -2) == 0)
59 typedef struct _enode
{
60 enode_type type
; /* polynomial or periodic or evector */
61 int size
; /* number of attached pointers */
62 int pos
; /* parameter position */
63 evalue arr
[1]; /* array of rational/pointer */
66 typedef struct _enumeration
{
68 Polyhedron
*ValidityDomain
; /* contraints on the parameters */
69 evalue EP
; /* dimension = combined space */
70 struct _enumeration
*next
; /* Ehrhart Polynomial, corresponding
71 to parameter values inside the
72 domain ValidityDomain below */
76 void evalue_set_si(evalue
*ev
, int n
, int d
);
77 void evalue_set(evalue
*ev
, Value n
, Value d
);
78 void evalue_copy(evalue
*dst
, const evalue
*src
);
79 evalue
*evalue_dup(const evalue
*e
);
80 evalue
*evalue_zero();
82 evalue
*evalue_var(int var
);
83 void evalue_free(evalue
*e
);
84 enode
*new_enode(enode_type type
,int size
,int pos
);
85 enode
*ecopy(enode
*e
);
86 int eequal(const evalue
*e1
, const evalue
*e2
);
87 void free_evalue_refs(evalue
*e
);
88 void print_evalue(FILE *DST
, const evalue
*e
, const char * const *pname
);
89 void print_enode(FILE *DST
, enode
*p
, const char *const *pname
);
90 void reduce_evalue (evalue
*e
);
91 void reduce_evalue_in_domain(evalue
*e
, Polyhedron
*D
);
92 void aep_evalue(evalue
*e
, int *ref
);
93 void addeliminatedparams_evalue(evalue
*e
,Matrix
*CT
);
94 void addeliminatedparams_enum(evalue
*e
, Matrix
*CT
, Polyhedron
*CEq
,
95 unsigned MaxRays
, unsigned nparam
);
96 void eadd(const evalue
*e1
, evalue
*res
);
97 void emul(const evalue
*e1
, evalue
*res
);
98 void emask(evalue
*mask
, evalue
*res
);
99 int in_domain(Polyhedron
*P
, Value
*list_args
);
100 double compute_evalue(const evalue
*e
, Value
*list_args
);
101 Value
*compute_poly(Enumeration
*en
,Value
*list_args
);
102 evalue
*evalue_eval(const evalue
*e
, Value
*values
);
103 void evalue_mod2table(evalue
*ev
, int nparam
);
104 void evalue_mod2relation(evalue
*e
);
105 void evalue_combine(evalue
*e
);
106 void evalue_range_reduction(evalue
*e
);
107 int evalue_range_reduction_in_domain(evalue
*e
, Polyhedron
*D
);
108 size_t evalue_size(evalue
*e
);
109 size_t domain_size(Polyhedron
*D
);
110 Enumeration
* partition2enumeration(evalue
*EP
);
111 void evalue_frac2floor(evalue
*e
);
112 int evalue_frac2floor_in_domain(evalue
*e
, Polyhedron
*D
);
113 void evalue_frac2floor2(evalue
*e
, int shift
);
114 int evalue_frac2floor_in_domain3(evalue
*e
, Polyhedron
*D
, int shift
);
115 void evalue_frac2polynomial(evalue
*e
, int sign
, unsigned MaxRays
);
116 evalue
*esum(evalue
*E
, int nvar
);
117 evalue
*evalue_sum(evalue
*E
, int nvar
, unsigned MaxRays
);
118 void eor(evalue
*e1
, evalue
*res
);
119 void evalue_denom(const evalue
*e
, Value
*d
);
120 void evalue_div(evalue
*e
, Value n
);
121 void evalue_mul(evalue
*e
, Value n
);
122 void evalue_mul_div(evalue
*e
, Value n
, Value d
);
123 void evalue_negate(evalue
*e
);
124 void evalue_add_constant(evalue
*e
, const Value cst
);
125 void evalue_split_domains_into_orthants(evalue
*e
, unsigned MaxRays
);
126 void evalue_split_periods(evalue
*e
, int max_periods
, unsigned int MaxRays
);
127 void evalue_extract_affine(const evalue
*e
, Value
*coeff
, Value
*cst
, Value
*d
);
128 evalue
*affine2evalue(Value
*coeff
, Value denom
, int nvar
);
129 void evalue_substitute(evalue
*e
, evalue
**subs
);
130 void evalue_backsubstitute(evalue
*e
, Matrix
*CP
, unsigned MaxRays
);
131 evalue
*evalue_polynomial(Vector
*c
, const evalue
* X
);
133 struct evalue_section
{
138 evalue
*evalue_from_section_array(struct evalue_section
*s
, int n
);
140 #if defined(__cplusplus)