[Bug math/552] New: Wrong return values for tgamma (0) and tgamma (-0)

[Bug math/552] New: Wrong return values for tgamma (0) and tgamma (-0)

[hjl@lucon.org]

[-- Warning: decoded text below may be mangled, UTF-8 assumed --]
[-- Attachment #1: Type: text/plain, Size: 846 bytes --]

According to IEEE Std 1003.1, 2004 Edition:

http://www.opengroup.org/onlinepubs/009695399/functions/tgamma.html

If x is ±0, a pole error shall occur, and tgamma(), tgammaf(), and tgammal()
shall return ±HUGE_VAL, ±HUGE_VALF, and ±HUGE_VALL, respectively.

But glibc returns NAN.

-- 
           Summary: Wrong return values for tgamma (0) and tgamma (-0)
           Product: glibc
           Version: unspecified
            Status: NEW
          Severity: normal
          Priority: P2
         Component: math
        AssignedTo: aj at suse dot de
        ReportedBy: hjl at lucon dot org
                CC: glibc-bugs at sources dot redhat dot com


http://sources.redhat.com/bugzilla/show_bug.cgi?id=552

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

[Bug math/552] Wrong return values for tgamma (0) and tgamma (-0)

[cvs-commit@gcc.gnu.org]

------- Additional Comments From cvs-commit at gcc dot gnu dot org  2004-11-21 11:42 -------
Subject: Bug 552

CVSROOT:	/cvs/glibc
Module name:	libc
Branch: 	fedora-branch
Changes by:	jakub@sources.redhat.com	2004-11-21 11:42:09

Modified files:
	math           : libm-test.inc 
	sysdeps/generic: w_tgamma.c w_tgammaf.c w_tgammal.c 
	sysdeps/ieee754: k_standard.c 
	sysdeps/ieee754/dbl-64: e_gamma_r.c 
	sysdeps/ieee754/flt-32: e_gammaf_r.c 
	sysdeps/ieee754/ldbl-128: e_gammal_r.c 
	sysdeps/ieee754/ldbl-96: e_gammal_r.c 

Log message:
	2004-11-20  Jakub Jelinek  <jakub@redhat.com>
	
	* sysdeps/ieee754/k_standard.c: Document code 50.
	(__kernel_standard) <case 50>: Avoid raising div-by-zero
	exception again.
	
	2004-11-19  H.J. Lu  <hongjiu.lu@intel.com>
	
	[BZ #552]
	* math/libm-test.inc (tgamma_test): Update tgamma (0) and
	tgamma (-0).
	
	* sysdeps/generic/w_tgamma.c (__tgamma): Properly handle
	|x| == 0.
	* sysdeps/generic/w_tgammaf.c (__tgammaf): Likewise.
	* sysdeps/generic/w_tgammal.c (__tgammal): Likewise.
	* sysdeps/ieee754/dbl-64/e_gamma_r.c (__ieee754_gamma_r): Likewise.
	* sysdeps/ieee754/flt-32/e_gammaf_r.c: Likewise.
	* sysdeps/ieee754/ldbl-128/e_gammal_r.c: Likewise.
	* sysdeps/ieee754/ldbl-96/e_gammal_r.c: Likewise.
	
	* sysdeps/ieee754/k_standard.c (__kernel_standard): Handle
	tgamma (0) and tgamma (-0).

Patches:
http://sources.redhat.com/cgi-bin/cvsweb.cgi/libc/math/libm-test.inc.diff?cvsroot=glibc&only_with_tag=fedora-branch&r1=1.61.2.1&r2=1.61.2.2
http://sources.redhat.com/cgi-bin/cvsweb.cgi/libc/sysdeps/generic/w_tgamma.c.diff?cvsroot=glibc&only_with_tag=fedora-branch&r1=1.1&r2=1.1.4.1
http://sources.redhat.com/cgi-bin/cvsweb.cgi/libc/sysdeps/generic/w_tgammaf.c.diff?cvsroot=glibc&only_with_tag=fedora-branch&r1=1.1&r2=1.1.4.1
http://sources.redhat.com/cgi-bin/cvsweb.cgi/libc/sysdeps/generic/w_tgammal.c.diff?cvsroot=glibc&only_with_tag=fedora-branch&r1=1.1&r2=1.1.4.1
http://sources.redhat.com/cgi-bin/cvsweb.cgi/libc/sysdeps/ieee754/k_standard.c.diff?cvsroot=glibc&only_with_tag=fedora-branch&r1=1.2&r2=1.2.4.1
http://sources.redhat.com/cgi-bin/cvsweb.cgi/libc/sysdeps/ieee754/dbl-64/e_gamma_r.c.diff?cvsroot=glibc&only_with_tag=fedora-branch&r1=1.4&r2=1.4.4.1
http://sources.redhat.com/cgi-bin/cvsweb.cgi/libc/sysdeps/ieee754/flt-32/e_gammaf_r.c.diff?cvsroot=glibc&only_with_tag=fedora-branch&r1=1.4&r2=1.4.4.1
http://sources.redhat.com/cgi-bin/cvsweb.cgi/libc/sysdeps/ieee754/ldbl-128/e_gammal_r.c.diff?cvsroot=glibc&only_with_tag=fedora-branch&r1=1.3&r2=1.3.2.1
http://sources.redhat.com/cgi-bin/cvsweb.cgi/libc/sysdeps/ieee754/ldbl-96/e_gammal_r.c.diff?cvsroot=glibc&only_with_tag=fedora-branch&r1=1.5&r2=1.5.2.1



-- 


http://sources.redhat.com/bugzilla/show_bug.cgi?id=552

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

[Bug math/552] Wrong return values for tgamma (0) and tgamma (-0)

[jakub@redhat.com]

------- Additional Comments From jakub at redhat dot com  2004-11-22 09:22 -------
The same change went into HEAD.

-- 
           What    |Removed                     |Added
----------------------------------------------------------------------------
             Status|NEW                         |RESOLVED
         Resolution|                            |FIXED
            Summary|Wrong return values for     |Wrong return values for
                   |tgamma (0) and tgamma (-0)  |tgamma (0) and tgamma (-0)


http://sources.redhat.com/bugzilla/show_bug.cgi?id=552

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