
    /*+-----------------------------------------------------------------**
     **                       OpenScop Library                          **
     **-----------------------------------------------------------------**
     **                             int.c                               **
     **-----------------------------------------------------------------**
     **                   First version: 18/07/2011                     **
     **-----------------------------------------------------------------**

 
 *****************************************************************************
 * OpenScop: Structures and formats for polyhedral tools to talk together    *
 *****************************************************************************
 *    ,___,,_,__,,__,,__,,__,,_,__,,_,__,,__,,___,_,__,,_,__,                *
 *    /   / /  //  //  //  // /   / /  //  //   / /  // /  /|,_,             *
 *   /   / /  //  //  //  // /   / /  //  //   / /  // /  / / /\             *
 *  |~~~|~|~~~|~~~|~~~|~~~|~|~~~|~|~~~|~~~|~~~|~|~~~|~|~~~|/_/  \            *
 *  | G |C| P | = | L | P |=| = |C| = | = | = |=| = |=| C |\  \ /\           *
 *  | R |l| o | = | e | l |=| = |a| = | = | = |=| = |=| L | \# \ /\          *
 *  | A |a| l | = | t | u |=| = |n| = | = | = |=| = |=| o | |\# \  \         *
 *  | P |n| l | = | s | t |=| = |d| = | = | = | |   |=| o | | \# \  \        *
 *  | H | | y |   | e | o | | = |l|   |   | = | |   | | G | |  \  \  \       *
 *  | I | |   |   | e |   | |   | |   |   |   | |   | |   | |   \  \  \      *
 *  | T | |   |   |   |   | |   | |   |   |   | |   | |   | |    \  \  \     *
 *  | E | |   |   |   |   | |   | |   |   |   | |   | |   | |     \  \  \    *
 *  | * |*| * | * | * | * |*| * |*| * | * | * |*| * |*| * | /      \* \  \   *
 *  | O |p| e | n | S | c |o| p |-| L | i | b |r| a |r| y |/        \  \ /   *
 *  '---'-'---'---'---'---'-'---'-'---'---'---'-'---'-'---'          '--'    *
 *                                                                           *
 * Copyright (C) 2008 University Paris-Sud 11 and INRIA                      *
 *                                                                           *
 * (3-clause BSD license)                                                    *
 * Redistribution and use in source  and binary forms, with or without       *
 * modification, are permitted provided that the following conditions        *
 * are met:                                                                  *
 *                                                                           *
 * 1. Redistributions of source code must retain the above copyright notice, *
 *    this list of conditions and the following disclaimer.                  *
 * 2. Redistributions in binary form must reproduce the above copyright      *
 *    notice, this list of conditions and the following disclaimer in the    *
 *    documentation and/or other materials provided with the distribution.   *
 * 3. The name of the author may not be used to endorse or promote products  *
 *    derived from this software without specific prior written permission.  *
 *                                                                           *
 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR      *
 * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES *
 * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.   *
 * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,          *
 * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT  *
 * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, *
 * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY     *
 * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT       *
 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF  *
 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.         *
 *                                                                           *
 * OpenScop Library, a library to manipulate OpenScop formats and data       *
 * structures. Written by:                                                   *
 * Cedric Bastoul     <Cedric.Bastoul@u-psud.fr> and                         *
 * Louis-Noel Pouchet <Louis-Noel.pouchet@inria.fr>                          *
 *                                                                           *
 *****************************************************************************/

#include <stdlib.h>
#include <stdio.h>
#ifdef OSL_GMP_IS_HERE
# include <gmp.h>
#endif

#include <osl/macros.h>
#include <osl/int.h>


/*+***************************************************************************
 *                                Basic Functions                            *
 *****************************************************************************/


/**
 * osl_int_dump_precision function:
 * this function prints in a human readable fashion the precision
 * corresponding to the "precision" parameter.
 * \param[in] file      The file where to print the precision.
 * \param[in] precision The precision to print.
 */
void osl_int_dump_precision(FILE * file, int precision) {

  switch (precision) {
    case OSL_PRECISION_SP:
      fprintf(file, "32 bits");
      break;
    case OSL_PRECISION_DP:
      fprintf(file, "64 bits");
      break;
#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      fprintf(file, "GMP");
      break;
#endif
    default:
      fprintf(file, "unknown precision %d", precision);
  }
}


int osl_int_sizeof(int precision) {
  switch (precision) {
    case OSL_PRECISION_SP:
      return sizeof(long int);

    case OSL_PRECISION_DP:
      return sizeof(long long int);

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return sizeof(mpz_t);
#endif

    default:
      OSL_error("unknown precision");
  }
}


