JehanneOS
/
newlib
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

newlib/libm/common: Don't re-convert float to bits in modf/modff 

These functions shared a pattern of re-converting the argument to bits
when returning +/-0. Skip that as the initial conversion still has the
sign bit.

Signed-off-by: Keith Packard <keithp@keithp.com>
This commit is contained in:
Keith Packard via Newlib 2020-03-25 17:18:20 -07:00 committed by Corinna Vinschen
parent 61cd34c1bf
commit 3439f3b0e9
2 changed files with 5 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
newlib/libm/common/s_modf.c
View File

@ -81,10 +81,8 @@ QUICKREF
} else { } else {
i = (0x000fffff)>>j0; i = (0x000fffff)>>j0;
if(((i0&i)|i1)==0) { /* x is integral */ if(((i0&i)|i1)==0) { /* x is integral */
__uint32_t high;
*iptr = x; *iptr = x;
GET_HIGH_WORD(high,x); INSERT_WORDS(x,i0&0x80000000,0); /* return +-0 */
INSERT_WORDS(x,high&0x80000000,0); /* return +-0 */
return x; return x;
} else { } else {
INSERT_WORDS(*iptr,i0&(~i),0); INSERT_WORDS(*iptr,i0&(~i),0);
@ -92,19 +90,15 @@ QUICKREF
} }
} }
} else if (j0>51) { /* no fraction part */ } else if (j0>51) { /* no fraction part */
__uint32_t high;
*iptr = x; *iptr = x;
if (__fpclassifyd(x) == FP_NAN) return *iptr = x+x; /* x is NaN, return NaN */ if (__fpclassifyd(x) == FP_NAN) return *iptr = x+x; /* x is NaN, return NaN */
GET_HIGH_WORD(high,x); INSERT_WORDS(x,i0&0x80000000,0); /* return +-0 */
INSERT_WORDS(x,high&0x80000000,0); /* return +-0 */
return x; return x;
} else { /* fraction part in low x */ } else { /* fraction part in low x */
i = ((__uint32_t)(0xffffffff))>>(j0-20); i = ((__uint32_t)(0xffffffff))>>(j0-20);
if((i1&i)==0) { /* x is integral */ if((i1&i)==0) { /* x is integral */
__uint32_t high;
*iptr = x; *iptr = x;
GET_HIGH_WORD(high,x); INSERT_WORDS(x,i0&0x80000000,0); /* return +-0 */
INSERT_WORDS(x,high&0x80000000,0); /* return +-0 */
return x; return x;
} else { } else {
INSERT_WORDS(*iptr,i0,i1&(~i)); INSERT_WORDS(*iptr,i0,i1&(~i));

8
newlib/libm/common/sf_modf.c
View File

@ -33,10 +33,8 @@
} else { } else {
i = (0x007fffff)>>j0; i = (0x007fffff)>>j0;
if((i0&i)==0) { /* x is integral */ if((i0&i)==0) { /* x is integral */
__uint32_t ix;
*iptr = x; *iptr = x;
GET_FLOAT_WORD(ix,x); SET_FLOAT_WORD(x,i0&0x80000000); /* return +-0 */
SET_FLOAT_WORD(x,ix&0x80000000); /* return +-0 */
return x; return x;
} else { } else {
SET_FLOAT_WORD(*iptr,i0&(~i)); SET_FLOAT_WORD(*iptr,i0&(~i));
@ -44,11 +42,9 @@
} }
} }
} else { /* no fraction part */ } else { /* no fraction part */
__uint32_t ix;
*iptr = x; *iptr = x;
if (__fpclassifyf(x) == FP_NAN) return *iptr = x+x; /* x is NaN, return NaN */ if (__fpclassifyf(x) == FP_NAN) return *iptr = x+x; /* x is NaN, return NaN */
GET_FLOAT_WORD(ix,x); SET_FLOAT_WORD(x,i0&0x80000000); /* return +-0 */
SET_FLOAT_WORD(x,ix&0x80000000); /* return +-0 */
return x; return x;
} }
} }