1/* 2 * Copyright 2008-2009 Katholieke Universiteit Leuven 3 * 4 * Use of this software is governed by the MIT license 5 * 6 * Written by Sven Verdoolaege, K.U.Leuven, Departement 7 * Computerwetenschappen, Celestijnenlaan 200A, B-3001 Leuven, Belgium 8 */ 9 10#ifndef ISL_INT_H 11#define ISL_INT_H 12 13#include <isl/hash.h> 14#include <string.h> 15#include <gmp.h> 16#if defined(__cplusplus) 17#include <iostream> 18#endif 19 20#if defined(__cplusplus) 21extern "C" { 22#endif 23 24#ifndef mp_get_memory_functions 25void mp_get_memory_functions( 26 void *(**alloc_func_ptr) (size_t), 27 void *(**realloc_func_ptr) (void *, size_t, size_t), 28 void (**free_func_ptr) (void *, size_t)); 29#endif 30 31/* isl_int is the basic integer type. It currently always corresponds 32 * to a gmp mpz_t, but in the future, different types such as long long 33 * or cln::cl_I will be supported. 34 */ 35typedef mpz_t isl_int; 36 37#define isl_int_init(i) mpz_init(i) 38#define isl_int_clear(i) mpz_clear(i) 39 40#define isl_int_set(r,i) mpz_set(r,i) 41#define isl_int_set_gmp(r,i) mpz_set(r,i) 42#define isl_int_set_si(r,i) mpz_set_si(r,i) 43#define isl_int_set_ui(r,i) mpz_set_ui(r,i) 44#define isl_int_get_gmp(i,g) mpz_set(g,i) 45#define isl_int_get_si(r) mpz_get_si(r) 46#define isl_int_get_ui(r) mpz_get_ui(r) 47#define isl_int_get_d(r) mpz_get_d(r) 48#define isl_int_get_str(r) mpz_get_str(0, 10, r) 49typedef void (*isl_int_print_gmp_free_t)(void *, size_t); 50#define isl_int_free_str(s) \ 51 do { \ 52 isl_int_print_gmp_free_t gmp_free; \ 53 mp_get_memory_functions(NULL, NULL, &gmp_free); \ 54 (*gmp_free)(s, strlen(s) + 1); \ 55 } while (0) 56#define isl_int_abs(r,i) mpz_abs(r,i) 57#define isl_int_neg(r,i) mpz_neg(r,i) 58#define isl_int_swap(i,j) mpz_swap(i,j) 59#define isl_int_swap_or_set(i,j) mpz_swap(i,j) 60#define isl_int_add_ui(r,i,j) mpz_add_ui(r,i,j) 61#define isl_int_sub_ui(r,i,j) mpz_sub_ui(r,i,j) 62 63#define isl_int_add(r,i,j) mpz_add(r,i,j) 64#define isl_int_sub(r,i,j) mpz_sub(r,i,j) 65#define isl_int_mul(r,i,j) mpz_mul(r,i,j) 66#define isl_int_mul_2exp(r,i,j) mpz_mul_2exp(r,i,j) 67#define isl_int_mul_ui(r,i,j) mpz_mul_ui(r,i,j) 68#define isl_int_pow_ui(r,i,j) mpz_pow_ui(r,i,j) 69#define isl_int_addmul(r,i,j) mpz_addmul(r,i,j) 70#define isl_int_submul(r,i,j) mpz_submul(r,i,j) 71 72#define isl_int_gcd(r,i,j) mpz_gcd(r,i,j) 73#define isl_int_lcm(r,i,j) mpz_lcm(r,i,j) 74#define isl_int_divexact(r,i,j) mpz_divexact(r,i,j) 75#define isl_int_divexact_ui(r,i,j) mpz_divexact_ui(r,i,j) 76#define isl_int_tdiv_q(r,i,j) mpz_tdiv_q(r,i,j) 77#define isl_int_cdiv_q(r,i,j) mpz_cdiv_q(r,i,j) 78#define isl_int_fdiv_q(r,i,j) mpz_fdiv_q(r,i,j) 79#define isl_int_fdiv_r(r,i,j) mpz_fdiv_r(r,i,j) 80#define isl_int_fdiv_q_ui(r,i,j) mpz_fdiv_q_ui(r,i,j) 81 82#define isl_int_read(r,s) mpz_set_str(r,s,10) 83#define isl_int_print(out,i,width) \ 84 do { \ 85 char *s; \ 86 s = mpz_get_str(0, 10, i); \ 87 fprintf(out, "%*s", width, s); \ 88 isl_int_free_str(s); \ 89 } while (0) 90 91#define isl_int_sgn(i) mpz_sgn(i) 92#define isl_int_cmp(i,j) mpz_cmp(i,j) 93#define isl_int_cmp_si(i,si) mpz_cmp_si(i,si) 94#define isl_int_eq(i,j) (mpz_cmp(i,j) == 0) 95#define isl_int_ne(i,j) (mpz_cmp(i,j) != 0) 96#define isl_int_lt(i,j) (mpz_cmp(i,j) < 0) 97#define isl_int_le(i,j) (mpz_cmp(i,j) <= 0) 98#define isl_int_gt(i,j) (mpz_cmp(i,j) > 0) 99#define isl_int_ge(i,j) (mpz_cmp(i,j) >= 0) 100#define isl_int_abs_eq(i,j) (mpz_cmpabs(i,j) == 0) 101#define isl_int_abs_ne(i,j) (mpz_cmpabs(i,j) != 0) 102#define isl_int_abs_lt(i,j) (mpz_cmpabs(i,j) < 0) 103#define isl_int_abs_gt(i,j) (mpz_cmpabs(i,j) > 0) 104#define isl_int_abs_ge(i,j) (mpz_cmpabs(i,j) >= 0) 105 106 107#define isl_int_is_zero(i) (isl_int_sgn(i) == 0) 108#define isl_int_is_one(i) (isl_int_cmp_si(i,1) == 0) 109#define isl_int_is_negone(i) (isl_int_cmp_si(i,-1) == 0) 110#define isl_int_is_pos(i) (isl_int_sgn(i) > 0) 111#define isl_int_is_neg(i) (isl_int_sgn(i) < 0) 112#define isl_int_is_nonpos(i) (isl_int_sgn(i) <= 0) 113#define isl_int_is_nonneg(i) (isl_int_sgn(i) >= 0) 114#define isl_int_is_divisible_by(i,j) mpz_divisible_p(i,j) 115 116uint32_t isl_gmp_hash(mpz_t v, uint32_t hash); 117#define isl_int_hash(v,h) isl_gmp_hash(v,h) 118 119#if defined(__cplusplus) 120} 121#endif 122 123#if defined(__cplusplus) 124extern "C" { typedef void (*isl_gmp_free_t)(void *, size_t); } 125 126static inline std::ostream &operator<<(std::ostream &os, isl_int i) 127{ 128 char *s; 129 s = mpz_get_str(0, 10, i); 130 os << s; 131 isl_int_free_str(s); 132 return os; 133} 134#endif 135 136#endif 137