huge error for pow

Paul Zimmermann <[email protected]>
Newsgroups gmane.comp.lib.newlib
Message-ID <[email protected]>
       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
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.