From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <glibc-bugs-return-9918-listarch-glibc-bugs=sources.redhat.com@sourceware.org>
Received: (qmail 26011 invoked by alias); 1 Jul 2009 21:07:31 -0000
Received: (qmail 25946 invoked by uid 48); 1 Jul 2009 21:07:18 -0000
Date: Wed, 01 Jul 2009 21:07:00 -0000
From: "pav at iki dot fi" <sourceware-bugzilla@sourceware.org>
To: glibc-bugs@sources.redhat.com
Message-ID: <20090701210718.10357.pav@iki.fi>
Reply-To: sourceware-bugzilla@sourceware.org
Subject: [Bug math/10357] New: casinh(x) loses precision near x=0
X-Bugzilla-Reason: CC
Mailing-List: contact glibc-bugs-help@sourceware.org; run by ezmlm
Precedence: bulk
List-Id: <glibc-bugs.sourceware.org>
List-Subscribe: <mailto:glibc-bugs-subscribe@sourceware.org>
List-Post: <mailto:glibc-bugs@sourceware.org>
List-Help: <mailto:glibc-bugs-help@sourceware.org>, <http://sourceware.org/lists.html#faqs>
Sender: glibc-bugs-owner@sourceware.org
X-SW-Source: 2009-07/txt/msg00003.txt.bz2

The `casinh(x)` implementation in glibc-2_10-branch loses precision for small x.

The following assertion should pass:

    assert(casinh(1e-20) == 9.9999999999999995e-21);

However, with Glibc 2.9 it does not: `casinh(1e-20) == 0` with it. In Glibc
sources (math/s_asinh.c), casinh is computed from y=log(x + sqrt(1 + x*x)). In
finite precision, this yields y=0 for |x| < floating-point epsilon.

-- 
           Summary: casinh(x) loses precision near x=0
           Product: glibc
           Version: 2.10
            Status: NEW
          Severity: minor
          Priority: P2
         Component: math
        AssignedTo: aj at suse dot de
        ReportedBy: pav at iki dot fi
                CC: glibc-bugs at sources dot redhat dot com


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

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