From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 3817 invoked by alias); 27 Nov 2013 21:42:11 -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 3773 invoked by uid 48); 27 Nov 2013 21:42:07 -0000 From: "jsm28 at gcc dot gnu.org" To: glibc-bugs@sourceware.org Subject: [Bug math/16271] New: dbl-64 e_sqrt.c not correctly rounding in non-default rounding modes Date: Wed, 27 Nov 2013 21:42: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.19 X-Bugzilla-Keywords: X-Bugzilla-Severity: normal X-Bugzilla-Who: jsm28 at gcc dot gnu.org X-Bugzilla-Status: NEW 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 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: 2013-11/txt/msg00299.txt.bz2 https://sourceware.org/bugzilla/show_bug.cgi?id=16271 Bug ID: 16271 Summary: dbl-64 e_sqrt.c not correctly rounding in non-default rounding modes Product: glibc Version: 2.19 Status: NEW Severity: normal Priority: P2 Component: math Assignee: unassigned at sourceware dot org Reporter: jsm28 at gcc dot gnu.org sysdeps/ieee754/dbl-64/e_sqrt.c does not produce correct results for all inputs in rounding modes other than round-to-nearest, as illustrated by (on ARM): Failure: Test: sqrt_towardzero (0xf.fffffp+124) Result: is: 1.84467435239537295360e+19 0x1.fffffeffffffc0000000p+63 should be: 1.84467435239537274880e+19 0x1.fffffeffffffb0000000p+63 difference: 2.04800000000000000000e+03 0x1.00000000000000000000p+11 ulp : 1.0000 max.ulp : 0.0000 Failure: Test: sqrt_towardzero (0xf.fffffp-4) Result: is: 9.99999970197677168215e-01 0x1.fffffeffffffc0000000p-1 should be: 9.99999970197677057193e-01 0x1.fffffeffffffb0000000p-1 difference: 1.11022302462515654042e-16 0x1.00000000000000000000p-53 ulp : 1.0000 max.ulp : 0.0000 Maximal error of `sqrt_towardzero' is : 1 ulp accepted: 0 ulp Failure: Test: sqrt_downward (0xf.fffffp+124) Result: is: 1.84467435239537295360e+19 0x1.fffffeffffffc0000000p+63 should be: 1.84467435239537274880e+19 0x1.fffffeffffffb0000000p+63 difference: 2.04800000000000000000e+03 0x1.00000000000000000000p+11 ulp : 1.0000 max.ulp : 0.0000 Failure: Test: sqrt_downward (0xf.fffffp-4) Result: is: 9.99999970197677168215e-01 0x1.fffffeffffffc0000000p-1 should be: 9.99999970197677057193e-01 0x1.fffffeffffffb0000000p-1 difference: 1.11022302462515654042e-16 0x1.00000000000000000000p-53 ulp : 1.0000 max.ulp : 0.0000 Maximal error of `sqrt_downward' is : 1 ulp accepted: 0 ulp Failure: Test: sqrt_upward (0x1.000002p+0) Result: is: 1.00000005960464299904e+00 0x1.000000ffffff80000000p+0 should be: 1.00000005960464322108e+00 0x1.000000ffffff90000000p+0 difference: 2.22044604925031308085e-16 0x1.00000000000000000000p-52 ulp : 1.0000 max.ulp : 0.0000 Maximal error of `sqrt_upward' is : 1 ulp accepted: 0 ulp -- You are receiving this mail because: You are on the CC list for the bug.