d0p1/tcc-stupidos
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Round mode of ll -> float conversion to nearest

Browse source 
Change rounding mode of long long to float conversion to nearest in
libtcc1.
This commit is contained in:
Thomas Preud'homme 2014-02-05 19:47:15 +08:00
parent d029507494
commit 88c9f1bb4e
1 changed files with 20 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

25
lib/armeabi.c
View file

@ -170,7 +170,7 @@ DEFINE__AEABI_D2XLZ(d2lz, 1)
#define DEFINE__AEABI_XL2F(name, with_sign) \ #define DEFINE__AEABI_XL2F(name, with_sign) \
unsigned __aeabi_ ## name(unsigned long long v) \ unsigned __aeabi_ ## name(unsigned long long v) \
{ \ { \
int s /* shift */, sign = 0; \ int s /* shift */, flb /* first lost bit */, sign = 0; \
unsigned p = 0 /* power */, ret; \ unsigned p = 0 /* power */, ret; \
double_unsigned_struct val; \ double_unsigned_struct val; \
\ \
@ -188,20 +188,28 @@ unsigned __aeabi_ ## name(unsigned long long v) \
if (s < FLOAT_FRAC_BITS) { \ if (s < FLOAT_FRAC_BITS) { \
ret <<= FLOAT_FRAC_BITS - s; \ ret <<= FLOAT_FRAC_BITS - s; \
ret |= val.low >> (32 - (FLOAT_FRAC_BITS - s)); \ ret |= val.low >> (32 - (FLOAT_FRAC_BITS - s)); \
} else \ flb = (val.low >> (32 - (FLOAT_FRAC_BITS - s - 1))) & 1; \
} else { \
flb = (ret >> (s - FLOAT_FRAC_BITS - 1)) & 1; \
ret >>= s - FLOAT_FRAC_BITS; \ ret >>= s - FLOAT_FRAC_BITS; \
} \
s += 32; \ s += 32; \
} else { \ } else { \
for (s = 31, p = 1 << 31; p && !(val.low & p); s--, p >>= 1); \ for (s = 31, p = 1 << 31; p && !(val.low & p); s--, p >>= 1); \
if (p) { \ if (p) { \
ret = val.low & (p - 1); \ ret = val.low & (p - 1); \
if (s <= FLOAT_FRAC_BITS) \ if (s <= FLOAT_FRAC_BITS) { \
ret <<= FLOAT_FRAC_BITS - s; \ ret <<= FLOAT_FRAC_BITS - s; \
else \ flb = 0; \
} else { \
flb = (ret >> (s - FLOAT_FRAC_BITS - 1)) & 1; \
ret >>= s - FLOAT_FRAC_BITS; \ ret >>= s - FLOAT_FRAC_BITS; \
} \
} else \ } else \
return 0; \ return 0; \
} \ } \
if (flb) \
ret++; \
\ \
/* fill exponent bits */ \ /* fill exponent bits */ \
ret |= (s + ONE_EXP(FLOAT)) << FLOAT_FRAC_BITS; \ ret |= (s + ONE_EXP(FLOAT)) << FLOAT_FRAC_BITS; \
@ -222,7 +230,7 @@ DEFINE__AEABI_XL2F(l2f, 1)
#define __AEABI_XL2D(name, with_sign) \ #define __AEABI_XL2D(name, with_sign) \
void __aeabi_ ## name(unsigned long long v) \ void __aeabi_ ## name(unsigned long long v) \
{ \ { \
int s, high_shift, sign = 0; \ int s /* shift */, high_shift, sign = 0; \
unsigned tmp, p = 0; \ unsigned tmp, p = 0; \
double_unsigned_struct val, ret; \ double_unsigned_struct val, ret; \
\ \
@ -248,6 +256,13 @@ void __aeabi_ ## name(unsigned long long v) \
ret.high = tmp >> high_shift; \ ret.high = tmp >> high_shift; \
ret.low = tmp << (32 - high_shift); \ ret.low = tmp << (32 - high_shift); \
ret.low |= val.low >> high_shift; \ ret.low |= val.low >> high_shift; \
if ((val.low >> (high_shift - 1)) & 1) { \
if (ret.low == UINT_MAX) { \
ret.high++; \
ret.low = 0; \
} else \
ret.low++; \
} \
} \ } \
s += 32; \ s += 32; \
} else { \ } else { \