isl_val_gmp.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  */
  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  */
  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  */
  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  */
  89 size_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  */
 113 int 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 }