C implementation of mod_1_1

Niels Möller nisse at lysator.liu.se
Tue Mar 1 13:40:39 CET 2011

nisse at lysator.liu.se (Niels Möller) writes:

> Could also be an error in the computations before mpq_get_d. I suspect
> the difference in exponent between "last" and "new" is unexpected (one
> could examine the printouts to see if its mpq_get_d gave a sensible
> result or not, I think, but I can't do that at the moment).

A quick check with pari's rational -> float conversion indicates that
"new" is a correct conversion of the rational A. 

So the bug is most likely trigged by the code adding EPS to A.


Niels Möller. PGP-encrypted email is preferred. Keyid C0B98E26.
Internet email is subject to wholesale government surveillance.

More information about the gmp-devel mailing list