void * osl_int_address(int precision, void * base, int offset) {
  switch (precision) {
    case OSL_PRECISION_SP:
      return (long int *)base + offset;

    case OSL_PRECISION_DP:
      return (long long int *)base + offset;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return (mpz_t *)base + offset;
#endif

    default:
      OSL_error("unknown precision");
  }
}


void osl_int_init(int precision, void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)value = 0;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)value = 0;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_init(*(mpz_t *)value);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


void * osl_int_malloc(int precision) {
  void * value;

  switch (precision) {
    case OSL_PRECISION_SP:
      value = malloc(sizeof(long int));
      break;

    case OSL_PRECISION_DP:
      value = malloc(sizeof(long long int));
      *(long long int *)value = 0;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      value = malloc(sizeof(mpz_t));
      break;
#endif

    default:
      OSL_error("unknown precision");
  }

  osl_int_init(precision, value, 0);
  return value;
}


/**
 * val1_base[val1_offset] = val2_base[val2_offset];
 */
void osl_int_assign(int precision,
                    void * val1_base, int val1_offset,
                    void * val2_base, int val2_offset) {
  void * val1 = osl_int_address(precision, val1_base, val1_offset);
  void * val2 = osl_int_address(precision, val2_base, val2_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)val1 = *(long int *)val2;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)val1 = *(long long int *)val2;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_set(*(mpz_t *)val1, *(mpz_t *)val2);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * value_base[value_offset] = i;
 */
void osl_int_set_si(int precision, void * value_base, int value_offset,
                    int i) {
  void * value = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)value = (long int)i;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)value = (long long int)i;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_set_si(*(mpz_t *)value, i);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * return value_base[value_offset];
 */
int osl_int_get_si(int precision, void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      return *(int *)value;

    case OSL_PRECISION_DP:
      return *(int *)value;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return mpz_get_si(*(mpz_t *)value);
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * value_base[value_offset] = i; // including initialization for GMP
 */
void osl_int_init_set_si(int precision,
                         void * value_base, int value_offset, int i) {
  void * value = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)value = (long int)i;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)value = (long long int)i;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_init_set_si(*(mpz_t *)value, i);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * value_base[value_offset] = 0; // Including cleaning for GMP
 */
void osl_int_clear(int precision, void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);
  
  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)value = 0;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)value = 0;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_clear(*(mpz_t *)value);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


void osl_int_free(int precision, void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);

  osl_int_clear(precision, value_base, value_offset);
  free(value);
}


/**
 * osl_int_print function:
 * this function displays an integer value into a file (file, possibly stdout).
 * \param file         The file where the integer has to be printed.
 * \param precision    The precision of the integer.
 * \param value_base   Address of the base integer value.
 * \param value_offset Offset in number of values from the base integer value.
 */
void osl_int_print(FILE * file, int precision,
                   void * value_base, int value_offset) {
  char string[OSL_MAX_STRING];
  
  osl_int_sprint(string, precision, value_base, value_offset);
  fprintf(file, "%s", string);
}


/**
 * osl_int_sprint function:
 * this function prints an integer value into a string, it uses the
 * OpenScop Library formats OSL_FMT_* to format the printing.
 * \param string       The string where the integer has to be printed.
 * \param precision    The precision of the integer.
 * \param value_base   Address of the base integer value.
 * \param value_offset Offset in number of values from the base integer value.
 */
void osl_int_sprint(char * string, int precision,
                    void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);
            
  switch (precision) {
    case OSL_PRECISION_SP:
      sprintf(string, OSL_FMT_SP, *(long int *)value);
      break;

    case OSL_PRECISION_DP:
      sprintf(string, OSL_FMT_DP, *(long long int *)value);
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP: {
      char * str;
      str = mpz_get_str(0, 10, *(mpz_t *)value); //TODO: 10 -> #define
      sprintf(string, OSL_FMT_MP, str);
      free(str);
      break;
    }
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * osl_int_sprint_txt function:
 * this function is similar to osl_int_sprintf but it prints the value
 * using OSL_TMT_TXT_* formats.
 * \see osl_int_sprintf
 */
void osl_int_sprint_txt(char * string, int precision,
                        void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);
            
  switch (precision) {
    case OSL_PRECISION_SP:
      sprintf(string, OSL_FMT_TXT_SP, *(long int *)value);
      break;

    case OSL_PRECISION_DP:
      sprintf(string, OSL_FMT_TXT_DP, *(long long int *)value);
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP: {
      char * str;
      str = mpz_get_str(0, 10, *(mpz_t *)value); //TODO: 10 -> #define
      sprintf(string, OSL_FMT_TXT_MP, str);
      free(str);
      break;
    }
#endif

    default:
      OSL_error("unknown precision");
  }
}


