Mercurial > hg > CbC > CbC_gcc

annotate libquadmath/quadmath-imp.h @ 158:494b0b89df80 default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
...
author Shinji KONO <kono@ie.u-ryukyu.ac.jp>
date Mon, 25 May 2020 18:13:55 +0900
parents 1830386684a0
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
1 /* GCC Quad-Precision Math Library
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
2 Copyright (C) 2010, 2011 Free Software Foundation, Inc.
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
3 Written by Francois-Xavier Coudert <fxcoudert@gcc.gnu.org>
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
4
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
5 This file is part of the libquadmath library.
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
6 Libquadmath is free software; you can redistribute it and/or
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
7 modify it under the terms of the GNU Library General Public
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
8 License as published by the Free Software Foundation; either
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
9 version 2 of the License, or (at your option) any later version.
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
10
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
11 Libquadmath is distributed in the hope that it will be useful,
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
14 Library General Public License for more details.
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
15
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
16 You should have received a copy of the GNU Library General Public
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
17 License along with libquadmath; see the file COPYING.LIB. If
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
18 not, write to the Free Software Foundation, Inc., 51 Franklin Street - Fifth Floor,
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
19 Boston, MA 02110-1301, USA. */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
20
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
21 #ifndef QUADMATH_IMP_H
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
22 #define QUADMATH_IMP_H
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
23
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
24 #include <errno.h>
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
25 #include <limits.h>
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
26 #include <stdbool.h>
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
27 #include <stdint.h>
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
28 #include <stdlib.h>
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
29 #include "quadmath.h"
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
30 #include "config.h"
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
31 #ifdef HAVE_FENV_H
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
32 # include <fenv.h>
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
33 #endif
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
34
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
35
111
kono
parents: 68
diff changeset
36 /* Under IEEE 754, an architecture may determine tininess of
kono
parents: 68
diff changeset
37 floating-point results either "before rounding" or "after
kono
parents: 68
diff changeset
38 rounding", but must do so in the same way for all operations
kono
parents: 68
diff changeset
39 returning binary results. Define TININESS_AFTER_ROUNDING to 1 for
kono
parents: 68
diff changeset
40 "after rounding" architectures, 0 for "before rounding"
kono
parents: 68
diff changeset
41 architectures. */
kono
parents: 68
diff changeset
42
kono
parents: 68
diff changeset
43 #define TININESS_AFTER_ROUNDING 1
kono
parents: 68
diff changeset
44
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
45 #define HIGH_ORDER_BIT_IS_SET_FOR_SNAN 0
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
46
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
47 #define FIX_FLT128_LONG_CONVERT_OVERFLOW 0
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
48 #define FIX_FLT128_LLONG_CONVERT_OVERFLOW 0
111
kono
parents: 68
diff changeset
49
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
50 /* Prototypes for internal functions. */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
51 extern int32_t __quadmath_rem_pio2q (__float128, __float128 *);
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
52 extern void __quadmath_kernel_sincosq (__float128, __float128, __float128 *,
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
53 __float128 *, int);
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
54 extern __float128 __quadmath_kernel_sinq (__float128, __float128, int);
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
55 extern __float128 __quadmath_kernel_cosq (__float128, __float128);
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
56 extern __float128 __quadmath_kernel_tanq (__float128, __float128, int);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
57 extern __float128 __quadmath_gamma_productq (__float128, __float128, int,
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
58 __float128 *);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
59 extern __float128 __quadmath_gammaq_r (__float128, int *);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
60 extern __float128 __quadmath_lgamma_negq (__float128, int *);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
61 extern __float128 __quadmath_lgamma_productq (__float128, __float128,
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
62 __float128, int);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
63 extern __float128 __quadmath_lgammaq_r (__float128, int *);
111
kono
parents: 68
diff changeset
64 extern __float128 __quadmath_x2y2m1q (__float128 x, __float128 y);
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
65 extern __complex128 __quadmath_kernel_casinhq (__complex128, int);
111
kono
parents: 68
diff changeset
66
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
67 static inline void
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
68 mul_splitq (__float128 *hi, __float128 *lo, __float128 x, __float128 y)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
69 {
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
70 /* Fast built-in fused multiply-add. */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
71 *hi = x * y;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
72 *lo = fmaq (x, y, -*hi);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
73 }
111
kono
parents: 68
diff changeset
74
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
75
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
76
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
77
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
78 /* Frankly, if you have __float128, you have 64-bit integers, right? */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
79 #ifndef UINT64_C
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
80 # error "No way!"
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
81 #endif
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
82
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
83
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
84 /* Main union type we use to manipulate the floating-point type. */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
85 typedef union
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
86 {
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
87 __float128 value;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
88
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
89 struct
111
kono
parents: 68
diff changeset
90 #ifdef __MINGW32__
kono
parents: 68
diff changeset
91 /* On mingw targets the ms-bitfields option is active by default.
kono
parents: 68
diff changeset
92 Therefore enforce gnu-bitfield style. */
kono
parents: 68
diff changeset
93 __attribute__ ((gcc_struct))
kono
parents: 68
diff changeset
94 #endif
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
95 {
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
96 #if __BYTE_ORDER__ == __ORDER_BIG_ENDIAN__
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
97 unsigned negative:1;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
98 unsigned exponent:15;
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
99 unsigned mantissa0:16;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
100 unsigned mantissa1:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
101 unsigned mantissa2:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
102 unsigned mantissa3:32;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
103 #else
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
104 unsigned mantissa3:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
105 unsigned mantissa2:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
106 unsigned mantissa1:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
107 unsigned mantissa0:16;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
108 unsigned exponent:15;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
109 unsigned negative:1;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
110 #endif
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
111 } ieee;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
112
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
113 struct
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
114 {
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
115 #if __BYTE_ORDER__ == __ORDER_BIG_ENDIAN__
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
116 uint64_t high;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
117 uint64_t low;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
118 #else
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
119 uint64_t low;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
120 uint64_t high;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
121 #endif
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
122 } words64;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
123
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
124 struct
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
125 {
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
126 #if __BYTE_ORDER__ == __ORDER_BIG_ENDIAN__
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
127 uint32_t w0;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
128 uint32_t w1;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
129 uint32_t w2;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
130 uint32_t w3;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
131 #else
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
132 uint32_t w3;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
133 uint32_t w2;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
134 uint32_t w1;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
135 uint32_t w0;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
136 #endif
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
137 } words32;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
138
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
139 struct
111
kono
parents: 68
diff changeset
140 #ifdef __MINGW32__
kono
parents: 68
diff changeset
141 /* Make sure we are using gnu-style bitfield handling. */
kono
parents: 68
diff changeset
142 __attribute__ ((gcc_struct))
kono
parents: 68
diff changeset
143 #endif
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
144 {
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
145 #if __BYTE_ORDER__ == __ORDER_BIG_ENDIAN__
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
146 unsigned negative:1;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
147 unsigned exponent:15;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
148 unsigned quiet_nan:1;
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
149 unsigned mantissa0:15;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
150 unsigned mantissa1:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
151 unsigned mantissa2:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
152 unsigned mantissa3:32;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
153 #else
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
154 unsigned mantissa3:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
155 unsigned mantissa2:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
156 unsigned mantissa1:32;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
157 unsigned mantissa0:15;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
158 unsigned quiet_nan:1;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
159 unsigned exponent:15;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
160 unsigned negative:1;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
161 #endif
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
162 } ieee_nan;
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
163
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
164 } ieee854_float128;
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
165
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
166
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
167 /* Get two 64 bit ints from a long double. */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
168 #define GET_FLT128_WORDS64(ix0,ix1,d) \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
169 do { \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
170 ieee854_float128 u; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
171 u.value = (d); \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
172 (ix0) = u.words64.high; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
173 (ix1) = u.words64.low; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
174 } while (0)
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
175
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
176 /* Set a long double from two 64 bit ints. */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
177 #define SET_FLT128_WORDS64(d,ix0,ix1) \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
178 do { \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
179 ieee854_float128 u; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
180 u.words64.high = (ix0); \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
181 u.words64.low = (ix1); \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
182 (d) = u.value; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
183 } while (0)
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
184
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
185 /* Get the more significant 64 bits of a long double mantissa. */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
186 #define GET_FLT128_MSW64(v,d) \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
187 do { \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
188 ieee854_float128 u; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
189 u.value = (d); \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
190 (v) = u.words64.high; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
191 } while (0)
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
192
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
193 /* Set the more significant 64 bits of a long double mantissa from an int. */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
194 #define SET_FLT128_MSW64(d,v) \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
195 do { \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
196 ieee854_float128 u; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
197 u.value = (d); \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
198 u.words64.high = (v); \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
199 (d) = u.value; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
200 } while (0)
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
201
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
202 /* Get the least significant 64 bits of a long double mantissa. */
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
203 #define GET_FLT128_LSW64(v,d) \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
204 do { \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
205 ieee854_float128 u; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
206 u.value = (d); \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
207 (v) = u.words64.low; \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
208 } while (0)
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
209
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
210
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
211 #define IEEE854_FLOAT128_BIAS 0x3fff
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
212
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
213 #define QUADFP_NAN 0
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
214 #define QUADFP_INFINITE 1
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
215 #define QUADFP_ZERO 2
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
216 #define QUADFP_SUBNORMAL 3
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
217 #define QUADFP_NORMAL 4
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
218 #define fpclassifyq(x) \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
219 __builtin_fpclassify (QUADFP_NAN, QUADFP_INFINITE, QUADFP_NORMAL, \
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
220 QUADFP_SUBNORMAL, QUADFP_ZERO, x)
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
221
111
kono
parents: 68
diff changeset
222 #ifndef math_opt_barrier
kono
parents: 68
diff changeset
223 # define math_opt_barrier(x) \
kono
parents: 68
diff changeset
224 ({ __typeof (x) __x = (x); __asm ("" : "+m" (__x)); __x; })
kono
parents: 68
diff changeset
225 # define math_force_eval(x) \
kono
parents: 68
diff changeset
226 ({ __typeof (x) __x = (x); __asm __volatile__ ("" : : "m" (__x)); })
68
561a7518be6b update gcc-4.6
Nobuyasu Oshiro <dimolto@cr.ie.u-ryukyu.ac.jp>
parents:
diff changeset
227 #endif
111
kono
parents: 68
diff changeset
228
kono
parents: 68
diff changeset
229 /* math_narrow_eval reduces its floating-point argument to the range
kono
parents: 68
diff changeset
230 and precision of its semantic type. (The original evaluation may
kono
parents: 68
diff changeset
231 still occur with excess range and precision, so the result may be
kono
parents: 68
diff changeset
232 affected by double rounding.) */
kono
parents: 68
diff changeset
233 #define math_narrow_eval(x) (x)
kono
parents: 68
diff changeset
234
kono
parents: 68
diff changeset
235 /* If X (which is not a NaN) is subnormal, force an underflow
kono
parents: 68
diff changeset
236 exception. */
kono
parents: 68
diff changeset
237 #define math_check_force_underflow(x) \
kono
parents: 68
diff changeset
238 do \
kono
parents: 68
diff changeset
239 { \
kono
parents: 68
diff changeset
240 __float128 force_underflow_tmp = (x); \
kono
parents: 68
diff changeset
241 if (fabsq (force_underflow_tmp) < FLT128_MIN) \
kono
parents: 68
diff changeset
242 { \
kono
parents: 68
diff changeset
243 __float128 force_underflow_tmp2 \
kono
parents: 68
diff changeset
244 = force_underflow_tmp * force_underflow_tmp; \
kono
parents: 68
diff changeset
245 math_force_eval (force_underflow_tmp2); \
kono
parents: 68
diff changeset
246 } \
kono
parents: 68
diff changeset
247 } \
kono
parents: 68
diff changeset
248 while (0)
kono
parents: 68
diff changeset
249 /* Likewise, but X is also known to be nonnegative. */
kono
parents: 68
diff changeset
250 #define math_check_force_underflow_nonneg(x) \
kono
parents: 68
diff changeset
251 do \
kono
parents: 68
diff changeset
252 { \
kono
parents: 68
diff changeset
253 __float128 force_underflow_tmp = (x); \
kono
parents: 68
diff changeset
254 if (force_underflow_tmp < FLT128_MIN) \
kono
parents: 68
diff changeset
255 { \
kono
parents: 68
diff changeset
256 __float128 force_underflow_tmp2 \
kono
parents: 68
diff changeset
257 = force_underflow_tmp * force_underflow_tmp; \
kono
parents: 68
diff changeset
258 math_force_eval (force_underflow_tmp2); \
kono
parents: 68
diff changeset
259 } \
kono
parents: 68
diff changeset
260 } \
kono
parents: 68
diff changeset
261 while (0)
kono
parents: 68
diff changeset
262
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
263 /* Likewise, for both real and imaginary parts of a complex
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
264 result. */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
265 #define math_check_force_underflow_complex(x) \
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
266 do \
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
267 { \
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
268 __typeof (x) force_underflow_complex_tmp = (x); \
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
269 math_check_force_underflow (__real__ force_underflow_complex_tmp); \
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
270 math_check_force_underflow (__imag__ force_underflow_complex_tmp); \
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
271 } \
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
272 while (0)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
273
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
274 #ifndef HAVE_FENV_H
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
275 # define feraiseexcept(arg) ((void) 0)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
276 typedef int fenv_t;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
277 # define feholdexcept(arg) ((void) 0)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
278 # define fesetround(arg) ((void) 0)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
279 # define feupdateenv(arg) ((void) (arg))
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
280 # define fesetenv(arg) ((void) (arg))
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
281 # define fetestexcept(arg) 0
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
282 # define feclearexcept(arg) ((void) 0)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
283 #else
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
284 # ifndef HAVE_FEHOLDEXCEPT
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
285 # define feholdexcept(arg) ((void) 0)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
286 # endif
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
287 # ifndef HAVE_FESETROUND
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
288 # define fesetround(arg) ((void) 0)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
289 # endif
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
290 # ifndef HAVE_FEUPDATEENV
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
291 # define feupdateenv(arg) ((void) (arg))
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
292 # endif
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
293 # ifndef HAVE_FESETENV
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
294 # define fesetenv(arg) ((void) (arg))
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
295 # endif
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
296 # ifndef HAVE_FETESTEXCEPT
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
297 # define fetestexcept(arg) 0
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
298 # endif
111
kono
parents: 68
diff changeset
299 #endif
145
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
300
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
301 #ifndef __glibc_likely
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
302 # define __glibc_likely(cond) __builtin_expect ((cond), 1)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
303 #endif
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
304
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
305 #ifndef __glibc_unlikely
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
306 # define __glibc_unlikely(cond) __builtin_expect ((cond), 0)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
307 #endif
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
308
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
309 #if defined HAVE_FENV_H && defined HAVE_FESETROUND && defined HAVE_FEUPDATEENV
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
310 struct rm_ctx
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
311 {
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
312 fenv_t env;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
313 bool updated_status;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
314 };
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
315
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
316 # define SET_RESTORE_ROUNDF128(RM) \
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
317 struct rm_ctx ctx __attribute__((cleanup (libc_feresetround_ctx))); \
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
318 libc_feholdsetround_ctx (&ctx, (RM))
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
319
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
320 static inline __attribute__ ((always_inline)) void
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
321 libc_feholdsetround_ctx (struct rm_ctx *ctx, int round)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
322 {
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
323 ctx->updated_status = false;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
324
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
325 /* Update rounding mode only if different. */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
326 if (__glibc_unlikely (round != fegetround ()))
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
327 {
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
328 ctx->updated_status = true;
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
329 fegetenv (&ctx->env);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
330 fesetround (round);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
331 }
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
332 }
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
333
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
334 static inline __attribute__ ((always_inline)) void
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
335 libc_feresetround_ctx (struct rm_ctx *ctx)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
336 {
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
337 /* Restore the rounding mode if updated. */
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
338 if (__glibc_unlikely (ctx->updated_status))
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
339 feupdateenv (&ctx->env);
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
340 }
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
341 #else
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
342 # define SET_RESTORE_ROUNDF128(RM) ((void) 0)
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
343 #endif
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
344
1830386684a0 gcc-9.2.0
anatofuz
parents: 111
diff changeset
345 #endif