xref: /openbmc/linux/arch/mips/math-emu/ieee754.c (revision 9d5a6349)
1 // SPDX-License-Identifier: GPL-2.0-only
2 /* ieee754 floating point arithmetic
3  * single and double precision
4  *
5  * BUGS
6  * not much dp done
7  * doesn't generate IEEE754_INEXACT
8  */
9 /*
10  * MIPS floating point support
11  * Copyright (C) 1994-2000 Algorithmics Ltd.
12  */
13 
14 #include <linux/compiler.h>
15 
16 #include "ieee754.h"
17 #include "ieee754sp.h"
18 #include "ieee754dp.h"
19 
20 /*
21  * Special constants
22  */
23 
24 /*
25  * Older GCC requires the inner braces for initialization of union ieee754dp's
26  * anonymous struct member.  Without an error will result.
27  */
28 #define xPCNST(s, b, m, ebias)						\
29 {									\
30 	{								\
31 		.sign	= (s),						\
32 		.bexp	= (b) + ebias,					\
33 		.mant	= (m)						\
34 	}								\
35 }
36 
37 #define DPCNST(s, b, m)							\
38 	xPCNST(s, b, m, DP_EBIAS)
39 
40 const union ieee754dp __ieee754dp_spcvals[] = {
41 	DPCNST(0, DP_EMIN - 1, 0x0000000000000ULL),	/* + zero   */
42 	DPCNST(1, DP_EMIN - 1, 0x0000000000000ULL),	/* - zero   */
43 	DPCNST(0, 0,	       0x0000000000000ULL),	/* + 1.0   */
44 	DPCNST(1, 0,	       0x0000000000000ULL),	/* - 1.0   */
45 	DPCNST(0, 3,           0x4000000000000ULL),	/* + 10.0   */
46 	DPCNST(1, 3,           0x4000000000000ULL),	/* - 10.0   */
47 	DPCNST(0, DP_EMAX + 1, 0x0000000000000ULL),	/* + infinity */
48 	DPCNST(1, DP_EMAX + 1, 0x0000000000000ULL),	/* - infinity */
49 	DPCNST(0, DP_EMAX + 1, 0x7FFFFFFFFFFFFULL),	/* + ind legacy qNaN */
50 	DPCNST(0, DP_EMAX + 1, 0x8000000000000ULL),	/* + indef 2008 qNaN */
51 	DPCNST(0, DP_EMAX,     0xFFFFFFFFFFFFFULL),	/* + max */
52 	DPCNST(1, DP_EMAX,     0xFFFFFFFFFFFFFULL),	/* - max */
53 	DPCNST(0, DP_EMIN,     0x0000000000000ULL),	/* + min normal */
54 	DPCNST(1, DP_EMIN,     0x0000000000000ULL),	/* - min normal */
55 	DPCNST(0, DP_EMIN - 1, 0x0000000000001ULL),	/* + min denormal */
56 	DPCNST(1, DP_EMIN - 1, 0x0000000000001ULL),	/* - min denormal */
57 	DPCNST(0, 31,          0x0000000000000ULL),	/* + 1.0e31 */
58 	DPCNST(0, 63,          0x0000000000000ULL),	/* + 1.0e63 */
59 };
60 
61 #define SPCNST(s, b, m)							\
62 	xPCNST(s, b, m, SP_EBIAS)
63 
64 const union ieee754sp __ieee754sp_spcvals[] = {
65 	SPCNST(0, SP_EMIN - 1, 0x000000),	/* + zero   */
66 	SPCNST(1, SP_EMIN - 1, 0x000000),	/* - zero   */
67 	SPCNST(0, 0,	       0x000000),	/* + 1.0   */
68 	SPCNST(1, 0,	       0x000000),	/* - 1.0   */
69 	SPCNST(0, 3,	       0x200000),	/* + 10.0   */
70 	SPCNST(1, 3,	       0x200000),	/* - 10.0   */
71 	SPCNST(0, SP_EMAX + 1, 0x000000),	/* + infinity */
72 	SPCNST(1, SP_EMAX + 1, 0x000000),	/* - infinity */
73 	SPCNST(0, SP_EMAX + 1, 0x3FFFFF),	/* + indef legacy quiet NaN */
74 	SPCNST(0, SP_EMAX + 1, 0x400000),	/* + indef 2008 quiet NaN */
75 	SPCNST(0, SP_EMAX,     0x7FFFFF),	/* + max normal */
76 	SPCNST(1, SP_EMAX,     0x7FFFFF),	/* - max normal */
77 	SPCNST(0, SP_EMIN,     0x000000),	/* + min normal */
78 	SPCNST(1, SP_EMIN,     0x000000),	/* - min normal */
79 	SPCNST(0, SP_EMIN - 1, 0x000001),	/* + min denormal */
80 	SPCNST(1, SP_EMIN - 1, 0x000001),	/* - min denormal */
81 	SPCNST(0, 31,	       0x000000),	/* + 1.0e31 */
82 	SPCNST(0, 63,	       0x000000),	/* + 1.0e63 */
83 };
84