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