^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 1) // SPDX-License-Identifier: GPL-2.0-only
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 2) /* IEEE754 floating point arithmetic
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 3) * single precision
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 4) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 5) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 6) * MIPS floating point support
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 7) * Copyright (C) 1994-2000 Algorithmics Ltd.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 8) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 9)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 10) #include "ieee754sp.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 11)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 12) union ieee754sp ieee754sp_fint(int x)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 13) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 14) unsigned int xm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 15) int xe;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 16) int xs;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 17)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 18) ieee754_clearcx();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 19)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 20) if (x == 0)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 21) return ieee754sp_zero(0);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 22) if (x == 1 || x == -1)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 23) return ieee754sp_one(x < 0);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 24) if (x == 10 || x == -10)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 25) return ieee754sp_ten(x < 0);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 26)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 27) xs = (x < 0);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 28) if (xs) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 29) if (x == (1 << 31))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 30) xm = ((unsigned) 1 << 31); /* max neg can't be safely negated */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 31) else
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 32) xm = -x;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 33) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 34) xm = x;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 35) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 36) xe = SP_FBITS + 3;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 37)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38) if (xm >> (SP_FBITS + 1 + 3)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39) /* shunt out overflow bits
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) while (xm >> (SP_FBITS + 1 + 3)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42) SPXSRSX1();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 43) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 44) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 45) /* normalize in grs extended single precision
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) while ((xm >> (SP_FBITS + 3)) == 0) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 48) xm <<= 1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 49) xe--;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 50) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 51) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 52) return ieee754sp_format(xs, xe, xm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 53) }