fabs(-0.0) returns -0.0
John Baldwin
jhb at freebsd.org
Tue Nov 14 15:58:06 UTC 2017
On Monday, November 13, 2017 10:52:10 PM Khilan Gudka wrote:
> Hi,
>
> The implementation of fabs(3) for MIPS (at least) appears to not be giving
> the right answer for when -0 is passed, as it returns -0 instead of +0. The
> implementation of fabs is:
>
> double fabs(double x) {
> if (x < 0)
> return -x;
> return x;
> }
>
> The if-test fails for -0 and thus it is returned. Is this a known issue? A
> simple fix would be to return x+0.0 instead of just x.
>
> The implication of this is that other functions which rely on fabs, such as
> hypot, also return -0. For example, hypot(-0.0, 0.0) returns -0 instead of
> +0.
The C version of fabs() for 32-bit arm has the same misimplementation btw.
All other architectures implement fabs in assembly in libc.
--
John Baldwin
More information about the freebsd-mips
mailing list