system/freebsd-src
system/freebsd-src
1
0
Fork 0
mirror of https://github.com/freebsd/freebsd-src synced 2024-10-05 08:00:30 +00:00
Code Issues Packages Projects Releases Wiki Activity

Simplified setiing up args for __kernel_rem_pio2(). We already have x

Browse source 
with a 24-bit fraction, so we don't need a loop to split it into up to
3 terms with 24-bit fractions.
This commit is contained in:
Bruce Evans 2005-11-23 03:03:09 +00:00
parent 33f8f56e09
commit 01231dd04c
Notes: svn2git 2020-12-20 02:59:44 +00:00 
svn path=/head/; revision=152707
1 changed files with 9 additions and 17 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

26
lib/msun/src/e_rem_pio2f.c
View file

@ -61,9 +61,10 @@ pio2_1t = 6.07710050650619224932e-11; /* 0x3DD0B461, 0x1A626331 */
int32_t __ieee754_rem_pio2f(float x, float *y)
{
double z,w,t,r,fn;
double tx[3],ty[2];
int32_t e0,i,nx,n,ix,hx;
double w,t,r,fn;
double tx[1],ty[2];
float z;
int32_t e0,n,ix,hx;
GET_FLOAT_WORD(hx,x);
ix = hx&0x7fffffff;
@ -85,20 +86,11 @@ pio2_1t = 6.07710050650619224932e-11; /* 0x3DD0B461, 0x1A626331 */
if(ix>=0x7f800000) { /* x is inf or NaN */
y[0]=y[1]=x-x; return 0;
}
/* set z = scalbn(|x|,ilogb(x)-23) */
z = x;
GET_HIGH_WORD(hx,z);
ix = hx&0x7fffffff;
e0 = (ix>>20)-1046; /* e0 = ilogb(z)-23; */
SET_HIGH_WORD(z, ix - ((int32_t)(e0<<20)));
for(i=0;i<2;i++) {
tx[i] = (double)((int32_t)(z));
z = (z-tx[i])*two24;
}
tx[2] = z;
nx = 3;
while(tx[nx-1]==zero) nx--; /* skip zero term */
n = __kernel_rem_pio2(tx,ty,e0,nx,1,two_over_pi);
/* set z = scalbn(|x|,ilogb(|x|)-23) */
e0 = (ix>>23)-150; /* e0 = ilogb(|x|)-23; */
SET_FLOAT_WORD(z, ix - ((int32_t)(e0<<23)));
tx[0] = z;
n = __kernel_rem_pio2(tx,ty,e0,1,1,two_over_pi);
y[0] = ty[0];
y[1] = ty[0] - y[0];
if(hx<0) {y[0] = -y[0]; y[1] = -y[1]; return -n;}