mpf_sub bug (essentially same as reldiff bug)

Kevin Ryde user42 at
Wed May 12 00:46:52 CEST 2004

abbott at writes:
>   mpf_sub(rd, a, b);   // should be about 2^(-129), but gives zero

Thanks, that's a bug.

I think this case of nearly equal values is meant to work, you can see
the loop in mpf/sub.c under ediff==0 which skips equal high limbs.
Have to check what's not happening.

All followups to gmp-bugs at please.

More information about the gmp-bugs mailing list