mpf_sub bug (essentially same as reldiff bug)

Kevin Ryde user42 at zip.com.au
Wed May 12 00:46:52 CEST 2004


abbott at module.dima.unige.it 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 swox.com please.


More information about the gmp-bugs mailing list