d/freebsd-dev
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Add a u suffix to the IEEEl2bits unions o_threshold and u_threshold,

Browse Source 
and use macros to access the e component of the unions.  This allows
the portions of the code in ld80 to be identical to the ld128 code.

Obtained from:	bde
This commit is contained in:
Steve Kargl 2013-06-03 18:07:04 +00:00
parent 4aa8c9453f
commit ad36b00fcb
Notes: svn2git 2020-12-20 02:59:44 +00:00 
svn path=/head/; revision=251328
1 changed files with 7 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
lib/msun/ld80/s_expl.c
View File

@ -61,9 +61,12 @@ static volatile const long double tiny = 0x1p-10000L;
static const union IEEEl2bits
/* log(2**16384 - 0.5) rounded towards zero: */
o_threshold = LD80C(0xb17217f7d1cf79ab, 13, 11356.5234062941439488L),
/* log(2**16384 - 0.5 + 1) rounded towards zero for expm1l() is the same: */
o_thresholdu = LD80C(0xb17217f7d1cf79ab, 13, 11356.5234062941439488L),
#define o_threshold (o_thresholdu.e)
/* log(2**(-16381-64-1)) rounded towards zero: */
u_threshold = LD80C(0xb21dfe7f09e2baa9, 13, -11399.4985314888605581L);
u_thresholdu = LD80C(0xb21dfe7f09e2baa9, 13, -11399.4985314888605581L);
#define u_threshold (u_thresholdu.e)
static const double
/*
@ -246,9 +249,9 @@ expl(long double x)
return (0.0L); /* x is -Inf */
return (x + x); /* x is +Inf, NaN or unsupported */
}
if (x > o_threshold.e)
if (x > o_threshold)
return (huge * huge);
if (x < u_threshold.e)
if (x < u_threshold)
return (tiny * tiny);
} else if (ix < BIAS - 66) { /* |x| < 0x1p-66 */
/* includes pseudo-denormals */