void osl_int_sread(char ** string, int precision,
                   void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);
  int nb_read = 0;

  switch (precision) {
    case OSL_PRECISION_SP:
      nb_read = sscanf(*string, OSL_FMT_TXT_SP, (long int *)value);
      if (nb_read == 0)
        OSL_error("failed to read an integer");
      break;

    case OSL_PRECISION_DP:
      nb_read = sscanf(*string, OSL_FMT_TXT_DP, (long long int *)value);
      if (nb_read == 0)
        OSL_error("failed to read an integer");
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP: {
      long long int tmp;
      nb_read = sscanf(*string, OSL_FMT_TXT_DP, &tmp);
      if (nb_read == 0)
        OSL_error("failed to read an integer");
      mpz_set_si(*(mpz_t *)value, tmp);
      break;
    }
#endif

    default:
      OSL_error("unknown precision");
  }

  // Update the position in the input string.
  *string = *string + nb_read;
}


/*+***************************************************************************
 *                            Arithmetic Operations                          *
 *****************************************************************************/


/**
 * result_base[result_offset] = value_base[value_offset] + 1;
 */
void osl_int_increment(int precision,
                       void * result_base, int result_offset,
                       void * value_base,  int value_offset) {
  void * result = osl_int_address(precision, result_base, result_offset);
  void * value  = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)result = *(long int *)value + (long int)1;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)result = *(long long int *)value + (long long int)1;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_add_ui(*(mpz_t *)result, *(mpz_t *)value, 1);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * result_base[result_offset] = value_base[value_offset] - 1;
 */
void osl_int_decrement(int precision,
                       void * result_base, int result_offset,
                       void * value_base,  int value_offset) {
  void * result = osl_int_address(precision, result_base, result_offset);
  void * value  = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)result = *(long int *)value - (long int)1;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)result = *(long long int *)value - (long long int)1;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP: {
      mpz_t one;
      mpz_init_set_si(one, 1);
      mpz_sub(*(mpz_t *)result, *(mpz_t *)value, one);
      mpz_clear(one);
      break;
    }
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * result_base[result_offset] = val1_base[val1_offset]+val2_base[val2_offset];
 */
void osl_int_add(int precision,
                 void * result_base, int result_offset,
                 void * val1_base,   int val1_offset,
                 void * val2_base,   int val2_offset) {
  void * result = osl_int_address(precision, result_base, result_offset);
  void * val1   = osl_int_address(precision, val1_base, val1_offset);
  void * val2   = osl_int_address(precision, val2_base, val2_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)result = *(long int *)val1 + *(long int *)val2;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)result = *(long long int *)val1 +
                                 *(long long int *)val2;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_add(*(mpz_t *)result, *(mpz_t *)val1, *(mpz_t *)val2);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * result_base[result_offset] = value_base[value_offset] + i;
 */
void osl_int_add_si(int precision,
                    void * result_base, int result_offset,
                    void * value_base,  int value_offset, int i) {
  void * result = osl_int_address(precision, result_base, result_offset);
  void * value  = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)result = *(long int *)value + (long int)i;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)result = *(long long int *)value + (long long int)i;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP: {
      mpz_t si;
      mpz_init_set_si(si, i);
      mpz_add(*(mpz_t *)result, *(mpz_t *)value, si);
      mpz_clear(si);
      break;
    }
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * result_base[result_offset] = val1_base[val1_offset]*val2_base[val2_offset];
 */
void osl_int_mul(int precision,
                 void * result_base, int result_offset,
                 void * val1_base,   int val1_offset,
                 void * val2_base,   int val2_offset) {
  void * result = osl_int_address(precision, result_base, result_offset);
  void * val1   = osl_int_address(precision, val1_base, val1_offset);
  void * val2   = osl_int_address(precision, val2_base, val2_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)result = *(long int *)val1 * *(long int *)val2;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)result = *(long long int *)val1 *
                                 *(long long int *)val2;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_mul(*(mpz_t *)result, *(mpz_t *)val1, *(mpz_t *)val2);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * result_base[result_offset] = value_base[value_offset] * i;
 */
void osl_int_mul_si(int precision,
                    void * result_base, int result_offset,
                    void * value_base,  int value_offset, int i) {
  void * result = osl_int_address(precision, result_base, result_offset);
  void * value  = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)result = *(long int *)value * (long int)i;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)result = *(long long int *)value * (long long int)i;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_mul_si(*(mpz_t *)result, *(mpz_t *)value, i);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * result_base[result_offset] = val1_base[val1_offset]-val2_base[val2_offset];
 */
