public inbox for glibc-bugs@sourceware.org
help / color / mirror / Atom feed
* [Bug math/10357] New: casinh(x) loses precision near x=0
@ 2009-07-01 21:07 pav at iki dot fi
  0 siblings, 0 replies; only message in thread
From: pav at iki dot fi @ 2009-07-01 21:07 UTC (permalink / raw)
  To: glibc-bugs

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.


^ permalink raw reply	[flat|nested] only message in thread
only message in thread, other threads:[~2009-07-01 21:07 UTC | newest]

Thread overview: (only message) (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2009-07-01 21:07 [Bug math/10357] New: casinh(x) loses precision near x=0 pav at iki dot fi

This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).