From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 14277 invoked by alias); 24 Jun 2015 20:31:06 -0000 Mailing-List: contact glibc-bugs-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Post: List-Help: , Sender: glibc-bugs-owner@sourceware.org Received: (qmail 13945 invoked by uid 48); 24 Jun 2015 20:31:02 -0000 From: "jsm28 at gcc dot gnu.org" To: glibc-bugs@sourceware.org Subject: [Bug math/18594] New: cexp, ccos, ccosh, csin, csinh spurious underflows Date: Wed, 24 Jun 2015 20:31:00 -0000 X-Bugzilla-Reason: CC X-Bugzilla-Type: new X-Bugzilla-Watch-Reason: None X-Bugzilla-Product: glibc X-Bugzilla-Component: math X-Bugzilla-Version: 2.21 X-Bugzilla-Keywords: X-Bugzilla-Severity: normal X-Bugzilla-Who: jsm28 at gcc dot gnu.org X-Bugzilla-Status: NEW X-Bugzilla-Resolution: X-Bugzilla-Priority: P2 X-Bugzilla-Assigned-To: unassigned at sourceware dot org X-Bugzilla-Target-Milestone: --- X-Bugzilla-Flags: X-Bugzilla-Changed-Fields: bug_id short_desc product version bug_status bug_severity priority component assigned_to reporter target_milestone Message-ID: Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: 7bit X-Bugzilla-URL: http://sourceware.org/bugzilla/ Auto-Submitted: auto-generated MIME-Version: 1.0 X-SW-Source: 2015-06/txt/msg00280.txt.bz2 https://sourceware.org/bugzilla/show_bug.cgi?id=18594 Bug ID: 18594 Summary: cexp, ccos, ccosh, csin, csinh spurious underflows Product: glibc Version: 2.21 Status: NEW Severity: normal Priority: P2 Component: math Assignee: unassigned at sourceware dot org Reporter: jsm28 at gcc dot gnu.org Target Milestone: --- cexp, ccos, ccosh, csin and csinh have spurious underflows in cases where they compute sin of the smallest normal, that produces an underflow exception (depending on which sin implementation is in use) but the final result does not underflow. ctan may also have such underflows, or they may be latent (the issue there is that e.g. ctan (DBL_MIN) should, rounded upwards, be the next double value above DBL_MIN, which under glibc's accuracy goals may not have an underflow exception, but the intermediate computation of sin (DBL_MIN) would legitimately underflow on before-rounding architectures). For example, for float on x86_64: Failure: Real part of: ccos_downward (0x4p-128 + 0x1p+0 i): Exception "Underflow" set Failure: Real part of: ccosh_downward (0x1p+0 + 0x4p-128 i): Exception "Underflow" set Failure: Real part of: cexp_downward (0x1p+0 + 0x4p-128 i): Exception "Underflow" set Failure: Real part of: csin_downward (-0x4p-128 + 0x1p+0 i): Exception "Underflow" set Failure: Real part of: csinh_downward (0x1p+0 + 0x4p-128 i): Exception "Underflow" set -- You are receiving this mail because: You are on the CC list for the bug.