include/cloog/int.h.in

   1 #ifndef CLOOG_INT_H
   2 #define CLOOG_INT_H
   3
   4 #undef CLOOG_INT_INT
   5 #undef CLOOG_INT_LONG
   6 #undef CLOOG_INT_LONG_LONG
   7 #undef CLOOG_INT_GMP
   8
   9 #include <assert.h>
  10 #include <stdio.h>
  11 #if defined(CLOOG_INT_GMP)
  12 #include <string.h>
  13 #include <gmp.h>
  14 #ifndef mp_get_memory_functions
  15 void mp_get_memory_functions(
  16                 void *(**alloc_func_ptr) (size_t),
  17                 void *(**realloc_func_ptr) (void *, size_t, size_t),
  18                 void (**free_func_ptr) (void *, size_t));
  19 #endif
  20 #endif
  21
  22 #if defined(__cplusplus)
  23 extern "C"
  24   {
  25 #endif
  26
  27 #if defined(CLOOG_INT_INT)
  28 typedef int             cloog_int_t;
  29 #define CLOOG_INT_FORMAT        "%d"
  30 #elif defined(CLOOG_INT_LONG)
  31 typedef long            cloog_int_t;
  32 #define CLOOG_INT_FORMAT        "%ld"
  33 #elif defined(CLOOG_INT_LONG_LONG)
  34 typedef long long       cloog_int_t;
  35 #define CLOOG_INT_FORMAT        "%lld"
  36 #elif defined(CLOOG_INT_GMP)
  37 typedef mpz_t           cloog_int_t;
  38 #else
  39 #error "No integer type defined"
  40 #endif
  41
  42 #if defined(CLOOG_INT_GMP)
  43
  44 #define cloog_int_init(i)       do {                                    \
  45                                         mpz_init(i);                    \
  46                                         cloog_int_leak_up();            \
  47                                 } while (0)
  48 #define cloog_int_clear(i)      do {                                    \
  49                                         mpz_clear(i);                   \
  50                                         cloog_int_leak_down();          \
  51                                 } while (0)
  52
  53 #define cloog_int_set(r,i)      mpz_set(r,i)
  54 #define cloog_int_set_si(r,i)   mpz_set_si(r,i)
  55 #define cloog_int_abs(r,i)      mpz_abs(r,i)
  56 #define cloog_int_neg(r,i)      mpz_neg(r,i)
  57 #define cloog_int_swap(i,j)     mpz_swap(i,j)
  58
  59 #define cloog_int_add(r,i,j)    mpz_add(r,i,j)
  60 #define cloog_int_add_ui(r,i,j) mpz_add_ui(r,i,j)
  61 #define cloog_int_sub(r,i,j)    mpz_sub(r,i,j)
  62 #define cloog_int_sub_ui(r,i,j) mpz_sub_ui(r,i,j)
  63 #define cloog_int_mul(r,i,j)    mpz_mul(r,i,j)
  64 #define cloog_int_addmul(r,i,j) mpz_addmul(r,i,j)
  65 #define cloog_int_divexact(r,i,j)       mpz_divexact(r,i,j)
  66 #define cloog_int_tdiv_q(r,i,j) mpz_tdiv_q(r,i,j)
  67 #define cloog_int_fdiv_q(r,i,j) mpz_fdiv_q(r,i,j)
  68 #define cloog_int_fdiv_r(r,i,j) mpz_fdiv_r(r,i,j)
  69 #define cloog_int_gcd(r,i,j)    mpz_gcd(r,i,j)
  70
  71 #define cloog_int_sgn(i)        mpz_sgn(i)
  72 #define cloog_int_cmp(i,j)      mpz_cmp(i,j)
  73 #define cloog_int_cmp_si(i,si)  mpz_cmp_si(i,si)
  74 #define cloog_int_eq(i,j)       (mpz_cmp(i,j) == 0)
  75 #define cloog_int_ne(i,j)       (mpz_cmp(i,j) != 0)
  76 #define cloog_int_gt(i,j)       (mpz_cmp(i,j) > 0)
  77
  78 #define cloog_int_is_divisible_by(i,j)  mpz_divisible_p(i,j)
  79
  80 typedef void (*cloog_int_print_gmp_free_t)(void *, size_t);
  81 #define cloog_int_print(out,i)                                          \
  82         do {                                                            \
  83                 char *s;                                                \
  84                 cloog_int_print_gmp_free_t gmp_free;                    \
  85                 s = mpz_get_str(0, 10, i);                              \
  86                 fprintf(out, "%s", s);                                  \
  87                 mp_get_memory_functions(NULL, NULL, &gmp_free);         \
  88                 (*gmp_free)(s, strlen(s)+1);                            \
  89         } while (0)
  90
  91 #else
  92
  93 #define cloog_int_init(i)       ((i) = 0)
  94 #define cloog_int_clear(i)      do { } while (0)
  95
  96 #define cloog_int_set(r,i)      ((r) = (i))
  97 #define cloog_int_set_si(r,i)   ((r) = (i))
  98 #define cloog_int_abs(r,i)      ((r) = (i) > 0 ? (i) : -(i))
  99 #define cloog_int_neg(r,i)      ((r) = -(i))
 100 #define cloog_int_swap(i,j)     do {                                    \
 101                                         cloog_int_t _t = i;             \
 102                                         i = j;                          \
 103                                         j = _t;                         \
 104                                 } while (0)
 105
 106 #define cloog_int_add(r,i,j)    ((r) = (i) + (j))
 107 #define cloog_int_add_ui(r,i,j) ((r) = (i) + (j))
 108 #define cloog_int_sub(r,i,j)    ((r) = (i) - (j))
 109 #define cloog_int_sub_ui(r,i,j) ((r) = (i) - (j))
 110 #define cloog_int_mul(r,i,j)    ((r) = (i) * (j))
 111 #define cloog_int_addmul(r,i,j) ((r) += (i) * (j))
 112 #define cloog_int_divexact(r,i,j)       ((r) = (i) / (j))
 113 #define cloog_int_tdiv_q(r,i,j) ((r) = (i) / (j))
 114 #define cloog_int_fdiv_q(r,i,j) do {                                    \
 115                                         assert((j) > 0);                \
 116                                         (r) = (i) >= 0 ? (i) / (j) :    \
 117                                             -((-(i)+(j)-1)/(j));        \
 118                                 } while (0)
 119 #define cloog_int_fdiv_r(r,i,j) do {                                    \
 120                                         assert((j) > 0);                \
 121                                         (r) = (i) >= 0 ? (i) % (j) :    \
 122                                             (j)-1 - ((-(i)+(j)-1)%(j)); \
 123                                 } while (0)
 124 cloog_int_t cloog_gcd(cloog_int_t a, cloog_int_t b);
 125 #define cloog_int_gcd(r,i,j)    (r) = cloog_gcd(i,j)
 126
 127 #define cloog_int_sgn(i)        ((i) > 0 ? 1 : (i) < 0 ? -1 : 0)
 128 #define cloog_int_cmp(i,j)      (i - j)
 129 #define cloog_int_cmp_si(i,si)  (i - si)
 130 #define cloog_int_eq(i,j)       ((i) == (j))
 131 #define cloog_int_ne(i,j)       ((i) != (j))
 132 #define cloog_int_gt(i,j)       ((i) > (j))
 133
 134 #define cloog_int_is_divisible_by(i,j)  ((i) % (j) == 0)
 135
 136 #define cloog_int_print(out,i)  fprintf(out, CLOOG_INT_FORMAT, i)
 137
 138 #endif
 139
 140 #define cloog_int_is_pos(i)     (cloog_int_sgn(i) > 0)
 141 #define cloog_int_is_neg(i)     (cloog_int_sgn(i) < 0)
 142 #define cloog_int_is_zero(i)    (cloog_int_sgn(i) == 0)
 143 #define cloog_int_is_one(i)     (cloog_int_cmp_si(i,1) == 0)
 144 #define cloog_int_is_neg_one(i) (cloog_int_cmp_si(i,-1) == 0)
 145 #define cloog_int_gt_si(i,si)   (cloog_int_cmp_si(i,si) > 0)
 146
 147 struct cloog_vec {
 148         unsigned        size;
 149         cloog_int_t     *p;
 150 };
 151
 152 struct cloog_vec *cloog_vec_alloc(unsigned size);
 153 void cloog_vec_free(struct cloog_vec *vec);
 154 void cloog_seq_neg(cloog_int_t *dst, cloog_int_t *src, unsigned len);
 155 void cloog_seq_combine(cloog_int_t *dst, cloog_int_t m1, cloog_int_t *src1,
 156                         cloog_int_t m2, cloog_int_t *src2, unsigned len);
 157
 158 #if defined(__cplusplus)
 159   }
 160 #endif
 161
 162 #endif