isl_test_int.c

Go to the documentation of this file.

/*
 * Copyright 2015 INRIA Paris-Rocquencourt
 *
 * Use of this software is governed by the MIT license
 *
 * Written by Michael Kruse, INRIA Paris-Rocquencourt,
 * Domaine de Voluceau, Rocquenqourt, B.P. 105,
 * 78153 Le Chesnay Cedex France
 */
 
#include <assert.h>
#include <stdio.h>
#include <isl_int.h>
 
#define ARRAY_SIZE(array) (sizeof(array)/sizeof(*array))
 
#ifdef USE_SMALL_INT_OPT
/* Test whether small and big representation of the same number have the same
 * hash.
 */
static void int_test_hash(isl_int val)
{
  uint32_t demotedhash, promotedhash;
  isl_int demoted, promoted;
 
  isl_int_init(demoted);
  isl_int_set(demoted, val);
 
  isl_int_init(promoted);
  isl_int_set(promoted, val);
 
  isl_sioimath_try_demote(demoted);
  isl_sioimath_promote(promoted);
 
  assert(isl_int_eq(demoted, promoted));
 
  demotedhash = isl_int_hash(demoted, 0);
  promotedhash = isl_int_hash(promoted, 0);
  assert(demotedhash == promotedhash);
 
  isl_int_clear(demoted);
  isl_int_clear(promoted);
}
 
struct {
  void (*fn)(isl_int);
  char *val;
} int_single_value_tests[] = {
  { &int_test_hash, "0" },
  { &int_test_hash, "1" },
  { &int_test_hash, "-1" },
  { &int_test_hash, "23" },
  { &int_test_hash, "-23" },
  { &int_test_hash, "107" },
  { &int_test_hash, "32768" },
  { &int_test_hash, "2147483647" },
  { &int_test_hash, "-2147483647" },
  { &int_test_hash, "2147483648" },
  { &int_test_hash, "-2147483648" },
};
 
static void int_test_single_value()
{
  int i;
 
  for (i = 0; i < ARRAY_SIZE(int_single_value_tests); i += 1) {
    isl_int val;
 
    isl_int_init(val);
    isl_int_read(val, int_single_value_tests[i].val);
 
    (*int_single_value_tests[i].fn)(val);
 
    isl_int_clear(val);
  }
}
 
static void invoke_alternate_representations_2args(char *arg1, char *arg2,
  void (*fn)(isl_int, isl_int))
{
  int j;
  isl_int int1, int2;
 
  isl_int_init(int1);
  isl_int_init(int2);
 
  for (j = 0; j < 4; ++j) {
    isl_int_read(int1, arg1);
    isl_int_read(int2, arg2);
 
    if (j & 1)
      isl_sioimath_promote(int1);
    else
      isl_sioimath_try_demote(int1);
 
    if (j & 2)
      isl_sioimath_promote(int2);
    else
      isl_sioimath_try_demote(int2);
 
    (*fn)(int1, int2);
  }
 
  isl_int_clear(int1);
  isl_int_clear(int2);
}
 
static void invoke_alternate_representations_3args(char *arg1, char *arg2,
  char *arg3, void (*fn)(isl_int, isl_int, isl_int))
{
  int j;
  isl_int int1, int2, int3;
 
  isl_int_init(int1);
  isl_int_init(int2);
  isl_int_init(int3);
 
  for (j = 0; j < 8; ++j) {
    isl_int_read(int1, arg1);
    isl_int_read(int2, arg2);
    isl_int_read(int3, arg3);
 
    if (j & 1)
      isl_sioimath_promote(int1);
    else
      isl_sioimath_try_demote(int1);
 
    if (j & 2)
      isl_sioimath_promote(int2);
    else
      isl_sioimath_try_demote(int2);
 
    if (j & 4)
      isl_sioimath_promote(int3);
    else
      isl_sioimath_try_demote(int3);
 
    (*fn)(int1, int2, int3);
  }
 
  isl_int_clear(int1);
  isl_int_clear(int2);
  isl_int_clear(int3);
}
#else  /* USE_SMALL_INT_OPT */
 
static void int_test_single_value()
{
}
 
static void invoke_alternate_representations_2args(char *arg1, char *arg2,
  void (*fn)(isl_int, isl_int))
{
  isl_int int1, int2;
 
  isl_int_init(int1);
  isl_int_init(int2);
 
  isl_int_read(int1, arg1);
  isl_int_read(int2, arg2);
 
  (*fn)(int1, int2);
 
  isl_int_clear(int1);
  isl_int_clear(int2);
}
 
