Again a mechanical conversion.
Signed-off-by: Alex Bennée <alex.ben...@linaro.org>
---
fpu/softfloat.c | 82 +++++++++++++++++++++++++++++++++++++++++++++++++
include/fpu/softfloat.h | 1 +
2 files changed, 83 insertions(+)
diff --git a/fpu/softfloat.c b/fpu/softfloat.c
index f7473f97e3..dc7f5f6d88 100644
--- a/fpu/softfloat.c
+++ b/fpu/softfloat.c
@@ -3532,6 +3532,88 @@ static void normalizeFloat16Subnormal(uint32_t aSig, int
*zExpPtr,
*zExpPtr = 1 - shiftCount;
}
+/*----------------------------------------------------------------------------
+| Rounds the half-precision floating-point value `a' to an integer,
+| and returns the result as a half-precision floating-point value. The
+| operation is performed according to the IEC/IEEE Standard for Binary
+| Floating-Point Arithmetic.
+*----------------------------------------------------------------------------*/
+
+float16 float16_round_to_int(float16 a, float_status *status)
+{
+ flag aSign;
+ int aExp;
+ uint16_t lastBitMask, roundBitsMask;
+ uint16_t z;
+ a = float16_squash_input_denormal(a, status);
+
+ aExp = extractFloat16Exp( a );
+ if ( 0x19 <= aExp ) {
+ if ( ( aExp == 0x1F ) && extractFloat16Frac( a ) ) {
+ return propagateFloat16NaN(a, a, status);
+ }
+ return a;
+ }
+ if ( aExp <= 0xE ) {
+ if ( (uint16_t) ( float16_val(a)<<1 ) == 0 ) return a;
+ status->float_exception_flags |= float_flag_inexact;
+ aSign = extractFloat16Sign( a );
+ switch (status->float_rounding_mode) {
+ case float_round_nearest_even:
+ if ( ( aExp == 0xE ) && extractFloat16Frac( a ) ) {
+ return packFloat16( aSign, 0xF, 0 );
+ }
+ break;
+ case float_round_ties_away:
+ if (aExp == 0xE) {
+ return packFloat16(aSign, 0xF, 0);
+ }
+ break;
+ case float_round_down:
+ return make_float16(aSign ? 0xBC00 : 0);
+ case float_round_up:
+ /* -0.0/1.0f */
+ return make_float16(aSign ? 0x8000 : 0x3C00);
+ }
+ return packFloat16( aSign, 0, 0 );
+ }
+ lastBitMask = 1;
+ lastBitMask <<= 0x19 - aExp;
+ roundBitsMask = lastBitMask - 1;
+ z = float16_val(a);
+ switch (status->float_rounding_mode) {
+ case float_round_nearest_even:
+ z += lastBitMask>>1;
+ if ((z & roundBitsMask) == 0) {
+ z &= ~lastBitMask;
+ }
+ break;
+ case float_round_ties_away:
+ z += lastBitMask >> 1;
+ break;
+ case float_round_to_zero:
+ break;
+ case float_round_up:
+ if (!extractFloat16Sign(make_float16(z))) {
+ z += roundBitsMask;
+ }
+ break;
+ case float_round_down:
+ if (extractFloat16Sign(make_float16(z))) {
+ z += roundBitsMask;
+ }
+ break;
+ default:
+ abort();
+ }
+ z &= ~ roundBitsMask;
+ if (z != float16_val(a)) {
+ status->float_exception_flags |= float_flag_inexact;
+ }
+ return make_float16(z);
+
+}
+
/*----------------------------------------------------------------------------
| Returns the result of adding the absolute values of the half-precision
| floating-point values `a' and `b'. If `zSign' is 1, the sum is negated
diff --git a/include/fpu/softfloat.h b/include/fpu/softfloat.h
index a7435e2a5b..856f67cf12 100644
--- a/include/fpu/softfloat.h
+++ b/include/fpu/softfloat.h
@@ -346,6 +346,7 @@ float64 float16_to_float64(float16 a, flag ieee,
float_status *status);
| Software half-precision operations.
*----------------------------------------------------------------------------*/
+float16 float16_round_to_int(float16, float_status *status);
float16 float16_add(float16, float16, float_status *status);
float16 float16_sub(float16, float16, float_status *status);
float16 float16_mul(float16, float16, float_status *status);
--
2.14.1
