huge error for pow
Tue Feb 7 21:48:14 GMT 2023
Do you think you can put together a patch?
-- Jeff J.
On Tue, Feb 7, 2023 at 9:27 AM Paul Zimmermann <Paul.Zimmermann@inria.fr>
> while testing Newlib 4.3.0, I noticed the following huge error:
> zimmerma@salade:~/svn/tbd/20/src/binary64$ VERBOSE=-v ./doitb.newlib pow
> Checking pow with newlib-184.108.40.20630120
> Using seed 1619945
> pow 0 -1 -0x1.647ff80007ff8p-576,-0x1.3d018267f12fp+48 [inf] inf inf
> libm gives inf
> mpfr gives -inf
> The exponent y in x^y is a large odd integer, thus since x is negative we
> should get a negative number as result.
> This issue was apparently already present in 3.3.0, but we only discovered
> it now.
> Best regards,
More information about the Newlib