[Bug math/6779] New: remainder() wrongly gives a domain error for (NaN, 0)

[Bug math/6779] New: remainder() wrongly gives a domain error for (NaN, 0)

[mtk dot manpages at gmail dot com]

The posix reaminder() spec says:

==
If x or y is NaN, a NaN shall be returned.

If x is infinite or y is 0 and the other is non-NaN, a domain error shall occur,
and either a NaN (if supported), or an implementation-defined value shall be
returned
==

On glibc 2.8, remainder(NaN, 0) gives a domain error (EDOM, FE_INVALID).  As per
the above spec, it should return a NaN without error.

-- 
           Summary: remainder() wrongly gives a domain error for (NaN, 0)
           Product: glibc
           Version: unspecified
            Status: NEW
          Severity: normal
          Priority: P2
         Component: math
        AssignedTo: aj at suse dot de
        ReportedBy: mtk dot manpages at gmail dot com
                CC: glibc-bugs at sources dot redhat dot com


http://sourceware.org/bugzilla/show_bug.cgi?id=6779

------- You are receiving this mail because: -------
You are on the CC list for the bug, or are watching someone who is.

[Bug math/6779] remainder() wrongly gives a domain error for (NaN, 0)

[mtk dot manpages at gmail dot com]

------- Additional Comments From mtk dot manpages at gmail dot com  2008-07-31 08:34 -------
Created an attachment (id=2844)
 --> (http://sourceware.org/bugzilla/attachment.cgi?id=2844&action=view)
test program

Sample run showing problem:

$ /tmp/mt_remainder nan 0
errno == EDOM
fetestexcept() says:  FE_INVALID
remainder(nan,0.00000000000000000e+00)=nan
EDOM FE_INVALID nan


-- 


http://sourceware.org/bugzilla/show_bug.cgi?id=6779

------- You are receiving this mail because: -------
You are on the CC list for the bug, or are watching someone who is.