Boost C++ Libraries: Ticket #8809: pow(cpp_dec_float_type(-1), cpp_dec_float_type(n)) != (+-)1 for n >= 0x8000000000000000u

attachment set

Wed, 10 Jul 2013 12:08:03 GMT

attachment → boost_multiprecision_power_demo2.cpp

Demonstration of the error

Wed, 10 Jul 2013 12:19:00 GMT

Curiously (or not ;-)) - my std library function pow(double, double) gives the same answers ... ... perhaps shifting the scope of this elsewhere?

status changed; cc set

John Maddock — Wed, 10 Jul 2013 17:44:00 GMT

cc pbristow@… e_float@… added
status new → assigned

This is actually quite tricky - and hard to deal with in a consistent (across different backends) way.

The issue is that these routines need to check for integer exponent when the base is < 1, but this is hard to do when the exponent becomes overly large. A related point is that once the exponent becomes sufficiently large it's no longer possible to determine whether it is even or odd (needed for the sign of the result) as the final digit has dropped off the end of the floating point type.

That said we should be able to fix this up to numeric_limits<uintmax_t>::max(), and this has uncovered some broken error handling too. I'm CC'ing Chris and Paul into this issue in case they have some ideas.

status changed; resolution set

John Maddock — Fri, 12 Jul 2013 16:13:36 GMT

status assigned → closed
resolution → fixed

(In [85008]) Fix mpfi pow function to correctly handle negative base with integer power. Fixed error handling to not return numeric_limits<>::quiet_NaN if there is no NaN supported. Fixed default pow implementation to handle integer args up to numeric_limits<uintmax_t>::max(). Fixed error handling in default pow implementation to return NaN for too large exponents (can't tell if they're integers or even or odd) and negative bases. Fixed array subscript bug in cpp_dec_float. Greatly increased pow testing. Fixes #8809.