[<prev] [next>] [<thread-prev] [thread-next>] [day] [month] [year] [list]
Message-ID: <20260730153746.GJ27423@brightrain.aerifal.cx>
Date: Thu, 30 Jul 2026 11:37:47 -0400
From: Rich Felker <dalias@...c.org>
To: Paul Zimmermann <Paul.Zimmermann@...ia.fr>
Cc: musl@...ts.openwall.com
Subject: Re: issue in acosh
On Thu, Jul 30, 2026 at 10:16:49AM +0200, Paul Zimmermann wrote:
> Hi,
>
> for x=-0x1.34e729fd08086p+21, Musl 1.2.6 yields 0x1.075fa39542c7fp+3
> instead of NaN:
>
> $ VERBOSE=-v ./doit.musl acosh 1000
> Checking acosh with musl-1.2.6
> Using seed 16199
> tested 1121438 numbers from extra file
> Using 20 thread(s)
> acosh 0 -1 -0x1.34e729fd08086p+21 [0] [inf] inf inf
> libm gives 0x1.075fa39542c7fp+3
> mpfr gives -nan
Thanks! I looked at the function briefly. Does the comment:
/* x < 1 domain error is handled in the called functions */
somehow fail to hold because of a rounding error in the intermediate
expression?
Rich
Powered by blists - more mailing lists
Confused about mailing lists and their use?
Read about mailing lists on Wikipedia
and check out these
guidelines on proper formatting of your messages.
