Orange Pi5 kernel

^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)  * double precision: common utilities
^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 "ieee754dp.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 11) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 12) union ieee754dp ieee754dp_flong(s64 x)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 13) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 14) 	u64 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 ieee754dp_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 ieee754dp_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 ieee754dp_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 == (1ULL << 63))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 30) 			xm = (1ULL << 63);	/* 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) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 37) 	/* normalize */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38) 	xe = DP_FBITS + 3;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39) 	if (xm >> (DP_FBITS + 1 + 3)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40) 		/* shunt out overflow bits */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) 		while (xm >> (DP_FBITS + 1 + 3)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42) 			XDPSRSX1();
^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 double precision */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) 		while ((xm >> (DP_FBITS + 3)) == 0) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) 			xm <<= 1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 48) 			xe--;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 49) 		}
^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 ieee754dp_format(xs, xe, xm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 53) }