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)  * 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) #include "ieee754dp.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 12) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 13) static inline union ieee754sp ieee754sp_nan_fdp(int xs, u64 xm)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 14) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 15) 	return buildsp(xs, SP_EMAX + 1 + SP_EBIAS,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 16) 		       xm >> (DP_FBITS - SP_FBITS));
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 17) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 18) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 19) union ieee754sp ieee754sp_fdp(union ieee754dp x)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 20) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 21) 	union ieee754sp y;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 22) 	u32 rm;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 23) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 24) 	COMPXDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 25) 	COMPYSP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 26) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 27) 	EXPLODEXDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 28) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 29) 	ieee754_clearcx();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 30) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 31) 	FLUSHXDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 32) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 33) 	switch (xc) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 34) 	case IEEE754_CLASS_SNAN:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 35) 		x = ieee754dp_nanxcpt(x);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 36) 		EXPLODEXDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 37) 		fallthrough;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38) 	case IEEE754_CLASS_QNAN:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39) 		y = ieee754sp_nan_fdp(xs, xm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40) 		if (!ieee754_csr.nan2008) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) 			EXPLODEYSP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42) 			if (!ieee754_class_nan(yc))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 43) 				y = ieee754sp_indef();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 44) 		}
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 45) 		return y;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) 	case IEEE754_CLASS_INF:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 48) 		return ieee754sp_inf(xs);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 49) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 50) 	case IEEE754_CLASS_ZERO:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 51) 		return ieee754sp_zero(xs);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 52) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 53) 	case IEEE754_CLASS_DNORM:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 54) 		/* can't possibly be sp representable */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 55) 		ieee754_setcx(IEEE754_UNDERFLOW);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 56) 		ieee754_setcx(IEEE754_INEXACT);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 57) 		if ((ieee754_csr.rm == FPU_CSR_RU && !xs) ||
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 58) 				(ieee754_csr.rm == FPU_CSR_RD && xs))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 59) 			return ieee754sp_mind(xs);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 60) 		return ieee754sp_zero(xs);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 61) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 62) 	case IEEE754_CLASS_NORM:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 63) 		break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 64) 	}
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 65) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 66) 	/*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 67) 	 * Convert from DP_FBITS to SP_FBITS+3 with sticky right shift.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 68) 	 */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 69) 	rm = (xm >> (DP_FBITS - (SP_FBITS + 3))) |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 70) 	     ((xm << (64 - (DP_FBITS - (SP_FBITS + 3)))) != 0);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 71) 
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 72) 	return ieee754sp_format(xs, xe, rm);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 73) }