[Bug math/18699] New: tilegx cproj() for various complex infinities does not yield infinity

[Bug math/18699] New: tilegx cproj() for various complex infinities does not yield infinity

[cmetcalf at ezchip dot com]

https://sourceware.org/bugzilla/show_bug.cgi?id=18699

            Bug ID: 18699
           Summary: tilegx cproj() for various complex infinities does not
                    yield infinity
           Product: glibc
           Version: unspecified
            Status: NEW
          Severity: normal
          Priority: P2
         Component: math
          Assignee: unassigned at sourceware dot org
          Reporter: cmetcalf at ezchip dot com
  Target Milestone: ---

The float and ifloat tests both fail with the appended output: cproj()
of various complex infinities should be yielding infinity, but instead
yields a very tiny positive number.

This is for glibc 2.22 (almost, as it hasn't fully frozen yet).

testing float (without inline functions)
Failure: Test: Real part of: cproj (inf + inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (inf - inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-inf + inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-inf - inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (inf + +0 i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (inf - 0 i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-inf + +0 i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-inf - 0 i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (+0 + inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (+0 - inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-0 + inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-0 - inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (inf + 1.0 i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (inf - 1.0 i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-inf + 1.0 i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-inf - 1.0 i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (1.0 + inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (1.0 - inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-1.0 + inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-1.0 - inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (inf + qNaN i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (inf - qNaN i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-inf + qNaN i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-inf - qNaN i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (qNaN + inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (qNaN - inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-qNaN + inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf
Failure: Test: Real part of: cproj (-qNaN - inf i)
Result:
 is:          3.57331108402828353086e-43   0x1.fe000000000000000000p-142
 should be:   inf   inf

Test suite completed:
  9998 test cases plus 9169 tests for exception flags and
    9169 tests for errno executed.
  28 errors occurred.

-- 
You are receiving this mail because:
You are on the CC list for the bug.

[Bug math/18699] tilegx cproj() for various complex infinities does not yield infinity

[cmetcalf at ezchip dot com]

https://sourceware.org/bugzilla/show_bug.cgi?id=18699

--- Comment #1 from Chris Metcalf <cmetcalf at ezchip dot com> ---
Joseph Myers writes, "This sounds like a compiler bug to me (whether in the
compilation of cprojf or something else involved in the test)."

-- 
You are receiving this mail because:
You are on the CC list for the bug.

[Bug math/18699] tilegx cproj() for various complex infinities does not yield infinity

[cmetcalf at ezchip dot com]

https://sourceware.org/bugzilla/show_bug.cgi?id=18699

Chris Metcalf <cmetcalf at ezchip dot com> changed:

           What    |Removed                     |Added
----------------------------------------------------------------------------
             Status|NEW                         |RESOLVED
         Resolution|---                         |FIXED

--- Comment #2 from Chris Metcalf <cmetcalf at ezchip dot com> ---
Fixed by:

commit fe8c2b33aed067282ecac0cc3fdff1feef88548b
Author: Wilco Dijkstra <wdijkstr@arm.com>
Date:   Fri Sep 18 20:51:52 2015 +0100

    Since we now inline isinf, isnan and isfinite in math.h, replace uses of
__isinf_ns(l/f)
    with isinf, and remove the unused inlines __isinf_ns(l/f), __isnan(f) and
__finite(f).

-- 
You are receiving this mail because:
You are on the CC list for the bug.

[Bug math/18699] tilegx cproj() for various complex infinities does not yield infinity

[cvs-commit at gcc dot gnu.org]

https://sourceware.org/bugzilla/show_bug.cgi?id=18699

--- Comment #3 from cvs-commit at gcc dot gnu.org <cvs-commit at gcc dot gnu.org> ---
This is an automated email from the git hooks/post-receive script. It was
generated because a ref change was pushed to the repository containing
the project "GNU C Library master sources".

The branch, master has been updated
       via  45ccb8346bba9b0f7a72a92d7f8b1c796dd496ba (commit)
      from  822095e8a4f757f02b25ac4b8bcfd8c4c8109eb7 (commit)

Those revisions listed above that are new to this repository have
not appeared on any other notification email; so we list those
revisions in full, below.

- Log -----------------------------------------------------------------
https://sourceware.org/git/gitweb.cgi?p=glibc.git;h=45ccb8346bba9b0f7a72a92d7f8b1c796dd496ba

commit 45ccb8346bba9b0f7a72a92d7f8b1c796dd496ba
Author: Chris Metcalf <cmetcalf@ezchip.com>
Date:   Thu Oct 22 13:36:53 2015 -0400

    Update NEWS to mention drive-by fix for bug 18699.

    Wilco fixed this tilegx bug in commit fe8c2b33aed0 ("Since we
    now inline isinf, isnan and isfinite in math.h, replace uses
    of __isinf_ns(l/f) with isinf, and remove the unused inlines
    __isinf_ns(l/f), __isnan(f) and __finite(f).")

    I verified that reverting math/s_cprojf.c back to using the
    sysdeps/ieee754/flt-32/math_private.h version of __isinf_nsf()
    instead of isinf() brings back the bug on tilegx.

-----------------------------------------------------------------------

Summary of changes:
 ChangeLog |    2 ++
 NEWS      |   15 ++++++++-------
 2 files changed, 10 insertions(+), 7 deletions(-)

-- 
You are receiving this mail because:
You are on the CC list for the bug.

[Bug math/18699] tilegx cproj() for various complex infinities does not yield infinity

[jsm28 at gcc dot gnu.org]

https://sourceware.org/bugzilla/show_bug.cgi?id=18699

Joseph Myers <jsm28 at gcc dot gnu.org> changed:

           What    |Removed                     |Added
----------------------------------------------------------------------------
   Target Milestone|---                         |2.23

-- 
You are receiving this mail because:
You are on the CC list for the bug.