libs/newlib-cygwin
libs/newlib-cygwin
4
0
Fork 0
mirror of git://sourceware.org/git/newlib-cygwin.git synced 2025-01-20 05:19:21 +08:00
Code Issues Projects Releases Wiki Activity

Mon Apr 17 12:46:00 2000 Marek Michalkiewicz <marekm@linux.org.pl>

Browse Source 
* libc/signal/signal.c (_signal_r) : Removed unused local variable temp.        * libc/stdio/findfp.c (std): Added declaration of flags and file.
        * libc/stdio/mktemp.c (_gettemp, _mkstemp_r, mkstemp): Added int
        return type.
        * libc/stdio/putchar.c (putchar): Added return statement.
        * libc/stdio/refill.c (lflush): Added correct parentheses.
        * libc/stdio/vfprintf.c (_VFPRINTF_R): Ditto.
        * libc/stdio/vfscanf.c (__svfscanf): Changed sprintf call which
        prints long value to use l qualifier.
        * libc/stdlib/dtoa.c (_dtoa_r): Added parentheses to remove warning
        messages and initialized local values: ilim, ilim1, and spec_case.
        * libc/stdlib/ecvtbuf.c (print_e): Removed unused variable dp.
        * libc/stdlib/mbctype.h (_issjis1, _issjis2): Added parentheses.
        * libc/stdlib/mprec.c: Ditto.
        * libc/stdlib/setenv_r.c: Ditto.
        * libc/stdlib/strtod.c: Ditto.
        * libc/stdlib/strtol.c: Ditto.
        * libc/stdlib/strtoul.c: Ditto.
        * libm/common/sf_expm1.c: Added curly braces to if else clauses.
        * libm/common/sf_log1p.c: Ditto.
        * libm/common/sf_scalbn.c: Ditto.
        * libm/math/ef_log.c: Ditto.
This commit is contained in:
Jeff Johnston 2000-04-17 17:10:18 +00:00
parent 08ce80eac2
commit c87be3e4d6
20 changed files with 74 additions and 45 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
newlib/ChangeLog
View File

@ -1,3 +1,28 @@
Mon Apr 17 12:46:00 2000 Marek Michalkiewicz <marekm@linux.org.pl>
* libc/signal/signal.c (_signal_r) : Removed unused local variable temp.
* libc/stdio/findfp.c (std): Added declaration of flags and file.
* libc/stdio/mktemp.c (_gettemp, _mkstemp_r, mkstemp): Added int
return type.
* libc/stdio/putchar.c (putchar): Added return statement.
* libc/stdio/refill.c (lflush): Added correct parentheses.
* libc/stdio/vfprintf.c (_VFPRINTF_R): Ditto.
* libc/stdio/vfscanf.c (__svfscanf): Changed sprintf call which
prints long value to use l qualifier.
* libc/stdlib/dtoa.c (_dtoa_r): Added parentheses to remove warning
messages and initialized local values: ilim, ilim1, and spec_case.
* libc/stdlib/ecvtbuf.c (print_e): Removed unused variable dp.
* libc/stdlib/mbctype.h (_issjis1, _issjis2): Added parentheses.
* libc/stdlib/mprec.c: Ditto.
* libc/stdlib/setenv_r.c: Ditto.
* libc/stdlib/strtod.c: Ditto.
* libc/stdlib/strtol.c: Ditto.
* libc/stdlib/strtoul.c: Ditto.
* libm/common/sf_expm1.c: Added curly braces to if else clauses.
* libm/common/sf_log1p.c: Ditto.
* libm/common/sf_scalbn.c: Ditto.
* libm/math/ef_log.c: Ditto.
Sun Apr 16 12:45:00 2000 Corinna Vinschen <corinna@vinschen.de> Sun Apr 16 12:45:00 2000 Corinna Vinschen <corinna@vinschen.de>
* libc/posix/execvp.c (execvp): Check path for * libc/posix/execvp.c (execvp): Check path for

2
newlib/libc/signal/signal.c
View File

