1 // SPDX-License-Identifier: GPL-2.0-only
2 /* IEEE754 floating point arithmetic
3  * double precision: common utilities
4  */
5 /*
6  * MIPS floating point support
7  * Copyright (C) 1994-2000 Algorithmics Ltd.
8  */
9 
10 #include "ieee754sp.h"
11 #include "ieee754dp.h"
12 
ieee754dp_nan_fsp(int xs,u64 xm)13 static inline union ieee754dp ieee754dp_nan_fsp(int xs, u64 xm)
14 {
15 	return builddp(xs, DP_EMAX + 1 + DP_EBIAS,
16 		       xm << (DP_FBITS - SP_FBITS));
17 }
18 
ieee754dp_fsp(union ieee754sp x)19 union ieee754dp ieee754dp_fsp(union ieee754sp x)
20 {
21 	COMPXSP;
22 
23 	EXPLODEXSP;
24 
25 	ieee754_clearcx();
26 
27 	FLUSHXSP;
28 
29 	switch (xc) {
30 	case IEEE754_CLASS_SNAN:
31 		return ieee754dp_nanxcpt(ieee754dp_nan_fsp(xs, xm));
32 
33 	case IEEE754_CLASS_QNAN:
34 		return ieee754dp_nan_fsp(xs, xm);
35 
36 	case IEEE754_CLASS_INF:
37 		return ieee754dp_inf(xs);
38 
39 	case IEEE754_CLASS_ZERO:
40 		return ieee754dp_zero(xs);
41 
42 	case IEEE754_CLASS_DNORM:
43 		/* normalize */
44 		while ((xm >> SP_FBITS) == 0) {
45 			xm <<= 1;
46 			xe--;
47 		}
48 		break;
49 
50 	case IEEE754_CLASS_NORM:
51 		break;
52 	}
53 
54 	/*
55 	 * Can't possibly overflow,underflow, or need rounding
56 	 */
57 
58 	/* drop the hidden bit */
59 	xm &= ~SP_HIDDEN_BIT;
60 
61 	return builddp(xs, xe + DP_EBIAS,
62 		       (u64) xm << (DP_FBITS - SP_FBITS));
63 }
64