1 /* Software floating-point emulation. 2 Definitions for IEEE Double Precision 3 Copyright (C) 1997, 1998, 1999, 2006, 2007, 2008, 2009, 2012 4 Free Software Foundation, Inc. 5 This file is part of the GNU C Library. 6 Contributed by Richard Henderson (rth@cygnus.com), 7 Jakub Jelinek (jj@ultra.linux.cz), 8 David S. Miller (davem@redhat.com) and 9 Peter Maydell (pmaydell@chiark.greenend.org.uk). 10 11 The GNU C Library is free software; you can redistribute it and/or 12 modify it under the terms of the GNU Lesser General Public 13 License as published by the Free Software Foundation; either 14 version 2.1 of the License, or (at your option) any later version. 15 16 In addition to the permissions in the GNU Lesser General Public 17 License, the Free Software Foundation gives you unlimited 18 permission to link the compiled version of this file into 19 combinations with other programs, and to distribute those 20 combinations without any restriction coming from the use of this 21 file. (The Lesser General Public License restrictions do apply in 22 other respects; for example, they cover modification of the file, 23 and distribution when not linked into a combine executable.) 24 25 The GNU C Library is distributed in the hope that it will be useful, 26 but WITHOUT ANY WARRANTY; without even the implied warranty of 27 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 28 Lesser General Public License for more details. 29 30 You should have received a copy of the GNU Lesser General Public 31 License along with the GNU C Library; if not, see 32 <http://www.gnu.org/licenses/>. */ 33 34 #if _FP_W_TYPE_SIZE < 32 35 #error "Here's a nickel kid. Go buy yourself a real computer." 36 #endif 37 38 #if _FP_W_TYPE_SIZE < 64 39 #define _FP_FRACTBITS_D (2 * _FP_W_TYPE_SIZE) 40 #else 41 #define _FP_FRACTBITS_D _FP_W_TYPE_SIZE 42 #endif 43 44 #define _FP_FRACBITS_D 53 45 #define _FP_FRACXBITS_D (_FP_FRACTBITS_D - _FP_FRACBITS_D) 46 #define _FP_WFRACBITS_D (_FP_WORKBITS + _FP_FRACBITS_D) 47 #define _FP_WFRACXBITS_D (_FP_FRACTBITS_D - _FP_WFRACBITS_D) 48 #define _FP_EXPBITS_D 11 49 #define _FP_EXPBIAS_D 1023 50 #define _FP_EXPMAX_D 2047 51 52 #define _FP_QNANBIT_D \ 53 ((_FP_W_TYPE)1 << (_FP_FRACBITS_D-2) % _FP_W_TYPE_SIZE) 54 #define _FP_QNANBIT_SH_D \ 55 ((_FP_W_TYPE)1 << (_FP_FRACBITS_D-2+_FP_WORKBITS) % _FP_W_TYPE_SIZE) 56 #define _FP_IMPLBIT_D \ 57 ((_FP_W_TYPE)1 << (_FP_FRACBITS_D-1) % _FP_W_TYPE_SIZE) 58 #define _FP_IMPLBIT_SH_D \ 59 ((_FP_W_TYPE)1 << (_FP_FRACBITS_D-1+_FP_WORKBITS) % _FP_W_TYPE_SIZE) 60 #define _FP_OVERFLOW_D \ 61 ((_FP_W_TYPE)1 << _FP_WFRACBITS_D % _FP_W_TYPE_SIZE) 62 63 typedef float DFtype __attribute__((mode(DF))); 64 65 #if _FP_W_TYPE_SIZE < 64 66 67 union _FP_UNION_D 68 { 69 DFtype flt; 70 struct _FP_STRUCT_LAYOUT { 71 #if __BYTE_ORDER == __BIG_ENDIAN 72 unsigned sign : 1; 73 unsigned exp : _FP_EXPBITS_D; 74 unsigned frac1 : _FP_FRACBITS_D - (_FP_IMPLBIT_D != 0) - _FP_W_TYPE_SIZE; 75 unsigned frac0 : _FP_W_TYPE_SIZE; 76 #else 77 unsigned frac0 : _FP_W_TYPE_SIZE; 78 unsigned frac1 : _FP_FRACBITS_D - (_FP_IMPLBIT_D != 0) - _FP_W_TYPE_SIZE; 79 unsigned exp : _FP_EXPBITS_D; 80 unsigned sign : 1; 81 #endif 82 } bits __attribute__((packed)); 83 }; 84 85 #define FP_DECL_D(X) _FP_DECL(2,X) 86 #define FP_UNPACK_RAW_D(X,val) _FP_UNPACK_RAW_2(D,X,val) 87 #define FP_UNPACK_RAW_DP(X,val) _FP_UNPACK_RAW_2_P(D,X,val) 88 #define FP_PACK_RAW_D(val,X) _FP_PACK_RAW_2(D,val,X) 89 #define FP_PACK_RAW_DP(val,X) \ 90 do { \ 91 if (!FP_INHIBIT_RESULTS) \ 92 _FP_PACK_RAW_2_P(D,val,X); \ 93 } while (0) 94 95 #define FP_UNPACK_D(X,val) \ 96 do { \ 97 _FP_UNPACK_RAW_2(D,X,val); \ 98 _FP_UNPACK_CANONICAL(D,2,X); \ 99 } while (0) 100 101 #define FP_UNPACK_DP(X,val) \ 102 do { \ 103 _FP_UNPACK_RAW_2_P(D,X,val); \ 104 _FP_UNPACK_CANONICAL(D,2,X); \ 105 } while (0) 106 107 #define FP_UNPACK_SEMIRAW_D(X,val) \ 108 do { \ 109 _FP_UNPACK_RAW_2(D,X,val); \ 110 _FP_UNPACK_SEMIRAW(D,2,X); \ 111 } while (0) 112 113 #define FP_UNPACK_SEMIRAW_DP(X,val) \ 114 do { \ 115 _FP_UNPACK_RAW_2_P(D,X,val); \ 116 _FP_UNPACK_SEMIRAW(D,2,X); \ 117 } while (0) 118 119 #define FP_PACK_D(val,X) \ 120 do { \ 121 _FP_PACK_CANONICAL(D,2,X); \ 122 _FP_PACK_RAW_2(D,val,X); \ 123 } while (0) 124 125 #define FP_PACK_DP(val,X) \ 126 do { \ 127 _FP_PACK_CANONICAL(D,2,X); \ 128 if (!FP_INHIBIT_RESULTS) \ 129 _FP_PACK_RAW_2_P(D,val,X); \ 130 } while (0) 131 132 #define FP_PACK_SEMIRAW_D(val,X) \ 133 do { \ 134 _FP_PACK_SEMIRAW(D,2,X); \ 135 _FP_PACK_RAW_2(D,val,X); \ 136 } while (0) 137 138 #define FP_PACK_SEMIRAW_DP(val,X) \ 139 do { \ 140 _FP_PACK_SEMIRAW(D,2,X); \ 141 if (!FP_INHIBIT_RESULTS) \ 142 _FP_PACK_RAW_2_P(D,val,X); \ 143 } while (0) 144 145 #define FP_ISSIGNAN_D(X) _FP_ISSIGNAN(D,2,X) 146 #define FP_NEG_D(R,X) _FP_NEG(D,2,R,X) 147 #define FP_ADD_D(R,X,Y) _FP_ADD(D,2,R,X,Y) 148 #define FP_SUB_D(R,X,Y) _FP_SUB(D,2,R,X,Y) 149 #define FP_MUL_D(R,X,Y) _FP_MUL(D,2,R,X,Y) 150 #define FP_DIV_D(R,X,Y) _FP_DIV(D,2,R,X,Y) 151 #define FP_SQRT_D(R,X) _FP_SQRT(D,2,R,X) 152 #define _FP_SQRT_MEAT_D(R,S,T,X,Q) _FP_SQRT_MEAT_2(R,S,T,X,Q) 153 154 #define FP_CMP_D(r,X,Y,un) _FP_CMP(D,2,r,X,Y,un) 155 #define FP_CMP_EQ_D(r,X,Y) _FP_CMP_EQ(D,2,r,X,Y) 156 #define FP_CMP_UNORD_D(r,X,Y) _FP_CMP_UNORD(D,2,r,X,Y) 157 158 #define FP_TO_INT_D(r,X,rsz,rsg) _FP_TO_INT(D,2,r,X,rsz,rsg) 159 #define FP_FROM_INT_D(X,r,rs,rt) _FP_FROM_INT(D,2,X,r,rs,rt) 160 161 #define _FP_FRAC_HIGH_D(X) _FP_FRAC_HIGH_2(X) 162 #define _FP_FRAC_HIGH_RAW_D(X) _FP_FRAC_HIGH_2(X) 163 164 #else 165 166 union _FP_UNION_D 167 { 168 DFtype flt; 169 struct _FP_STRUCT_LAYOUT { 170 #if __BYTE_ORDER == __BIG_ENDIAN 171 unsigned sign : 1; 172 unsigned exp : _FP_EXPBITS_D; 173 _FP_W_TYPE frac : _FP_FRACBITS_D - (_FP_IMPLBIT_D != 0); 174 #else 175 _FP_W_TYPE frac : _FP_FRACBITS_D - (_FP_IMPLBIT_D != 0); 176 unsigned exp : _FP_EXPBITS_D; 177 unsigned sign : 1; 178 #endif 179 } bits __attribute__((packed)); 180 }; 181 182 #define FP_DECL_D(X) _FP_DECL(1,X) 183 #define FP_UNPACK_RAW_D(X,val) _FP_UNPACK_RAW_1(D,X,val) 184 #define FP_UNPACK_RAW_DP(X,val) _FP_UNPACK_RAW_1_P(D,X,val) 185 #define FP_PACK_RAW_D(val,X) _FP_PACK_RAW_1(D,val,X) 186 #define FP_PACK_RAW_DP(val,X) \ 187 do { \ 188 if (!FP_INHIBIT_RESULTS) \ 189 _FP_PACK_RAW_1_P(D,val,X); \ 190 } while (0) 191 192 #define FP_UNPACK_D(X,val) \ 193 do { \ 194 _FP_UNPACK_RAW_1(D,X,val); \ 195 _FP_UNPACK_CANONICAL(D,1,X); \ 196 } while (0) 197 198 #define FP_UNPACK_DP(X,val) \ 199 do { \ 200 _FP_UNPACK_RAW_1_P(D,X,val); \ 201 _FP_UNPACK_CANONICAL(D,1,X); \ 202 } while (0) 203 204 #define FP_UNPACK_SEMIRAW_D(X,val) \ 205 do { \ 206 _FP_UNPACK_RAW_1(D,X,val); \ 207 _FP_UNPACK_SEMIRAW(D,1,X); \ 208 } while (0) 209 210 #define FP_UNPACK_SEMIRAW_DP(X,val) \ 211 do { \ 212 _FP_UNPACK_RAW_1_P(D,X,val); \ 213 _FP_UNPACK_SEMIRAW(D,1,X); \ 214 } while (0) 215 216 #define FP_PACK_D(val,X) \ 217 do { \ 218 _FP_PACK_CANONICAL(D,1,X); \ 219 _FP_PACK_RAW_1(D,val,X); \ 220 } while (0) 221 222 #define FP_PACK_DP(val,X) \ 223 do { \ 224 _FP_PACK_CANONICAL(D,1,X); \ 225 if (!FP_INHIBIT_RESULTS) \ 226 _FP_PACK_RAW_1_P(D,val,X); \ 227 } while (0) 228 229 #define FP_PACK_SEMIRAW_D(val,X) \ 230 do { \ 231 _FP_PACK_SEMIRAW(D,1,X); \ 232 _FP_PACK_RAW_1(D,val,X); \ 233 } while (0) 234 235 #define FP_PACK_SEMIRAW_DP(val,X) \ 236 do { \ 237 _FP_PACK_SEMIRAW(D,1,X); \ 238 if (!FP_INHIBIT_RESULTS) \ 239 _FP_PACK_RAW_1_P(D,val,X); \ 240 } while (0) 241 242 #define FP_ISSIGNAN_D(X) _FP_ISSIGNAN(D,1,X) 243 #define FP_NEG_D(R,X) _FP_NEG(D,1,R,X) 244 #define FP_ADD_D(R,X,Y) _FP_ADD(D,1,R,X,Y) 245 #define FP_SUB_D(R,X,Y) _FP_SUB(D,1,R,X,Y) 246 #define FP_MUL_D(R,X,Y) _FP_MUL(D,1,R,X,Y) 247 #define FP_DIV_D(R,X,Y) _FP_DIV(D,1,R,X,Y) 248 #define FP_SQRT_D(R,X) _FP_SQRT(D,1,R,X) 249 #define _FP_SQRT_MEAT_D(R,S,T,X,Q) _FP_SQRT_MEAT_1(R,S,T,X,Q) 250 251 /* The implementation of _FP_MUL_D and _FP_DIV_D should be chosen by 252 the target machine. */ 253 254 #define FP_CMP_D(r,X,Y,un) _FP_CMP(D,1,r,X,Y,un) 255 #define FP_CMP_EQ_D(r,X,Y) _FP_CMP_EQ(D,1,r,X,Y) 256 #define FP_CMP_UNORD_D(r,X,Y) _FP_CMP_UNORD(D,1,r,X,Y) 257 258 #define FP_TO_INT_D(r,X,rsz,rsg) _FP_TO_INT(D,1,r,X,rsz,rsg) 259 #define FP_FROM_INT_D(X,r,rs,rt) _FP_FROM_INT(D,1,X,r,rs,rt) 260 261 #define _FP_FRAC_HIGH_D(X) _FP_FRAC_HIGH_1(X) 262 #define _FP_FRAC_HIGH_RAW_D(X) _FP_FRAC_HIGH_1(X) 263 264 #endif /* W_TYPE_SIZE < 64 */ 265