On Mon, 18 Jun 2012, Zimmermann Paul wrote:
I am contributing the following function:long double mpz_get_ld (mpz_t z) { double h, l; mpz_t tmp; h = mpz_get_d (z); mpz_init (tmp); mpz_set_d (tmp, h); mpz_sub (tmp, tmp, z); l = mpz_get_d (tmp); mpz_clear (tmp); return (long double) h + (long double) l; }
If the function mpz_get_ld is provided, I expect it to work for the whole (exponent) range of long double value, not just the range of double. Should we first do something about large z, then apply your method and fix the result with ldexpl?
-- Marc Glisse _______________________________________________ gmp-devel mailing list gmp-devel@gmplib.org http://gmplib.org/mailman/listinfo/gmp-devel
