]>
Commit | Line | Data |
---|---|---|
302bb1a3 | 1 | /* _Float128 overrides for building ldbl-128 as _Float128. |
dff8da6b | 2 | Copyright (C) 2017-2024 Free Software Foundation, Inc. |
302bb1a3 PM |
3 | This file is part of the GNU C Library. |
4 | ||
5 | The GNU C Library is free software; you can redistribute it and/or | |
6 | modify it under the terms of the GNU Lesser General Public | |
7 | License as published by the Free Software Foundation; either | |
8 | version 2.1 of the License, or (at your option) any later version. | |
9 | ||
10 | The GNU C Library is distributed in the hope that it will be useful, | |
11 | but WITHOUT ANY WARRANTY; without even the implied warranty of | |
12 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
13 | Lesser General Public License for more details. | |
14 | ||
15 | You should have received a copy of the GNU Lesser General Public | |
16 | License along with the GNU C Library; if not, see | |
5a82c748 | 17 | <https://www.gnu.org/licenses/>. */ |
302bb1a3 PM |
18 | |
19 | /* This must be included before the function renames below. */ | |
32bf1d09 | 20 | #include <gmp.h> |
302bb1a3 PM |
21 | #include <math.h> |
22 | #undef HUGE_VALL | |
23 | #define HUGE_VALL HUGE_VAL_F128 | |
24 | #include <math/mul_splitl.h> | |
25 | ||
db7a548d JM |
26 | /* This must be included before the renames of types and macros from |
27 | it. */ | |
28 | #include <ieee754.h> | |
29 | ||
302bb1a3 PM |
30 | /* Renames derived from math_private.h. */ |
31 | #include <math_private.h> | |
70e2ba33 | 32 | #include <fenv_private.h> |
302bb1a3 PM |
33 | #include <ieee754_float128.h> |
34 | #define ieee854_long_double_shape_type ieee854_float128_shape_type | |
35 | #define ieee854_long_double ieee854_float128 | |
36 | ||
37 | #undef GET_LDOUBLE_LSW64 | |
38 | #undef GET_LDOUBLE_MSW64 | |
39 | #undef GET_LDOUBLE_WORDS64 | |
40 | #undef SET_LDOUBLE_LSW64 | |
41 | #undef SET_LDOUBLE_MSW64 | |
42 | #undef SET_LDOUBLE_WORDS64 | |
43 | #define GET_LDOUBLE_LSW64(x,y) GET_FLOAT128_LSW64 (x, y) | |
44 | #define GET_LDOUBLE_MSW64(x,y) GET_FLOAT128_MSW64 (x, y) | |
45 | #define GET_LDOUBLE_WORDS64(x,y,z) GET_FLOAT128_WORDS64 (x, y, z) | |
46 | #define SET_LDOUBLE_LSW64(x,y) SET_FLOAT128_LSW64 (x, y) | |
47 | #define SET_LDOUBLE_MSW64(x,y) SET_FLOAT128_MSW64 (x, y) | |
48 | #define SET_LDOUBLE_WORDS64(x,y,z) SET_FLOAT128_WORDS64 (x, y, z) | |
49 | ||
50 | #undef IEEE854_LONG_DOUBLE_BIAS | |
51 | #define IEEE854_LONG_DOUBLE_BIAS IEEE854_FLOAT128_BIAS | |
52 | ||
53 | #ifdef SET_RESTORE_ROUNDF128 | |
54 | # undef SET_RESTORE_ROUNDL | |
f20079d2 | 55 | # define SET_RESTORE_ROUNDL(RM) SET_RESTORE_ROUNDF128 (RM) |
302bb1a3 PM |
56 | #endif |
57 | ||
63716ab2 JM |
58 | #ifdef libc_feholdexcept_setroundf128 |
59 | # undef libc_feholdexcept_setroundl | |
60 | # define libc_feholdexcept_setroundl(ENV, RM) \ | |
61 | libc_feholdexcept_setroundf128 (ENV, RM) | |
62 | #endif | |
63 | ||
64 | #ifdef libc_feupdateenv_testf128 | |
65 | # undef libc_feupdateenv_testl | |
66 | # define libc_feupdateenv_testl(ENV, EX) libc_feupdateenv_testf128 (ENV, EX) | |
67 | #endif | |
302bb1a3 PM |
68 | |
69 | /* misc macros from the header below. */ | |
70 | #include <fix-fp-int-convert-overflow.h> | |
71 | #undef FIX_LDBL_LONG_CONVERT_OVERFLOW | |
72 | #undef FIX_LDBL_LLONG_CONVERT_OVERFLOW | |
73 | #define FIX_LDBL_LONG_CONVERT_OVERFLOW FIX_FLT128_LONG_CONVERT_OVERFLOW | |
74 | #define FIX_LDBL_LLONG_CONVERT_OVERFLOW FIX_FLT128_LLONG_CONVERT_OVERFLOW | |
75 | ||
76 | ||
77 | /* float.h constants. */ | |
78 | #include <float.h> | |
79 | #undef LDBL_DIG | |
80 | #undef LDBL_EPSILON | |
81 | #undef LDBL_MANT_DIG | |
82 | #undef LDBL_MAX | |
83 | #undef LDBL_MAX_10_EXP | |
84 | #undef LDBL_MAX_EXP | |
85 | #undef LDBL_MIN | |
86 | #undef LDBL_MIN_10_EXP | |
87 | #undef LDBL_MIN_EXP | |
88 | #undef LDBL_TRUE_MIN | |
89 | #define LDBL_DIG FLT128_DIG | |
90 | #define LDBL_EPSILON FLT128_EPSILON | |
91 | #define LDBL_MANT_DIG FLT128_MANT_DIG | |
92 | #define LDBL_MAX FLT128_MAX | |
93 | #define LDBL_MAX_10_EXP FLT128_MAX_10_EXP | |
94 | #define LDBL_MAX_EXP FLT128_MAX_EXP | |
95 | #define LDBL_MIN FLT128_MIN | |
96 | #define LDBL_MIN_10_EXP FLT128_MIN_10_EXP | |
97 | #define LDBL_MIN_EXP FLT128_MIN_EXP | |
98 | #define LDBL_TRUE_MIN FLT128_TRUE_MIN | |
99 | ||
100 | ||
101 | /* math.h GNU constants. */ | |
102 | #undef M_El | |
103 | #undef M_LOG2El | |
104 | #undef M_LOG10El | |
105 | #undef M_LN2l | |
106 | #undef M_LN10l | |
107 | #undef M_PIl | |
108 | #undef M_PI_2l | |
109 | #undef M_PI_4l | |
110 | #undef M_1_PIl | |
111 | #undef M_2_PIl | |
112 | #undef M_2_SQRTPIl | |
113 | #undef M_SQRT2l | |
114 | #undef M_SQRT1_2l | |
115 | #define M_El M_Ef128 | |
116 | #define M_LOG2El M_LOG2Ef128 | |
117 | #define M_LOG10El M_LOG10Ef128 | |
118 | #define M_LN2l M_LN2f128 | |
119 | #define M_LN10l M_LN10f128 | |
120 | #define M_PIl M_PIf128 | |
121 | #define M_PI_2l M_PI_2f128 | |
122 | #define M_PI_4l M_PI_4f128 | |
123 | #define M_1_PIl M_1_PIf128 | |
124 | #define M_2_PIl M_2_PIf128 | |
125 | #define M_2_SQRTPIl M_2_SQRTPIf128 | |
126 | #define M_SQRT2l M_SQRT2f128 | |
127 | #define M_SQRT1_2l M_SQRT1_2f128 | |
128 | ||
129 | ||
fd3b4e7c JM |
130 | #include <libm-alias-ldouble.h> |
131 | #include <libm-alias-float128.h> | |
132 | #undef libm_alias_ldouble_r | |
133 | #define libm_alias_ldouble_r(from, to, r) libm_alias_float128_r (from, to, r) | |
134 | ||
135 | ||
63716ab2 JM |
136 | #include <math/math-narrow.h> |
137 | #undef libm_alias_float_ldouble | |
138 | #define libm_alias_float_ldouble(func) libm_alias_float32_float128 (func) | |
139 | #undef libm_alias_double_ldouble | |
140 | #define libm_alias_double_ldouble(func) libm_alias_float64_float128 (func) | |
abd38358 JM |
141 | #undef libm_alias_ldouble_narrow |
142 | #define libm_alias_ldouble_narrow(from, to) \ | |
143 | libm_alias_float128_narrow (from, to) | |
63716ab2 | 144 | |
ae3577f6 SL |
145 | #include <math-use-builtins.h> |
146 | #undef USE_NEARBYINTL_BUILTIN | |
147 | #define USE_NEARBYINTL_BUILTIN USE_NEARBYINTF128_BUILTIN | |
a2a9b004 SL |
148 | #undef USE_RINTL_BUILTIN |
149 | #define USE_RINTL_BUILTIN USE_RINTF128_BUILTIN | |
6c1b6a5e SL |
150 | #undef USE_FLOORL_BUILTIN |
151 | #define USE_FLOORL_BUILTIN USE_FLOORF128_BUILTIN | |
62560ee8 SL |
152 | #undef USE_CEILL_BUILTIN |
153 | #define USE_CEILL_BUILTIN USE_CEILF128_BUILTIN | |
1ac9c1cf SL |
154 | #undef USE_TRUNCL_BUILTIN |
155 | #define USE_TRUNCL_BUILTIN USE_TRUNCF128_BUILTIN | |
f82996f8 SL |
156 | #undef USE_ROUNDL_BUILTIN |
157 | #define USE_ROUNDL_BUILTIN USE_ROUNDF128_BUILTIN | |
eb906620 STH |
158 | #undef USE_ROUNDEVENL_BUILTIN |
159 | #define USE_ROUNDEVENL_BUILTIN USE_ROUNDEVENF128_BUILTIN | |
f818afdd SL |
160 | #undef USE_COPYSIGNL_BUILTIN |
161 | #define USE_COPYSIGNL_BUILTIN USE_COPYSIGNF128_BUILTIN | |
628d90c5 VG |
162 | #undef USE_FMAL_BUILTIN |
163 | #define USE_FMAL_BUILTIN USE_FMAF128_BUILTIN | |
2b23ab1f XT |
164 | #undef USE_LRINTL_BUILTIN |
165 | #define USE_LRINTL_BUILTIN USE_LRINTF128_BUILTIN | |
a1981ecb XT |
166 | #undef USE_LLRINTL_BUILTIN |
167 | #define USE_LLRINTL_BUILTIN USE_LLRINTF128_BUILTIN | |
2e2485ce XT |
168 | #undef USE_LOGBL_BUILTIN |
169 | #define USE_LOGBL_BUILTIN USE_LOGBF128_BUILTIN | |
63716ab2 | 170 | |
302bb1a3 PM |
171 | /* IEEE function renames. */ |
172 | #define __ieee754_acoshl __ieee754_acoshf128 | |
173 | #define __ieee754_acosl __ieee754_acosf128 | |
174 | #define __ieee754_asinhl __ieee754_asinhf128 | |
175 | #define __ieee754_asinl __ieee754_asinf128 | |
176 | #define __ieee754_atan2l __ieee754_atan2f128 | |
177 | #define __ieee754_atanhl __ieee754_atanhf128 | |
178 | #define __ieee754_coshl __ieee754_coshf128 | |
179 | #define __ieee754_cosl __ieee754_cosf128 | |
180 | #define __ieee754_exp10l __ieee754_exp10f128 | |
181 | #define __ieee754_exp2l __ieee754_exp2f128 | |
182 | #define __ieee754_expl __ieee754_expf128 | |
183 | #define __ieee754_fmodl __ieee754_fmodf128 | |
184 | #define __ieee754_gammal_r __ieee754_gammaf128_r | |
185 | #define __ieee754_hypotl __ieee754_hypotf128 | |
186 | #define __ieee754_ilogbl __ieee754_ilogbf128 | |
187 | #define __ieee754_j0l __ieee754_j0f128 | |
188 | #define __ieee754_j1l __ieee754_j1f128 | |
189 | #define __ieee754_jnl __ieee754_jnf128 | |
190 | #define __ieee754_lgammal_r __ieee754_lgammaf128_r | |
191 | #define __ieee754_log10l __ieee754_log10f128 | |
192 | #define __ieee754_log2l __ieee754_log2f128 | |
193 | #define __ieee754_logl __ieee754_logf128 | |
194 | #define __ieee754_powl __ieee754_powf128 | |
195 | #define __ieee754_rem_pio2l __ieee754_rem_pio2f128 | |
196 | #define __ieee754_remainderl __ieee754_remainderf128 | |
197 | #define __ieee754_sinhl __ieee754_sinhf128 | |
198 | #define __ieee754_sqrtl __ieee754_sqrtf128 | |
199 | #define __ieee754_y0l __ieee754_y0f128 | |
200 | #define __ieee754_y1l __ieee754_y1f128 | |
201 | #define __ieee754_ynl __ieee754_ynf128 | |
202 | ||
203 | ||
204 | /* finite math entry points. */ | |
205 | #define __acoshl_finite __acoshf128_finite | |
206 | #define __acosl_finite __acosf128_finite | |
207 | #define __asinl_finite __asinf128_finite | |
208 | #define __atan2l_finite __atan2f128_finite | |
209 | #define __atanhl_finite __atanhf128_finite | |
210 | #define __coshl_finite __coshf128_finite | |
211 | #define __cosl_finite __cosf128_finite | |
212 | #define __exp10l_finite __exp10f128_finite | |
213 | #define __exp2l_finite __exp2f128_finite | |
214 | #define __expl_finite __expf128_finite | |
215 | #define __fmodl_finite __fmodf128_finite | |
216 | #define __hypotl_finite __hypotf128_finite | |
217 | #define __ilogbl_finite __ilogbf128_finite | |
218 | #define __j0l_finite __j0f128_finite | |
219 | #define __j1l_finite __j1f128_finite | |
220 | #define __jnl_finite __jnf128_finite | |
221 | #define __lgammal_r_finite __lgammaf128_r_finite | |
222 | #define __log10l_finite __log10f128_finite | |
223 | #define __log2l_finite __log2f128_finite | |
224 | #define __logl_finite __logf128_finite | |
225 | #define __powl_finite __powf128_finite | |
226 | #define __remainderl_finite __remainderf128_finite | |
227 | #define __sinhl_finite __sinhf128_finite | |
228 | #define __y0l_finite __y0f128_finite | |
229 | #define __y1l_finite __y1f128_finite | |
230 | #define __ynl_finite __ynf128_finite | |
231 | ||
232 | ||
233 | /* internal function names. */ | |
234 | #define __asinhl __asinhf128 | |
235 | #define __atanl __atanf128 | |
236 | #define __cbrtl __cbrtf128 | |
237 | #define __ceill __ceilf128 | |
238 | #define __copysignl __copysignf128 | |
239 | #define __cosl __cosf128 | |
240 | #define __erfcl __erfcf128 | |
241 | #define __erfl __erff128 | |
242 | #define __expl __expf128 | |
243 | #define __expm1l __expm1f128 | |
244 | #define __fabsl __fabsf128 | |
245 | #define __fdiml __fdimf128 | |
246 | #define __finitel __finitef128 | |
247 | #define __floorl __floorf128 | |
248 | #define __fmal __fmaf128 | |
249 | #define __fmaxl __fmaxf128 | |
250 | #define __fminl __fminf128 | |
251 | #define __fpclassifyl __fpclassifyf128 | |
252 | #define __frexpl __frexpf128 | |
253 | #define __gammal_r_finite __gammaf128_r_finite | |
01f28812 | 254 | #define __getpayloadl __getpayloadf128 |
302bb1a3 PM |
255 | #define __isinfl __isinff128 |
256 | #define __isnanl __isnanf128 | |
257 | #define __issignalingl __issignalingf128 | |
258 | #define __ldexpl __ldexpf128 | |
259 | #define __llrintl __llrintf128 | |
260 | #define __llroundl __llroundf128 | |
261 | #define __log1pl __log1pf128 | |
262 | #define __logbl __logbf128 | |
263 | #define __logl __logf128 | |
264 | #define __lrintl __lrintf128 | |
265 | #define __lroundl __lroundf128 | |
266 | #define __modfl __modff128 | |
267 | #define __nearbyintl __nearbyintf128 | |
268 | #define __nextafterl __nextafterf128 | |
269 | #define __nextdownl __nextdownf128 | |
270 | #define __nextupl __nextupf128 | |
271 | #define __remquol __remquof128 | |
272 | #define __rintl __rintf128 | |
01f28812 | 273 | #define __roundevenl __roundevenf128 |
302bb1a3 PM |
274 | #define __roundl __roundf128 |
275 | #define __scalblnl __scalblnf128 | |
276 | #define __scalbnl __scalbnf128 | |
277 | #define __signbitl __signbitf128 | |
278 | #define __sincosl __sincosf128 | |
279 | #define __sinl __sinf128 | |
280 | #define __sqrtl __sqrtf128 | |
281 | #define __tanhl __tanhf128 | |
282 | #define __tanl __tanf128 | |
01f28812 | 283 | #define __totalorderl __totalorderf128 |
42760d76 | 284 | #define __totalorder_compatl __totalorder_compatf128 |
01f28812 | 285 | #define __totalordermagl __totalordermagf128 |
42760d76 | 286 | #define __totalordermag_compatl __totalordermag_compatf128 |
302bb1a3 PM |
287 | #define __truncl __truncf128 |
288 | #define __x2y2m1l __x2y2m1f128 | |
289 | ||
d8742dd8 JM |
290 | #define __faddl __f32addf128 |
291 | #define __daddl __f64addf128 | |
632a6cbe JM |
292 | #define __fdivl __f32divf128 |
293 | #define __ddivl __f64divf128 | |
b3f27d81 JM |
294 | #define __ffmal __f32fmaf128 |
295 | #define __dfmal __f64fmaf128 | |
69a01461 JM |
296 | #define __fmull __f32mulf128 |
297 | #define __dmull __f64mulf128 | |
abd38358 JM |
298 | #define __fsqrtl __f32sqrtf128 |
299 | #define __dsqrtl __f64sqrtf128 | |
8d3f9e85 JM |
300 | #define __fsubl __f32subf128 |
301 | #define __dsubl __f64subf128 | |
d8742dd8 | 302 | |
220622dd WD |
303 | /* Used on __finite compat alias. */ |
304 | #define __acosl __acosf128 | |
305 | #define __acoshl __acoshf128 | |
306 | #define __asinl __asinf128 | |
307 | #define __atan2l __atan2f128 | |
308 | #define __atanhl __atanhf128 | |
309 | #define __coshl __coshf128 | |
310 | #define __exp10l __exp10f128 | |
311 | #define __expl __expf128 | |
312 | #define __fmodl __fmodf128 | |
313 | #define __gammal_r __gammaf128_r | |
314 | #define __hypotl __hypotf128 | |
315 | #define __j0l __j0f128 | |
316 | #define __j1l __j1f128 | |
317 | #define __jnl __jnf128 | |
318 | #define __lgammal_r __lgammaf128_r | |
319 | #define __log10l __log10f128 | |
320 | #define __log2l __log2f128 | |
321 | #define __logl __logf128 | |
322 | #define __powl __powf128 | |
323 | #define __remainderl __remainderf128 | |
324 | #define __sinhl __sinhf128 | |
325 | #define __y0l __y0f128 | |
326 | #define __y1l __y1f128 | |
327 | #define __ynl __ynf128 | |
328 | ||
302bb1a3 PM |
329 | /* __nexttowardf128 is not _Float128 API. */ |
330 | #define __nexttowardl __nexttowardf128_do_not_use | |
331 | #define nexttowardl nexttowardf128_do_not_use | |
332 | ||
333 | ||
334 | /* public entry points. */ | |
335 | #define asinhl asinhf128 | |
336 | #define atanl atanf128 | |
337 | #define cbrtl cbrtf128 | |
338 | #define ceill ceilf128 | |
339 | #define copysignl copysignf128 | |
340 | #define cosl cosf128 | |
341 | #define erfcl erfcf128 | |
342 | #define erfl erff128 | |
343 | #define expl expf128 | |
344 | #define expm1l expm1f128 | |
345 | #define fabsl fabsf128 | |
346 | #define fdiml fdimf128 | |
347 | #define finitel finitef128_do_not_use | |
348 | #define floorl floorf128 | |
349 | #define fmal fmaf128 | |
350 | #define fmaxl fmaxf128 | |
351 | #define fminl fminf128 | |
352 | #define frexpl frexpf128 | |
353 | #define getpayloadl getpayloadf128 | |
354 | #define isinfl isinff128_do_not_use | |
355 | #define isnanl isnanf128_do_not_use | |
356 | #define ldexpl ldexpf128 | |
357 | #define llrintl llrintf128 | |
358 | #define llroundl llroundf128 | |
359 | #define log1pl log1pf128 | |
360 | #define logbl logbf128 | |
bb014f50 | 361 | #define logp1l logp1f128 |
302bb1a3 PM |
362 | #define logl logf128 |
363 | #define lrintl lrintf128 | |
364 | #define lroundl lroundf128 | |
365 | #define modfl modff128 | |
366 | #define nanl nanf128 | |
367 | #define nearbyintl nearbyintf128 | |
368 | #define nextafterl nextafterf128 | |
369 | #define nextdownl nextdownf128 | |
370 | #define nextupl nextupf128 | |
371 | #define remquol remquof128 | |
372 | #define rintl rintf128 | |
373 | #define roundevenl roundevenf128 | |
374 | #define roundl roundf128 | |
375 | #define scalbnl scalbnf128 | |
376 | #define sincosl sincosf128 | |
377 | #define sinl sinf128 | |
378 | #define sqrtl sqrtf128 | |
379 | #define tanhl tanhf128 | |
380 | #define tanl tanf128 | |
381 | #define totalorderl totalorderf128 | |
382 | #define totalordermagl totalordermagf128 | |
383 | #define truncl truncf128 | |
384 | ||
385 | ||
386 | /* misc internal renames. */ | |
387 | #define __builtin_fmal __builtin_fmaf128 | |
388 | #define __expl_table __expf128_table | |
389 | #define __gamma_productl __gamma_productf128 | |
390 | #define __kernel_cosl __kernel_cosf128 | |
391 | #define __kernel_rem_pio2l __kernel_rem_pio2f128 | |
392 | #define __kernel_sincosl __kernel_sincosf128 | |
393 | #define __kernel_sinl __kernel_sinf128 | |
394 | #define __kernel_tanl __kernel_tanf128 | |
395 | #define __lgamma_negl __lgamma_negf128 | |
396 | #define __lgamma_productl __lgamma_productf128 | |
32bf1d09 | 397 | #define __mpn_extract_long_double __mpn_extract_float128 |
302bb1a3 PM |
398 | #define __sincosl_table __sincosf128_table |
399 | #define mul_splitl mul_splitf128 | |
400 | ||
401 | /* Builtin renames. */ | |
402 | #define __builtin_copysignl __builtin_copysignf128 | |
403 | #define __builtin_signbitl __builtin_signbit | |
ae3577f6 | 404 | #define __builtin_nearbyintl __builtin_nearbyintf128 |
a2a9b004 | 405 | #define __builtin_rintl __builtin_rintf128 |
6c1b6a5e | 406 | #define __builtin_floorl __builtin_floorf128 |
62560ee8 | 407 | #define __builtin_ceill __builtin_ceilf128 |
1ac9c1cf | 408 | #define __builtin_truncl __builtin_truncf128 |
f82996f8 | 409 | #define __builtin_roundl __builtin_roundf128 |
8f27dc1a JM |
410 | #define __builtin_roundevenl __builtin_roundevenf128 |
411 | #define __builtin_lrintl __builtin_lrintf128 | |
412 | #define __builtin_llrintl __builtin_llrintf128 | |
f818afdd | 413 | #define __builtin_copysignl __builtin_copysignf128 |
8f27dc1a | 414 | #define __builtin_fabsl __builtin_fabsf128 |
302bb1a3 PM |
415 | |
416 | /* Get the constant suffix from bits/floatn-compat.h. */ | |
417 | #define L(x) __f128 (x) | |
418 | ||
419 | static inline void | |
420 | mul_splitf128 (_Float128 *hi, _Float128 *lo, _Float128 x, _Float128 y) | |
421 | { | |
422 | #ifdef __FP_FAST_FMAF128 | |
423 | /* Fast built-in fused multiply-add. */ | |
424 | *hi = x * y; | |
425 | *lo = __builtin_fmal (x, y, -*hi); | |
426 | #else | |
427 | /* Apply Dekker's algorithm. */ | |
428 | *hi = x * y; | |
429 | # define C ((1LL << (FLT128_MANT_DIG + 1) / 2) + 1) | |
430 | _Float128 x1 = x * C; | |
431 | _Float128 y1 = y * C; | |
432 | # undef C | |
433 | x1 = (x - x1) + x1; | |
434 | y1 = (y - y1) + y1; | |
435 | _Float128 x2 = x - x1; | |
436 | _Float128 y2 = y - y1; | |
437 | *lo = (((x1 * y1 - *hi) + x1 * y2) + x2 * y1) + x2 * y2; | |
438 | #endif | |
439 | } |