static void invoke_alternate_representations_3args(char *arg1, char *arg2,
  char *arg3, void (*fn)(isl_int, isl_int, isl_int))
{
  isl_int int1, int2, int3;
 
  isl_int_init(int1);
  isl_int_init(int2);
  isl_int_init(int3);
 
  isl_int_read(int1, arg1);
  isl_int_read(int2, arg2);
  isl_int_read(int3, arg3);
 
  (*fn)(int1, int2, int3);
 
  isl_int_clear(int1);
  isl_int_clear(int2);
  isl_int_clear(int3);
}
#endif /* USE_SMALL_INT_OPT */
 
static void int_test_neg(isl_int expected, isl_int arg)
{
  isl_int result;
  isl_int_init(result);
 
  isl_int_neg(result, arg);
  assert(isl_int_eq(result, expected));
 
  isl_int_neg(result, expected);
  assert(isl_int_eq(result, arg));
 
  isl_int_clear(result);
}
 
static void int_test_abs(isl_int expected, isl_int arg)
{
  isl_int result;
  isl_int_init(result);
 
  isl_int_abs(result, arg);
  assert(isl_int_eq(result, expected));
 
  isl_int_clear(result);
}
 
struct {
  void (*fn)(isl_int, isl_int);
  char *expected, *arg;
} int_unary_tests[] = {
  { &int_test_neg, "0", "0" },
  { &int_test_neg, "-1", "1" },
  { &int_test_neg, "-2147483647", "2147483647" },
  { &int_test_neg, "-2147483648", "2147483648" },
  { &int_test_neg, "-9223372036854775807", "9223372036854775807" },
  { &int_test_neg, "-9223372036854775808", "9223372036854775808" },
 
  { &int_test_abs, "0", "0" },
  { &int_test_abs, "1", "1" },
  { &int_test_abs, "1", "-1" },
  { &int_test_abs, "2147483647", "2147483647" },
  { &int_test_abs, "2147483648", "-2147483648" },
  { &int_test_abs, "9223372036854775807", "9223372036854775807" },
  { &int_test_abs, "9223372036854775808", "-9223372036854775808" },
};
 
