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#include <isl_config.h> 17 18#ifndef mp_get_memory_functions 19void mp_get_memory_functions( 20 void *(**alloc_func_ptr) (size_t), 21 void *(**realloc_func_ptr) (void *, size_t, size_t), 22 void (**free_func_ptr) (void *, size_t)); 23#endif 24 25/* isl_int is the basic integer type. It currently always corresponds 26 * to a gmp mpz_t, but in the future, different types such as long long 27 * or cln::cl_I will be supported. 28 */ 29typedef mpz_t isl_int; 30 31#define isl_int_init(i) mpz_init(i) 32#define isl_int_clear(i) mpz_clear(i) 33 34#define isl_int_set(r,i) mpz_set(r,i) 35#define isl_int_set_gmp(r,i) mpz_set(r,i) 36#define isl_int_set_si(r,i) mpz_set_si(r,i) 37#define isl_int_set_ui(r,i) mpz_set_ui(r,i) 38#define isl_int_get_gmp(i,g) mpz_set(g,i) 39#define isl_int_fits_slong(r) mpz_fits_slong_p(r) 40#define isl_int_get_si(r) mpz_get_si(r) 41#define isl_int_fits_ulong(r) mpz_fits_ulong_p(r) 42#define isl_int_get_ui(r) mpz_get_ui(r) 43#define isl_int_get_d(r) mpz_get_d(r) 44#define isl_int_get_str(r) mpz_get_str(0, 10, r) 45typedef void (*isl_int_print_gmp_free_t)(void *, size_t); 46#define isl_int_free_str(s) \ 47 do { \ 48 isl_int_print_gmp_free_t gmp_free; \ 49 mp_get_memory_functions(NULL, NULL, &gmp_free); \ 50 (*gmp_free)(s, strlen(s) + 1); \ 51 } while (0) 52#define isl_int_abs(r,i) mpz_abs(r,i) 53#define isl_int_neg(r,i) mpz_neg(r,i) 54#define isl_int_swap(i,j) mpz_swap(i,j) 55#define isl_int_swap_or_set(i,j) mpz_swap(i,j) 56#define isl_int_add_ui(r,i,j) mpz_add_ui(r,i,j) 57#define isl_int_sub_ui(r,i,j) mpz_sub_ui(r,i,j) 58 59#define isl_int_add(r,i,j) mpz_add(r,i,j) 60#define isl_int_sub(r,i,j) mpz_sub(r,i,j) 61#define isl_int_mul(r,i,j) mpz_mul(r,i,j) 62#define isl_int_mul_2exp(r,i,j) mpz_mul_2exp(r,i,j) 63#define isl_int_mul_si(r,i,j) mpz_mul_si(r,i,j) 64#define isl_int_mul_ui(r,i,j) mpz_mul_ui(r,i,j) 65#define isl_int_pow_ui(r,i,j) mpz_pow_ui(r,i,j) 66#define isl_int_addmul(r,i,j) mpz_addmul(r,i,j) 67#define isl_int_addmul_ui(r,i,j) mpz_addmul_ui(r,i,j) 68#define isl_int_submul(r,i,j) mpz_submul(r,i,j) 69#define isl_int_submul_ui(r,i,j) mpz_submul_ui(r,i,j) 70 71#define isl_int_gcd(r,i,j) mpz_gcd(r,i,j) 72#ifdef GMP_NORMALIZE_GCDEXT 73void isl_gmp_gcdext(mpz_t G, mpz_t S, mpz_t T, mpz_t A, mpz_t B); 74#define isl_int_gcdext(g,x,y,i,j) isl_gmp_gcdext(g,x,y,i,j) 75#else 76#define isl_int_gcdext(g,x,y,i,j) mpz_gcdext(g,x,y,i,j) 77#endif 78#define isl_int_lcm(r,i,j) mpz_lcm(r,i,j) 79#define isl_int_divexact(r,i,j) mpz_divexact(r,i,j) 80#define isl_int_divexact_ui(r,i,j) mpz_divexact_ui(r,i,j) 81#define isl_int_tdiv_q(r,i,j) mpz_tdiv_q(r,i,j) 82#define isl_int_cdiv_q(r,i,j) mpz_cdiv_q(r,i,j) 83#define isl_int_fdiv_q(r,i,j) mpz_fdiv_q(r,i,j) 84#define isl_int_fdiv_r(r,i,j) mpz_fdiv_r(r,i,j) 85#define isl_int_fdiv_q_ui(r,i,j) mpz_fdiv_q_ui(r,i,j) 86 87#define isl_int_read(r,s) mpz_set_str(r,s,10) 88#define isl_int_print(out,i,width) \ 89 do { \ 90 char *s; \ 91 s = mpz_get_str(0, 10, i); \ 92 fprintf(out, "%*s", width, s); \ 93 isl_int_free_str(s); \ 94 } while (0) 95 96#define isl_int_sgn(i) mpz_sgn(i) 97#define isl_int_cmp(i,j) mpz_cmp(i,j) 98#define isl_int_cmp_si(i,si) mpz_cmp_si(i,si) 99#define isl_int_eq(i,j) (mpz_cmp(i,j) == 0) 100#define isl_int_ne(i,j) (mpz_cmp(i,j) != 0) 101#define isl_int_lt(i,j) (mpz_cmp(i,j) < 0) 102#define isl_int_le(i,j) (mpz_cmp(i,j) <= 0) 103#define isl_int_gt(i,j) (mpz_cmp(i,j) > 0) 104#define isl_int_ge(i,j) (mpz_cmp(i,j) >= 0) 105#define isl_int_abs_eq(i,j) (mpz_cmpabs(i,j) == 0) 106#define isl_int_abs_ne(i,j) (mpz_cmpabs(i,j) != 0) 107#define isl_int_abs_lt(i,j) (mpz_cmpabs(i,j) < 0) 108#define isl_int_abs_gt(i,j) (mpz_cmpabs(i,j) > 0) 109#define isl_int_abs_ge(i,j) (mpz_cmpabs(i,j) >= 0) 110 111 112#define isl_int_is_zero(i) (isl_int_sgn(i) == 0) 113#define isl_int_is_one(i) (isl_int_cmp_si(i,1) == 0) 114#define isl_int_is_negone(i) (isl_int_cmp_si(i,-1) == 0) 115#define isl_int_is_pos(i) (isl_int_sgn(i) > 0) 116#define isl_int_is_neg(i) (isl_int_sgn(i) < 0) 117#define isl_int_is_nonpos(i) (isl_int_sgn(i) <= 0) 118#define isl_int_is_nonneg(i) (isl_int_sgn(i) >= 0) 119#define isl_int_is_divisible_by(i,j) mpz_divisible_p(i,j) 120 121uint32_t isl_gmp_hash(mpz_t v, uint32_t hash); 122#define isl_int_hash(v,h) isl_gmp_hash(v,h) 123 124#endif 125