mpf_sub bug (essentially same as reldiff bug)
Giulio Malinverno
giulio.malinverno at republic-spaceworks.org
Wed May 12 10:22:09 CEST 2004
If we use the same precision for all quantities:
changing
mpf_init2(rd, 8);
to
mpf_init2(rd, 256);
we obtain the correct result (both in C and C++ sources)
mantissa=100000000000000000000000000000000000000146936793852785938496
diff is approx -0.5 times 2^(-128)
a: 0.1e1
b:
0.1000000000000000000000000000000000000001469367938527859384960920671527807097273e1
rd:
-0.1469367938527859384960920671527807097273331945965109401885939632848021574318409e-38
Best regards
Giulio Malinverno
--
==================================================
Necesse est quod sub certo numero omnia creata comprehendantur
==================================================
More information about the gmp-bugs
mailing list