huge error for pow
Paul Zimmermann
Paul.Zimmermann@inria.fr
Tue Feb 7 14:27:42 GMT 2023
Hi,
while testing Newlib 4.3.0, I noticed the following huge error:
zimmerma@salade:~/svn/tbd/20/src/binary64$ VERBOSE=-v ./doitb.newlib pow 1000
Checking pow with newlib-4.3.0.20230120
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,
Paul
More information about the Newlib
mailing list