void osl_int_sub(int precision,
                 void * result_base, int result_offset,
                 void * val1_base,   int val1_offset,
                 void * val2_base,   int val2_offset) {
  void * result = osl_int_address(precision, result_base, result_offset);
  void * val1   = osl_int_address(precision, val1_base, val1_offset);
  void * val2   = osl_int_address(precision, val2_base, val2_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)result = *(long int *)val1 - *(long int *)val2;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)result = *(long long int *)val1 -
                                 *(long long int *)val2;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_sub(*(mpz_t *)result, *(mpz_t *)val1, *(mpz_t *)val2);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * result_base[result_offset] = -value_base[value_offset];
 */
void osl_int_oppose(int precision,
                    void * result_base, int result_offset,
                    void * value_base,  int value_offset) {
  void * result = osl_int_address(precision, result_base, result_offset);
  void * value  = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      *(long int *)result = -*(long int *)value;
      break;

    case OSL_PRECISION_DP:
      *(long long int *)result = -*(long long int *)value;
      break;

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      mpz_neg(*(mpz_t *)result, *(mpz_t *)value);
      break;
#endif

    default:
      OSL_error("unknown precision");
  }
}


/*+***************************************************************************
 *                            Conditional Operations                         *
 *****************************************************************************/


/**
 * (val1_base[val1_offset] == val2_base[val2_offset])
 */
int osl_int_eq(int precision,
               void * val1_base, int val1_offset,
               void * val2_base, int val2_offset) {
  void * val1 = osl_int_address(precision, val1_base, val1_offset);
  void * val2 = osl_int_address(precision, val2_base, val2_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      return (*(long int *)val1 == *(long int *)val2);

    case OSL_PRECISION_DP:
      return (*(long long int *)val1 == *(long long int *)val2);

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return (mpz_cmp(*(mpz_t *)val1, *(mpz_t *)val2) == 0);
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * (val1_base[val1_offset] != val2_base[val2_offset])
 */
int osl_int_ne(int precision,
               void * val1_base, int val1_offset,
               void * val2_base, int val2_offset) {
  return !osl_int_eq(precision,
                          val1_base, val1_offset,
                          val2_base, val2_offset);
}


/**
 * (value_base[value_offset] > 0)
 */
int osl_int_pos(int precision, void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      return (*(long int *)value > 0);

    case OSL_PRECISION_DP:
      return (*(long long int *)value > 0);

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return (mpz_sgn(*(mpz_t *)value) > 0);
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * (value_base[value_offset] < 0)
 */
int osl_int_neg(int precision, void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      return (*(long int *)value < 0);

    case OSL_PRECISION_DP:
      return (*(long long int *)value < 0);

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return (mpz_sgn(*(mpz_t *)value) < 0);
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * (value_base[value_offset] == 0)
 */
int osl_int_zero(int precision, void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      return (*(long int *)value == 0);

    case OSL_PRECISION_DP:
      return (*(long long int *)value == 0);

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return (mpz_sgn(*(mpz_t *)value) == 0);
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * (value_base[value_offset] == 1)
 */
int osl_int_one(int precision, void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      return (*(long int *)value == (long int)1);

    case OSL_PRECISION_DP:
      return (*(long long int *)value == (long long int)1);

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return (mpz_cmp_si(*(mpz_t *)value, 1) == 0);
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * (value_base[value_offset] == -1)
 */
int osl_int_mone(int precision, void * value_base, int value_offset) {
  void * value = osl_int_address(precision, value_base, value_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      return (*(long int *)value == (long int)-1);

    case OSL_PRECISION_DP:
      return (*(long long int *)value == (long long int)-1);

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return (mpz_cmp_si(*(mpz_t *)value, -1) == 0);
#endif

    default:
      OSL_error("unknown precision");
  }
}


/**
 * ((val1_base[val1_offset] % val2_base[val2_offset]) == 0)
 */
int osl_int_divisible(int precision,
                      void * val1_base, int val1_offset,
                      void * val2_base, int val2_offset) {
  void * val1 = osl_int_address(precision, val1_base, val1_offset);
  void * val2 = osl_int_address(precision, val2_base, val2_offset);

  switch (precision) {
    case OSL_PRECISION_SP:
      return ((*(long int *)val1 % *(long int *)val2) == 0);

    case OSL_PRECISION_DP:
      return ((*(long long int *)val1 % *(long long int *)val2) == 0);

#ifdef OSL_GMP_IS_HERE
    case OSL_PRECISION_MP:
      return mpz_divisible_p(*(mpz_t *)val1, *(mpz_t *)val2);
#endif

    default:
      OSL_error("unknown precision");
  }
}