src/mmx_gmp.cpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : mmx_gmp.cpp
* DESCRIPTION: Additional code for interface to gmp and mpfr
* COPYRIGHT  : (C) 2003  Joris van der Hoeven
*******************************************************************************
* This software falls under the GNU general public license and comes WITHOUT
* ANY WARRANTY WHATSOEVER. See the file $TEXMACS_PATH/LICENSE for more details.
* If you don't have this file, write to the Free Software Foundation, Inc.,
* 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
******************************************************************************/

#include <string.h>
#include <numerix/rational.hpp>
#include <basix/evaluator.hpp>
#include <basix/literal.hpp>
namespace mmx {

string
as_string (const integer& i) {
  char* s= mpz_get_str (NULL, 10, *i);
  string r= s;
  mmx_free (s, strlen (s) + 1);
  /* FIXME: for upcoming version of GMP
  void (*free_func) (void *, size_t);
  mp_get_memory_functions (NULL, NULL, &free_func);
  free_func ((void*) s, strlen (s) + 1);
  */
  return r;
}

generic
construct (const integer& i) {
  return construct (as<generic> (i));
}

generic
construct (const rational& x) {
  return construct (as<generic> (x));
}

syntactic
flatten (const integer& i) {
  return syntactic (as_string (i));
}

syntactic
flatten (const rational& x) {
  if (denominator (x) == 1) return flatten (numerator (x));
  return flatten (numerator (x)) / flatten (denominator (x));
}

void
mpz_binary_write (const port& p, const mpz_t arg) {
  int n= arg->_mp_size;
  nat a= (nat) (n<0? -n: n);
  binary_write<int> (p, n);
  for (nat i=0; i<a; i++)
    binary_write<mp_limb_t> (p, arg->_mp_d[i]);
}

void
mpz_binary_read (const port& p, mpz_t dest) {
  int n= binary_read<int> (p);
  nat a= (nat) (n<0? -n: n);
  (void) _mpz_realloc (dest, a);
  for (nat i=0; i<a; i++)
    dest->_mp_d[i]= binary_read<mp_limb_t> (p);
  dest->_mp_size= n;
}

void
numeric_to_string (const integer& x, string& s) {
  s= as_string (x);
}

bool
string_to_numeric (const string& s, integer& val) {
  val= integer (s);
  return is_integer_string (s);
}

bool
is_literal_integer (const literal& l) {
  return is_integer_string (*l);
}

integer
as_literal_integer (const literal& l) {
  return integer (*l);
}

} // namespace mmx

---