Re: [PATCH] Update sparc ULPs.

[David Miller <davem at davemloft dot net>]

From: "Joseph S. Myers" <joseph@codesourcery.com>
Date: Sat, 17 Nov 2012 22:23:47 +0000

> On Sat, 17 Nov 2012, David Miller wrote:
> 
>> diff --git a/sysdeps/ieee754/ldbl-128/e_j0l.c b/sysdeps/ieee754/ldbl-128/e_j0l.c
>> index 112a8f3..7ee4731 100644
>> --- a/sysdeps/ieee754/ldbl-128/e_j0l.c
>> +++ b/sysdeps/ieee754/ldbl-128/e_j0l.c
>> @@ -831,6 +831,8 @@ long double
>>        return -HUGE_VALL + x;
>>      }
>>    xx = fabsl (x);
>> +  if (xx < 0x2p-27)
>> +    return TWOOPI * __ieee754_logl (x);
> 
> The logic in the double version is:
> 
>  *         Note: For tiny x, U/V = u0 and j0(x)~1, hence
>  *              y0(tiny) = u0 + (2/pi)*ln(tiny), (choose tiny<2**-27)

That's all fine and well, but I should note that the float versions of
these routines use the same exact value of tiny as the double ones.