1 /* GCC Quad-Precision Math Library
2 Copyright (C) 2010 Free Software Foundation, Inc.
3 Written by Francois-Xavier Coudert <fxcoudert@gcc.gnu.org>
5 This file is part of the libiberty library.
6 Libiberty is free software; you can redistribute it and/or
7 modify it under the terms of the GNU Library General Public
8 License as published by the Free Software Foundation; either
9 version 2 of the License, or (at your option) any later version.
11 Libiberty is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 Library General Public License for more details.
16 You should have received a copy of the GNU Library General Public
17 License along with libiberty; see the file COPYING.LIB. If
18 not, write to the Free Software Foundation, Inc., 51 Franklin Street - Fifth Floor,
19 Boston, MA 02110-1301, USA. */
21 #ifndef QUADMATH_IMP_H
22 #define QUADMATH_IMP_H
30 // Prototypes for internal functions
31 extern int32_t __quadmath_rem_pio2q (__float128, __float128 *);
32 extern void __quadmath_kernel_sincosq (__float128, __float128, __float128 *,
34 extern __float128 __quadmath_kernel_sinq (__float128, __float128, int);
35 extern __float128 __quadmath_kernel_cosq (__float128, __float128);
39 // Frankly, if you have __float128, you have 64-bit integers, right?
45 // If we don't have macros to know endianess, assume little endian
46 #if !defined(__BIG_ENDIAN__) && !defined(__LITTLE_ENDIAN__)
47 # define __LITTLE_ENDIAN__ 1
51 // Main union type we use to manipulate the floating-point type
61 uint64_t mant_high:48;
66 uint64_t mant_high:48;
104 unsigned exponent:15;
105 unsigned quiet_nan:1;
106 uint64_t mant_high:47;
107 uint64_t mant_low:64;
109 #if __LITTLE_ENDIAN__
110 uint64_t mant_low:64;
111 uint64_t mant_high:47;
112 unsigned quiet_nan:1;
113 unsigned exponent:15;
121 /* Get two 64 bit ints from a long double. */
122 #define GET_FLT128_WORDS64(ix0,ix1,d) \
124 ieee854_float128 u; \
126 (ix0) = u.words64.high; \
127 (ix1) = u.words64.low; \
130 /* Set a long double from two 64 bit ints. */
131 #define SET_FLT128_WORDS64(d,ix0,ix1) \
133 ieee854_float128 u; \
134 u.words64.high = (ix0); \
135 u.words64.low = (ix1); \
139 /* Get the more significant 64 bits of a long double mantissa. */
140 #define GET_FLT128_MSW64(v,d) \
142 ieee854_float128 u; \
144 (v) = u.words64.high; \
147 /* Set the more significant 64 bits of a long double mantissa from an int. */
148 #define SET_FLT128_MSW64(d,v) \
150 ieee854_float128 u; \
152 u.words64.high = (v); \
156 /* Get the least significant 64 bits of a long double mantissa. */
157 #define GET_FLT128_LSW64(v,d) \
159 ieee854_float128 u; \
161 (v) = u.words64.low; \
165 #define IEEE854_FLOAT128_BIAS 0x3fff
168 #define QUADFP_INFINITE 1
169 #define QUADFP_ZERO 2
170 #define QUADFP_SUBNORMAL 3
171 #define QUADFP_NORMAL 4
172 #define fpclassifyq(x) \
173 __builtin_fpclassify (QUADFP_NAN, QUADFP_INFINITE, QUADFP_NORMAL, \
174 QUADFP_SUBNORMAL, QUADFP_ZERO, x)