Re: [PACTH,libm] hypothl(x) mishandles subnormal numbers.

2021-02-09 Thread Steve Kargl
On Wed, Feb 10, 2021 at 12:26:29AM +0100, Dimitry Andric wrote: > > > On 10 Feb 2021, at 00:15, Steve Kargl > > wrote: > > > > This patch fixes the issue. t1 is used to scale the subnormal > > numbers. It is generated by scaling the exponent, but that > > only works if t1 is 1 not 0. > > > >

Re: [PACTH,libm] hypothl(x) mishandles subnormal numbers.

2021-02-09 Thread Dimitry Andric
> On 10 Feb 2021, at 00:15, Steve Kargl > wrote: > > On Sat, Feb 06, 2021 at 10:32:33PM +0100, Dimitry Andric wrote: >> On 6 Feb 2021, at 22:04, Steve Kargl >> wrote: >>> >>> On Sat, Feb 06, 2021 at 12:39:29PM -0800, Steve Kargl wrote: I've long forgotten by freebsd bugzilla password.

[PACTH,libm] hypothl(x) mishandles subnormal numbers.

2021-02-09 Thread Steve Kargl
On Sat, Feb 06, 2021 at 10:32:33PM +0100, Dimitry Andric wrote: > On 6 Feb 2021, at 22:04, Steve Kargl > wrote: > > > > On Sat, Feb 06, 2021 at 12:39:29PM -0800, Steve Kargl wrote: > >> I've long forgotten by freebsd bugzilla password. > >> So, if someone would like to submit a bug report, > >>

Re: hypothl(x) mishandles subnormal numbers.

2021-02-06 Thread Steve Kargl
On Sat, Feb 06, 2021 at 10:32:33PM +0100, Dimitry Andric wrote: > On 6 Feb 2021, at 22:04, Steve Kargl > wrote: > > > > On Sat, Feb 06, 2021 at 12:39:29PM -0800, Steve Kargl wrote: > >> I've long forgotten by freebsd bugzilla password. > >> So, if someone would like to submit a bug report, > >>

Re: hypothl(x) mishandles subnormal numbers.

2021-02-06 Thread Dimitry Andric
On 6 Feb 2021, at 22:04, Steve Kargl wrote: > > On Sat, Feb 06, 2021 at 12:39:29PM -0800, Steve Kargl wrote: >> I've long forgotten by freebsd bugzilla password. >> So, if someone would like to submit a bug report, >> here's a test program. >> > Forgot to include that issue was identified from >

Re: hypothl(x) mishandles subnormal numbers.

2021-02-06 Thread Steve Kargl
On Sat, Feb 06, 2021 at 12:39:29PM -0800, Steve Kargl wrote: > I've long forgotten by freebsd bugzilla password. > So, if someone would like to submit a bug report, > here's a test program. > Forgot to include that issue was identified from a bug report in the OpenLibm bug mailing list. https:/

hypothl(x) mishandles subnormal numbers.

2021-02-06 Thread Steve Kargl
I've long forgotten by freebsd bugzilla password. So, if someone would like to submit a bug report, here's a test program. % cc -o z -O2 -fno-builtin z.c -lm % ./z Via scaling and sqrtl x=2.853684e-4932 y=1.444012e-4922 a=1.444012e-4922 x=0x1.b2933cafa0bb8p-16383 y=0x1.fp-16351 a=0x