C implementation of mod_1_1
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