File: /usr/src/linux/arch/ppc/math-emu/double.h
1 /*
2 * BK Id: SCCS/s.double.h 1.5 05/17/01 18:14:22 cort
3 */
4 /*
5 * Definitions for IEEE Double Precision
6 */
7
8 #if _FP_W_TYPE_SIZE < 32
9 #error "Here's a nickel kid. Go buy yourself a real computer."
10 #endif
11
12 #if _FP_W_TYPE_SIZE < 64
13 #define _FP_FRACTBITS_D (2 * _FP_W_TYPE_SIZE)
14 #else
15 #define _FP_FRACTBITS_D _FP_W_TYPE_SIZE
16 #endif
17
18 #define _FP_FRACBITS_D 53
19 #define _FP_FRACXBITS_D (_FP_FRACTBITS_D - _FP_FRACBITS_D)
20 #define _FP_WFRACBITS_D (_FP_WORKBITS + _FP_FRACBITS_D)
21 #define _FP_WFRACXBITS_D (_FP_FRACTBITS_D - _FP_WFRACBITS_D)
22 #define _FP_EXPBITS_D 11
23 #define _FP_EXPBIAS_D 1023
24 #define _FP_EXPMAX_D 2047
25
26 #define _FP_QNANBIT_D \
27 ((_FP_W_TYPE)1 << ((_FP_FRACBITS_D-2) % _FP_W_TYPE_SIZE))
28 #define _FP_IMPLBIT_D \
29 ((_FP_W_TYPE)1 << ((_FP_FRACBITS_D-1) % _FP_W_TYPE_SIZE))
30 #define _FP_OVERFLOW_D \
31 ((_FP_W_TYPE)1 << (_FP_WFRACBITS_D % _FP_W_TYPE_SIZE))
32
33 #if _FP_W_TYPE_SIZE < 64
34
35 union _FP_UNION_D
36 {
37 double flt;
38 struct {
39 #if __BYTE_ORDER == __BIG_ENDIAN
40 unsigned sign : 1;
41 unsigned exp : _FP_EXPBITS_D;
42 unsigned frac1 : _FP_FRACBITS_D - (_FP_IMPLBIT_D != 0) - _FP_W_TYPE_SIZE;
43 unsigned frac0 : _FP_W_TYPE_SIZE;
44 #else
45 unsigned frac0 : _FP_W_TYPE_SIZE;
46 unsigned frac1 : _FP_FRACBITS_D - (_FP_IMPLBIT_D != 0) - _FP_W_TYPE_SIZE;
47 unsigned exp : _FP_EXPBITS_D;
48 unsigned sign : 1;
49 #endif
50 } bits __attribute__((packed));
51 };
52
53 #define FP_DECL_D(X) _FP_DECL(2,X)
54 #define FP_UNPACK_RAW_D(X,val) _FP_UNPACK_RAW_2(D,X,val)
55 #define FP_PACK_RAW_D(val,X) _FP_PACK_RAW_2(D,val,X)
56
57 #define FP_UNPACK_D(X,val) \
58 do { \
59 _FP_UNPACK_RAW_2(D,X,val); \
60 _FP_UNPACK_CANONICAL(D,2,X); \
61 } while (0)
62
63 #define FP_PACK_D(val,X) \
64 do { \
65 _FP_PACK_CANONICAL(D,2,X); \
66 _FP_PACK_RAW_2(D,val,X); \
67 } while (0)
68
69 #define FP_NEG_D(R,X) _FP_NEG(D,2,R,X)
70 #define FP_ADD_D(R,X,Y) _FP_ADD(D,2,R,X,Y)
71 #define FP_SUB_D(R,X,Y) _FP_SUB(D,2,R,X,Y)
72 #define FP_MUL_D(R,X,Y) _FP_MUL(D,2,R,X,Y)
73 #define FP_DIV_D(R,X,Y) _FP_DIV(D,2,R,X,Y)
74 #define FP_SQRT_D(R,X) _FP_SQRT(D,2,R,X)
75
76 #define FP_CMP_D(r,X,Y,un) _FP_CMP(D,2,r,X,Y,un)
77 #define FP_CMP_EQ_D(r,X,Y) _FP_CMP_EQ(D,2,r,X,Y)
78
79 #define FP_TO_INT_D(r,X,rsz,rsg) _FP_TO_INT(D,2,r,X,rsz,rsg)
80 #define FP_FROM_INT_D(X,r,rs,rt) _FP_FROM_INT(D,2,X,r,rs,rt)
81
82 #else
83
84 union _FP_UNION_D
85 {
86 double flt;
87 struct {
88 #if __BYTE_ORDER == __BIG_ENDIAN
89 unsigned sign : 1;
90 unsigned exp : _FP_EXPBITS_D;
91 unsigned long frac : _FP_FRACBITS_D - (_FP_IMPLBIT_D != 0);
92 #else
93 unsigned long frac : _FP_FRACBITS_D - (_FP_IMPLBIT_D != 0);
94 unsigned exp : _FP_EXPBITS_D;
95 unsigned sign : 1;
96 #endif
97 } bits __attribute__((packed));
98 };
99
100 #define FP_DECL_D(X) _FP_DECL(1,X)
101 #define FP_UNPACK_RAW_D(X,val) _FP_UNPACK_RAW_1(D,X,val)
102 #define FP_PACK_RAW_D(val,X) _FP_PACK_RAW_1(D,val,X)
103
104 #define FP_UNPACK_D(X,val) \
105 do { \
106 _FP_UNPACK_RAW_1(D,X,val); \
107 _FP_UNPACK_CANONICAL(D,1,X); \
108 } while (0)
109
110 #define FP_PACK_D(val,X) \
111 do { \
112 _FP_PACK_CANONICAL(D,1,X); \
113 _FP_PACK_RAW_1(D,val,X); \
114 } while (0)
115
116 #define FP_NEG_D(R,X) _FP_NEG(D,1,R,X)
117 #define FP_ADD_D(R,X,Y) _FP_ADD(D,1,R,X,Y)
118 #define FP_SUB_D(R,X,Y) _FP_SUB(D,1,R,X,Y)
119 #define FP_MUL_D(R,X,Y) _FP_MUL(D,1,R,X,Y)
120 #define FP_DIV_D(R,X,Y) _FP_DIV(D,1,R,X,Y)
121 #define FP_SQRT_D(R,X) _FP_SQRT(D,1,R,X)
122
123 /* The implementation of _FP_MUL_D and _FP_DIV_D should be chosen by
124 the target machine. */
125
126 #define FP_CMP_D(r,X,Y,un) _FP_CMP(D,1,r,X,Y,un)
127 #define FP_CMP_EQ_D(r,X,Y) _FP_CMP_EQ(D,1,r,X,Y)
128
129 #define FP_TO_INT_D(r,X,rsz,rsg) _FP_TO_INT(D,1,r,X,rsz,rsg)
130 #define FP_FROM_INT_D(X,r,rs,rt) _FP_FROM_INT(D,1,X,r,rs,rt)
131
132 #endif /* W_TYPE_SIZE < 64 */
133