1 #include <string.h>
2 #include <isl/val_gmp.h>
3 #include <isl_val_private.h>
4
5 /* Return a reference to an isl_val representing the integer "z".
6 */
isl_val_int_from_gmp(isl_ctx * ctx,mpz_t z)7 __isl_give isl_val *isl_val_int_from_gmp(isl_ctx *ctx, mpz_t z)
8 {
9 isl_val *v;
10
11 v = isl_val_alloc(ctx);
12 if (!v)
13 return NULL;
14
15 isl_int_set(v->n, z);
16 isl_int_set_si(v->d, 1);
17
18 return v;
19 }
20
21 /* Return a reference to an isl_val representing the rational value "n"/"d".
22 */
isl_val_from_gmp(isl_ctx * ctx,const mpz_t n,const mpz_t d)23 __isl_give isl_val *isl_val_from_gmp(isl_ctx *ctx, const mpz_t n, const mpz_t d)
24 {
25 isl_val *v;
26
27 v = isl_val_alloc(ctx);
28 if (!v)
29 return NULL;
30
31 isl_int_set(v->n, n);
32 isl_int_set(v->d, d);
33
34 return isl_val_normalize(v);
35 }
36
37 /* Extract the numerator of a rational value "v" in "z".
38 *
39 * If "v" is not a rational value, then the result is undefined.
40 */
isl_val_get_num_gmp(__isl_keep isl_val * v,mpz_t z)41 int isl_val_get_num_gmp(__isl_keep isl_val *v, mpz_t z)
42 {
43 if (!v)
44 return -1;
45 if (!isl_val_is_rat(v))
46 isl_die(isl_val_get_ctx(v), isl_error_invalid,
47 "expecting rational value", return -1);
48 mpz_set(z, v->n);
49 return 0;
50 }
51
52 /* Extract the denominator of a rational value "v" in "z".
53 *
54 * If "v" is not a rational value, then the result is undefined.
55 */
isl_val_get_den_gmp(__isl_keep isl_val * v,mpz_t z)56 int isl_val_get_den_gmp(__isl_keep isl_val *v, mpz_t z)
57 {
58 if (!v)
59 return -1;
60 if (!isl_val_is_rat(v))
61 isl_die(isl_val_get_ctx(v), isl_error_invalid,
62 "expecting rational value", return -1);
63 mpz_set(z, v->d);
64 return 0;
65 }
66
67 /* Return a reference to an isl_val representing the unsigned
68 * integer value stored in the "n" chunks of size "size" at "chunks".
69 * The least significant chunk is assumed to be stored first.
70 */
isl_val_int_from_chunks(isl_ctx * ctx,size_t n,size_t size,const void * chunks)71 __isl_give isl_val *isl_val_int_from_chunks(isl_ctx *ctx, size_t n,
72 size_t size, const void *chunks)
73 {
74 isl_val *v;
75
76 v = isl_val_alloc(ctx);
77 if (!v)
78 return NULL;
79
80 mpz_import(v->n, n, -1, size, 0, 0, chunks);
81 isl_int_set_si(v->d, 1);
82
83 return v;
84 }
85
86 /* Return the number of chunks of size "size" required to
87 * store the absolute value of the numerator of "v".
88 */
isl_val_n_abs_num_chunks(__isl_keep isl_val * v,size_t size)89 isl_size isl_val_n_abs_num_chunks(__isl_keep isl_val *v, size_t size)
90 {
91 if (!v)
92 return isl_size_error;
93
94 if (!isl_val_is_rat(v))
95 isl_die(isl_val_get_ctx(v), isl_error_invalid,
96 "expecting rational value", return isl_size_error);
97
98 size *= 8;
99 return (mpz_sizeinbase(v->n, 2) + size - 1) / size;
100 }
101
102 /* Store a representation of the absolute value of the numerator of "v"
103 * in terms of chunks of size "size" at "chunks".
104 * The least significant chunk is stored first.
105 * The number of chunks in the result can be obtained by calling
106 * isl_val_n_abs_num_chunks. The user is responsible for allocating
107 * enough memory to store the results.
108 *
109 * In the special case of a zero value, isl_val_n_abs_num_chunks will
110 * return one, while mpz_export will not fill in any chunks. We therefore
111 * do it ourselves.
112 */
isl_val_get_abs_num_chunks(__isl_keep isl_val * v,size_t size,void * chunks)113 isl_stat isl_val_get_abs_num_chunks(__isl_keep isl_val *v, size_t size,
114 void *chunks)
115 {
116 if (!v || !chunks)
117 return isl_stat_error;
118
119 if (!isl_val_is_rat(v))
120 isl_die(isl_val_get_ctx(v), isl_error_invalid,
121 "expecting rational value", return isl_stat_error);
122
123 mpz_export(chunks, NULL, -1, size, 0, 0, v->n);
124 if (isl_val_is_zero(v))
125 memset(chunks, 0, size);
126
127 return isl_stat_ok;
128 }
129