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