129 lines · c
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 */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 */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 */41int 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 */56int 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 unsigned68 * integer value stored in the "n" chunks of size "size" at "chunks".69 * The least significant chunk is assumed to be stored first.70 */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 to87 * store the absolute value of the numerator of "v".88 */89isl_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 calling106 * isl_val_n_abs_num_chunks. The user is responsible for allocating107 * enough memory to store the results.108 *109 * In the special case of a zero value, isl_val_n_abs_num_chunks will110 * return one, while mpz_export will not fill in any chunks. We therefore111 * do it ourselves.112 */113isl_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