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 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 */ 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 */ 89size_t isl_val_n_abs_num_chunks(__isl_keep isl_val *v, size_t size) 90{ 91 if (!v) 92 return 0; 93 94 if (!isl_val_is_rat(v)) 95 isl_die(isl_val_get_ctx(v), isl_error_invalid, 96 "expecting rational value", return 0); 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 */ 113int isl_val_get_abs_num_chunks(__isl_keep isl_val *v, size_t size, 114 void *chunks) 115{ 116 if (!v || !chunks) 117 return -1; 118 119 if (!isl_val_is_rat(v)) 120 isl_die(isl_val_get_ctx(v), isl_error_invalid, 121 "expecting rational value", return -1); 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 0; 128} 129