isl_val_sioimath.c
1.85 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
#include <isl_val_private.h>
/* Return a reference to an isl_val representing the unsigned
* integer value stored in the "n" chunks of size "size" at "chunks".
* The least significant chunk is assumed to be stored first.
*/
__isl_give isl_val *isl_val_int_from_chunks(isl_ctx *ctx, size_t n,
size_t size, const void *chunks)
{
isl_val *v;
v = isl_val_alloc(ctx);
if (!v)
return NULL;
impz_import(isl_sioimath_reinit_big(v->n), n, -1, size, 0, 0, chunks);
isl_sioimath_try_demote(v->n);
isl_int_set_si(v->d, 1);
return v;
}
/* Store a representation of the absolute value of the numerator of "v"
* in terms of chunks of size "size" at "chunks".
* The least significant chunk is stored first.
* The number of chunks in the result can be obtained by calling
* isl_val_n_abs_num_chunks. The user is responsible for allocating
* enough memory to store the results.
*
* In the special case of a zero value, isl_val_n_abs_num_chunks will
* return one, while impz_export will not fill in any chunks. We therefore
* do it ourselves.
*/
int isl_val_get_abs_num_chunks(__isl_keep isl_val *v, size_t size,
void *chunks)
{
isl_sioimath_scratchspace_t scratch;
if (!v || !chunks)
return -1;
if (!isl_val_is_rat(v))
isl_die(isl_val_get_ctx(v), isl_error_invalid,
"expecting rational value", return -1);
impz_export(chunks, NULL, -1, size, 0, 0,
isl_sioimath_bigarg_src(*v->n, &scratch));
if (isl_val_is_zero(v))
memset(chunks, 0, size);
return 0;
}
/* Return the number of chunks of size "size" required to
* store the absolute value of the numerator of "v".
*/
size_t isl_val_n_abs_num_chunks(__isl_keep isl_val *v, size_t size)
{
if (!v)
return 0;
if (!isl_val_is_rat(v))
isl_die(isl_val_get_ctx(v), isl_error_invalid,
"expecting rational value", return 0);
size *= 8;
return (isl_sioimath_sizeinbase(*v->n, 2) + size - 1) / size;
}