static void int_test_divexact(isl_int expected, isl_int lhs, isl_int rhs)
{
  isl_int result;
  unsigned long rhsulong;
 
  if (isl_int_sgn(rhs) == 0)
    return;
 
  isl_int_init(result);
 
  isl_int_divexact(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_tdiv_q(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_fdiv_q(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_cdiv_q(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  if (isl_int_fits_ulong(rhs)) {
    rhsulong = isl_int_get_ui(rhs);
 
    isl_int_divexact_ui(result, lhs, rhsulong);
    assert(isl_int_eq(expected, result));
 
    isl_int_fdiv_q_ui(result, lhs, rhsulong);
    assert(isl_int_eq(expected, result));
 
    isl_int_cdiv_q_ui(result, lhs, rhsulong);
    assert(isl_int_eq(expected, result));
  }
 
  isl_int_clear(result);
}
 
static void int_test_mul(isl_int expected, isl_int lhs, isl_int rhs)
{
  isl_int result;
  isl_int_init(result);
 
  isl_int_mul(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  if (isl_int_fits_ulong(rhs)) {
    unsigned long rhsulong = isl_int_get_ui(rhs);
 
    isl_int_mul_ui(result, lhs, rhsulong);
    assert(isl_int_eq(expected, result));
  }
 
  if (isl_int_fits_slong(rhs)) {
    unsigned long rhsslong = isl_int_get_si(rhs);
 
    isl_int_mul_si(result, lhs, rhsslong);
    assert(isl_int_eq(expected, result));
  }
 
  isl_int_clear(result);
}
 
/* Use a triple that satisfies 'product = factor1 * factor2' to check the
 * operations mul, divexact, tdiv, fdiv and cdiv.
 */
static void int_test_product(isl_int product, isl_int factor1, isl_int factor2)
{
  int_test_divexact(factor1, product, factor2);
  int_test_divexact(factor2, product, factor1);
 
  int_test_mul(product, factor1, factor2);
  int_test_mul(product, factor2, factor1);
}
 
static void int_test_add(isl_int expected, isl_int lhs, isl_int rhs)
{
  isl_int result;
  isl_int_init(result);
 
  isl_int_add(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_clear(result);
}
 
static void int_test_sub(isl_int expected, isl_int lhs, isl_int rhs)
{
  isl_int result;
  isl_int_init(result);
 
  isl_int_sub(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_clear(result);
}
 
/* Use a triple that satisfies 'sum = term1 + term2' to check the operations add
 * and sub.
 */
static void int_test_sum(isl_int sum, isl_int term1, isl_int term2)
{
  int_test_sub(term1, sum, term2);
  int_test_sub(term2, sum, term1);
 
  int_test_add(sum, term1, term2);
  int_test_add(sum, term2, term1);
}
 
static void int_test_fdiv(isl_int expected, isl_int lhs, isl_int rhs)
{
  unsigned long rhsulong;
  isl_int result;
  isl_int_init(result);
 
  isl_int_fdiv_q(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  if (isl_int_fits_ulong(rhs)) {
    rhsulong = isl_int_get_ui(rhs);
 
    isl_int_fdiv_q_ui(result, lhs, rhsulong);
    assert(isl_int_eq(expected, result));
  }
 
  isl_int_clear(result);
}
 
static void int_test_cdiv(isl_int expected, isl_int lhs, isl_int rhs)
{
  unsigned long rhsulong;
  isl_int result;
  isl_int_init(result);
 
  isl_int_cdiv_q(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  if (isl_int_fits_ulong(rhs)) {
    rhsulong = isl_int_get_ui(rhs);
 
    isl_int_cdiv_q_ui(result, lhs, rhsulong);
    assert(isl_int_eq(expected, result));
  }
 
  isl_int_clear(result);
}
 
static void int_test_tdiv(isl_int expected, isl_int lhs, isl_int rhs)
{
  isl_int result;
  isl_int_init(result);
 
  isl_int_tdiv_q(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_clear(result);
}
 
static void int_test_fdiv_r(isl_int expected, isl_int lhs, isl_int rhs)
{
  isl_int result;
  isl_int_init(result);
 
  isl_int_fdiv_r(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_clear(result);
}
 
static void int_test_gcd(isl_int expected, isl_int lhs, isl_int rhs)
{
  isl_int result;
  isl_int_init(result);
 
  isl_int_gcd(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_gcd(result, rhs, lhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_clear(result);
}
 
static void int_test_lcm(isl_int expected, isl_int lhs, isl_int rhs)
{
  isl_int result;
  isl_int_init(result);
 
  isl_int_lcm(result, lhs, rhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_lcm(result, rhs, lhs);
  assert(isl_int_eq(expected, result));
 
  isl_int_clear(result);
}
 
static int sgn(int val)
{
  if (val > 0)
    return 1;
  if (val < 0)
    return -1;
  return 0;
}
 
static void int_test_cmp(int exp, isl_int lhs, isl_int rhs)
{
  long rhslong;
 
  assert(exp == sgn(isl_int_cmp(lhs, rhs)));
 
  if (isl_int_fits_slong(rhs)) {
    rhslong = isl_int_get_si(rhs);
    assert(exp == sgn(isl_int_cmp_si(lhs, rhslong)));
  }
}
 
/* Test the comparison relations over two numbers.
 * expected is the sign (1, 0 or -1) of 'lhs - rhs'.
 */
static void int_test_cmps(isl_int expected, isl_int lhs, isl_int rhs)
{
  int exp;
  isl_int diff;
 
  exp = isl_int_get_si(expected);
 
  isl_int_init(diff);
  isl_int_sub(diff, lhs, rhs);
  assert(exp == isl_int_sgn(diff));
  isl_int_clear(diff);
 
  int_test_cmp(exp, lhs, rhs);
  int_test_cmp(-exp, rhs, lhs);
}
 
static void int_test_abs_cmp(isl_int expected, isl_int lhs, isl_int rhs)
{
  int exp;
 
  exp = isl_int_get_si(expected);
  assert(exp == sgn(isl_int_abs_cmp(lhs, rhs)));
  assert(-exp == sgn(isl_int_abs_cmp(rhs, lhs)));
}
 
/* If "expected" is equal to 1, then check that "rhs" divides "lhs".
 * If "expected" is equal to 0, then check that "rhs" does not divide "lhs".
 */
static void int_test_divisible(isl_int expected, isl_int lhs, isl_int rhs)
{
  int exp;
 
  exp = isl_int_get_si(expected);
  assert(isl_int_is_divisible_by(lhs, rhs) == exp);
}
 
struct {
  void (*fn)(isl_int, isl_int, isl_int);
  char *expected, *lhs, *rhs;
} int_binary_tests[] = {
  { &int_test_sum, "0", "0", "0" },
  { &int_test_sum, "1", "1", "0" },
  { &int_test_sum, "2", "1", "1" },
  { &int_test_sum, "-1", "0", "-1" },
  { &int_test_sum, "-2", "-1", "-1" },
 
  { &int_test_sum, "2147483647", "1073741823", "1073741824" },
  { &int_test_sum, "-2147483648", "-1073741824", "-1073741824" },
 
  { &int_test_sum, "2147483648", "2147483647", "1" },
  { &int_test_sum, "-2147483648", "-2147483647", "-1" },
 
  { &int_test_product, "0", "0", "0" },
  { &int_test_product, "0", "0", "1" },
  { &int_test_product, "1", "1", "1" },
 
  { &int_test_product, "6", "2", "3" },
  { &int_test_product, "-6", "2", "-3" },
  { &int_test_product, "-6", "-2", "3" },
  { &int_test_product, "6", "-2", "-3" },
 
  { &int_test_product, "2147483648", "65536", "32768" },
  { &int_test_product, "-2147483648", "65536", "-32768" },
 
  { &int_test_product,
    "4611686014132420609", "2147483647", "2147483647" },
  { &int_test_product,
    "-4611686014132420609", "-2147483647", "2147483647" },
 
  { &int_test_product,
    "4611686016279904256", "2147483647", "2147483648" },
  { &int_test_product,
    "-4611686016279904256", "-2147483647", "2147483648" },
  { &int_test_product,
    "-4611686016279904256", "2147483647", "-2147483648" },
  { &int_test_product,
    "4611686016279904256", "-2147483647", "-2147483648" },
 
  { &int_test_product, "85070591730234615847396907784232501249",
    "9223372036854775807", "9223372036854775807" },
  { &int_test_product, "-85070591730234615847396907784232501249",
    "-9223372036854775807", "9223372036854775807" },
 
  { &int_test_product, "85070591730234615856620279821087277056",
    "9223372036854775807", "9223372036854775808" },
  { &int_test_product, "-85070591730234615856620279821087277056",
    "-9223372036854775807", "9223372036854775808" },
  { &int_test_product, "-85070591730234615856620279821087277056",
    "9223372036854775807", "-9223372036854775808" },
  { &int_test_product, "85070591730234615856620279821087277056",
    "-9223372036854775807", "-9223372036854775808" },
 
  { &int_test_product, "340282366920938463426481119284349108225",
    "18446744073709551615", "18446744073709551615" },
  { &int_test_product, "-340282366920938463426481119284349108225",
    "-18446744073709551615", "18446744073709551615" },
 
  { &int_test_product, "340282366920938463444927863358058659840",
    "18446744073709551615", "18446744073709551616" },
  { &int_test_product, "-340282366920938463444927863358058659840",
    "-18446744073709551615", "18446744073709551616" },
  { &int_test_product, "-340282366920938463444927863358058659840",
    "18446744073709551615", "-18446744073709551616" },
  { &int_test_product, "340282366920938463444927863358058659840",
    "-18446744073709551615", "-18446744073709551616" },
 
  { &int_test_fdiv, "0", "1", "2" },
  { &int_test_fdiv_r, "1", "1", "3" },
  { &int_test_fdiv, "-1", "-1", "2" },
  { &int_test_fdiv_r, "2", "-1", "3" },
  { &int_test_fdiv, "-1", "1", "-2" },
  { &int_test_fdiv_r, "-2", "1", "-3" },
  { &int_test_fdiv, "0", "-1", "-2" },
  { &int_test_fdiv_r, "-1", "-1", "-3" },
 
  { &int_test_cdiv, "1", "1", "2" },
  { &int_test_cdiv, "0", "-1", "2" },
  { &int_test_cdiv, "0", "1", "-2" },
  { &int_test_cdiv, "1", "-1", "-2" },
 
  { &int_test_cdiv, "1073741824", "2147483647", "2" },
  { &int_test_cdiv, "1073741824", "2147483648", "2" },
  { &int_test_cdiv, "-1073741824", "-2147483648", "2" },
  { &int_test_cdiv, "-1073741823", "-2147483647", "2" },
 
  { &int_test_tdiv, "0", "1", "2" },
  { &int_test_tdiv, "0", "-1", "2" },
  { &int_test_tdiv, "0", "1", "-2" },
  { &int_test_tdiv, "0", "-1", "-2" },
 
  { &int_test_gcd, "0", "0", "0" },
  { &int_test_lcm, "0", "0", "0" },
  { &int_test_gcd, "7", "0", "7" },
  { &int_test_lcm, "0", "0", "7" },
  { &int_test_gcd, "1", "1", "1" },
  { &int_test_lcm, "1", "1", "1" },
  { &int_test_gcd, "1", "1", "-1" },
  { &int_test_lcm, "1", "1", "-1" },
  { &int_test_gcd, "1", "-1", "-1" },
  { &int_test_lcm, "1", "-1", "-1" },
  { &int_test_gcd, "3", "6", "9" },
  { &int_test_lcm, "18", "6", "9" },
  { &int_test_gcd, "1", "14", "2147483647" },
  { &int_test_lcm, "15032385529", "7", "2147483647" },
  { &int_test_gcd, "2", "6", "-2147483648" },
  { &int_test_lcm, "6442450944", "6", "-2147483648" },
  { &int_test_gcd, "1", "6", "9223372036854775807" },
  { &int_test_lcm, "55340232221128654842", "6", "9223372036854775807" },
  { &int_test_gcd, "2", "6", "-9223372036854775808" },
  { &int_test_lcm, "27670116110564327424", "6", "-9223372036854775808" },
  { &int_test_gcd, "1", "18446744073709551616", "18446744073709551615" },
  { &int_test_lcm, "340282366920938463444927863358058659840",
    "18446744073709551616", "18446744073709551615" },
 
  { &int_test_cmps, "0", "0", "0" },
  { &int_test_abs_cmp, "0", "0", "0" },
  { &int_test_cmps, "1", "1", "0" },
  { &int_test_abs_cmp, "1", "1", "0" },
  { &int_test_cmps, "-1", "-1", "0" },
  { &int_test_abs_cmp, "1", "-1", "0" },
  { &int_test_cmps, "-1", "-1", "1" },
  { &int_test_abs_cmp, "0", "-1", "1" },
 
  { &int_test_cmps, "-1", "5", "2147483647" },
  { &int_test_abs_cmp, "-1", "5", "2147483647" },
  { &int_test_cmps, "1", "5", "-2147483648" },
  { &int_test_abs_cmp, "-1", "5", "-2147483648" },
  { &int_test_cmps, "-1", "5", "9223372036854775807" },
  { &int_test_abs_cmp, "-1", "5", "9223372036854775807" },
  { &int_test_cmps, "1", "5", "-9223372036854775809" },
  { &int_test_abs_cmp, "-1", "5", "-9223372036854775809" },
 
  { &int_test_divisible, "1", "0", "0" },
  { &int_test_divisible, "0", "1", "0" },
  { &int_test_divisible, "0", "2", "0" },
  { &int_test_divisible, "0", "2147483647", "0" },
  { &int_test_divisible, "0", "9223372036854775807", "0" },
  { &int_test_divisible, "1", "0", "1" },
  { &int_test_divisible, "1", "1", "1" },
  { &int_test_divisible, "1", "2", "1" },
  { &int_test_divisible, "1", "2147483647", "1" },
  { &int_test_divisible, "1", "9223372036854775807", "1" },
  { &int_test_divisible, "1", "0", "2" },
  { &int_test_divisible, "0", "1", "2" },
  { &int_test_divisible, "1", "2", "2" },
  { &int_test_divisible, "0", "2147483647", "2" },
  { &int_test_divisible, "0", "9223372036854775807", "2" },
};
 
/* Tests the isl_int_* function to give the expected results. Tests are
 * grouped by the number of arguments they take.
 *
 * If small integer optimization is enabled, we also test whether the results
 * are the same in small and big representation.
 */
int main()
{
  int i;
 
  int_test_single_value();
 
  for (i = 0; i < ARRAY_SIZE(int_unary_tests); i += 1) {
    invoke_alternate_representations_2args(
        int_unary_tests[i].expected, int_unary_tests[i].arg,
        int_unary_tests[i].fn);
  }
 
  for (i = 0; i < ARRAY_SIZE(int_binary_tests); i += 1) {
    invoke_alternate_representations_3args(
        int_binary_tests[i].expected, int_binary_tests[i].lhs,
        int_binary_tests[i].rhs, int_binary_tests[i].fn);
  }
 
  return 0;
}