@ -144,7 +144,7 @@ _DEFUN (_signal_r, (ptr, sig, func),
int sig _AND int sig _AND
_sig_func_ptr func) _sig_func_ptr func)
{ {
_sig_func_ptr old_func, *temp; _sig_func_ptr old_func;
if (sig < 0 || sig >= NSIG) if (sig < 0 || sig >= NSIG)
{ {

2
newlib/libc/stdio/findfp.c
View File

@ -26,6 +26,8 @@
static void static void
std (ptr, flags, file, data) std (ptr, flags, file, data)
FILE *ptr; FILE *ptr;
int flags;
int file;
struct _reent *data; struct _reent *data;
{ {
ptr->_p = 0; ptr->_p = 0;

4
newlib/libc/stdio/mktemp.c
View File

@ -102,7 +102,7 @@ Supporting OS subroutines required: <<getpid>>, <<open>>, <<stat>>.
#include <ctype.h> #include <ctype.h>
#include <reent.h> #include <reent.h>
static static int
_DEFUN (_gettemp, (ptr, path, doopen), _DEFUN (_gettemp, (ptr, path, doopen),
struct _reent *ptr _AND struct _reent *ptr _AND
char *path _AND char *path _AND
@ -182,6 +182,7 @@ _DEFUN (_gettemp, (ptr, path, doopen),
/*NOTREACHED*/ /*NOTREACHED*/
} }
int
_DEFUN (_mkstemp_r, (ptr, path), _DEFUN (_mkstemp_r, (ptr, path),
struct _reent *ptr _AND struct _reent *ptr _AND
char *path) char *path)
@ -201,6 +202,7 @@ _DEFUN (_mktemp_r, (ptr, path),
#ifndef _REENT_ONLY #ifndef _REENT_ONLY
int
_DEFUN (mkstemp, (path), _DEFUN (mkstemp, (path),
char *path) char *path)
{ {

2
newlib/libc/stdio/putchar.c
View File

@ -88,7 +88,7 @@ putchar (c)
{ {
/* CHECK_INIT is (eventually) called by __swbuf. */ /* CHECK_INIT is (eventually) called by __swbuf. */
_putchar_r (_REENT, c); return _putchar_r (_REENT, c);
} }
#endif #endif

2
newlib/libc/stdio/refill.c
View File

@ -25,7 +25,7 @@ static int
lflush (fp) lflush (fp)
FILE *fp; FILE *fp;
{ {
if ((fp->_flags & (__SLBF | __SWR)) == __SLBF | __SWR) if ((fp->_flags & (__SLBF | __SWR)) == (__SLBF | __SWR))
return fflush (fp); return fflush (fp);
return 0; return 0;
} }

2
newlib/libc/stdio/vfprintf.c
View File

@ -824,7 +824,7 @@ number: if ((dprec = prec) >= 0)
ox[0] = *cp++; ox[0] = *cp++;
ox[1] = '.'; ox[1] = '.';
PRINT(ox, 2); PRINT(ox, 2);
if (_double || flags & ALT == 0) { if (_double || (flags & ALT) == 0) {
PRINT(cp, ndig-1); PRINT(cp, ndig-1);
} else /* 0.[0..] */ } else /* 0.[0..] */
/* __dtoa irregularity */ /* __dtoa irregularity */

2
newlib/libc/stdio/vfscanf.c
View File

@ -825,7 +825,7 @@ __svfscanf (fp, fmt0, ap)
truncate some trailing digits to make room. */ truncate some trailing digits to make room. */
if (exp_start >= buf + sizeof (buf) - MAX_LONG_LEN) if (exp_start >= buf + sizeof (buf) - MAX_LONG_LEN)
exp_start = buf + sizeof (buf) - MAX_LONG_LEN - 1; exp_start = buf + sizeof (buf) - MAX_LONG_LEN - 1;
sprintf (exp_start, "e%d", new_exp); sprintf (exp_start, "e%ld", new_exp);
} }
res = atof (buf); res = atof (buf);
if (flags & LONG) if (flags & LONG)

29
newlib/libc/stdlib/dtoa.c
View File

@ -291,7 +291,7 @@ _DEFUN (_dtoa_r,
#ifdef Sudden_Underflow #ifdef Sudden_Underflow
i = (int) (word0 (d) >> Exp_shift1 & (Exp_mask >> Exp_shift1)); i = (int) (word0 (d) >> Exp_shift1 & (Exp_mask >> Exp_shift1));
#else #else
if (i = (int) (word0 (d) >> Exp_shift1 & (Exp_mask >> Exp_shift1))) if ((i = (int) (word0 (d) >> Exp_shift1 & (Exp_mask >> Exp_shift1))) != 0)
{ {
#endif #endif
d2.d = d.d; d2.d = d.d;
@ -337,8 +337,8 @@ _DEFUN (_dtoa_r,
/* d is denormalized */ /* d is denormalized */
i = bbits + be + (Bias + (P - 1) - 1); i = bbits + be + (Bias + (P - 1) - 1);
x = i > 32 ? word0 (d) << 64 - i | word1 (d) >> i - 32 x = (i > 32) ? (word0 (d) << (64 - i)) | (word1 (d) >> (i - 32))
: word1 (d) << 32 - i; : (word1 (d) << (32 - i));
d2.d = x; d2.d = x;
word0 (d2) -= 31 * Exp_msk1; /* adjust exponent */ word0 (d2) -= 31 * Exp_msk1; /* adjust exponent */
i -= (Bias + (P - 1) - 1) + 1; i -= (Bias + (P - 1) - 1) + 1;
@ -388,11 +388,11 @@ _DEFUN (_dtoa_r,
try_quick = 0; try_quick = 0;
} }
leftright = 1; leftright = 1;
ilim = ilim1 = -1;
switch (mode) switch (mode)
{ {
case 0: case 0:
case 1: case 1:
ilim = ilim1 = -1;
i = 18; i = 18;
ndigits = 0; ndigits = 0;
break; break;
@ -449,7 +449,7 @@ _DEFUN (_dtoa_r,
} }
d.d /= ds; d.d /= ds;
} }
else if (j1 = -k) else if ((j1 = -k) != 0)
{ {
d.d *= tens[j1 & 0xf]; d.d *= tens[j1 & 0xf];
for (j = j1 >> 4; j; j >>= 1, i++) for (j = j1 >> 4; j; j >>= 1, i++)
@ -564,7 +564,7 @@ _DEFUN (_dtoa_r,
if (i == ilim) if (i == ilim)
{ {
d.d += d.d; d.d += d.d;
if (d.d > ds || d.d == ds && L & 1) if ((d.d > ds) || ((d.d == ds) && (L & 1)))
{ {
bump_up: bump_up:
while (*--s == '9') while (*--s == '9')
@ -640,7 +640,7 @@ _DEFUN (_dtoa_r,
Bfree (ptr, b); Bfree (ptr, b);
b = b1; b = b1;
} }
if (j = b5 - m5) if ((j = b5 - m5) != 0)
b = pow5mult (ptr, b, j); b = pow5mult (ptr, b, j);
} }
else else
@ -652,6 +652,7 @@ _DEFUN (_dtoa_r,
/* Check for special case that d is a normalized power of 2. */ /* Check for special case that d is a normalized power of 2. */
spec_case = 0;
if (mode < 2) if (mode < 2)
{ {
if (!word1 (d) && !(word0 (d) & Bndry_mask) if (!word1 (d) && !(word0 (d) & Bndry_mask)
@ -665,8 +666,6 @@ _DEFUN (_dtoa_r,
s2 += Log2P; s2 += Log2P;
spec_case = 1; spec_case = 1;
} }
else
spec_case = 0;
} }
/* Arrange for convenient computation of quotients: /* Arrange for convenient computation of quotients:
@ -678,10 +677,10 @@ _DEFUN (_dtoa_r,
*/ */
#ifdef Pack_32 #ifdef Pack_32
if (i = ((s5 ? 32 - hi0bits (S->_x[S->_wds - 1]) : 1) + s2) & 0x1f) if ((i = ((s5 ? 32 - hi0bits (S->_x[S->_wds - 1]) : 1) + s2) & 0x1f) != 0)
i = 32 - i; i = 32 - i;
#else #else
if (i = ((s5 ? 32 - hi0bits (S->_x[S->_wds - 1]) : 1) + s2) & 0xf) if ((i = ((s5 ? 32 - hi0bits (S->_x[S->_wds - 1]) : 1) + s2) & 0xf) != 0)
i = 16 - i; i = 16 - i;
#endif #endif
if (i > 4) if (i > 4)
@ -765,17 +764,17 @@ _DEFUN (_dtoa_r,
goto ret; goto ret;
} }
#endif #endif
if (j < 0 || j == 0 && !mode if ((j < 0) || ((j == 0) && !mode
#ifndef ROUND_BIASED #ifndef ROUND_BIASED
&& !(word1 (d) & 1) && !(word1 (d) & 1)
#endif #endif
) ))
{ {
if (j1 > 0) if (j1 > 0)
{ {
b = lshift (ptr, b, 1); b = lshift (ptr, b, 1);
j1 = cmp (b, S); j1 = cmp (b, S);
if ((j1 > 0 || j1 == 0 && dig & 1) if (((j1 > 0) || ((j1 == 0) && (dig & 1)))
&& dig++ == '9') && dig++ == '9')
goto round_9_up; goto round_9_up;
} }
@ -819,7 +818,7 @@ _DEFUN (_dtoa_r,
b = lshift (ptr, b, 1); b = lshift (ptr, b, 1);
j = cmp (b, S); j = cmp (b, S);
if (j > 0 || j == 0 && dig & 1) if ((j > 0) || ((j == 0) && (dig & 1)))
{ {
roundoff: roundoff:
while (*--s == '9') while (*--s == '9')

1
newlib/libc/stdlib/ecvtbuf.c
View File

@ -151,7 +151,6 @@ _DEFUN (print_e, (ptr, buf, invalue, width, type, dot),
char type _AND char type _AND
int dot) int dot)
{ {
int dp;
int sign; int sign;
char *end; char *end;
char *p; char *p;

4
newlib/libc/stdlib/mbctype.h
View File

@ -12,8 +12,8 @@ int _EXFUN(_issjis2, (int c));
int _EXFUN(_iseucjp, (int c)); int _EXFUN(_iseucjp, (int c));
int _EXFUN(_isjis, (int c)); int _EXFUN(_isjis, (int c));
#define _issjis1(c) ((c) >= 0x81 && (c) <= 0x9f || (c) >= 0xe0 && (c) <= 0xef) #define _issjis1(c) (((c) >= 0x81 && (c) <= 0x9f) || ((c) >= 0xe0 && (c) <= 0xef))
#define _issjis2(c) ((c) >= 0x40 && (c) <= 0x7e || (c) >= 0x80 && (c) <= 0xfc) #define _issjis2(c) (((c) >= 0x40 && (c) <= 0x7e) || ((c) >= 0x80 && (c) <= 0xfc))
#define _iseucjp(c) ((c) >= 0xa1 && (c) <= 0xfe) #define _iseucjp(c) ((c) >= 0xa1 && (c) <= 0xfe)
#define _isjis(c) ((c) >= 0x21 && (c) <= 0x7e) #define _isjis(c) ((c) >= 0x21 && (c) <= 0x7e)

18
newlib/libc/stdlib/mprec.c
View File

@ -107,7 +107,7 @@ _DEFUN (Balloc, (ptr, k), struct _reent *ptr _AND int k)
} }
} }
if (rv = ptr->_freelist[k]) if ((rv = ptr->_freelist[k]) != 0)
{ {
ptr->_freelist[k] = rv->_next; ptr->_freelist[k] = rv->_next;
} }
@ -354,7 +354,7 @@ _DEFUN (mult, (ptr, a, b), struct _reent * ptr _AND _Bigint * a _AND _Bigint * b
#ifdef Pack_32 #ifdef Pack_32
for (; xb < xbe; xb++, xc0++) for (; xb < xbe; xb++, xc0++)
{ {
if (y = *xb & 0xffff) if ((y = *xb & 0xffff) != 0)
{ {
x = xa; x = xa;
xc = xc0; xc = xc0;
@ -370,7 +370,7 @@ _DEFUN (mult, (ptr, a, b), struct _reent * ptr _AND _Bigint * a _AND _Bigint * b
while (x < xae); while (x < xae);
*xc = carry; *xc = carry;
} }
if (y = *xb >> 16) if ((y = *xb >> 16) != 0)
{ {
x = xa; x = xa;
xc = xc0; xc = xc0;
@ -420,7 +420,7 @@ _DEFUN (pow5mult,
int i; int i;
static _CONST int p05[3] = {5, 25, 125}; static _CONST int p05[3] = {5, 25, 125};
if (i = k & 3) if ((i = k & 3) != 0)
b = multadd (ptr, b, p05[i - 1], 0); b = multadd (ptr, b, p05[i - 1], 0);
if (!(k >>= 2)) if (!(k >>= 2))
@ -484,7 +484,7 @@ _DEFUN (lshift, (ptr, b, k), struct _reent * ptr _AND _Bigint * b _AND int k)
z = *x++ >> k1; z = *x++ >> k1;
} }
while (x < xe); while (x < xe);
if (*x1 = z) if ((*x1 = z) != 0)
++n1; ++n1;
} }
#else #else
@ -697,17 +697,17 @@ _DEFUN (b2d, (a, e),
#ifdef Pack_32 #ifdef Pack_32
if (k < Ebits) if (k < Ebits)
{ {
d0 = Exp_1 | y >> Ebits - k; d0 = Exp_1 | y >> (Ebits - k);
w = xa > xa0 ? *--xa : 0; w = xa > xa0 ? *--xa : 0;
#ifndef _DOUBLE_IS_32BITS #ifndef _DOUBLE_IS_32BITS
d1 = y << (32 - Ebits) + k | w >> Ebits - k; d1 = y << ((32 - Ebits) + k) | w >> (Ebits - k);
#endif #endif
goto ret_d; goto ret_d;
} }
z = xa > xa0 ? *--xa : 0; z = xa > xa0 ? *--xa : 0;
if (k -= Ebits) if (k -= Ebits)
{ {
d0 = Exp_1 | y << k | z >> 32 - k; d0 = Exp_1 | y << k | z >> (32 - k);
y = xa > xa0 ? *--xa : 0; y = xa > xa0 ? *--xa : 0;
#ifndef _DOUBLE_IS_32BITS #ifndef _DOUBLE_IS_32BITS
d1 = z << k | y >> 32 - k; d1 = z << k | y >> 32 - k;
@ -789,7 +789,7 @@ _DEFUN (d2b,
z |= Exp_msk11; z |= Exp_msk11;
#endif #endif
#else #else
if (de = (int) (d0 >> Exp_shift)) if ((de = (int) (d0 >> Exp_shift)) != 0)
z |= Exp_msk1; z |= Exp_msk1;
#endif #endif
#ifdef Pack_32 #ifdef Pack_32

4
newlib/libc/stdlib/setenv_r.c
View File

@ -64,7 +64,7 @@ _DEFUN (_setenv_r, (reent_ptr, name, value, rewrite),
} }
if (strlen (C) >= l_value) if (strlen (C) >= l_value)
{ /* old larger; copy over */ { /* old larger; copy over */
while (*C++ = *value++); while ((*C++ = *value++) != 0);
ENV_UNLOCK; ENV_UNLOCK;
return 0; return 0;
} }
@ -108,7 +108,7 @@ _DEFUN (_setenv_r, (reent_ptr, name, value, rewrite),
return -1; return -1;
} }
for (C = environ[offset]; (*C = *name++) && *C != '='; ++C); for (C = environ[offset]; (*C = *name++) && *C != '='; ++C);
for (*C++ = '='; *C++ = *value++;); for (*C++ = '='; (*C++ = *value++) != 0;);
ENV_UNLOCK; ENV_UNLOCK;

4
newlib/libc/stdlib/strtod.c
View File

@ -316,7 +316,7 @@ dig_done:
if (e1 > 0) if (e1 > 0)
{ {
if (i = e1 & 15) if ((i = e1 & 15) != 0)
rv.d *= tens[i]; rv.d *= tens[i];
if (e1 &= ~15) if (e1 &= ~15)
{ {
@ -373,7 +373,7 @@ dig_done:
else if (e1 < 0) else if (e1 < 0)
{ {
e1 = -e1; e1 = -e1;
if (i = e1 & 15) if ((i = e1 & 15) != 0)
rv.d /= tens[i]; rv.d /= tens[i];
if (e1 &= ~15) if (e1 &= ~15)
{ {

2
newlib/libc/stdlib/strtol.c
View File

@ -194,7 +194,7 @@ _DEFUN (_strtol_r, (rptr, nptr, endptr, base),
break; break;
if (c >= base) if (c >= base)
break; break;
if (any < 0 || acc > cutoff || acc == cutoff && c > cutlim) if (any < 0 || acc > cutoff || (acc == cutoff && c > cutlim))
any = -1; any = -1;
else { else {
any = 1; any = 1;

2
newlib/libc/stdlib/strtoul.c
View File

@ -174,7 +174,7 @@ _DEFUN (_strtoul_r, (rptr, nptr, endptr, base),
break; break;
if (c >= base) if (c >= base)
break; break;
if (any < 0 || acc > cutoff || acc == cutoff && c > cutlim) if (any < 0 || acc > cutoff || (acc == cutoff && c > cutlim))
any = -1; any = -1;
else { else {
any = 1; any = 1;

3
newlib/libm/common/sf_expm1.c
View File

@ -102,9 +102,10 @@ Q5 = -2.0109921195e-07; /* 0xb457edbb */
e = (x*(e-c)-c); e = (x*(e-c)-c);
e -= hxs; e -= hxs;
if(k== -1) return (float)0.5*(x-e)-(float)0.5; if(k== -1) return (float)0.5*(x-e)-(float)0.5;
if(k==1) if(k==1) {
if(x < (float)-0.25) return -(float)2.0*(e-(x+(float)0.5)); if(x < (float)-0.25) return -(float)2.0*(e-(x+(float)0.5));
else return one+(float)2.0*(x-e); else return one+(float)2.0*(x-e);
}
if (k <= -2 || k>56) { /* suffice to return exp(x)-1 */ if (k <= -2 || k>56) { /* suffice to return exp(x)-1 */
__int32_t i; __int32_t i;
y = one-(e-x); y = one-(e-x);

4
newlib/libm/common/sf_log1p.c
View File

@ -93,8 +93,8 @@ static float zero = 0.0;
} }
hfsq=(float)0.5*f*f; hfsq=(float)0.5*f*f;
if(hu==0) { /* |f| < 2**-20 */ if(hu==0) { /* |f| < 2**-20 */
if(f==zero) if(k==0) return zero; if(f==zero) { if(k==0) return zero;
else {c += k*ln2_lo; return k*ln2_hi+c;} else {c += k*ln2_lo; return k*ln2_hi+c;}}
R = hfsq*((float)1.0-(float)0.66666666666666666*f); R = hfsq*((float)1.0-(float)0.66666666666666666*f);
if(k==0) return f-R; else if(k==0) return f-R; else
return k*ln2_hi-((R-(k*ln2_lo+c))-f); return k*ln2_hi-((R-(k*ln2_lo+c))-f);

3
newlib/libm/common/sf_scalbn.c
View File

@ -54,10 +54,11 @@ tiny = 1.0e-30;
if (k > 0xfe) return huge*copysignf(huge,x); /* overflow */ if (k > 0xfe) return huge*copysignf(huge,x); /* overflow */
if (k > 0) /* normal result */ if (k > 0) /* normal result */
{SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); return x;} {SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); return x;}
if (k <= -25) if (k <= -25) {
if (n > OVERFLOW_INT) /* in case integer overflow in n+k */ if (n > OVERFLOW_INT) /* in case integer overflow in n+k */
return huge*copysignf(huge,x); /*overflow*/ return huge*copysignf(huge,x); /*overflow*/
else return tiny*copysignf(tiny,x); /*underflow*/ else return tiny*copysignf(tiny,x); /*underflow*/
}
k += 25; /* subnormal result */ k += 25; /* subnormal result */
SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23));
return x*twom25; return x*twom25;

4
newlib/libm/math/ef_log.c
View File

@ -65,8 +65,8 @@ static float zero = 0.0;
k += (i>>23); k += (i>>23);
f = x-(float)1.0; f = x-(float)1.0;
if((0x007fffff&(15+ix))<16) { /* |f| < 2**-20 */ if((0x007fffff&(15+ix))<16) { /* |f| < 2**-20 */
if(f==zero) if(k==0) return zero; else {dk=(float)k; if(f==zero) { if(k==0) return zero; else {dk=(float)k;
return dk*ln2_hi+dk*ln2_lo;} return dk*ln2_hi+dk*ln2_lo;}}
R = f*f*((float)0.5-(float)0.33333333333333333*f); R = f*f*((float)0.5-(float)0.33333333333333333*f);
if(k==0) return f-R; else {dk=(float)k; if(k==0) return f-R; else {dk=(float)k;
return dk*ln2_hi-((R-dk*ln2_lo)-f);} return dk*ln2_hi-((R-dk*ln2_lo)-f);}